15 svar
235 visningar
naytte behöver inte mer hjälp
naytte Online 8369 – Moderator
Postad: 11 aug 19:46 Redigerad: 13 aug 15:33

Komplex analys - lösa integral med residykalkyl

Hej!

Jag sitter med b)-uppgiften nedan från en gammal tentamen i komplex analys och skulle behöva lite vägledning och även hjälp att se om mina logiska slutledningnar är korrekta.

Vi inleder med att beteckna den sökta integralen med II och att definiera en komplex funktion ff enligt

fz=11+z2N\displaystyle f\left(z\right)=\frac{1}{\left(1+z^2\right)^N}

samt två kurvor i planet enligt

γ1R:zt=t,t∈-R,Rγ2R:zθ=Reiθ,θ∈0,π\displaystyle \left\{\begin{array}{ll}\gamma_1\left(R\right):z\left(t\right)=t,& t\in\left[-R,R\right]\\gamma_2\left(R\right): z\left(\theta\right)=Re^{i\theta},& \theta \in\left[0,\pi\right] \end{array}\right.

Notera speciellt att

I=limR→∞∫γ1Rfzdz=:limR→∞IR\displaystyle I=\lim_{R\to \infty} \int_{\gamma_1\left(R\right)}f\left(z\right)dz=:\lim_{R\to\infty} I_R

Eftersom ΓR:=γ1R∪γ2R\Gamma\left(R\right) :=\gamma_1\left(R\right) \cup \gamma_2\left(R\right) är en enkel, positivt orienterad och sluten kurva, inom och på vilken ff är analytisk, förutom för tillräckligt stora RR vid en uppräknelig mängd isolerade singulariteter z0,…,znz_0,\dots, z_n, gäller det enligt Cauchys residysats att

∮ΓRfzdz=2πi∑j=1nReszj\displaystyle \oint_{\Gamma\left(R\right)}f\left(z\right)dz=2\pi i\sum_{j=1}^{n}\mathrm{Res}\left(z_j\right)

För tillräckligt stora RR ser vi att ff har en pol av ordning NN inom ΓR\Gamma\left(R\right), nämligen vid z=iz=i. Residyn i denna punkt ges av

Resi=limz→i1N-1!dN-1dzN-1z-iNfz=limz→i1N-1!dN-1dzN-1z+i-N\displaystyle \mathrm{Res}\left(i\right)= \lim_{z\to i}\frac{1}{\left(N-1\right)!}\frac{d ^{N-1}}{{dz}^{N-1}}\left(\left(z-i\right)^N f\left(z\right)\right)=\lim_{z\to i}\frac{1}{\left(N-1\right)!}\frac{d ^{N-1}}{{dz}^{N-1}}\left(z+i\right)^{-N}

Vi inför en variabel k∈N>0k\in \mathbb{N}_{\gt 0} och härleder ett allmänt uttryck för derivatan för kk. Sedan kan vi stoppa in k=N-1k=N-1. Vi ser att

dkdzkz+i-N=-1kN+k-1!N-1!z+i-N+k\displaystyle \frac{d ^{k}}{{dz}^{k}}\left(z+i\right)^{-N}=\frac{\left(-1\right)^k\left(N+k-1\right)!}{\left(N-1\right)!}\left(z+i\right)^{-\left(N+k\right)}

Så då k=N-1k=N-1 har vi

dN-1dzN-1z+i-N=-1N-12N-2!N-1!z+i-2N-1\displaystyle \frac{d ^{N-1}}{{dz}^{N-1}}\left(z+i\right)^{-N}=\frac{\left(-1\right)^{N-1}\left(2N-2\right)!}{\left(N-1\right)!}\left(z+i\right)^{-\left(2N-1\right)}

Gränsvärdet utvärderas nu trivialt:

limz→i1N-1!dN-1dzN-1z-iNfz=-1N-12N-2!22N-1N-1!2i-(2N-1)\displaystyle \lim_{z\to i}\frac{1}{\left(N-1\right)!}\frac{d ^{N-1}}{{dz}^{N-1}}\left(\left(z-i\right)^N f\left(z\right)\right)=\frac{\left(-1\right)^{N-1}\left(2N-2\right)!}{2^{2N-1}\left(\left(N-1\right)!\right)^2}i^{-(2N-1)}

Om vi återgår till Residysatsen har vi nu alltså

Resi=2πi-1N-12N-2!22N-1N-1!2i-(2N-1)=π-1N-12N-2!22N-2N-1!2i2N-2=π2N-2!22N-2N-1!2\displaystyle \mathrm{Res}\left(i\right) = 2\pi i\frac{\left(-1\right)^{N-1}\left(2N-2\right)!}{2^{2N-1}\left(\left(N-1\right)!\right)^2}i^{-(2N-1)}=\pi \frac{\left(-1\right)^{N-1}\left(2N-2\right)!}{2^{2N-2}\left(\left(N-1\right)!\right)^2i^{2N-2}}=\frac{\pi\left(2N-2\right)!}{2^{2N-2}\left(\left(N-1\right)!\right)^2}

Vi har då alltså:

I=limR→∞∫γ1Rfzdz=π2N-2!22N-2N-1!2-limR→∞∫γ2Rfzdz\displaystyle I=\lim_{R\to\infty}\int_{\gamma_1\left(R\right)}f\left(z\right)dz=\frac{\pi\left(2N-2\right)!}{2^{2N-2}\left(\left(N-1\right)!\right)^2}-\lim_{R\to\infty}\int_{\gamma_2\left(R\right)}f\left(z\right)dz

Nu gäller det bara att lyckas beräkna den andra integralen, men här är jag lite fast. Helst skulle jag vilja använda typ MLML-olikheten på något sätt men jag är osäker på de tekniska detaljerna.

LuMa07 769
Postad: 11 aug 20:14 Redigerad: 11 aug 20:37

Sätt in parametriseringen av kurvan γ2\gamma_2 och sedan är det "bara" triangelolikhet som behövs för att visa att gränsvärdet är noll.

∫γ2dz(1+z2)N=∫0πiR eiθ(1+R2ei2θ)N dθ\displaystyle \int_{\gamma_2} \frac{dz}{(1+z^2)^N} = \int_0^\pi \frac{iR\,e^{i\theta}}{(1+R^2 e^{i2\theta})^N}\,d\theta.

Triangelolikheten ger att |1+R2ei2θ|≥R2|ei2θ|-1=R2-1| 1+R^2 e^{i2\theta} | \ge R^2 |e^{i2\theta}| - 1 = R^2-1.

Visa spoiler

Därmed blir

∫0πiR eiθ(1+R2ei2θ)N dθ≤ ∫0πiR eiθ(1+R2ei2θ)N dθ≤∫0πR(R2-1)N dθ→0\displaystyle \left| \int_0^\pi \frac{iR\,e^{i\theta}}{(1+R^2 e^{i2\theta})^N}\,d\theta \right| \le  \int_0^\pi \left| \frac{iR\,e^{i\theta}}{(1+R^2 e^{i2\theta})^N} \right| \,d\theta \le \int_0^\pi \frac{R}{(R^2-1)^N} \,d\theta \to 0 då R→∞R\to\infty

naytte Online 8369 – Moderator
Postad: 12 aug 01:14 Redigerad: 12 aug 01:15

Jag hänger inte riktigt med på vad som händer här. Jag köper olikheten med beloppen. Men sedan hänger jag inte riktigt med. Är tanken att

∫0πiReiθ1+R2e2iθNdθ=∫0πiReiθ1+R2e2iθNdθ=∫0πR1+R2e2iθNdθ=∫0πR1+R2e2iθNdθ\displaystyle \int_{0}^{\pi}\left|\frac{iRe^{i\theta}}{\left(1+R^2e^{2i\theta}\right)^N}\right|d\theta =\int_{0}^{\pi}\frac{\left| iRe^{i\theta}\right|}{\left| \left(1+R^2e^{2i\theta}\right)^N\right |}d\theta=\int_{0}^{\pi}\frac{R}{\left|\left(1+R^2e^{2i\theta}\right)^N\right |}d\theta =\int_{0}^{\pi}\frac{R}{\left|1+R^2e^{2i\theta}\right |^N}d\theta

Och eftersom triangelolikheten ger att 1+R2e2iθ≥R2-1\left| 1+R^2e^{2i\theta}\right| \ge R^2-1, så måste

R1+R2e2iθN≤RR2-1N\displaystyle \frac{R}{\left|1+R^2e^{2i\theta}\right |^N} \le \frac{R}{\left(R^2-1\right)^N}

varför

∫0πR1+R2e2iθNdθ≤∫0πRR2-1Ndθ\displaystyle \int_{0}^{\pi}\frac{R}{\left|1+R^2e^{2i\theta}\right |^N} d\theta \le \int_{0}^{\pi}\frac{R}{\left(R^2-1\right)^N}d\theta

och eftersom något större eller lika stort än beloppet av vår integral går mot noll, måste beloppet av vår integral gå mot noll, och då också alltså själva integralen (beloppet kan enbart gå mot noll om integralen går mot noll)?

naytte Online 8369 – Moderator
Postad: 12 aug 01:42 Redigerad: 12 aug 01:42

Hur som helst, om vi vet att

limR→∞∫γ2Rfzdz=0\displaystyle \lim_{R\to\infty}\int_{\gamma_2\left(R\right)}f\left(z\right)dz=0

då har vi att

I=π2N-2!22N-2N-1!2\displaystyle I= \frac{\pi\left(2N-2\right)!}{2^{2N-2}\left(\left(N-1\right)!\right)^2}

Speciellt för N=2N=2 (den första deluppgiften) har vi att

I=π2\displaystyle I= \frac{\pi}{2}

Vilket är rätt svar.

Hur som helst undrar jag om mitt resonemang och mina slutledningar är logiskt korrekta? Jag har försökt motivera varje steg så noggrant som möjligt och explicit förklara varför t.ex. residysatsen är tillämpbar här.

LuMa07 769
Postad: 12 aug 05:37 Redigerad: 12 aug 05:53
naytte skrev:

...

och eftersom något större eller lika stort än beloppet av vår integral går mot noll, måste beloppet av vår integral gå mot noll, och då också alltså själva integralen (beloppet kan enbart gå mot noll om integralen går mot noll)?

Exakt.

Det gäller även för komplexvärda gränsvärden att

limR→∞I2(R)=0\displaystyle \lim_{R\to \infty} {I_2(R)} = 0 om och endast om limR→∞I2(R)=0\displaystyle \lim_{R\to \infty} \left| I_2(R)\right| = 0

Sedan är det bara den vanliga instängningssatsen för reellvärda gränsvärden:

0≤∫γ2(R)f(z) dz≤∫0πR(R2-1)N dθ=π R(R2-1)N⏟→0\displaystyle 0 \le \left | \int_{\gamma_2(R)} {f(z)}\,dz \right| \le \int_0^\pi \dfrac{R}{(R^2 - 1)^N}\,d\theta = \underbrace{\dfrac{\pi\,R}{(R^2 - 1)^N}}_{\to 0}

så

limR→∞0≤limR→∞∫γ2(R)f(z) dz≤limR→∞π R(R2-1)N \displaystyle \lim_{R\to\infty} 0 \le \lim_{R\to \infty}\left | \int_{\gamma_2(R)} {f(z)}\,dz \right| \le \lim_{R\to\infty} \dfrac{\pi\,R}{(R^2 - 1)^N}.


Btw: Det gick också att skaffa sig en ML-olikhet som dock också baseras på triangelolikheten:

|1+z2|≥|z2|-|1|=|z|2-1 ⇒|1(1+z2)N|=1|1+z2|N≤1(|z|2-1)N=1(R2-1)N=:M {|1 + z^2|} \ge {|z^2|} - {|1|} = {|z|^2} - 1  \implies |\dfrac{1}{(1+z^2)^N}| = \dfrac{1}{|1+z^2|^N} \le \dfrac{1}{(|z|^2 - 1)^N} = \dfrac{1}{(R^2-1)^N} =: M

(OBS: Olikhetshantering som ovan är korrekt endast under förutsättningen att |z|>1|z| > 1. Vi vill att R→∞R\to\infty, så det är inget problem egentligen.)

Halvcirkelns längd är då L:=π RL := \pi\,R, så ∫γ2(R)f(z) dz≤M·L=1(R2-1)N·πR\displaystyle \left | \int_{\gamma_2(R)} {f(z)}\,dz \right| \le M \cdot L = \dfrac{1}{(R^2-1)^N} \cdot \pi R.

Okej, tack!

En fråga till. Borde inte triangelolikheten ge följande?

1+R2e2iθ≤1+R2e2iθ=1+R2\displaystyle \left|1+R^2e^{2i\theta} \right| \le \left| 1\right| +\left| R^2e^{2i\theta}\right|=1+R^2

Jag förstår inte riktigt hur man kommer fram till din olikhet. Men hur som helst borde väl olikheten ovan ge samma svar ändå?

D4NIEL 3600
Postad: 12 aug 20:26 Redigerad: 12 aug 20:34

Det stämmer att |z1+z2|≤|z1|+|z2||z_1+z_2|\leq |z_1|+|z_2|. Men därifrån kan du härleda två till, sätt till exempel z2z_2 till -z2-z_2 och så har du 

|z1-z2|≤|z1|+|z2||z_1-z_2|\leq |z_1|+|z_2|

Och än viktigare, (till exempel genom att rita upp triangeln och studera möjliga sidlängder)

|z1+z2|≥|z1|-|z2|\displaystyle |z_1+z_2|\geq \left| | z_1|-|z_2|\right|

Tänk också på att du vill boxa in värdet åt rätt håll, mindre nämnare ger större kvot. Alltså ska du använda den tredje varianten.

Och så tips. I er kurs ingår förmodligen hjälpsatser så ni inte behöver genomföra ett ML-resonemang.

Till exempel finns det en färdig sats som säger att  om polynomen P(x)P(x) och Q(x)≠0Q(x)\neq 0 skiljer sig åt med grad två eller mer så är

∫-∞∞P(x)Q(x) dx=2πi∑resP(z)Q(z)\displaystyle \int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}\,dx=2\pi i \sum_{\mathrm{res}}\frac{P(z)}{Q(z)} för alla poler i övre halvplan (u.h.p).

Det är naturligtvis bra att öva ML, men på en tenta kan det vara smart att använda hjälpsatser för att spara tid och minska kognitiv last (om de ingår i kursen såklart)

naytte Online 8369 – Moderator
Postad: 12 aug 20:37 Redigerad: 12 aug 20:37

Tack för svar!

Jag är med på den andra olikheten, men hur kommer den tredje från triangelolikheten?

Dessutom, vad står res\text{res} för under summatecknet? Är de zz som kvantifieras över i summan polerna i ÖHP?

D4NIEL 3600
Postad: 12 aug 20:45 Redigerad: 12 aug 20:49

Ja, det är residyerna vid alla poler (summan av dem alltså)  i övre halvplanet  (upper half plane). Så satsen säger att har du kvoten av två polynom som skiljer sig åt med grad 2 eller mer så går det att summera residyerna vid polerna för att få fram ett värde på den reella integralen. Jag hoppas och tror att den är med i er kurs (den var med på min tid juh!)

Du kan bevisa olikheten geometriskt. Den säger egentligen bara att längden av summan av två vektorer är minst skillnaden mellan deras längder.

Vill du hellre bevisa det algebraiskt kan du sätta z1=(z1+z2)-z2z_1=(z_1+z_2)-z_2 i den andra olikheten olikheten. Sen subtraherar du |z2||z_2| osv

naytte Online 8369 – Moderator
Postad: 12 aug 20:48 Redigerad: 12 aug 20:48

Okej, jag ska testa att bevisa olikheten!

Vad gäller satsen som du hänvisar till, innebär det alltså man helt kan gå runt att använda residysatsen som jag gjorde i detta fall?

D4NIEL 3600
Postad: 12 aug 20:54

Nja, du måste fortfarande beräkna Res[f(z),i]\mathrm{Res}[f(z),i] precis som du gjorde. Men du slipper visa att integralen runt halvcirkeln blir noll.

naytte Online 8369 – Moderator
Postad: 12 aug 20:58 Redigerad: 12 aug 21:16

Jaha, jag missförstod vad satsen sade! Det är alltså

∫-∞∞P(x)Q(x)dx=2πi∑j=1nReszj\displaystyle \int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}dx=2\pi i\sum_{j=1}^n \mathrm{Res}\left(z_j\right)

där z1,…,znz_1,\dots, z_n är singulariteterna till integranden i det övre halvplanet?

I så fall tror jag också att jag förstår varför man har det kravet på graden. Jag antar att det man egentligen gör i beviset är att man väljer halvcirkeln (samma som jag) som sin kontur, och sedan visar man att bidraget från cirkelbågen går mot noll då radien på halvcirkeln går mot oändligheten? Det vill säga, det är ett specialfall av residysatsen? Om vi använder samma notation som jag gjorde i ursprungsinlägget visar man alltså att

limR→∞∫γ2RP(z)Q(z)dz=0\displaystyle \lim_{R\to\infty} \int_{\gamma_2\left(R\right)}\frac{P(z)}{Q(z)}dz=0

och kvar blir

∫-∞∞P(x)Q(x)dx=2πi∑jReszj\displaystyle \int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}dx=2\pi i\sum_{j}\mathrm{Res}\left(z_j\right)

där residyerna förstås ligger i övre halvplanet eftersom det var så vi valde halvcirkeln.

D4NIEL 3600
Postad: 12 aug 22:09

Just det! Och eftersom polynomen är godtyckliga kan man sedan använda den ganska brett. Det finns fler satser för mer komplicerade funktioner, men man kommer ganska långt med polynom.

naytte Online 8369 – Moderator
Postad: 12 aug 22:20 Redigerad: 12 aug 22:20

Ja, det kan säkert bli användbart!

Eftersom många av övningsfrågorna på äldre tentor är av detta slag kan det säkert vara en mycket användbar sats. Råkar du veta vad den heter (så att jag kan se om den ingår)?

Tack för hjälpen, förresten!

LuMa07 769
Postad: 13 aug 09:26
D4NIEL skrev:...

∫-∞∞P(x)Q(x) dx=2πi∑resP(z)Q(z)\displaystyle \int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}\,dx=2\pi i \sum_{\mathrm{res}}\frac{P(z)}{Q(z)} för alla poler i övre halvplan (u.h.p).

....

Exakt detta påstående är en övning i kursboken: Exercise 13 i Chapter 9 i Beck-Marchesi-Pixton-Sabalkas bok https://matthbeck.github.io/papers/complex.pdf (och a-uppgiften från tentan är Exercise 14)

Där ser man!

Det är faktiskt den boken vi har som kursbok, så borde det ju inte vara några problem för mig att använda satsen. Så länge jag hänvisar till "enligt lösningen till övning 13 i kapitel 9 i kursboken...".

Svara
Close