Impulzní odezva nehomogenního lineárního diferenciálního operátoru
Tento článek pojednává o klasickém přístupu k funkcím Greena. Pro moderní diskusi viz
zásadní řešení .
v matematika , a Greenova funkce je impulsní odezva z nehomogenní lineární operátor diferenciálu definované v doméně se zadanými počátečními podmínkami nebo okrajovými podmínkami.
To znamená, že pokud L je tedy operátor lineárního diferenciálu
Greenova funkce G je řešení rovnice LG = δ , kde δ je Diracova delta funkce ; řešení problému počáteční hodnoty Ly = F je konvoluce (G * F ), kde G je funkce Green. Skrz princip superpozice , vzhledem k tomu, lineární obyčejná diferenciální rovnice (ÓDA), L (solution) = source, one can first resolve L (zelená) = δ s , pro každého s a uvědomit si to, protože zdroj je součet delta funkce , řešením je také součet Greenových funkcí, podle linearity L .
Funkce Greena jsou pojmenovány po Britech matematik George Green , kteří tento koncept poprvé vyvinuli ve 30. letech 19. století. V moderní studii lineárního parciální diferenciální rovnice , Greenovy funkce jsou studovány převážně z hlediska zásadní řešení namísto.
Pod teorie mnoha těl , termín je také používán v fyzika , konkrétně v kvantová teorie pole , aerodynamika , aeroakustika , elektrodynamika , seismologie a statistická teorie pole , odkazovat na různé typy korelační funkce , dokonce i ty, které neodpovídají matematické definici. V kvantové teorii pole mají Greenovy funkce role propagátoři .
Definice a použití Funkce Green, G (x, s ) , a operátor lineárního diferenciálu L = L ( X ) { displaystyle operatorname {L} = operatorname {L} (x)} jednající na distribuce přes podmnožinu Euklidovský prostor R n { displaystyle mathbb {R} ^ {n}} , v určitém okamžiku s , je jakékoli řešení
L G ( X , s ) = δ ( s − X ) , { displaystyle operatorname {L} , G (x, s) = delta (s-x) ,,} (1)
kde δ je Diracova delta funkce . Tuto vlastnost Greenovy funkce lze využít k řešení diferenciálních rovnic formy
L u ( X ) = F ( X ) . { displaystyle operatorname {L} , u (x) = f (x) ~.} (2)
Pokud jádro z L je netriviální, pak funkce Green není jedinečná. V praxi však nějaká kombinace symetrie , okrajové podmínky a / nebo jiná externě stanovená kritéria poskytnou jedinečnou Greenovu funkci. Greenovy funkce lze kategorizovat podle typu splněných okrajových podmínek, podle a Greenovo číslo funkce . Greenovy funkce obecně také jsou distribuce , ne nutně funkce skutečné proměnné.
Greenovy funkce jsou také užitečnými nástroji při řešení vlnové rovnice a difúzní rovnice . v kvantová mechanika , Greenova funkce Hamiltonian je klíčový koncept s důležitými vazbami na koncept hustota stavů .
Greenova funkce používaná ve fyzice je obvykle definována opačným znaménkem. To znamená
L G ( X , s ) = δ ( X − s ) . { displaystyle operatorname {L} , G (x, s) = delta (x-s) ~.} Tato definice významně nemění žádnou z vlastností Greenovy funkce kvůli rovnoměrnosti delta funkce Dirac.
Pokud je operátor překlad neměnný , tedy kdy L { displaystyle operatorname {L}} má konstantní koeficienty s ohledem na X , pak lze Greenovu funkci považovat za a konvoluční jádro , to znamená,
G ( X , s ) = G ( X − s ) . { displaystyle G (x, s) = G (x-s) ~.} V tomto případě je funkce greenu stejná jako impulsní odezva lineární časově invariantní systémová teorie .
Motivace Volně řečeno, pokud taková funkce G lze najít pro operátora L { displaystyle operatorname {L}} , pak, když vynásobíme rovnici (1) pro Greenovu funkci s F (s ) , a poté integrovat s ohledem na s , získáváme,
∫ L G ( X , s ) F ( s ) d s = ∫ δ ( X − s ) F ( s ) d s = F ( X ) . { displaystyle int operatorname {L} , G (x, s) , f (s) , ds = int delta (xs) , f (s) , ds = f (x) ~ .} Protože operátor L = L ( X ) { displaystyle operatorname {L} = operatorname {L} (x)} je lineární a působí pouze na proměnnou X (a ne o proměnné integrace s ), jeden může vzít operátora L { displaystyle operatorname {L}} mimo integraci, výnos
L ( ∫ G ( X , s ) F ( s ) d s ) = F ( X ) . { displaystyle operatorname {L} , left ( int G (x, s) , f (s) , ds right) = f (x) ~.} Tohle znamená tamto
u ( X ) = ∫ G ( X , s ) F ( s ) d s { displaystyle u (x) = int G (x, s) , f (s) , ds} (3)
je řešení rovnice L u ( X ) = F ( X ) . { displaystyle operatorname {L} u (x) = f (x) ~.}
Lze tedy získat funkci u (X ) prostřednictvím znalosti funkce Greena v rovnici (1) a zdrojového výrazu na pravé straně v rovnici (2). Tento proces závisí na linearitě operátora L { displaystyle operatorname {L}} .
Jinými slovy, řešení rovnice (2), u (X ) , lze určit integrací uvedenou v rovnici (3). Ačkoli F (X ) je známo, tuto integraci nelze provést, pokud G je také známo. Problém nyní spočívá v nalezení funkce Greena G který splňuje rovnici (1). Z tohoto důvodu se funkce Zelené také někdy nazývá zásadní řešení přidružené k operátorovi L { displaystyle operatorname {L}} .
Ne každý operátor L { displaystyle operatorname {L}} připouští funkci Greena. Greenovu funkci lze také považovat za a pravý inverzní z L { displaystyle operatorname {L}} . Kromě obtíží při hledání funkce Greena pro konkrétního operátora může být celkem obtížné vyhodnotit integrál v rovnici (3). Metoda však poskytuje teoreticky přesný výsledek.
To lze považovat za expanzi F podle a Diracova delta funkce základ (projektování F přes δ ( X − s ) { displaystyle delta (x-s) ,} ; a superpozice řešení na každém z nich projekce . Taková integrální rovnice je známá jako a Fredholmova integrální rovnice , jehož studium představuje Fredholmova teorie .
Greenovy funkce pro řešení nehomogenních problémů okrajových hodnot Primárním využitím Greenových funkcí v matematice je řešení nehomogenních problémy s hraniční hodnotou . V moderní teoretická fyzika , Greenovy funkce se také obvykle používají jako propagátoři v Feynmanovy diagramy ; termín Greenova funkce se často dále používá pro všechny korelační funkce .
Rámec Nechat L { displaystyle operatorname {L}} být Sturm – Liouville operátor, lineární diferenciální operátor formuláře
L = d d X [ p ( X ) d d X ] + q ( X ) { displaystyle operatorname {L} = { dfrac {d} {dx}} doleva [p (x) { dfrac {d} {dx}} doprava] + q (x)} a nechte D → { displaystyle { vec { operatorname {D}}}} být vektorovou hodnotou okrajové podmínky operátor
D → u = [ α 1 u ′ ( 0 ) + β 1 u ( 0 ) α 2 u ′ ( ℓ ) + β 2 u ( ℓ ) ] . { displaystyle { vec { operatorname {D}}} , u = left [{ begin {matrix} alpha _ {1} u '(0) + beta _ {1} u (0) alpha _ {2} u '( ell) + beta _ {2} u ( ell) end {matrix}} right] ~.} Nechat F ( X ) { displaystyle f (x)} být spojitá funkce v [ 0 , ℓ ] . { displaystyle [0, ell] ~.} Dále předpokládejme, že problém
L u = F D → u = 0 → { displaystyle { begin {zarovnáno} operatorname {L} , u & = f { vec { operatorname {D}}} , u & = { vec {0}} end {zarovnáno}}} je „běžné“, tj. jediné řešení pro F ( X ) = 0 { displaystyle f (x) = 0} pro všechny X je u ( X ) = 0 { displaystyle u (x) = 0} .[A]
Teorém Existuje jedno a jediné řešení u ( X ) { displaystyle u (x)} to uspokojuje
L u = F D → u = 0 → { displaystyle { begin {zarovnáno} operatorname {L} , u & = f { vec { operatorname {D}}} , u & = { vec {0}} end {zarovnáno}}} a je to dáno
u ( X ) = ∫ 0 ℓ F ( s ) G ( X , s ) d s , { displaystyle u (x) = int _ {0} ^ { ell} f (s) , G (x, s) , ds ~,} kde G ( X , s ) { displaystyle G (x, s)} je funkce Greena splňující následující podmínky:
G ( X , s ) { displaystyle G (x, s)} je nepřetržitý v X { displaystyle x} a s { displaystyle s} .Pro X ≠ s { displaystyle x neq s ~} , L G ( X , s ) = 0 { displaystyle quad operatorname {L} , G (x, s) = 0 ~} . Pro s ≠ 0 { displaystyle s neq 0 ~} , D → G ( X , s ) = 0 → { displaystyle quad { vec { operatorname {D}}} , G (x, s) = { vec {0}} ~} . Derivát "skok": G ′ ( s 0 + , s ) − G ′ ( s 0 − , s ) = 1 / p ( s ) { displaystyle quad G '(s_ {0 +}, s) -G' (s_ {0 -}, s) = 1 / p (s) ~} .Symetrie: G ( X , s ) = G ( s , X ) { displaystyle quad G (x, s) = G (s, x) ~} . Pokročilé a retardované Greenovy funkce Někdy lze Greenovu funkci rozdělit na součet dvou funkcí. Jeden s proměnnou pozitivní (+) a druhý s proměnnou negativní (-). Jedná se o pokročilé a retardované Greenovy funkce, a když zkoumaná rovnice závisí na čase, je to jedna z částí kauzální a další anti-kauzální. U těchto problémů je obvykle důležitá kauzální část. Toto jsou často řešení nehomogenní elektromagnetická vlnová rovnice .
Hledání funkcí Greena Jednotky I když to jednoznačně neopravuje formu, jakou bude mít Greenova funkce, provedením a rozměrová analýza najít jednotky, které Greenova funkce musí mít, je důležitá kontrola zdravého rozumu u jakékoli Greenovy funkce nalezené jinými prostředky. Rychlé prozkoumání definující rovnice,
L G ( X , s ) = δ ( X − s ) , { displaystyle LG (x, s) = delta (x-s),} ukazuje, že jednotky G { displaystyle G} závisí nejen na jednotkách L { displaystyle L} ale také na počtu a jednotkách prostoru, ve kterém jsou polohové vektory X { displaystyle x} a s { displaystyle s} jsou prvky. To vede ke vztahu:
[ [ G ] ] = [ [ L ] ] − 1 [ [ d X ] ] − 1 , { displaystyle [[G]] = [[L]] ^ {- 1} [[ mathrm {d} x]] ^ {- 1},} kde [ [ G ] ] { displaystyle [[G]]} je definován jako, "fyzické jednotky G { displaystyle G} ", a d X { displaystyle mathrm {d} x} je objemový prvek prostoru (nebo vesmírný čas ).
Například pokud L = ∂ t 2 { displaystyle L = částečné _ {t} ^ {2}} a čas je potom jedinou proměnnou:
[ [ L ] ] = [ [ t i m E ] ] − 2 , { displaystyle [[L]] = [[ mathrm {time}]] ^ {- 2},} [ [ d X ] ] = [ [ t i m E ] ] , A n d { displaystyle [[ mathrm {d} x]] = [[ mathrm {čas}]], mathrm {a}} [ [ G ] ] = [ [ t i m E ] ] . { displaystyle [[G]] = [[ mathrm {time}]].} Li L = ◻ = 1 C 2 ∂ t 2 − ∇ 2 { displaystyle L = square = { frac {1} {c ^ {2}}} částečný _ {t} ^ {2} - nabla ^ {2}} , operátor d'Alembert , a prostor má poté 3 rozměry:
[ [ L ] ] = [ [ l E n G t h ] ] − 2 , { displaystyle [[L]] = [[ mathrm {délka}]] ^ {- 2},} [ [ d X ] ] = [ [ t i m E ] ] [ [ l E n G t h ] ] 3 , A n d { displaystyle [[ mathrm {d} x]] = [[ mathrm {čas}]] [[ mathrm {délka}]] ^ {3}, mathrm {a}} [ [ G ] ] = [ [ t i m E ] ] − 1 [ [ l E n G t h ] ] − 1 . { displaystyle [[G]] = [[ mathrm {time}]] ^ {- 1} [[ mathrm {length}]] ^ {- 1}.} Expanze vlastních čísel Pokud operátor diferenciálu L připouští soubor vlastní vektory Ψn (X ) (tj. sada funkcí Ψn a skaláry λ n takhle L Ψn = λ n Ψn ), který je kompletní, pak je možné z těchto vlastních vektorů a vlastní čísla .
„Kompletní“ znamená, že sada funkcí { Ψn } splňuje následující vztah úplnosti ,
δ ( X − X ′ ) = ∑ n = 0 ∞ Ψ n † ( X ) Ψ n ( X ′ ) . { displaystyle delta (x-x ') = součet _ {n = 0} ^ { infty} Psi _ {n} ^ { dagger} (x) Psi _ {n} (x'). } Pak platí
G ( X , X ′ ) = ∑ n = 0 ∞ Ψ n † ( X ) Ψ n ( X ′ ) λ n , { displaystyle G (x, x ') = součet _ {n = 0} ^ { infty} { dfrac { Psi _ {n} ^ { dagger} (x) Psi _ {n} (x ')} { lambda _ {n}}},}
kde † { displaystyle dýka} představuje komplexní konjugaci.
Použití operátora L na každou stranu této rovnice vznikne vztah úplnosti, který se předpokládal.
Obecná studie o funkci Green napsaná ve výše uvedené formě a její vztah k funkční prostory tvořený vlastními vektory, je známý jako Fredholmova teorie .
Existuje několik dalších metod pro nalezení funkcí Greena, včetně metoda obrázků , oddělení proměnných , a Laplaceovy transformace (Cole 2011).
Kombinace Greenových funkcí Pokud je operátor diferenciálu L { displaystyle L} lze započítat jako L = L 1 L 2 { displaystyle L = L_ {1} L_ {2}} pak Greenova funkce L { displaystyle L} lze sestavit z funkcí Green pro L 1 { displaystyle L_ {1}} a L 2 { displaystyle L_ {2}} :
G ( X , s ) = ∫ G 2 ( X , s 1 ) G 1 ( s 1 , s ) d s 1 . { displaystyle G (x, s) = int G_ {2} (x, s_ {1}) , G_ {1} (s_ {1}, s) , mathrm {d} s_ {1}. } Výše uvedená totožnost vyplývá okamžitě z převzetí G ( X , s ) { displaystyle G (x, s)} být představou inverze správného operátoru L { displaystyle L} , obdobně jako pro invertibilní lineární operátor C { displaystyle C} , definován C = ( A B ) − 1 = B − 1 A − 1 { displaystyle C = (AB) ^ {- 1} = B ^ {- 1} A ^ {- 1}} , je reprezentován svými maticovými prvky C i , j { displaystyle C_ {i, j}} .
Následuje další identita pro diferenciální operátory, které jsou skalárními polynomy derivace, L = P N ( ∂ X ) { displaystyle L = P_ {N} ( částečné _ {x})} . The základní věta o algebře , v kombinaci s tím, že ∂ X { displaystyle částečné _ {x}} dojíždí sama se sebou , zaručuje, že polynom lze zapracovat L { displaystyle L} ve formě:
L = ∏ i = 1 N ( ∂ X − z i ) , { displaystyle L = prod _ {i = 1} ^ {N} ( částečné _ {x} -z_ {i}),} kde z i { displaystyle z_ {i}} jsou nuly P N ( z ) { displaystyle P_ {N} (z)} . Užívání Fourierova transformace z L G ( X , s ) = δ ( X − s ) { displaystyle LG (x, s) = delta (x-s)} s ohledem na oba X { displaystyle x} a s { displaystyle s} dává:
G ^ ( k X , k s ) = δ ( k X − k s ) ∏ i = 1 N ( i k X − z i ) . { displaystyle { widehat {G}} (k_ {x}, k_ {s}) = { frac { delta (k_ {x} -k_ {s})} { prod _ {i = 1} ^ {N} (ik_ {x} -z_ {i})}}.} Zlomek lze poté rozdělit na součet pomocí a Rozklad částečné frakce před Fourierovou transformací zpět na X { displaystyle x} a s { displaystyle s} prostor. Tento proces poskytuje identity, které souvisejí s integrály Greenových funkcí a jejich součty. Například pokud L = ( ∂ X + y ) ( ∂ X + α ) 2 { displaystyle L = ( částečné _ {x} + gamma) ( částečné _ {x} + alfa) ^ {2}} pak jeden formulář pro jeho funkci Green je:
G ( X , s ) = 1 ( α − y ) 2 Θ ( X − s ) E − y ( X − s ) − 1 ( α − y ) 2 Θ ( X − s ) E − α ( X − s ) + 1 y − α Θ ( X − s ) ( X − s ) E − α ( X − s ) = ∫ Θ ( X − s 1 ) ( X − s 1 ) E − α ( X − s 1 ) Θ ( s 1 − s ) E − y ( s 1 − s ) d s 1 . { displaystyle { begin {zarovnáno} G (x, s) & = { frac {1} {( alpha - gama) ^ {2}}} Theta (xs) operatorname {e} ^ {- gamma (xs)} - { frac {1} {( alpha - gamma) ^ {2}}} Theta (xs) operatorname {e} ^ {- alpha (xs)} + { frac {1} { gamma - alpha}} Theta (xs) , (xs) operatorname {e} ^ {- alpha (xs)} [5pt] & = int Theta (x-s_ {1}) (x-s_ {1}) operatorname {e} ^ {- alpha (x-s_ {1})} Theta (s_ {1} -s) operatorname {e} ^ {- gamma (s_ {1} -s)} , mathrm {d} s_ {1}. end {zarovnáno}}} I když je uvedený příklad analyticky použitelný, ilustruje proces, který funguje, když integrál není triviální (například když ∇ 2 { displaystyle nabla ^ {2}} je operátor v polynomu).
Tabulka Greenových funkcí Následující tabulka poskytuje přehled funkcí Greenových často se objevujících diferenciálních operátorů, kde r = X 2 + y 2 + z 2 { displaystyle textstyle r = { sqrt {x ^ {2} + y ^ {2} + z ^ {2}}}} , ρ = X 2 + y 2 { displaystyle textstyle rho = { sqrt {x ^ {2} + y ^ {2}}}} , Θ ( t ) { displaystyle textstyle Theta (t)} je Funkce Heaviside step , J ν ( z ) { displaystyle textstyle J _ { nu} (z)} je Besselova funkce , Já ν ( z ) { displaystyle textstyle I _ { nu} (z)} je upravená Besselova funkce prvního druhu , a K. ν ( z ) { displaystyle textstyle K _ { nu} (z)} je upravená Besselova funkce druhého druhu .[1] Kde čas (t ) se objeví v prvním sloupci, je uvedena pokročilá (kauzální) Greenova funkce.
Diferenciální operátor L Greenova funkce G Příklad aplikace ∂ t n + 1 { displaystyle částečné _ {t} ^ {n + 1}} t n n ! Θ ( t ) { displaystyle { frac {t ^ {n}} {n!}} Theta (t)} ∂ t + y { displaystyle částečné _ {t} + gamma} Θ ( t ) E − y t { displaystyle Theta (t) mathrm {e} ^ {- gama t}} ( ∂ t + y ) 2 { displaystyle left ( partial _ {t} + gamma right) ^ {2}} Θ ( t ) t E − y t { displaystyle Theta (t) t mathrm {e} ^ {- gama t}} ∂ t 2 + 2 y ∂ t + ω 0 2 { displaystyle částečné _ {t} ^ {2} +2 gamma částečné _ {t} + omega _ {0} ^ {2}} Θ ( t ) E − y t hřích ( ω t ) ω { displaystyle Theta (t) mathrm {e} ^ {- gamma t} ~ { frac { sin ( omega t)} { omega}}} s ω = ω 0 2 − y 2 { displaystyle omega = { sqrt { omega _ {0} ^ {2} - gamma ^ {2}}}} 1D tlumený harmonický oscilátor 2D Laplaceův operátor Δ 2D = ∂ X 2 + ∂ y 2 { displaystyle Delta _ { text {2D}} = částečné _ {x} ^ {2} + částečné _ {y} ^ {2}} 1 2 π ln ρ { displaystyle { frac {1} {2 pi}} ln rho} s ρ = X 2 + y 2 { displaystyle rho = { sqrt {x ^ {2} + y ^ {2}}}} 2D Poissonova rovnice Operátor 3D Laplace Δ 3D = ∂ X 2 + ∂ y 2 + ∂ z 2 { displaystyle Delta _ { text {3D}} = částečné _ {x} ^ {2} + částečné _ {y} ^ {2} + částečné _ {z} ^ {2}} − 1 4 π r { displaystyle { frac {-1} {4 pi r}}} s r = X 2 + y 2 + z 2 { displaystyle r = { sqrt {x ^ {2} + y ^ {2} + z ^ {2}}}} Poissonova rovnice Provozovatel Helmholtz Δ 3D + k 2 { displaystyle Delta _ { text {3D}} + k ^ {2}} − E − i k r 4 π r = i k 32 π r { displaystyle { frac {- mathrm {e} ^ {- ikr}} {4 pi r}} = i { sqrt { frac {k} {32 pi r}}}} H 1 / 2 ( 2 ) ( k r ) { displaystyle H_ {1/2} ^ {(2)} (kr)} = i k 4 π { displaystyle = i { frac {k} {4 pi}} ,} h 0 ( 2 ) ( k r ) { displaystyle h_ {0} ^ {(2)} (kr)} stacionární 3D Schrödingerova rovnice pro volná částice Δ − k 2 { displaystyle Delta -k ^ {2}} v n { displaystyle n} rozměry − ( 2 π ) − n / 2 ( k r ) n / 2 − 1 K. n / 2 − 1 ( k r ) { displaystyle - (2 pi) ^ {- n / 2} left ({ frac {k} {r}} right) ^ {n / 2-1} K_ {n / 2-1} (kr )} Yukawa potenciál , Feynmanův propagátor ∂ t 2 − C 2 ∂ X 2 { displaystyle částečné _ {t} ^ {2} -c ^ {2} částečné _ {x} ^ {2}} 1 2 C Θ ( t − | X / C | ) { displaystyle { frac {1} {2c}} Theta (t- | x / c |)} 1D vlnová rovnice ∂ t 2 − C 2 Δ 2D { displaystyle částečný _ {t} ^ {2} -c ^ {2} , Delta _ { text {2D}}} 1 2 π C C 2 t 2 − ρ 2 Θ ( t − ρ / C ) { displaystyle { frac {1} {2 pi c { sqrt {c ^ {2} t ^ {2} - rho ^ {2}}}}} Theta (t- rho / c)} 2D vlnová rovnice D'Alembertův operátor ◻ = 1 C 2 ∂ t 2 − Δ 3D { displaystyle square = { frac {1} {c ^ {2}}} částečný _ {t} ^ {2} - Delta _ { text {3D}}} δ ( t − r C ) 4 π r { displaystyle { frac { delta (t - { frac {r} {c}})} {4 pi r}}} 3D vlnová rovnice ∂ t − k ∂ X 2 { displaystyle částečné _ {t} -k částečné _ {x} ^ {2}} Θ ( t ) ( 1 4 π k t ) 1 / 2 E − X 2 / 4 k t { displaystyle Theta (t) vlevo ({ frac {1} {4 pi kt}} vpravo) ^ {1/2} mathrm {e} ^ {- x ^ {2} / 4kt}} 1D difúze ∂ t − k Δ 2D { displaystyle částečné _ {t} -k , Delta _ { text {2D}}} Θ ( t ) ( 1 4 π k t ) E − ρ 2 / 4 k t { displaystyle Theta (t) vlevo ({ frac {1} {4 pi kt}} vpravo) mathrm {e} ^ {- rho ^ {2} / 4kt}} 2D difúze ∂ t − k Δ 3D { displaystyle částečné _ {t} -k , Delta _ { text {3D}}} Θ ( t ) ( 1 4 π k t ) 3 / 2 E − r 2 / 4 k t { displaystyle Theta (t) vlevo ({ frac {1} {4 pi kt}} vpravo) ^ {3/2} mathrm {e} ^ {- r ^ {2} / 4kt}} 3D difúze 1 C 2 ∂ t 2 − ∂ X 2 + μ 2 { displaystyle { frac {1} {c ^ {2}}} částečné _ {t} ^ {2} - částečné _ {x} ^ {2} + mu ^ {2}} 1 2 [ ( 1 − hřích μ C t ) ( δ ( C t − X ) + δ ( C t + X ) ) + μ Θ ( C t − | X | ) J 0 ( μ u ) ] , u = C 2 t 2 − X 2 { displaystyle { frac {1} {2}} doleva [ doleva (1- sin { mu ct} doprava) ( delta (ct-x) + delta (ct + x)) + mu Theta (ct- | x |) J_ {0} ( mu u) right], , u = { sqrt {c ^ {2} t ^ {2} -x ^ {2}}}} 1D Klein-Gordonova rovnice 1 C 2 ∂ t 2 − Δ 2D + μ 2 { displaystyle { frac {1} {c ^ {2}}} částečný _ {t} ^ {2} - Delta _ { text {2D}} + mu ^ {2}} 1 4 π [ ( 1 + cos ( μ C t ) ) δ ( C t − ρ ) ρ + μ 2 Θ ( C t − ρ ) upřímně ( μ u ) ] , u = C 2 t 2 − ρ 2 { displaystyle { frac {1} {4 pi}} vlevo [(1+ cos ( mu ct)) { frac { delta (ct- rho)} { rho}} + mu ^ {2} Theta (ct- rho) operatorname {sinc} ( mu u) right], , u = { sqrt {c ^ {2} t ^ {2} - rho ^ {2 }}}} 2D Klein-Gordonova rovnice ◻ + μ 2 { displaystyle square + mu ^ {2}} 1 4 π [ δ ( t − r C ) r + μ C Θ ( C t − r ) J 1 ( μ u ) u ] , u = C 2 t 2 − r 2 { displaystyle { frac {1} {4 pi}} vlevo [{ frac { delta vlevo (t - { frac {r} {c}} vpravo)} {r}} + mu c Theta (ct-r) { frac {J_ {1} left ( mu u right)} {u}} right], , u = { sqrt {c ^ {2} t ^ { 2} -r ^ {2}}}} 3D Klein-Gordonova rovnice ∂ t 2 + 2 y ∂ t − C 2 ∂ X 2 { displaystyle částečné _ {t} ^ {2} +2 gamma částečné _ {t} -c ^ {2} částečné _ {x} ^ {2}} 1 2 E − y t [ δ ( C t − X ) + δ ( C t + X ) + Θ ( C t − | X | ) ( y C Já 0 ( y u C ) + y t u Já 1 ( y u C ) ) ] , u = C 2 t 2 − X 2 { displaystyle { frac {1} {2}} e ^ {- gamma t} left [ delta (ct-x) + delta (ct + x) + theta (ct- | x |) left ({ frac { gamma} {c}} I_ {0} left ({ frac { gamma u} {c}} right) + { frac { gamma t} {u}} I_ { 1} left ({ frac { gamma u} {c}} right) right) right], , u = { sqrt {c ^ {2} t ^ {2} -x ^ {2 }}}} telegrafní rovnice ∂ t 2 + 2 y ∂ t − C 2 Δ 2D { displaystyle částečné _ {t} ^ {2} +2 gamma částečné _ {t} -c ^ {2} , Delta _ { text {2D}}} E − y t 4 π [ ( 1 + E − y t + 3 y t ) δ ( C t − ρ ) ρ + Θ ( C t − ρ ) ( y sinh ( y u C ) C u + 3 y t hovadina ( y u C ) u 2 − 3 C t sinh ( y u C ) u 3 ) ] , u = C 2 t 2 − ρ 2 { displaystyle { frac {e ^ {- gamma t}} {4 pi}} left [(1 + e ^ {- gamma t} +3 gamma t) { frac { delta (ct - rho)} { rho}} + Theta (ct- rho) left ({ frac { gamma sinh left ({ frac { gamma u} {c}} right)} { cu}} + { frac {3 gamma t cosh left ({ frac { gamma u} {c}} right)} {u ^ {2}}} - { frac {3ct sinh left ({ frac { gamma u} {c}} right)} {u ^ {3}}} right) right], , u = { sqrt {c ^ {2} t ^ {2 } - rho ^ {2}}}} 2D relativistické vedení tepla ∂ t 2 + 2 y ∂ t − C 2 Δ 3D { displaystyle částečné _ {t} ^ {2} +2 gamma částečné _ {t} -c ^ {2} , Delta _ { text {3D}}} E − y t 20 π [ ( 8 − 3 E − y t + 2 y t + 4 y 2 t 2 ) δ ( C t − r ) r 2 + y 2 C Θ ( C t − r ) ( 1 C u Já 1 ( y u C ) + 4 t u 2 Já 2 ( y u C ) ) ] , u = C 2 t 2 − r 2 { displaystyle { frac {e ^ {- gamma t}} {20 pi}} left [ left (8-3e ^ {- gamma t} +2 gamma t + 4 gamma ^ {2 } t ^ {2} right) { frac { delta (ct-r)} {r ^ {2}}} + { frac { gamma ^ {2}} {c}} Theta (ct- r) left ({ frac {1} {cu}} I_ {1} left ({ frac { gamma u} {c}} right) + { frac {4t} {u ^ {2} }} I_ {2} left ({ frac { gamma u} {c}} right) right) right], , u = { sqrt {c ^ {2} t ^ {2} - r ^ {2}}}} 3D relativistické vedení tepla
Greenovy funkce pro Laplacian Greenovy funkce pro lineární diferenciální operátory zahrnující Laplacian lze snadno použít pomocí druhého z Greenovy identity .
Chcete-li odvodit Greenovu větu, začněte s věta o divergenci (jinak známé jako Gaussova věta ),
∫ PROTI ∇ ⋅ A → d PROTI = ∫ S A → ⋅ d σ ^ . { displaystyle int _ {V} nabla cdot { vec {A}} dV = int _ {S} { vec {A}} cdot d { widehat { sigma}} ~.} Nechat A → = φ ∇ ψ − ψ ∇ φ { displaystyle { vec {A}} = varphi , nabla psi - psi , nabla varphi} a nahradit Gaussovým zákonem.
Vypočítat ∇ ⋅ A → { displaystyle nabla cdot { vec {A}}} a použít pravidlo produktu pro operátora,,
∇ ⋅ A → = ∇ ⋅ ( φ ∇ ψ − ψ ∇ φ ) = ( ∇ φ ) ⋅ ( ∇ ψ ) + φ ∇ 2 ψ − ( ∇ φ ) ⋅ ( ∇ ψ ) − ψ ∇ 2 φ = φ ∇ 2 ψ − ψ ∇ 2 φ . { displaystyle { begin {zarovnáno} nabla cdot { vec {A}} & = nabla cdot ( varphi , nabla psi ; - ; psi , nabla varphi) & = ( nabla varphi) cdot ( nabla psi) ; + ; varphi , nabla ^ {2} psi ; - ; ( nabla varphi) cdot ( nabla psi) ; - ; psi nabla ^ {2} varphi & = varphi , nabla ^ {2} psi ; - ; psi , nabla ^ {2} varphi. end {zarovnáno}}} Zapojením do věty o divergenci vznikne Greenova věta ,
∫ PROTI ( φ ∇ 2 ψ − ψ ∇ 2 φ ) d PROTI = ∫ S ( φ ∇ ψ − ψ ∇ φ ) ⋅ d σ ^ . { displaystyle int _ {V} ( varphi , nabla ^ {2} psi - psi , nabla ^ {2} varphi) , dV = int _ {S} ( varphi , nabla psi - psi nabla , varphi) cdot d { widehat { sigma}}.} Předpokládejme, že operátor lineárního diferenciálu L je Laplacian , ∇², a že existuje Greenova funkce G pro Laplaciany. Definující vlastnost funkce Greena stále platí,
L G ( X , X ′ ) = ∇ 2 G ( X , X ′ ) = δ ( X − X ′ ) . { displaystyle LG (x, x ') = nabla ^ {2} G (x, x') = delta (x-x ').} Nechat ψ = G { displaystyle psi = G} v Greenově druhé identitě, viz Greenovy identity . Pak,
∫ PROTI [ φ ( X ′ ) δ ( X − X ′ ) − G ( X , X ′ ) ∇ ′ 2 φ ( X ′ ) ] d 3 X ′ = ∫ S [ φ ( X ′ ) ∇ ′ G ( X , X ′ ) − G ( X , X ′ ) ∇ ′ φ ( X ′ ) ] ⋅ d σ ^ ′ . { displaystyle int _ {V} vlevo [ varphi (x ') delta (x-x') - G (x, x ') , { nabla'} ^ {2} , varphi ( x ') right] d ^ {3} x' = int _ {S} left [ varphi (x ') , { nabla'} G (x, x ') - G (x, x ') , { nabla'} varphi (x ') right] cdot d { widehat { sigma}}'.} Pomocí tohoto výrazu je možné vyřešit Laplaceova rovnice ∇2 φ (X ) = 0 nebo Poissonova rovnice ∇2 φ (X ) = −ρ (X ), s výhradou buď Neumann nebo Dirichlet okrajové podmínky. Jinými slovy, můžeme vyřešit pro φ (X ) všude uvnitř svazku, kde (1) je hodnota φ (X ) je uveden na ohraničujícím povrchu objemu (Dirichletovy okrajové podmínky), nebo (2) normální derivace φ (X ) je uveden na ohraničující ploše (okrajové podmínky Neumanna).
Předpokládejme, že problém je vyřešit φ (X ) uvnitř regionu. Pak integrál
∫ PROTI φ ( X ′ ) δ ( X − X ′ ) d 3 X ′ { displaystyle int limity _ {V} varphi (x ') delta (x-x') , d ^ {3} x '} redukuje na jednoduše φ (X ) vzhledem k určující vlastnosti Diracova delta funkce a máme
φ ( X ) = − ∫ PROTI G ( X , X ′ ) ρ ( X ′ ) d 3 X ′ + ∫ S [ φ ( X ′ ) ∇ ′ G ( X , X ′ ) − G ( X , X ′ ) ∇ ′ φ ( X ′ ) ] ⋅ d σ ^ ′ . { displaystyle varphi (x) = - int _ {V} G (x, x ') rho (x') d ^ {3} x '+ int _ {S} doleva [ varphi ( x ') , nabla' G (x, x ') - G (x, x') , nabla ' varphi (x') doprava] cdot d { widehat { sigma}} '. } Tato forma vyjadřuje známou vlastnost harmonické funkce , že pokud je hodnota nebo normální derivace známa na ohraničující ploše, pak je hodnota funkce uvnitř objemu známá všude .
v elektrostatika , φ (X ) je interpretován jako elektrický potenciál , ρ (X ) tak jako elektrický náboj hustota a normální derivát ∇ φ ( X ′ ) ⋅ d σ ^ ′ { displaystyle nabla varphi (x ') cdot d { widehat { sigma}}'}} jako normální složka elektrického pole.
Pokud je problém vyřešit problém s hraniční hodnotou Dirichlet, měla by být zvolena Greenova funkce tak, že G (X ,X ′) Zmizí, když buď X nebo X ′ Je na ohraničující ploše. Tedy pouze jeden ze dvou termínů v EU povrchový integrál Zůstává. Pokud je problémem vyřešit problém s hraniční hodnotou Neumanna, zvolí se Greenova funkce tak, že její normální derivace zmizí na ohraničující ploše, což se zdá být nejlogičtější volbou. (Viz Jackson J.D. klasická elektrodynamika, strana 39). Aplikace Gaussovy věty na diferenciální rovnici definující Greenovu funkci se však získá
∫ S ∇ ′ G ( X , X ′ ) ⋅ d σ ^ ′ = ∫ PROTI ∇ ′ 2 G ( X , X ′ ) d 3 X ′ = ∫ PROTI δ ( X − X ′ ) d 3 X ′ = 1 , { displaystyle int _ {S} nabla 'G (x, x') cdot d { widehat { sigma}} '= int _ {V} nabla' ^ {2} G (x, x ') d ^ {3} x' = int _ {V} delta (x-x ') d ^ {3} x' = 1 ~,} což znamená normální derivaci G (X ,X ′) Nemůže zmizet na povrchu, protože se musí integrovat do 1 na povrchu. (Opět viz Jackson J.D. klasická elektrodynamika, strana 39 pro tento a následující argument).
Nejjednodušší forma, kterou může mít normální derivace, je konstanta, konkrétně 1 /S , kde S je povrchová plocha povrchu. Povrchový člen v řešení se stává
∫ S φ ( X ′ ) ∇ ′ G ( X , X ′ ) ⋅ d σ ^ ′ = ⟨ φ ⟩ S { displaystyle int _ {S} varphi (x ') , nabla' G (x, x ') cdot d { widehat { sigma}}' = langle varphi rangle _ {S} } kde ⟨ φ ⟩ S { displaystyle langle varphi rangle _ {S}} je průměrná hodnota potenciálu na povrchu. Toto číslo není obecně známé, ale je často nedůležité, protože cílem je často získat elektrické pole dané spíše gradientem potenciálu než samotným potenciálem.
Bez okrajových podmínek je Greenova funkce pro Laplacian (Greenova funkce pro Laplacovu rovnici se třemi proměnnými ) je
G ( X , X ′ ) = − 1 4 π | X − X ′ | . { displaystyle G (x, x ') = - { dfrac {1} {4 pi | x-x' |}}.} Předpokládejme, že ohraničující povrch zhasne do nekonečna a zapojením tohoto výrazu pro funkci Green nakonec získá standardní výraz pro elektrický potenciál, pokud jde o hustotu elektrického náboje jako
φ ( X ) = ∫ PROTI ρ ( X ′ ) 4 π ε | X − X ′ | d 3 X ′ . { displaystyle varphi (x) = int _ {V} { dfrac { rho (x ')} {4 pi varepsilon | x-x' |}} , d ^ {3} x '~ .}
Příklad Příklad. Najděte funkci Zelená pro následující problém, jehož Greenovo číslo funkce je X11:
L u = u ″ + k 2 u = F ( X ) u ( 0 ) = 0 , u ( π 2 k ) = 0. { displaystyle { begin {zarovnáno} Lu & = u '' + k ^ {2} u = f (x) u (0) & = 0, quad u left ({ tfrac { pi} { 2k}} vpravo) = 0. End {zarovnáno}}} První krok: Greenova funkce pro lineárního operátora je definována jako řešení
G ″ ( X , s ) + k 2 G ( X , s ) = δ ( X − s ) . { displaystyle g '' (x, s) + k ^ {2} g (x, s) = delta (x-s).} Li X ≠ s { displaystyle x neq s} , pak funkce delta dává nulu a obecné řešení je
G ( X , s ) = C 1 cos k X + C 2 hřích k X . { displaystyle g (x, s) = c_ {1} cos kx + c_ {2} sin kx.} Pro X < s { displaystyle x , okrajová podmínka v X = 0 { displaystyle x = 0} naznačuje
G ( 0 , s ) = C 1 ⋅ 1 + C 2 ⋅ 0 = 0 , C 1 = 0 { displaystyle g (0, s) = c_ {1} cdot 1 + c_ {2} cdot 0 = 0, quad c_ {1} = 0} -li X < s { displaystyle x a s ≠ π 2 k { displaystyle s neq { tfrac { pi} {2k}}} .
Pro X > s { displaystyle x> s} , okrajová podmínka v X = π 2 k { displaystyle x = { tfrac { pi} {2k}}} naznačuje
G ( π 2 k , s ) = C 3 ⋅ 0 + C 4 ⋅ 1 = 0 , C 4 = 0 { displaystyle g left ({ tfrac { pi} {2k}}, s right) = c_ {3} cdot 0 + c_ {4} cdot 1 = 0, quad c_ {4} = 0 } Rovnice G ( 0 , s ) = 0 { displaystyle g (0, s) = 0} je z podobných důvodů přeskočen.
Souhrn dosavadních výsledků:
G ( X , s ) = { C 2 hřích k X , pro X < s , C 3 cos k X , pro s < X . { displaystyle g (x, s) = { begin {cases} c_ {2} sin kx, & { text {for}} x Druhý krok: Dalším úkolem je určit C 2 { displaystyle c_ {2}} a C 3 { displaystyle c_ {3}} .
Zajištění kontinuity funkce Green v X = s { displaystyle x = s} naznačuje
C 2 hřích k s = C 3 cos k s { displaystyle c_ {2} sin ks = c_ {3} cos ks} Lze zajistit správnou diskontinuitu v první derivaci integrací definující diferenciální rovnice z X = s − ε { displaystyle x = s- varepsilon} na X = s + ε { displaystyle x = s + varepsilon} a brát limit jako ε { displaystyle varepsilon} jde na nulu:
C 3 ⋅ ( − k hřích k s ) − C 2 ⋅ ( k cos k s ) = 1 { displaystyle c_ {3} cdot (-k sin ks) -c_ {2} cdot (k cos ks) = 1} Lze vyřešit dvě rovnice (ne) spojitosti C 2 { displaystyle c_ {2}} a C 3 { displaystyle c_ {3}} získat
C 2 = − cos k s k ; C 3 = − hřích k s k { displaystyle c_ {2} = - { frac { cos ks} {k}} quad; quad c_ {3} = - { frac { sin ks} {k}}} Greenova funkce pro tento problém tedy je:
G ( X , s ) = { − cos k s k hřích k X , X < s , − hřích k s k cos k X , s < X . { displaystyle g (x, s) = { begin {cases} - { frac { cos ks} {k}} sin kx, & x Další příklady G ( X , y , X 0 , y 0 ) = 1 2 π [ ln ( X − X 0 ) 2 + ( y − y 0 ) 2 − ln ( X + X 0 ) 2 + ( y − y 0 ) 2 + ln ( X − X 0 ) 2 + ( y + y 0 ) 2 − ln ( X + X 0 ) 2 + ( y + y 0 ) 2 ] . { displaystyle { begin {aligned} G (x, y, x_ {0}, y_ {0}) = { dfrac {1} {2 pi}} & left [ ln { sqrt {(x -x_ {0}) ^ {2} + (y-y_ {0}) ^ {2}}} - ln { sqrt {(x + x_ {0}) ^ {2} + (y-y_ { 0}) ^ {2}}} vpravo. [5 bodů] & vlevo. {} + Ln { sqrt {(x-x_ {0}) ^ {2} + (y + y_ {0} ) ^ {2}}} - ln { sqrt {(x + x_ {0}) ^ {2} + (y + y_ {0}) ^ {2}}} , vpravo]. End { zarovnaný}}} Nechat A < X < b { displaystyle a a všechny tři jsou prvky reálných čísel. Pak pro jakoukoli funkci od reálných k reálným, F ( X ) { displaystyle f (x)} , s n { displaystyle n} th derivace, která je integrovatelná v daném intervalu [ A , b ] { displaystyle [a, b]} : F ( X ) = ∑ m = 0 n − 1 ( X − A ) m m ! [ d m F d X m ] X = A + ∫ A b [ ( X − s ) n − 1 ( n − 1 ) ! Θ ( X − s ) ] [ d n F d X n ] X = s d s . { displaystyle { begin {aligned} f (x) & = sum _ {m = 0} ^ {n-1} { frac {(xa) ^ {m}} {m!}} left [{ frac { mathrm {d} ^ {m} f} { mathrm {d} x ^ {m}}} vpravo] _ {x = a} + int _ {a} ^ {b} vlevo [ { frac {(xs) ^ {n-1}} {(n-1)!}} Theta (xs) right] left [{ frac { mathrm {d} ^ {n} f} { mathrm {d} x ^ {n}}} right] _ {x = s} mathrm {d} s end {zarovnáno}} ~.} Greenova funkce ve výše uvedené rovnici, G ( X , s ) = ( X − s ) n − 1 ( n − 1 ) ! Θ ( X − s ) { displaystyle G (x, s) = { frac {(x-s) ^ {n-1}} {(n-1)!}} Theta (x-s)} , není jedinečný. Jak je upravena rovnice, pokud G ( X − s ) { displaystyle g (x-s)} je přidán do G ( X , s ) { displaystyle G (x, s)} , kde G ( X ) { displaystyle g (x)} splňuje d n G d X n = 0 { displaystyle { frac { mathrm {d} ^ {n} g} { mathrm {d} x ^ {n}}} = 0} pro všechny X ∈ [ A , b ] { displaystyle x v [a, b]} (například, G ( X ) = − X / 2 { displaystyle g (x) = - x / 2} s n = 2 { displaystyle n = 2} )? Rovněž porovnejte výše uvedenou rovnici s formou a Taylor série se středem na X = A { displaystyle x = a} . Viz také ^ V odborném žargonu „běžný“ znamená, že pouze triviální řešení ( u ( X ) = 0 { displaystyle u (x) = 0} ) existuje pro homogenní problém ( F ( X ) = 0 { displaystyle f (x) = 0} ). Reference ^ několik příkladů převzatých od Schulze, Hermanna: Physik mit Bleistift. Frankfurt nad Mohanem: Deutsch, 2001. ISBN 3-8171-1661-6 (Němec) Bayin, S. S. (2006). Matematické metody ve vědě a inženýrství . Wiley. Kapitoly 18 a 19. Eyges, Leonard (1972). Klasické elektromagnetické pole . New York, NY: Dover Publications. ISBN 0-486-63947-9 . Kapitola 5 obsahuje velmi čitelný popis použití Greenových funkcí k řešení problémů s hraničními hodnotami v elektrostatice. Polyanin, A.D .; Zaitsev, V.F. (2003). Příručka přesných řešení pro obyčejné diferenciální rovnice (2. vyd.). Boca Raton, FL: Chapman & Hall / CRC Press. ISBN 1-58488-297-2 . Polyanin, A.D. (2002). Příručka lineárních parciálních diferenciálních rovnic pro inženýry a vědce . Boca Raton, FL: Chapman & Hall / CRC Press. ISBN 1-58488-299-9 . Mathews, Jon; Walker, Robert L. (1970). Matematické metody fyziky (2. vyd.). New York: W. A. Benjamin. ISBN 0-8053-7002-1 . Folland, G.B. Fourierova analýza a její aplikace . Matematická řada. Wadsworth a Brooks / Cole.Cole, K.D .; Beck, J.V .; Haji-Sheikh, A .; Litkouhi, B. (2011). "Metody pro získání funkcí Greena". Vedení tepla pomocí Greenových funkcí . Taylor a Francis. 101–148. ISBN 978-1-4398-1354-6 . Zelená, G (1828). Esej o aplikaci matematické analýzy na teorie elektřiny a magnetismu . Nottingham, Anglie: T. Wheelhouse. stránky 10-12 . Faryad a M .; Lakhtakia, A. (2018). Dyadická zelená funkce nekonečného prostoru v elektromagnetismu . London, UK / San Rafael, CA: IoP Science (UK) / Morgan and Claypool (USA). Bibcode :2018idgf.book ..... F . externí odkazy Kontrolní úřad