2.2.3. Rozwiązanie dla równania niejednorodnego Niejednorodnym odpowiednikiem równania (2.3) w jednym kierunku jest | |
(2.32) | |
W podręcznikach ([L2]–Marcinkowska, [L3]–Krzyżański) znajdujemy rozwiązanie tego równania, które przy przyjętych przez nas oznaczeniach ma postać | |
F(x, τ) = g1(x, τ) + g2(x, τ) , | (2.33) |
gdzie g1(x, τ) i g2(x, τ) są rozwiązaniami pochodzącymi od warunków początkowych (2.16-2.28) i wymuszenia. | |
(2.34) | |
(2.35) | |
Mimo że jest to poprawne rozwiązanie,
to nie możemy go wykorzystać dla
rozwiązania równań Maxwella głównie z powodu przyjętych w tych podręcznikach założeń potrzebnych
do przeprowadzenia dowodu jego poprawności i stanowiących istotne ograniczenia w jego zastosowaniu. Elementem decydującym o nieprzydatności jest przyjęcie, że wymuszenie f (x, τ) pojawia się w chwili τ = 0 i dopiero od tego momentu jest określone rozwiązanie. Przyjęcie takich założeń jest zasadne przy analizie równania na przykład wspomnianej już struny, ale przy równaniach Maxwella interesuje nas rozwiązanie przy wymuszeniu poprzedzającym bez ograniczeń czasu dla analizowanego punktu. W podręczniku [L2]–Marcinkowska (strona 96) jest dodatkowy wymóg ograniczonego i zamkniętego obszaru, co dodatkowo ogranicza możliwość powoływania się na to rozwiązanie. Należy zauważyć, że przy równaniach Maxwella nieistotne staje się w pewnym zakresie rozwiązanie pochodzące od warunków początkowych ponieważ to co jest warunkami początkowymi w rozwiązaniu równań Maxwella jest falą pochodzącą od wymuszenia w innym punkcie przestrzeni i tym samym już będącą elementem rozwiązania. Z powodu wymienionych ograniczeń koniecznym było opracowanie własnego, wolnego od ograniczeń rozwiązania. Przedstawione rozwiązanie jest łatwe do zrozumienia, bo oparte na elementarnej wiedzy o równaniach różniczkowych bez odwoływania się do twierdzeń wymagających dodatkowej wiedzy. Po zamianie zmiennych w równaniu (2.32) mamy | |
(2.36) | |
i otrzymujemy formalne rozwiązanie | |
(2.37) | |
ale wyprowadzenie z niego praktycznych wzorów nie jest łatwe i dlatego wykorzystamy łatwiejszą drogę. Rysunek (Rys.2) wyjaśnia znaczenie zmiennych i wykonywanych operacji. |
x = x0
τ = τ0
τ = τC
O(u0,v0)= O(x0, τ0)
A(uA,v0)
C(x0,τc)
B(u0,vB)
v
vC
vB
u
uC
uA
| |
Rys.2 | |
Na tym rysunku – u, v, x, τ oznaczają odpowiednie osie współrzędnych, – Fo(uo , vo) jest wartością funkcji w punkcie (uo , vo), czyli (xo , τo) oznaczonym na rysunku O(xo , τo) = O(uo , vo), – Fo(uc= xo + τc , vc= xo - τc) jest wartością w punkcie (xo , τc) oznaczonym na rysunku C(xo , τc) , którą możemy uważać za wartość początkową funkcji F(). – Kolorem czerwonym zaznaczono falę biegnącą w kierunku rosnących wartości x i odpowiednio kolorem niebieskim falę falę biegnącą w kierunku malejących wartosci x Całka przedstawia zmianę wartości funkcji F(x, τ) przy przejściu z punktu C(xo , τc) do punktu O(xo , τo), czyli zmianę wartości funkcji F(x, τ) w punkcie xo w czasie od τc do τo , przy czym τc ≤ τo . Dla naszego rozwiązania wykorzystamy dla przejścia z punktu C do punktu O na rysunku (Rys.2) zależność | |
(2.38) | |
Ponieważ będziemy wykorzystywać różne i częste zamiany funkcji oraz zmiennych, to aby czytelnik nie tracił czasu na weryfikację tych zmian przypomnimy podstawowe zależności: | |
u = x + τ , v = x - τ | (2.39) |
(2.40) | |
u - v = 2 ⋅ τ , u + v = 2 ⋅ x; | (2.41) |
F(x, τ) = Fo(u, v) | (2.42) |
(2.43) | |
Fo(x, τ) = Fo(x + τ, x - τ) | (2.44a) |
Fo(u, v) = Fo(u, 2 ⋅ x - u) = Fo(2 ⋅ x - v, v) | (2.44b) |
Fo(u, v) = Fo(u, u - 2 ⋅ τ) = Fo(v + 2 ⋅ τ, v) | (2.44c) |
Konsekwencją tych zależności są następne; | |
(2.45) | |
(2.46) | |
(2.47) | |
(2.48) | |
(2.49) | |
(2.50) | |
(2.51) | |
(2.52) | |
(2.53) | |
(2.54) | |
Zależność (2.50) możemy wykorzystać w (2.38) i wtedy otrzymujemy | |
(2.55) | |
Składniki po prawej stronie znaku równości we wzorze (2.55) możemy łatwo przy stałym τ określić względem wybranego punktu xo na osi x z zależności | |
(2.56) | |
(2.57) | |
co po skorzystaniu z (2.54) daje | |
(2.58) | |
(2.59) | |
Na podstawie (2.40) mamy | |
(2.60) | |
czyli | |
(2.61) | |
(2.62) | |
i w konsekwencji | |
(2.63) | |
(2.64) | |
Wstawiając te wielkości do (2.55) otrzymujemy ostatecznie | |
(2.65) | |
i jest to odpowiednik wzoru (2.35) otrzymany tym razem bez żadnych
założeń ograniczających czas τc , który możemy tym razem wydłużać nawet do -∞. Stała F(xo, τ) jest wynikiem całki funkcji f(x, τc) dla czasu przed τc . Przy wyprowadzaniu wzoru (2.35) zakładano, że funkcja f(x, τc) = 0 dla τ < τc i wtedy całka z tej funkcji, czyli nasza wartość F(xo , τc) , jest też równa zero. W rozwiązaniu równań Maxwella nie występuje potrzeba uwzględniania rozwiązań pochodzących od warunków początkowych F(x, τc) oraz ∂F(x, τc) /∂τ , które można dodać dla innych przypadków niż równania Maxwella, ustalając zależność tych warunków od źródeł f(x, τ) . Zmienną τo wprowadzono dla zademonstrowania identyczności wzorów (2.35) i (2.65), ale ponieważ oznacza ona w naszym przypadku czas aktualny, to przyjmując τo = 0 możemy, nie tracąc ogólności, uczynić krótsze i bardziej przejrzyste wzory. Wzór (2.65) w tej sytuacji wygląda następująco | |
(2.66) | |
gdzie jeszcze możemy pozbyć się stałej
F(xo,
τc)
, wydłużając czas τc
do minus nieskończoności i mamy wtedy | |
(2.67) | |
(2.68) | |
Stałe całkowania (2.61) i (2.62)
we wzorach (2.63) i (2.64) są jednakowe i we wzorze (2.65) występują z przeciwnymi znakami, co powoduje, że znoszą się wzajemnie. Wartości ∂Fo(u, vo) / ∂u oraz ∂Fo(uo , v) / ∂v na rysunku (Rys.2) odpowiadają wartościom tych pochodnych w punktach leżących na przecięciu linii odpowiadającej czasowi τ oraz linii odpowiadających zmiennym uo oraz vo , które to zmienne oznaczają odpowiednio falę w stronę malejących wartości x oraz falę w stronę rosnących wartości x . Możemy te fale na wzór równania jednorodnego rozdzielić i nazwać je jak na rysunku (Rys.2) Fv(u, vo) i Fu(uo , v) . Mamy zatem | |
Fo(uo, vo) = Fu(uo , vo) + Fv(uo, vo) | (2.69) |
(2.70) | |
(2.71) | |
We wzorze (2.70) całkowanie względem τ ,
odbywa się po prostej vo
a to oznacza, że v = x - τ = vo = xo - τo = xo , czyli τ = xo - x i dτ = dx . Dla odróżnienia x pochodzącego od zastąpienia τ (na podstawie dτ = dx ) przypiszemy mu apostrof x' i wzór (2.70) możemy zapisać | |
(2.72) | |
Podobnie we wzorze (2.71) całkowanie odbywa się po prostej
uo , co oznacza, że u = x + τ = uo = xo + τo = xo , czyli τ = xo - x i dτ = - dx i wzór (2.71) ma postać | |
(2.73) | |
Wzory stają się bardziej czytelne, gdy skorzystamy z jeszcze jednej zmiennej pomocniczej r . | |
(2.74) | |
(2.75) | |
Wzory (2.69) , (2.74) i (2.75)
są rozwiązaniami ogólnej postaci (2.32) jednokierunkowego, niejednorodnego równania Maxwella. |