25 svar
262 visningar
naytte behöver inte mer hjälp
naytte 7419 – Moderator
Postad: 16 okt 18:58 Redigerad: 16 okt 19:00

Komplex Analys - lös reell integral med komplex integral

Halloj!

Jag studerar tentauppgiften nedan

Jag tror att det är meningen att man ska göra om detta till en kurvintegral i \mathbb{C}. Denna integral är formellt ett gränsvärde:

I=limR-RR1x4+6x2+5dx=:limRIR\displaystyle I=\lim_{R \to \infty} \int_{-R}^{R}\frac{1}{x^4+6x^2+5}dx=:\lim_{R \to \infty} I_R

Sedan kanske man kan tänka att den slutna kurvan vi vill integrera över i \mathbb{C} kan vara dem som parametriseras de två delkurvorna

γ1:z(t):=t,t[-R,R]\displaystyle \gamma_1: z(t) := t, t\in[-R, R]

γ2:z2(θ):=Reiθ,θ[0,π]\displaystyle \gamma_2: z_2(\theta) := Re^{i\theta}, \theta\in[0,\pi]

Då har vi

IR+något annat=γ1γ21z4+6z2+5dz\displaystyle I_R+\text{något annat}=\oint_{\gamma_1\cup \gamma_2}\frac{1}{z^4+6z ^2+5}dz

Enligt residysatsen vet vi att

γ1γ21z4+6z2+5:=f(z)dz=2πik=1nResf,ak\displaystyle \oint_{\gamma_1\cup \gamma_2}\underbrace{\frac{1}{z^4+6z ^2+5}}_{:=f(z)}dz=2\pi i\sum_{k=1}^{n}\mathrm{Res}\left(f,a_k\right)

där aka_k är de punkter i γ1γ2\gamma_1 \cup \gamma_2 där funktionen inte är holomorf. Hur ska man ta sig vidare härifrån? Jag försöker lära mig genom att lösa uppgifter eftersom jag finner kurslitteraturen väldigt svårsmält så jag ursäkter om detta är "enkelt".

AlexMu 940
Postad: 16 okt 19:07

Vi kan börja med att beräkna kurvintegralen och därmed residyerna, vart är funktionen inte holomorf?

naytte 7419 – Moderator
Postad: 16 okt 20:00

Där vi har division med noll. Alltså de punkter som satisfierar z4+6z2+5=0z^4+6z^2+5=0. Vi ser snabbt att z=iz=i är en sådan punkt så vi kan faktorisera med polynomdivision:

z4+6z2+5=(z-i)(z3+iz2+5z+5i)z^4+6z^2+5=(z-i)(z^3+iz^2+5z+5i)

Men hur ska man faktorisera resten? Gäller det bara att helt enkelt gissa på rötter?

AlexMu 940
Postad: 16 okt 20:20

I det här fallet beter sig ff ungefär som en andragradare, vi har ju inga zz-termer eller z3z^3-termer, så vi kan använda kvadratkomplettering/pq-formeln.

Med din gissade rot kan vi också använda satsen som säger att om ww löser en polynomekvation med reella koefficienter så löser w¯\bar w den också

naytte 7419 – Moderator
Postad: 16 okt 20:22

Ah, såklart! Vi gör en substitution z2=tz^2=t. Återkommer.

naytte 7419 – Moderator
Postad: 16 okt 20:25 Redigerad: 16 okt 20:25

Jag erhåller rötterna z1,2=±i5z_{1,2}=\pm i\sqrt{5} och z3,4=±iz_{3,4}=\pm i om jag har tänkt rätt.

naytte 7419 – Moderator
Postad: 16 okt 20:26 Redigerad: 16 okt 20:27

Men nu kommer vi till det steget jag har lite svårt att begripa. Hur ska man komma fram till Laurentserieutvecklingen runt dessa punkter? Måste man sitta och manuellt ta fram koefficienterna? Vi behöver ju koefficienten för termen av grad -1-1 för att beräkna residyerna (denna koefficient kanske till och med ÄR residyn?).

AlexMu 940
Postad: 16 okt 20:31

Den koefficienten ÄR residyn, ja. Hur tar du fram koefficienterna i nuläget? Det finns några ganska trevliga sätt att bestämma dem. 

naytte 7419 – Moderator
Postad: 16 okt 20:35 Redigerad: 16 okt 20:35

Det enda sättet jag är säker på är att kötta med en formel på formelbladet, men det hinner man inte under en tentamen. Jag vet att det finns lite olika regler beroende på vilken typ av pol man diskuterar. Dessa fyra poler skulle väl klassas som enkla poler och om jag minns rätt från föreläsningarna ges residyn i en enkel pol z=cz=c av

limzcz-cfz\lim_{z\to c}\left(z-c\right)f\left(z\right)

Är polerna i den här uppgiften "enkla"?

AlexMu 940
Postad: 16 okt 20:46 Redigerad: 16 okt 20:55

Dessa poler är enkla, ja. Om z0z_0 är en pol till ff gäller det även att z0z_0 är en rot till 1/f1/f. En pols "ordning/multiplicitet" är densamma som multipliciteten hos roten till 1/f1/f. Eftersom alla rötter du har tagit fram till z4+6z2+5z^4 + 6z^2 + 5 har multiplicitet 1 så gäller det även att de polerna är enkla/har ordning 1. 

Om vi istället hade funktionen 1z2\frac{1}{z^2} skulle polen z=0z=0 vara av ordning 2, eftersom roten till det polynomet har multiplicitet 2. 


Iden bakom formeln du skrev är att i punkten cc (enkel pol till ff) så ges laurentutvecklingen av

fz=b1z-c+a0+a1z-c+\displaystyle f\left(z\right) = \frac{b_1}{z-c} + a_0 + a_1\left(z-c\right) + \dots

För att lösa ut residyn som jag döpt här till b1b_1, kan vi multiplicera leden med z-cz-c:

fzz-c=b1+a0z-c+a1z-c2+\displaystyle f\left(z\right)\left(z-c\right)=b_1 + a_0\left(z-c\right)+a_1\left(z-c\right)^2 + \dots

Sedan när zcz\to c får vi bara b1b_1 kvar i högerled!

Det enda sättet jag är säker på är att kötta med en formel på formelbladet,

Vi får väl inga formelblad?

naytte 7419 – Moderator
Postad: 17 okt 12:31

Nej, jag menade bara att jag kunde en formel.

Anto 351
Postad: 17 okt 13:52 Redigerad: 17 okt 13:53

Det finns ju fler Residyregler. En säger exempelvis att du kan behålla täljaren, derivera nämnaren och sedan sätta in singulariteterna. Laurent är alltså inte nödvändigt.

naytte 7419 – Moderator
Postad: 17 okt 22:31 Redigerad: 17 okt 22:31

Så blir då t.ex. residyn i punkten z=i5z=i\sqrt{5}:

Resf,z=i5=1(2i5)(i5-i)(i5+i)\displaystyle \mathrm{Res}\left(f,z=i\sqrt{5}\right)=\frac{1}{(2i\sqrt{5})(i\sqrt5-i)(i\sqrt5+i)}

?

AlexMu 940
Postad: 18 okt 10:04

Det ser rätt ut!

Tomten 2001
Postad: 18 okt 17:26

Noterar. Din kurva omsluter (för tillräckligt stora R)även z=i. Dess residu måste du också ta med för att få det sökta värdet på integralen.

naytte 7419 – Moderator
Postad: 18 okt 17:52 Redigerad: 18 okt 18:09

Okej, nu har jag kommit fram till följande:

IR=π5-145-t=0t=πiReitR4e4it+6R2e2it+5dt\displaystyle I_R=\frac{\pi\left(\sqrt5-1\right)}{4\sqrt 5}-\int_{t=0}^{t=\pi}\frac{iRe^{it}}{R^4e^{4it}+6R^2e^{2it}+5}dt

Om jag nu vil ha II, då måste jag ju tillämpa ett gränsvärde då RR\to \infty. Då jag gör detta, gäller då följande?

I=π5-145-t=0t=πlimRiReitR4e4it+6R2e2it+5dt\displaystyle I=\frac{\pi\left(\sqrt5-1\right)}{4\sqrt 5}-\int_{t=0}^{t=\pi}\left(\lim_{R \to \infty} \frac{iRe^{it}}{R^4e^{4it}+6R^2e^{2it}+5}\right)dt

Annars är jag lite osäker på hur man ska beräkna integralen. Frågan är alltså om den komplexa integralen har egenskapen att man kan "flytta in" gränsvärdet?

AlexMu 940
Postad: 18 okt 18:21 Redigerad: 18 okt 18:21

Det finns en trevlig sats, ML-inequality, för att beräkna dessa. Den säger att om |f(z)|M|f(z)| \leq M för alla zz på kurvan γ\gamma, så är

γfzdzML\displaystyle \left|\int_{\gamma}f\left(z\right)dz\right|\leq M L, där

LL är längden på kurvan γ\gamma

Jag vet inte om man kan flytta in gränsvärdet. 

naytte 7419 – Moderator
Postad: 18 okt 18:24 Redigerad: 18 okt 18:25

Längden på γ2\gamma_2 är ju L=πRL=\pi R, och då RR\to \infty har vi ju LL\to \infty, så jag är inte riktigt med på hur vi ska använda det.

Däremot kanske man kan triangelolikheta integranden (nämnaren) på något sätt.

AlexMu 940
Postad: 18 okt 18:27

Iden är att hitta ett MM som går snabbare mot 00 än LL mot \infty. Funktionen i integrandet beter sig ungefär som R-3R^{-3}, så om man finner ett MM kommer det säkert ha den storleksordningen. 

naytte 7419 – Moderator
Postad: 18 okt 18:42 Redigerad: 18 okt 22:40

M=R-2M=R^{-2}

För tillräckligt stora RR har vi då att MM är större än integranden. Funkar det? Då beter sig LL ungefär som RR medan MM beter sig som R-2R^{-2} så hela gränsvärdet går mot noll.

AlexMu 940
Postad: 19 okt 01:17

Visst! Vi har exempelvis att 

iReitR4e4it+6R2e2it+5=RR4e4it+6R2e2it+5=R-2R2e2it+6+5R-2e-2it\displaystyle \left|\frac{iRe^{it}}{R^{4}e^{4it}+6R^{2}e^{2it}+5}\right| =\frac{R}{\left|R^{4}e^{4it}+6R^{2}e^{2it}+5\right|}= \frac{R^{-2}}{\left|R^{2}e^{2it}+6+5R^{-2}e^{-2it}\right|}

och nämnaren är rätt uppenbart större än 11 för stora RR, vilket ger att integrandet är mindre än R-2R^{-2}

D4NIEL 3345
Postad: 19 okt 01:28 Redigerad: 19 okt 01:45

Det finns faktiskt en sats som säger att linjeintegralen av f(z)f(z) utmed γ2\gamma_2 som alltså utgörs av halvcirkeln ReiθRe^{i\theta}, 0θπ0\leq \theta\leq \pi, R>R0R>R_0 är

limRγ2fzdz=0\displaystyle \lim_{R\to \infty} \int_{\gamma_2} f\left(z\right)\,\mathrm{d}z=0

om du kan hitta k>1k>1, R0R_0 och λ\lambda så att för alla R>R0R>R_0 och zγ2z\in \gamma_2 gäller

fzλRk\displaystyle \left|f\left(z\right)\right|\leq \frac{\lambda }{R^k}

Detta visar man rättframt med ML-olikheten genom valet M=λRkM=\frac{\lambda}{R^k} och längden av kurvan γ2,L=πR\gamma_2,\, L=\pi R som ni redan påpekat.

En för studenter särskilt användbar variant av satsen med samma bevisidé gäller polynom och ingår förmodligen i er kurs. Den kan lite löst formuleras ungefär så här:

Låt polynomet Q(x)Q(x) vara av grad NN och sakna nollställen utmed x-axeln. Låt vidare polynomet P(x)P(x) vara av högst grad N-2N-2. Då gäller

-P(x)Q(x)dx=2πiresP(z)Q(z),  ö.h.p\displaystyle \int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}\,\mathrm{d}x=2\pi i \sum_{\mathrm{res}}\frac{P(z)}{Q(z)},\quad \mathrm{ö.h.p}

På en tenta kan du i de flesta fall hänvisa till någon av satserna, där den senare är mest bekväm. Det finns också en snarlik variant som kallas "Jordans lemma".

naytte 7419 – Moderator
Postad: 19 okt 13:51 Redigerad: 19 okt 13:52

Spännande! I den sista satsen du diskuterar som gäller polynom, vad innebär detta?

res\displaystyle \sum_{\text{res}}

Res står ju för residy antar jag, men residyerna inom vilket område?

D4NIEL 3345
Postad: 19 okt 14:11

övre halv planet (ö.h.p), dvs summa res för alla poler i övre halvplanet :)

Men leta upp motsvarande sats i er litteratur och använd samma formuleringar, krav och notation som föreläsaren förväntar sig. Det är förmodligen inte tänkt att ni ska behöva föra ett ML-resonemang för varje integral.

naytte 7419 – Moderator
Postad: 19 okt 14:48 Redigerad: 19 okt 14:48

Ska bli! Tack för tipset!

Men så för att försöka fullständiga uppgiften:

Vi har för tillräckligt stora RR att

ft=R|R4e4it+6R2e2it+5|1R3M:=1R2\displaystyle \left|f\left(t\right)\right|=\frac{R}{|R^4e^{4it}+6R^2e^{2it}+5|}\sim\frac{1}{R^3}\le M:=\frac{1}{R^2}

Vi har då alltså att

γ2fzdzML\displaystyle \left|\int_{\gamma_2}f\left(z\right)dz\right|\le ML

för tillräckligt stora RR. Om vi tillämpar ett gränsvärde ser vi att 

limRγ2fzdz=0\displaystyle \lim_{R \to \infty} \left|\int_{\gamma_2}f\left(z\right)dz\right|=0

Om beloppet av integralen går mot noll måste integralen i sig ge noll då RR\to \infty. Så slutligen har vi

I=π5-145\displaystyle I=\frac{\pi\left(\sqrt5-1\right)}{4\sqrt5}

Svara
Close