Definite Integration through Trigonometric Substitution on the Complex Plane
📂Complex AnaylsisDefinite Integration through Trigonometric Substitution on the Complex Plane
Theorem
∫02πf(cosθ,sinθ)dθ=∫Cf(z)dz=2πi∑Resf(z)
Description
Integrating real functions that are difficult to integrate can often be more straightforward through complex analysis. Among these, let’s look into integration techniques for integrands made up of trigonometric functions. The basic strategy involves changing the integration range to z(θ)=eiθ,0<θ<2π and taking the necessary parts from there. If needed, a slight manipulation to turn it into the form of −π<θ<π is also perfectly fine.
Viewing f as a function of cosθ and sinθ allows for the substitution of the trigonometric functions with sinz=2i1(z−z1) and cosz=21(z+z1), respectively. This is essentially using Euler’s formula z=eiθ=cosθ+isinθ in reverse to transform the integrand from a trigonometric function to a polynomial one. Once the integrand takes the form of a rational function, it becomes easy to find residues, so we can find the value using the residue theorem and then take the real or imaginary part as needed.
Example
As an example, let’s consider when I:=∫0π1−2acosθ+a2cos2θdθ,∣a∣<1, and find the value of I.
Solution
Since 1−2acosθ+a2cos2θ is an even function,
∫0π1−2acosθ+a2cos2θdθ=21∫−ππ1−2acosθ+a2cos2θdθ
substituting gives us dz=ieiθdθ⟺dθ=iz1dz, so
I==21∫C1−2a21(z+z1)+a2cos2θiz1dz2i1∫Cz−a(z2+1)+a2zcos2θdz
therefore, I becomes the real part of the complex integral 21∫Cz−a(z2+1)+a2z−iz2dz. Factoring the denominator simplifies calculations,
I=Re2a1∫C(z−a1)(z−a)iz2dz
by the Residue theorem,
2a1∫C(z−a1)(z−a)iz2dz=2a2πi∑Res(z−a1)(z−a)iz2
∣a∣<1, thus a is a singularity inside the unit circle C, while a1 lies outside C and does not need to be considered. Calculating the residue at the simple pole a gives us
Res(z−a1)(z−a)iz2=a−a1ia2=a2−1ia3
therefore
2a1∫C(z−a1)(z−a)iz2dz=aπia2−1ia3=1−a2πa2
Since I was the real part of 1−a2πa2, we get I=1−a2πa2.
This solution might seem overly complex and difficult at first glance. However, attempting the calculations without using complex analysis will reveal how much easier and more convenient this method actually is.
Furthermore, upon examining the unique shape of the example’s denominator, if we consider two vectors p and q, then ∣p−q∣ would be the length of the difference between the two vectors.
∣p−q∣2===(p−q)⋅(p−q)∣p∣2−2∣p∣∣q∣cosθ+∣q∣2∣p∣2(1−2∣p∣∣q∣cosθ+∣p∣2∣q∣2)
Setting it as a:=∣p∣∣q∣ gives us the shape we observed. Such a shape is frequently seen in physics, and questioning whether ‘integrals of this form are possible’ is undoubtedly important.