Fale/Fale
| Fale | |||
|
||||
[edytuj] Układy fizyczne, a jego drgania swobodne
Ruchem drgającym harmonicznym nazywamy drgania, które zależą od częstotliwości kołowej drgań własnych ω i przesunięcia fazowego φ. Drgania te opisujemy względem czasu "t" wzorem:
![]() |
(1.1)
|
Ruch pokazany przy pomocy równości (1.1), są to ruchy o jednym stopniu swobody. Okresem drgań harmonicznych nazywamy odwrotność częstości ν, którą możemy przepisać przy pomocy definicji częstotliwości kołowej:
![]() |
(1.2)
|
[edytuj] Siła kierujące i bezwładność ciała dla drgań o jednym stopniu swobody
Ruch harmoniczny (1.1) jest rezultatem działania dwóch przeciwstawnych czynników, tzn. siły kierującej i bezwładności. Siła kierująca w ruchu harmonicznym usiłuje przywrócić układ do stanu równowagi, a bezwładność jest cechą układu, która przeciwstawia się zmianom ruchu układu.
[edytuj] Drgania tłumione
Ruch harmoniczny opisywany wzorem (1.1) nie jest ruchem tłumionym, bo jego amplituda nie zmienia się w miarę upływu czasu. W ruchu tłumionym, amplituda drgań maleje z czasem. Przykładem drgań gasnących są drgania opisane równością:
![]() |
(1.3)
|
[edytuj] Wahadło matematyczne
Wahadłem matematycznym nazywamy układ, w którym ciało o masie M jest zawieszone na nieważkiej nici lub pręcie o długości l. Siłą kierującą w ruchu wahadła jest składowa siły grawitacji prostopadła do nici, jej wartość określa wzór F=Mg sinψ, to z drugiej zasady dynamiki Newtona w tym przypadku wynika:
![]() |
(1.4)
|
Funkcję sinus możemy rozłożyć w szereg Taylora, którego zapis dla naszej funkcji jest następujący:
![]() |
(1.5)
|
Dla małych kątów ψ, można pominąć wszystkie wyrazy szeregu poza pierwszym, wówczas równość (1.4) możemy przestawić (1.5):
![]() |
(1.6)
|
Rozwiązaniem równania (1.6) są funkcje harmoniczne w postaci (1.1) (co udowodnimy poniżej), dla której druga pochodna tego właśnie przypuszczalnego rozwiązania naszej funkcji jest:
![]() |
(1.7)
|
Jeśli funkcję (1.7) wraz (1.1) podstawimy do wzoru (1.6), w ten sposób otrzymamy tożsamość, co dowodzi o harmoniczności ruchu masy M w wahadle matematycznym.
[edytuj] Ciężarek o masie M przyczepiony do dwóch sprężynek
Rozważmy dwie sprężynki przeczepione do ścianek i połączone ze sobą poprzez masę M. Sprężynki mają znikomą masę. Całkowita siła działająca na masę M jest wyrażona:
![]() |
(1.8)
|
Ruch ciężarka możemy opisać przy pomocy drugiego prawa dynamiki Newtona, który piszemy:
![]() |
(1.9)
|
Ruch ciała M będziemy opisywać funkcją ψ(t), którą definiujemy jako starą funkcję z(t) minus stałą "a", zatem w (1.9) możemy dokonać zamiany zmiennych według:
![]() |
(1.10)
|
Równanie (1.9) przy pomocy podstawienia (1.10), którym występuje poszukiwana funkcja ψ(t) jest zależna od czasu, dla równania, która jest równaniem oscylatora harmonicznego, możemy zapisać w postaci:
![]() |
(1.11)
|
Kwadrat częstości kołowej ruchu harmonicznego ciała, a także rozwiązanie ruchu (1.9) na podstawie (1.11) i (1.10), przestawiamy:
|
|
[edytuj] Drgania poprzeczne ciężarka przyczepionego do sprężyn
Rozważmy ciężarek o masie M przyczepiony do dwóch nieważkich sprężyn. Współczynnik sprężystości sprężyn wynosi K. W położeniu równowagi każda ze sprężyn zostaje napięta do długości a. Opisywany układ może wykonywać dwa rodzaje drgań poprzecznych i jedno drganie podłużne. W stanie równowagi siła działająca na ciało z jednej ze sprężyn jest opisywana:
![]() |
(1.14)
|
Opisu drgań podłużnych nie będziemy powtarzać, ponieważ jest zamieszczony w poprzednim rozdziale. W drganiach poprzecznych siła, z jaką pojedyncza sprężyna działa na ciężarek M, jest opisana wzorem:
![]() |
(1.15)
|
Z drugiej zasady dynamiki Newtona ruch opisujący drgania poprzeczne jest opisywany przez równanie różniczkowe:
![]() |
(1.16)
|
Należy zauważyć, że długość l występująca po prawej stronie wzoru końcowego (1.16) jest funkcją w ogólności od zmiennej x.
[edytuj] Mała długość swobodna sprężyny jako pierwsze przybliżenie
Równanie ruchu ciała o masie M (1.16) w przybliżeniu małej długości sprężyny, tzn. stosunek długości sprężyny swobodnej a0 przez długość sprężyny l jest w naszym przypadku tutaj rozważanym wielkością pomijalną. Równanie ruchu (1.16) przybliżeniu tutaj omawianym jest w postaci:
![]() |
(1.17)
|
Wzór (1.17) uzyskaliśmy na podstawie dokładnego rozwiązania (1.16), w której definicja kwadratu częstotliwości kołowej dla przybliżenia małej długości sprężyny przy zachodzącym (1.14) jest:
![]() |
(1.18)
|
Rozwiązaniem równania ruchu (1.17) jest to równanie zależne od częstotliwości kołowej ω, które jest opisane wzorem (1.1). Widzimy, ze zgodnie ze wzorami (1.18) i (1.12), że częstość drgań poprzecznych w naszym przybliżeniu jest równa częstości drgań podłużnych w ruchu drgającym ciała o masie M.
[edytuj] Przybliżenie małych drgań
Aby rozwiązać równanie różniczkowe (1.16) stosowaliśmy przybliżenia małej długości swobodnej sprężyny, wtedy stosunek długości swobodnej sprężyny a0 do długości aktualnej sprężyny jest względnie mały, których stosunek w takiej sytuacji możemy pominąć. Jeśli tego stosunku nie można pominąć, to możemy powiedzieć, że dla małego ε, względem którego będziemy rozkładali odwrotność długości chwilowej sprężyny w szereg Taylora, który to podstawimy do równania różniczkowego (1.16), i pominiemy wyrazy wyższe niż liniowe, w ten sposób dojdziemy do równania różniczkowego, które jest równaniem różniczkowym ruchu harmonicznego, to co chcieliśmy otrzymać. Zatem napiszmy teraz definicję kwadratu chwilowej długości sprężyny "l" i kwadratu wielkości ε:
|
|
Odwrotność długości sprężyny liczymy rozkładając tą odwrotność w szereg Taylora względem zmiennej x, która jest położeniem ciała o masie M:
![]() |
(1.21)
|
Wzór (1.21) podstawiamy do (1.16) do odjemnika w nawiasie, w ten sposób możemy otrzymać równość zależną od zmiennej x, którego to równanie jest równaniem oscylatora harmonicznego:
![]() ![]() |
(1.22)
|
Równanie ruchu (1.22) możemy przepisać dokonując odpowiednich przekształceń w nawiasie w końcowej wspomnianej tożsamości:
![]() |
(1.23)
|
W równaniu ruchu masy M zapisaną w punkcie (1.23) możemy pominąć wyższe wyrazy niż liniowe przy przybliżeniu małych drgań, by potem wykorzystać wzór (1.14):
![]() |
(1.24)
|
Kwadrat częstotliwości kołowej ruchu drgającego drgań poprzecznych wynika z równania różniczkowego (1.24):
![]() |
(1.25)
|
Rozwiązanie równania ruchu (1.24) przy definicji kwadratu częstotliwości kołowej (1.25) jest równaniem zapisane w punkcie (1.1). Możemy wyznaczyć iloraz częstotliwość drgań podłużnych (1.12) przez częstotliwość kołową drgań poprzecznych (1.25) przy oznaczeniu (1.14) w następującej postaci:
![]() |
(1.26)
|
Patrząc na wzór (1.26), gdy długość własna sprężyny a0 jest względnie mała w porównaniu z długością sprężyny w stanie równowagi "a", to częstotliwości kołowe drgań podłużnych i poprzecznych są sobie w przybliżeniu równe.
[edytuj] Obwód LC, czyli solenoid z dwoma kondensatorami
Równanie opisujących układ LC z dwoma kondensatorami o pojemności C i z cewką od indukcji elektrycznej L przestawiamy:
![]() |
(1.27)
|
Ładunek Q1 jest powiązany z ładunkiem Q2, a także ten drugi ładunek jest powiązany z natężeniem prądu elektrycznego płynącego w obwodzie:
|
|
Równanie (1.27) możemy tak napisać przy wykorzystaniu tożsamości (1.28), które to możemy dalej zróżniczkować obustronnie względem czasu przy wykorzystaniu (1.29), wtedy otrzymamy równość różniczkową drugiego stopnia:
![]() |
(1.30)
|
Końcowe równanie (1.30) jest równaniem oscylatora harmonicznego, którego kwadrat częstotliwości kołowej przestawiamy:
![]() |
(1.31)
|
Natężenie prądu elektrycznego, który jest rozwiązaniem równania różniczkowego (1.30) przestawiamy podobnym wzorem do (1.1), jest to równanie w postaci:
![]() |
(1.32)
|
[edytuj] Zasada superpozycji i liniowość równań ruchu
Omawiane tutaj były przypadki, w których siła kierunkowa była wprost proporcjonalna do ψ(t), natomiast niezależała od wyższych potęg ψ, tzn. ψ2,ψ3,.... Równanie niezawierające wyższych potęg ψ niż pierwszego stopnia nazywamy równaniem liniowym. Gdy w naszym równaniu występują wyższe potęgi niż pierwszego stopnia, tak powstałe równanie nazywamy nieliniowym. Natomiast, gdy równanie nie zawiera wyrazów z niezależnych od ψ nazywamy równaniem jednorodnym.
[edytuj] Wprowadzenie do równań liniowych jednorodnych
Równania liniowe mają tę właściwość, że suma dwóch rozwiązań liniowych też spełnia równanie liniowe i jednorodne, co jest wiadomo. Udowodnimy czy też jest dla równań różniczkowych nieliniowych, że jeśli mamy dwa rozwiązania poniższego równania różniczkowego, to suma ich ma spełniać jest też jego rozwiązaniem, jak wykażemy co jest jedynie możliwe, gdy mamy równanie różniczkowe liniowe. A oto mamy równość różniczkową w ogólności zależną od wyższych potęg niż pierwszego stopnia:
![]() |
(1.33)
|
Weźmy dwa rozwiązania ψ1 i ψ2, które spełniają równanie różniczkowe (1.33):
|
|
Sumę rozwiązań ψa i ψb też spełnia równanie różniczkowe (1.33) jak i składowe tejże sumy:
![]() |
(1.36)
|
Równanie (1.34) możemy dodać do (1.35) do siebie i porównać ze wzorem (1.36), w ten sposób otrzymujemy tożsamości:
|
|
||||
|
|
Równanie (1.37) i (1.38) są tożsamościami, tzn. które są zawsze spełnione, a równości (1.39) i (1.40) są zawsze spełnione tylko dla α=0 i β=0, zatem co zachodzi tylko dla rozwiązania liniowego przy postawionych warunkach powyżej, to superpozycja rozwiązań równania różniczkowego liniowego jest też rozwiązaniem tego samego równania różniczkowego. Co kończy dowód.
[edytuj] Warunki początkowe a superpozycja rozwiązań
Załóżmy, że mamy dwa rozwiązania jako szczególne rozwiązania w chwilach szczególnych ψ1(0) i ψ2(0), które odpowiadają warunkom brzegowym rozwiązań ψ1(t) i ψ2(t), a rozwiązania równania będące superpozycją naszych rozwiązań, czyli ψ1(t)+ψ2(t), odpowiadają warunki brzegowe w postaci: ψ1(0)+ψ2(0).
[edytuj] Niejednorodne równanie liniowe
Wprowadźmy do równania liniowego siłę wymuszającą F(t) zależną od czasu, wtedy po tej operacji otrzymujemy równanie różniczkowe:
![]() |
(1.41)
|
Załóżmy, że mamy dwa równania różniczkowe (1.41), każde odpowiada innej sile wymuszającej, to superpozycja dwóch rozwiązań spełnia równanie różniczkowe odpowiadające sumie dwóch sił wymuszających w równaniu (1.41), tzn. F1(t)+F2(t).
[edytuj] Wahadło sferyczne
W celu zilustrowania zasady superpozycji rozpatrzmy wahadło sferyczne zawieszone na nieważkiej nici, która posiada możliwość ruchu w dowolnym kierunku. Położenie równowagi znajduje się w punkcie x=0 i y=0, co równania ruchu takiego wahadła piszemy:
|
|
Powyższe dwa równania są rozdzielone, tzn. każde równanie zależy od innych zmiennych i ich rozwiązania są niezależne i są opisywane wzorami tylko od zmiennej t i są to dwa rozwiązania:
|
|
Kwadrat częstotliwości kołowej występujących w rozwiązaniach (1.44) i (1.45) wynikających z równań (1.42) i (1.43), zapisujemy wzorem (1.6).
[edytuj] Układy o dwóch stopniach swobody i jego drgania swobodne
Opis układów w ogólnej konfiguracji wymaga znajomości dwóch zmiennych, tzn. ψa i ψb, co oznacza np. ruch wahadła wokół płaszczyzn wzajemnie prostopadłych kierunkach w płaszczyźnie drgań, a w przypadku układów LC, co oznacza, że ψa i ψb oznaczają ładunki na dwóch kondensatorach lub też natężenia prądów elektrycznych. Najbardziej ogólny ruch harmoniczny jest superpozycja dwóch drgań harmonicznych, których jak wiemy też jest w pewnym senie ruchem harmonicznych. Składowe tych rozważanych drgań harmonicznych nazywamy drganiami normalnymi lub drganiami własnymi. Można tak dobrać warunki początkowe, by układ wykonywał jedno lub drugie drganie harmoniczne. Załóżmy, że układ wykonuje drgania z częstotliwościami kołowymi ω1 i ω2. Przesunięcia fazowe dla obu drgań normalnych są sobie równe. Dla obu tychże drgań wychylenia od stanu równowagi są opisywane przez:
|
|
Dla drugiej postaci drgań normalnych, zachodzą podobne związki jak w punktach (1.46) I (1.47), tylko jedynkę zamieniamy na dwójkę. W ogólnej konfiguracji wychylenia od położenia równowagi dla ψa i ψb zapisujemy jako superpozycja drgań harmonicznych o częstotliwościach ω1 i ω2:
![]() |
(1.48)
|
a dla drugiej zmiennej, tzn. ψb drgania opisujemy:
![]() |
(1.48)
|
[edytuj] Wahadło sferyczne jako układ drgający o dwóch stopniach swobody
Dla wahadła sferycznego na podstawie równań różniczkowy opisujący te nasze drgania , tzn. równań (1.42) i (1.43) dla dwóch niezależnych drgań, tzn. w niezależnych zmiennych, dla którego kwadrat częstotliwości kołowej jest opisywany wzorem (1.6), których to wychylenia od stanu równowagi dla tych rozważanych zmiennych są przestawione jako:
|
|
[edytuj] Przybliżenie harmoniczne dla dwuwymiarowego układu współrzędnych
Rysunek obok przestawia ruch ciała o masie M przyczepionych do czterech sprężyn o stałych proporcjonalnych K1 i K2, które są przyczepione do odpowiednich ścianek. Masy sprężynek są w przybliżeniu zerowe w porównaniu z masą kulki M. Kulka może wykonywać dwa niezależne drgania wzdłuż osi iksowej i igrekowej lub w postaci złożenia tychże drgań. We wzorach opisujących ruch masy M pomijamy człony x2/a2, y2/a2 i xy/a2. W ten sposób możemy napisać przybliżenie klasyczne małych drgań. Jeśli dokonamy odpowiedniego małego przesunięcia ciała M w kierunku osi iksowej, to wtedy z dokładnością do wyrazów wyżej podanych otrzymamy wzór na siłę iksową i igrekową w postaci:
|
|
Podobnie postępujemy, gdy dokonamy małego przesunięcia wzdłuż osi igrekowej. Ogólne równania ruchu opisujące drgania masy M są to równania ruchu zapisane w przybliżeniu małej długości sprężyny, które możemy napisać je jako dwa równania opisujące ruch ciała M na płaszczyźnie xy, których rozwiązania są superpozycją drgań wzdłuż osi iksowej i igrekowej:
|
|
Rozwiązania równania (1.53), a także kwadrat częstotliwości kołowej możemy zapisać w postaci poniższych tożsamości:
|
|
Wychylenie od punktu początkowego względem osi igrekowej jako rozwiązanie równania różniczkowego (1.54) i kwadrat jej częstotliwości kołowej piszemy:
|
|
[edytuj] Współrzędne normalne w układzie obróconym o pewien kąt
Albowiem równania (1.53) i (1.54) opisują ruchy harmoniczne, ale w ogólności to nie ma ogólnych rozwiązań w postaci rozwiązań (1.48) i (1.49). W naszym przypadku opisujące ruch harmoniczny. Aby łatwiej można było opisywać ruch harmoniczny kulki wprowadza się współrzędne x'. y', które są obliczone względem starych współrzędnych obrócone o kąt α, których są kombinacjami liniowymi współrzędnych starych x i y. Są to współrzędne, których równania ruchu masy M opisywane są względem jednej współrzędnej, a każde równanie różniczkowe zależy od innej zmiennej. W ogólnym przypadku wybór współrzędnych normalnych nie jest taki łatwy. Zazwyczaj jest tak, że równania ruchu zawierają współrzędne x i y wziętej razem. W ogólnym przypadku prosta transformacja nie wystarcza i należy dokonać ogólnej transformacji liniowej w układzie równań różniczkowych opisujących nasze ruchy, by w ten sposób każde równanie było opisywane we współrzędnych normalnych, tzn. każde równanie zależało tylko od jednej zmiennej.
[edytuj] Rozwiązywanie równań różniczkowych ruchu
Załóżmy, że mamy układ fizyczny, która drga, który jest opisywany przy pomocy równań różniczkowych zapisane we współrzędnych x i y w postaci:
|
|
Przyjmijmy, że rozwiązaniami równań ruchu (1.59) i (1.60) są to rozwiązania w postaci funkcji x i y, które są opisywane w postaci funkcji zależnych od czasu i są zależne od takiego samego przesunięcia fazowego.
|
|
Jeśli układ wykonuje drgania normalne wzdłuż jednej osi, to wtedy równania (1.59) i (1.60) są opisywane równaniami dla każdego ruchu z osobna:
|
|
W ogólnym przypadku, gdy nie mamy doczynienia ze współrzędnymi normalnymi, to wtedy należy rozwiązania (1.61) i (1.62) podstawić do (1.59) i (1.64), w ten sposób otrzymać następujący układ równań:
|
|
Z równań (1.62) i (1.66) możemy wyznaczyć stosunek x/y wedle następujących sposobów:
|
|
Warunkiem koniecznym by rozwiązać równania różniczkowe ruchu (1.59) i (1.60), jest równość obu stron (1.67) i (1.68):
![]() |
(1.69)
|
Równanie (1.69) jest równaniem kwadratowym względem zmiennych ω2, którego rozwiązaniem są dwa kwadraty częstotliwości ω12 i ω22. Stosunek y/z względem równości (1.67) dla jednej i drugiej postaci drgań normalnych nazywamy stosunkiem w postaci:
|
|
Rozwiązaniami równania różniczkowego falowego (1.59) i (1.60) dla współrzędnej x jest superpozycją rozwiązań równania falowego dla częstotliwości ω1 i ω2, co jest to rozwiązaniem w postaci:
![]() |
(1.72)
|
A rozwiązaniami y(y) według (1.71) są to rozwiązania przestawione w postaci równania:
![]() ![]() |
(1.73)
|
[edytuj] Drgania podłużne, układ sprężynek i ciężarków
Rozpatrzmy układ trzech sprężynek o masie zaniedbywanej, w których dwie sprężynki są połączone ze ściankami, a sprężynka środkowa jest połączona z ciężarkami o masie M. Rozważmy sytuację, gdy dwa ciężarki są przesunięte o to samo przesunięcie, wtedy sprężynka środkowa nie jest napięta, Siła z jaką sprężynka działa na pierwsze ciało jest określona wzorem F1=-Kψa, a na drugie ciało wzorem F2=-Kψb, przy założeniu ψa=ψb, i pamiętając że sprężynka środkowa nie oddziaływuje w takim przypadku z ciężarkami, równanie ruchu jest opisane:
|
|
Kwadrat częstotliwości kołowej z jaką drga ten układ dwóch mas jest określony przez wzór:
![]() |
(1.76)
|
Rozważmy teraz sytuację, gdy ψa=-ψb, wtedy siła działająca na pierwszy ciężarek, której wartość jest równa potrojonej wartości współczynnika sprężystości i przesunięcia od stanu równowagi ciężarka pierwszego, której kierunek tej siły jest przeciwny niż kierunek osi zet:
![]() |
(1.77)
|
A całkowita siła działająca na drugi ciężarek, której wartość jest równa potrojonej wartości współczynnika sprężystości i przesunięcia od stanu równowagi ciężarka pierwszego, kierunek tej siły jest przeciwny niż kierunek osi zet:
![]() |
(1.78)
|
Równania ruchu dla ciężarka pierwszego i drugiego są określone przez wzory wynikających z drugiego prawa Newtona:
|
|
Kwadrat częstotliwości drgań swobodnych drgań normalnych, który dla drgań dwóch ciężarków opisywanych przez równania różniczkowe (1.79) i (1.80), przestawiamy:
![]() |
(1.81)
|
Na podstawie przypadków szczególnych jakie dokonaliśmy w powyższych rozważaniach dowiadujemy się, że układ może drgać na dwa sposoby, której kwadrat częstotliwości jest opisany wzorem (1.76) lub (1.81).
Rozważmy teraz przypadek ogólny, w którym opisywać będziemy ruch dwóch tych ciężarków, które są opisywane równaniami wynikających z drugiego prawa Newtona:
|
|
Równania (1.82) i (1.83) możemy do siebie dodać i je odjąć, w ten sposób otrzymać następne dwa równania zależne od zmiennych ψa+ψb i ψa-ψb:
|
|
Na podstawie równości (1.84) otrzymujemy pierwsze równanie na kwadrat częstotliwości kołowej drgań (1.76), a także z równania (1.85) wynika drugi kwadrat częstotliwości kołowej drgań normalnych (1.81). Z równań (1.84) i (1.85) wnioskujemy, że można napisać następujące dwa równania:
|
|
Całkowite rozwiązanie ruchu opisujące przez równanie (1.86) i równanie (1.87) są opisywane dla Ψa:
![]() |
(1.88)
|
A także rozwiązanie opisujące ruch ciała drugiego jest opisywane:
![]() |
(1.89)
|
[edytuj] Drgania poprzeczne, układ sprężynek i ciężarków
Rozważmy teraz sobie konfigurację sprężynek wykonujących drgania poprzeczne. Rozpatrzmy teraz konfigurację, w którym środkowa sprężynka nie jest wcale naprężona, wtedy równanie ruchu pierwszej masy przestawiamy:
![]() ![]() |
(1.90)
|
Równanie ruchu drugiego ciała jest bardzo podobne do (1.90). To równanie rozwiązuje się tak samo jak (1.16) pamiętając przy tym o (1.15), zatem kwadrat częstotliwości kołowej opisujący pierwsze drganie normalne dla ψa=ψb liczymy według:
![]() |
(1.91)
|
Następną postacią drgań normalnych jest to drganie, którego dla tego ruchu zachodzi ψA=-ψb, kąt pomiędzy sprężynką pierwszą a poziomem jest dwukrotnie mniejszy niż kąt pomiędzy sprężynką drugą a poziomem, co będziemy korzystali w wyprowadzeniu równania ruchu dla drugiej postaci drgań normalnych:
![]() |
(1.92)
|
Równanie ruchu dla drugiej sprężynki jest bardzo podobne do (1.92). Równanie (1.92) rozwiązuje się bardzo podobnie jak równanie (1.16) pamiętając przy okazji o (1.15), zatem kwadrat częstotliwości drgań normalnych możemy przestawić:
![]() |
(1.93)
|
[edytuj] Dwa sprzężone ze sobą obwody LC
Układ dwóch sprzężonych ze sobą obwodów LC nazywamy układ składający się z trzech kondensatorów o pojemności C równoległych względem siebie i dwóch cewek o indukcyjności własnej L też równoległych wobec siebie. Napiszmy teraz równania różniczkowe według drugiego prawa Kirchhoffa, które rządzą dwoma oczkami z osobna, które to powiązane są ze sobą względem prądów płynących przez cewki Ia i Ib, i przy pomocy ładunków znajdujących się na poszczególnych kondensatorach Q1, Q2 i Q3:
|
|
Możemy z różniczkować równania (1.94) i (1.95) względem czasu i w ten sposób otrzymujemy stąd dwie dalsze tożsamości zawierające prądy Ia, Ib, a także pochodnych ładunków znajdujących się na kondensatorach Q1, Q2, Q3 policzonych względem czasu:
|
|
Z pierwszego prawa Kirchhoffa możemy otrzymać trzy tożsamości, które wiążą ładunki na kondensatorach z natężeniami prądów płynących w poszczególnych obwodach:
|
|
|
Tożsamości (1.98), (1.99) i (1.100) możemy podstawić do wzorów (1.96) i (1.97), w ten sposób otrzymujemy dwie równości różniczkowe zależne od zmiennych Ia i Ib i od pozostałych parametrów:
|
|
Równania (1.101) i (1.102) możemy do siebie dodać i odjąć, w ten sposób otrzymać dwa równania zależne od zmiennych Ia+Ib oraz od Ia-Ib z osobna, z których możemy wyznaczyć te wielkości w postaci ich drgań z pewnymi częstotliwościami kołowymi określonych kolejno przez ω1 i ω2, zatem te dwa równania piszemy jako:
|
|
Z równań (1.103) i (1.104) możemy otrzymać dwie wspomniane wcześniej częstotliwości kołowe drgań normalnych, której to kwadraty:
|
|
Układ drga z jedną częstotliwości kołowych, którego kwadraty są opisywane wzorami (1.105), (1.106), gdy dla pierwszej mamy Ia=Ib, a dla drugiej Ia=-IB. Ogólne rozwiązania na natężenia prądów Ia i Ib są zapisane jako:
|
|
[edytuj] Dudnienia
Jest wiele przykładów układów drgających z dwiema częstotliwościami podstawowymi ω1 i ω2, przykładem takich drgań mogą być drgania o dwóch stopniach swobody. Przykładem układów o dwóch stopniach swobody są drgania spowodowane przez dwie siły wymuszające drgających z dwoma częstotliwościami. Aby zobaczyć co to są dudnienia, rozpatrzmy dwa drania o takich samych amplitudach i różnych częstotliwościach podstawowych, które są drania opisane przez:
|
|
Całkowite wychylenie od położenia równowagi określamy jako sumę wychyleń (1.109) i (1.110) i są one zależne od częstotliwości kołowych ω1 i ω2:
![]() |
(1.111)
|
[edytuj] Modulacje
Średnią częstotliwość kołową ωśr, która jest sumą częstotliwości kołowych ω1 i ω2 podzielonych przez dwa, oraz częstotliwość kołową modulacji ωmod, która jest różnicą dwóch częstotliwości kołowych podanych wcześniej i podzielonych przez dwa:
|
|
Częstotliwości drgań kołowych drgań normalnych możemy napisać w zależności od częstotliwości drgań średnich i częstotliwości drgań modulacji:
|
|
Jeszcze raz patrząc na wzór (1.111) przy takich samych amplitudach A i B, które tutaj oznaczymy przez A, a także na wzory (1.114) i (1.115), wtedy możemy dojść do wniosku:
![]() ![]() |
(1.116)
|
- gdzie amplituda modulacyjna określamy według wzoru (1.116):
![]() |
(1.117)
|
Wzór (1.116) jest wzorem ogólnym, ale największe korzyści daje nam, gdy częstotliwości ω1 i ω2 są do siebie zbliżone, to wtedy mamy na pewno ω1≈ω2 i ωmod<<ωsr, a więc amplituda praktycznie się nie zmienia i wtedy układ wykonuje drania prawie harmoniczne o częstości ωsr. Natomiast gdy ωmod=0, to wtedy mamy ωsr=ω1=ω2, jeśli te częstotliwości kołowe różnią się niewiele, to taki ruch nazywamy prawie harmonicznym lub prawie monochromatycznym o częstotliwości ωsr i o amplitudzie prawie się nie zmieniającej.
[edytuj] Dudnienia wywołane przez dwa kamertony
Rozważmy dwa kamertony, które wywołują fale dźwiękowe, które powodują zmiany ciśnienia w pobliżu błony bębenkowej, której jako ψ jest różnicą ciśnień błony bębenkowej pomiędzy jej warstwą zewnętrzną a wewnętrzną. Jeśli mamy dwa kamertony, to ta różnica na naszej błonie jest równa sumie różnic ciśnień ψa+ψb, czyli jest superpozycją dwóch fal pochodzących od kamertonów. Gdy różnica częstotliwości ν różnią się o więcej niż 6%, to nasz mózg, wolą opis drgań opisany wzorem (1.111). Oznacza to, że fale o dwóch częstotliwościach są słyszane jako dwa odrębne dźwięki. Jeśli natomiast mamy
, to słyszane dwa dźwięki oddzielona są wielką tercją. Jeśli natomiast
, wtedy dźwięk ν2 słyszany jest jako półtonu wyższy od ν1. Jeśli natomiast dźwięki różnią o mniej niż 10Hz, to nasze ucho ma kłopoty z rozróżnianiem ich jako odrębne tony. W tym wypadku dźwięki są słyszane jako dźwięk o częstotliwości νsr o wolno zmiennej amplitudzie, który jest napisany zależnością (1.117) dla różnicy ciśnień ψ(t) danym wzorem (1.116).
[edytuj] Detektor mocy superpozycji fali dźwiękowej
Według wzoru (1.117) amplituda fali dźwiękowej oscyluje z częstością kołową modulacji ωmod, co za każdym razem ωmodt zwiększa się o 2π, w której amplituda modulacji Amod wykonuje pełny cykl modulacji (mamy tu na myśl powolne zmiany zmiany z częstotliwością modulacji ωmod). W tym czasie amplituda Amod dwa razy przyjmuje wartość równą zero, w której ludzkie ucho nic nie słyszy, a także przyjmuje dwa razy wartość bezwzględną równą wartości maksymalnej, tzn. raz maksymalną wartość wziętej z minusem, a za drugim razem z plusem, czyli w tym czasie amplituda przyjmuje dwa różne znaki o takiej samej wartości bezwzględnej amplitudy. Ludzkie ucho nie rozróżnia z jakim znakiem amplitudy mamy do czynienia, ale rozróżnia jedynie kwadrat amplitudy fali dźwiękowej, który ma dwa maksima w jednym trakcie okresu.
Częstością dudnień nazywamy częstość powtarzania się maksymalnego kwadratu amplitudy fali modulowanej, jest ona pisana wzorem:
![]() |
(1.118)
|
Aby przekonać się o tym, że kwadrat amplitudy modulowanej, czyli kwadrat z (1.117) drga z częstotliwością kołową (1.118), to napiszmy następujące przekształcenia wykorzystując znane związki z trygonometrii, wtedy:
![]() |
(1.119)
|
Widzimy na podstawie obliczeń (1.119), że kwadrat amplitudy oscyluje z podwojoną częstotliwością modulowaną (1.118), co było do wykazania.
[edytuj] Dudnienia słabo sprzężonych ze sobą wahadeł
Weźmy układ dwóch wahadeł słabo sprzężone ze sobą, czyli dwie masy M zawieszone na nieważkich niciach, oraz połączone ze sobą przy pomocy sprężynki o stałej sprężystości K. Rozwiązując równania ruchu dla tych mas M, dowiemy się, że układ może drgać z dwiema częstotliwościami kołowymi, które są opisywane w zależności od stałej sprężystości K, przyspieszenia grawitacyjnego g i masy M. Mając równania ruchu możemy określić częstotliwości dudnień lub częstotliwość modulowaną, a także częstotliwość średnią. Napiszmy teraz równania ruchu obu kulek w sposób:
|
|
Różniczkowe równania ruchu dla obu kulek M, tzn. (1.120) I (1.121) możemy dodać i odjąć od siebie, w ten sposób otrzymać wnioski zależne od zmiennych ψa+ψb (pierwsza równość) i od ψa-ψb (druga równość), by z tych równań wyznaczyć wzory na wspomniane zmienne przy pomocy definicji częstotliwości ω1 i ω2:
|
|
Patrząc na wzory (1.123) i (1.124), w ten sposób otrzymujemy formuły na kwadraty wspomnianych częstotliwości kołowych drgań harmonicznych, których przedstawienie jest zależne od przyspieszenia ziemskiego "g", długości wahadła matematycznego "l" i stałej sprężystości K i masy obu kulek M:
|
|
W celu określenia jak wyglądają dudnienia napiszmy, czemu są równe ψa I ψb dla osobnych drgań harmonicznych, to napiszmy równania na przemieszczenie ciała "a" dla ψa i ciała "b" dla ψb przy definicji częstotliwości kołowych ω1 i ω2, w takim wypadku piszemy je:
![]() |
(1.126)
|
![]() |
(1.127)
|
Najsilniejszy efekt obserwuje się, gdy obie amplitudy są jednakowego znaku, tzn. A1=A2=A. Jeśli któraś z amplitud jest równe zero lub w przybliżeniu równe zero, to mamy do czynienia z drganiami harmonicznymi z drgających z jedną częstotliwością kołową, tzn. ω1 i ω2. Drgania dwóch mas przy jednakowych amplitudach możemy zapisać:
|
|
Wykorzystamy wzory (1.114) i (1.115), by potem je użyć do (1.128) (1.129), by dalej przeprowadzić końcowe obliczenia na drgania obu mas. Dla pierwszej masy mamy według wzoru na przemieszczenia ciała "a", czyli ψa (1.128):
![]() |
(1.129)
|
Dla drugiej masy mamy według wzoru na przemieszczenia ciała "b", czyli ψb (1.129):
![]() |
(1.129)
|
Energię kulki pierwszej określamy w zależności od częstotliwości modulacyjnej ωmod i ωsr, które określamy względem częstotliwości drgań podstawowych, którego piszemy przy pomocy kwadratu funkcji kosinus z iloczynu częstotliwości modulowanej i czasu, i częstotliwości średniej:
![]() |
(1.130)
|
A energię drugiej kulki podobnie jak dla pierwszej kulki określamy przez:
![]() |
(1.131)
|
Określmy teraz wielkość zwaną E, która jest sumą całkowitych energii kulki pierwszej (1.130) i drugiej (1.131), którego ta energia jest wielkością stałą jak można udowodnić, tzn. tylko zależy od częstotliwości średniej obu drgań podstawowych:
![]() |
(1.132)
|
Całkowite energie kulki pierwszej (1.130) i drugiej (1.131), wykorzystując przy tym fakt (1.132), są napisane jako z definicji kwadratu sinusa i kosinusa kąta połówkowego:
|
|
A różnica energii obu kulek, tzn. Ea i Eb na podstawie (1.133) i (1.134) przestawiamy jako funkcję wprost proporcjonalną do kosinusa z argumentu (ω1-ω2)t, i wyrażamy ją:
![]() |
(1.135)
|
[edytuj] Układy fizyczne, a jego drgania o wielu stopniach swobody
Rozważmy układ o N stopniach swobody, której N może być nieskończone. Przy N skończonym układ ma n częstotliwości drgań normalnych, dla całego układu istnieje jednakowe przesunięcie kątowe. Dla przykładu weźmy stosunek amplitud o następujących stosunkach:

To stosunek przesunięć drgań od położenia równowagi względem drgań podstawowych piszemy:

Jeśli układ fizyczny ma dużo elementów, to odległość pomiędzy nimi powinna być bardzo mała, praktycznie można powiedzieć, że ich odległość dąży do zera przy liczbie elementów dążących do nieskończoności, wtedy układ zachowuje się jak układ ciągły. Przesunięcia od położenia równowagi zapisujemy jako funkcja położenia (x,y,z) i czasu "t", którego zapis matematyczny jest ψ(x,y,z,t), która zastępuje przemieszczenia ψa(t), ψb(t), ψc(t),,.., itd. dla osobnych elementów układu. Takie zachowanie się funkcji ψ nazywamy falami.
[edytuj] Drgania poprzeczne struny dyskretnej
Weźmy sobie ciała od 1 do N połączonych sprężynkami o znikomych masach, które spełniają prawo Hooka.
Wraz rosnącymi numerami drgań postaci drgań swobodnych rośnie częstość kołowa drgań, a także rośnie siła kierująca powodująca te drgania. Dla postaci N=2 drgań swobodnych ma jeden węzeł. Na rysunku powyżej zaznaczono układy od jednego stopnia swobody do N-tego stopnia swobody. Dla N stopni swobody mamy N-1 węzłów. Kształt sprężymy ma charakter zygzakowaty dla największej postaci drgań, dla której każdy odcinek łączący masy przecina oś równowagi.
[edytuj] Drgania poprzeczne struny ciągłej
Załóżmy, że mamy N ciężarków połączonych ze sobą sprężynkami, których liczba N jest bardzo duża, powiedźmy, że wynosi N=1000000, których każdy z ciężarków wykonują drgania ψa, φb, φc, itd.. , wtedy układ zachowuje się prawie jak ciągły. Zatem całkowite wychylenie od stanu równowagi w zależności od położenia ciężarków, a także od czasu, jest wyrażone:
![]() |
(2.1)
|
[edytuj] Drgania podłużne i poprzeczne dla płaszczyźnie
Wybierzmy sobie drgania jednomiarowej struny dla drgań poprzecznych, której przemieszczenie dla naszego przypadku określamy:
![]() |
(2.2)
|
Drgania podłużne są wykonywane wzdłuż osi zetowej, a nas interesują drgania poprzeczne, co na podstawie wzoru (2.2) możemy przepisać wzór na postać tychże drgań:
![]() |
(2.3)
|
[edytuj] Polaryzacja struny ciągłej
Dokonajmy dalszych uproszczeń, i powiedźmy, że drania są wykonywane wzdłuż osi iksowej. Gęstość takiej struny jest równa ρ. A masa na odcinku Δx wyrażamy poprzez wzór:
![]() |
(2.4)
|
Całkowita siła działająca na wycinek struny przy pomocy siły kierującej T wyrażamy przez:
![]() |
(2.5)
|
Chcemy wyrazić siłę Fx poprzez pochodne cząstkową wychylenia od stanu równowagi względem położenia zetowego małego odcinka struny przy działających siłach na jej prawy i lewy koniec argumenty, który zapis nachylenia w tymże punkcie jest:
, który jest równy tangensowi nachylenia struny względem osi z dla chwili t dla wycinka struny o położeniu z. Siłę (2.5) możemy przestawić przy pomocy wzoru:
![]() ![]() |
(2.6)
|
Zdefiniujmy funkcję f(z) zależną od z jako pochodną cząstkową wielkości ψ(z,t) względem położenia z, którego definicja:
![]() |
(2.7)
|
Rozłóżmy funkcję (2.7) względem położenia wokół punktu z1, wykorzystując szereg Taylora, by potem wyciąć z niego wyrazy więcej niż liniowe:
![]() |
(2.8)
|
Biorąc Δz=z2-z1, możemy napisać różnice funkcji f(z) w puntach z2 i z1, którego różnica jest wprost proporcjonalna do drugiej pochodnej funkcji ψ względem położenia przy stałej proporcjonalności równej Δz:
![]() |
(2.9)
|
Wniosek wypływający z tożsamości (2.9) i wzoru (2.6) na definicję siły Fx jest taki:
![]() |
(2.10)
|
Zastosujmy teraz drugie prawo Newtona dla wycinka drgającej struny przy definicji siły Fx (2.10), który po skróceniu przez Δz otrzymujemy tożsamość robiąc te czynności pokolei:
![]() |
(2.11)
|
[edytuj] Fale stojące drgającej struny
Załóżmy, że struna drga z jedną częstością kołową, wtedy mamy do czynienia z drganiami harmonicznymi normalnymi. Równanie wychylenia struny od położenia równowagi określamy przez:
![]() |
(2.12)
|
Druga pochodna cząstkowa wyrażenia (2.12) względem czasu określamy jako wielkość wprost proporcjonalną do kwadratu częstotliwości kołowej drgań i do funkcji A(z):
![]() |
(2.13)
|
A druga pochodna funkcji Ψ(z,t) (2.12) względem położenia z wycinka struny ciągłej piszemy przez wzór poniżej, która jest zależna od drugiej pochodnej zupełnej wielkości A(z) względem czasu:
![]() |
(2.14)
|
Biorąc równanie różniczkowe falowe (2.11) i po podstawieniu do niego wzoru (2.13) i (2.14) i skróceniu przez kosinus, którego argumentami są te same określone wyrażenia, co uzasadnia takie działanie:
![]() |
(2.15)
|
Rozwiązując równanie (2.15) dla oscylatora harmonicznego, których ogólna postać tych oscylacji jest kombinacją funkcji sinus i cosinus względem wyrazenia 2πz/λ, wtedy:
![]() |
(2.16)
|
A także możemy napisać wniosek wynikający z (2.16), którą to możemy napisać jako druga pochodna zupełna względem położenia "z", która jest wprost proporcjonalna do funkcji A(z) i do odwrotności kwadratu długości fali kreślonej i to wszystko wziętej razem z minusem:
![]() |
(2.17)
|
Możemy porównać obie strony wzorów (2.15) z (2.17), a właściwie jej prawe strony, a lewe strony są takie same, do siebie, w ten sposób mamy tożsamość na prędkość fali, który jest pierwiastkiem z ilorazu siły kierującej T0 i gęstości liniowej masy dla tej opisywanej tutaj struny:
![]() |
(2.18)
|
[edytuj] Fale stojące
Całkowita funkcja falowa opisujące fale jest funkcją przedstawiona na podstawie wzorów (2.12) i (2.16), jest tutaj przedstawiona w jego pełnej postaci dla fali stojącej:
![]() |
(2.19)
|
[edytuj] Warunki brzegowe
Rozwiązanie na funkcje stojące (2.19) nie zawiera w sobie warunków brzegowych, bo ono jest ogólnym rozwiązaniem, w taki sposób, by dla z=0 był jeden koniec strony a drugi koniec znajduje się dla z=L, wtedy drugi współczynnik zetowy stojący w nawiasie (2.19) jest równy zero dla zerowego wychylenia lewego końca, bo ten koniec jest zaczepiony, wtedy:
![]() |
(2.20)
|
Patrząc na tożsamość (2.20) dowiadujemy się, że stała B jest zawsze równa zero, co było udowodnienia. Równanie fali stojącej (2.19) możemy zapisać, na podstawie dowodu (2.20), że stała B=0, jako:
![]() |
(2.21)
|
Aby dla prawego końca struny wychylenie dla fali stojącej (2.20) jest równe zero, czyli dla z=L, wtedy funkcja sinus powinna być równa zero, co matematycznie:
![]() |
(2.22)
|
Wykorzystując z własności funkcji sinus we równaniu (2.22), w ten sposób otrzymujemy, że L powinno być całkowitą wielokrotnością z liczby λ/2, co można udowodnić przeprowadzając odpowiednie obliczenia:
![]() |
(2.23)
|
W obliczeniach (2.23) pominęliśmy przypadek k=0, ponieważ wymagało to by zerową długość struny, lub gdy długość fali drgań była nieskończona. Dla k=1 warunek brzegowy wygląda
, a z warunku (2.23) i z definicji na stałą λ1 dalsze wzory na długość fali λk w zależności od k>1 zapisujemy:
![]() |
(2.24)
|
[edytuj] Stosunki częstotliwości kołowych w drganiach harmonicznych
Częstotliwości ν poszczególnych drgań możemy napisać w zależności częstotliwości podstawowej ν1 względem (2.24) w postaci wzorów:
![]() |
(2.25)
|
[edytuj] Kątowa liczba falowa
Liczbę falowa nazywamy liczbę falową σ, która jest odwrotnością długości fali, a kątową liczbę falową nazywamy obiekt, która powstaje po pomnożeniu jego przez liczbę 2π, wtedy;
|
|
Równanie fali (2.20) możemy przestawić przy pomocy (2.26) lub (2.27) na kilka równoważnych sposobów:
![]() ![]() = Asin ωtsin kz |
(2.28)
|
Patrząc na wzór (2.28), który przestawia drgania struny z obiema końcami nieruchomymi, tzn. dla z=0 dla naszej funkcji warunek nieruchomego lewego końca jest automatycznie spełnione, a także dla z=L warunek prawego nieruchomego końca jest spełnione, gdy kątowa liczba falowa drgań spełnia warunki:
![]() |
(2.29)
|
lub w przestawieniu liczby falowej σ:
![]() |
(2.30)
|
[edytuj] Związek pomiędzy długością fali a jej częstotliwością
Patrząc na wzór (2.17) otrzymujemy wzór na związek dyspersyjny liczby falowej i częstotliwości ν:
![]() |
(2.31)
|
Równość (2.31) mnożymy obustronnie przez liczbę 2π i po wykorzystaniu definicji częstotliwości kołowej i kątowej liczby falowej:
![]() |
(2.32)
|
[edytuj] Struny fortepianowe i jego związek dyspersyjny
Związek dyspersyjny dla struny fortepianowej, czyli związku
w zależności od liczby kątowej przestawiamy:
![]() |
(2.33)
|
Związek (2.33) jest bardzo podobny do związku (3.22), tylko jest ta różnica, że jeśli struna jest idealnie giętka, to wtedy stała α jest w rzeczywistości jest równa zero, ale struna nie jest taka tak naprawdę, tzn. dla której w ogólności dla struny nie idealnie giętkiej nie zachodzą już zależności λ1=2L, λ2=1/2λ1, λ3=1/3λ1, itd. Częstości dla naszej struny nie tworzą tutaj ciągu harmonicznego ν2=2ν1, ν3=3ν1, itd., ponieważ powyższy związek dyspersyjny nie daje takiej zależności, ten ciąg otrzymujemy w granicy α dążącego do zera, tzn., wtedy gdy zachodzi związek λν=const.
[edytuj] Analiza Fouriera ogólnych drgań struny ciągłej
Ogólne wychylenie wycinka ciągłej strony przestawiamy jako sumę drgań prostych harmonicznych:
![]() |
(2.34)
|
[edytuj] Analiza Fouriera struny drgań jako złożenie drgań harmonicznych
Każdą funkcję F(z) możemy rozłożyć w szereg Fouriera względem zmiennej znając współczynniki An i Bn naszego szeregu dla naszej funkcji, której postać tego szeregu jest:
![]() ![]() |
(2.35)
|
Następnym krokiem jest wyznaczenie współczynników Fouriera B0 i Bn, a także An, które wyznaczymy pamiętając przy tym, że zachodzi z2=z1+λ:
|
|
![]() |
(2.38)
|
![]() |
(2.39)
|
Funkcję F(z) możemy obustronnie przecałkować względem zmiennej w przedziale (z1,z2), wykorzystując przy tym związek (2.36) oraz (2.37), w ten sposób otrzymać związek zawierający w sobie stałą B0:
![]() |
(2.40)
|
Wykorzystując fakty (2.38) i (2.39), a także fakty w postaci całek (2.36) i (2.37), wtedy możemy napisać całki, napisane w przedziale (z1,z2), w których funkcją podcałkową jest iloczyn funkcji F(z) i funkcji sinus lub kosinus z argumentu mk1z, całkowalną względem zmiennej "z", którego to przepisy w postaci już policzonych całek zapisujemy jako:
|
|
[edytuj] Współczynniki Fouriera
Patrząc na wzory (2.40) i (2.41) oraz (2.42), wtedy możemy napisać wzory na współczynniki szeregu Fouriera (2.35) dla funkcji F(z), tzn. wzory na współczynniki Am, Bn:
![]() |
(2.43)
|
|
|
[edytuj] Fala prostokątna i jego analiza Fouriera
Rozpatrzmy teraz falę prostokątną, która jest falą periodyczną i patrząc na rysunek obok dowiadujemy się, że ona jest nieciągła dla z=0 i z=L. Wiec nie należy oczekiwać, że szereg Fouriera w pewnym przybliżeniu przestawi w miarę dokładnie falę prostokątna. Korzystając z (2.41) dowiadujemy się, że współczynniki Bm są równe zero, a współczynniki Am istnieją tylko dla nieparzystych m. Aby to wykazać policzymy kolejno współczynniki na nasze współczynniki szeregu Fouriera, zatem wyznaczmy czemu jest równy na współczynnik B0, którego definicja jest w punkcie (2.41):
![]() |
(2.46)
|
i dalej policzymy współczynniki Bm, którego definicja jest w punkcie (2.43):
![]() ![]() ![]() ![]() ![]() ![]() |
(2.47)
|
i dalej policzymy współczynniki Am, którego definicja jest w punkcie (2.44):
![]() ![]() ![]() ![]() ![]() |
(2.48)
|
Jak widzimy, że współczynniki Am (2.48), B0 (2.46), Bm (2.47) są takie jak przypuszczaliśmy. Analiza Fouriera pokazuje, że impuls prostokątny możemy rozłożyć na fale harmoniczne, wykorzystując przy tym wzór na współczynnik B0 (2.46), na współczynnik Bm (2.47) i na samym końcu na współczynnik Am (2.48), wtedy zapisujemy nasz impuls w postaci szeregu Fouriera:
![]() ![]() ![]() |
(2.49)
|
[edytuj] Fourierowska analiza okresowej funkcji zależnej od czasu
Funkcje okresowa zależna od czasu spełnia taki warunek, że dodanie do funkcji dla argumentu t1 okresu T1, otrzymujemy tą samą funkcję bez tej dokonanej operacji, co matematycznie zapisujemy:
![]() |
(2.50)
|
Funkcję zależną od czasu możemy rozłożyć w szereg Fouriera względem czasu "t" do postaci:
![]() |
(2.51)
|
Mając funkcję okresową F(t) możemy rozłożyć tą funkcję w szereg Fouriera mając definicję współczynników Fouriera o postaciach:
![]() |
(2.52)
|
|
|
[edytuj] Drgania struny niejednorodnej
Zakładamy, że siła T0 działa w stanie równowagi, która jest zależna od zmiennej zetowej, którą piszemy wzorem (2.6). Obierzmy sobie funkcję f(z), którego postać jest w zależności od tej siły i pierwszej pochodnej wychylenia od stanu równowagi dla drgań poprzecznych:
![]() |
(2.55)
|
Funkcję możemy rozłożyć w szereg Taylora (2.8) funkcję f(z) (2.55), w ten sposób otrzymujemy, że siła Fx działająca na na bardzo mały wycinek struny niejednorodnej jest:
![]() |
(2.56)
|
Drugie prawo dynamiki Newtona dla naszego przypadku niejednorodnej struny, który piszemy wzorem poniżej, którą dzielimy przez długość małego wycinka struny Δz i przez gęstość tegoż wycinka ρ(z), wtedy:
![]() |
(2.57)
|
Niech fala będzie opisywana wzorem zależnym od amplitudy A(z) i części periodycznej cos(ωt+φ), piszemy:
![]() |
(2.58)
|
Dla drgań normalnych struny możemy napisać zależności jako drugą pochodną przemieszczenie względem czasu i pierwszą pochodną tego samego przemieszczenia względem położenia "z", która jak się okazuje jest wprost proporcjonalna do pierwszej pochodnej zupełnej funkcji A(z) względem położenia małego wycinka struny względem położenia, zatem te wielkości:
|
|
Mając fakty (2.59) i (2.60), co one możemy podstawić do wzoru (2.57), w ten sposób otrzymać równość różniczkową z pochodnymi zupełnymi zależnych od współrzędnej zetowej:
![]() |
(2.61)
|
[edytuj] Drgania układu dyskretnego o N stopniach swobody
Rozpatrzmy układ w stanie równowagi składający się z N kulek o masach M połączonych sprężynkami napiętej w taki sposób, w których każda taka sprężynka działa na każdą kulkę z osobną z siłą T0, układ składający się z sprężynek o pewnej stałej sprężystości i kulek o masach M ma długość L=(N+1)a.
[edytuj] Równanie ruchu struny dyskretnej
Wzór na siłę działająca na układ dyskretnym piszemy podobnie jak według wzoru (2.6), tylko tam zamiast pochodnej występują różnice, zatem równanie ruchu w naszym przypadku piszemy:
![]() |
(2.62)
|
Obierzmy sobie ogólną postać wychylenia n-tej cząstki od położenia równowagi w postaci:
![]() |
(2.63)
|
Druga pochodna zupełna wyrażenia (2.63) względem czasu, która jest wprost proporcjonalna do amplitudy An i do kwadratu częstotliwości kołowej, i dalej mówiąc od funkcji kosinus, to wszystko wziętej z minusem, piszemy w postaci:
![]() |
(2.64)
|
Równanie (2.62) możemy w taki sposób napisać po wykorzystaniu (2.63) i (2.64), które to piszemy w zależności częstotliwości kołowej drgań, a także od amplitud zależnych od n, n+1 i n-1, który to końcowy rozważany wzór piszemy w postaci:
![]() |
(2.65)
|
Rozpatrzmy tak jak w przypadku struny ciągłej, którego to dla struny na obu jego końcach naszej struny nie ma wychylenia od stanu równowagi, zatem amplitudę An możemy napisać w postaci:
![]() |
(2.66)
|
Zatem dzięki definicji amplitudy An według (2.66) wówczas możemy napisać następujące wzory na amplitudy An+1 i An-1, co na samym końcu możemy policzyć różnicę amplitud dla n+1 i n-1, wtedy według wspomnianego wzoru możemy napisać trzy dalsze wzory:
![]() |
(2.67)
|
![]() |
(2.68)
|
![]() |
(2.69)
|
Możemy wykorzystać końcowy związek (2.69) do (2.65), w ten sposób otrzymać wzór w ogólności dla niezerowego w An, który jest zależny od częstotliwości kołowej drgań i innych parametrów:
![]() |
(2.70)
|
We związku (2.70) możemy wyznaczyć kwadrat częstotliwości kołowej i w ten sposób otrzymać zależność od kątowej liczby falowej "k" lub też w zależności do liczby falowej λ:
![]() |
(2.71)
|
[edytuj] Ogólna postać drgań i warunki brzegowe
Ogólna postać drgań An możemy przestawić w postaci wzoru jako kombinację funkcji sinus i kosinus z argumentu kna:
![]() |
(2.72)
|
Na końcu lewym wychylenie od położenia równowagi jest zerowe i wykorzystując fakt (2.72) dowiadujemy się, że B=0, w ten sposób otrzymujemy równość:
![]() |
(2.73)
|
Aby prawy koniec struny dyskretnej miał wychylenie dyskretne zerowe, to powinien zachodzić związek:
![]() |
(2.74)
|
Na rysunku obok zaznaczono zależność częstotliwości ω(k) od kątowej liczby falowej, którą tą zależność piszemy na podstawie (2.71) w postaci:
![]() |
(2.75)
|
[edytuj] Przejścia układu dyskretnego do granicy fal długich
Weźmy sobie fale o długości fali o wiele mniejszych niż zetowa odległość pomiędzy ciężarkami, tzn. a<<λ, co jest równoważne 2πa/λ<<1, wtedy możemy rozłożyć funkcję sinus w szereg Taylora w postaci:
![]() |
(2.76)
|
Do wzoru (2.75) możemy wykorzystać rozwinięcie funkcji sin (2.76), by w ten sposób napisać:
![]() |
(2.77)
|
[edytuj] Drgania podłużne układu mas połączonych ze sobą sprężynkami
Równanie ruchu dla masy ściśle określonej zapisujemy względem wychylenia od stanu równowagi poszczególnych sprężynek ψn, którego to przepis:
![]() |
(2.78)
|
równanie tego typu (2.78) było już rozwiązywane, zatem podajmy związek dyspersyjny łączący częstotliwość kołową z kątową liczbą falową, którego postać:
![]() |
(2.79)
|
Wychylenia poszczególnych mas dla drgań podłużnych możemy zapisać pod postacią funkcji:
![]() |
(2.80)
|
[edytuj] Układ LC wzajemnie sprzężonych ze sobą cewek i kondensatorów
Patrząc na rysunek b) obok możemy napisać drugie prawo Kirchhoffa w postaci:
![]() |
(2.81)
|
Równanie (2.81) możemy zróżniczkować obustronnie względem czasu i wykorzystywać będziemy pierwsze prawo Kirchhoffa:
![]() |
(2.82)
|
Co (2.82) możemy zapisać równoważnie w postaci:
![]() |
(2.82)
|
Równanie (2.82) jest równaniem, którego typ poznaliśmy wcześniej, którego związek dyspersyjny częstotliwości kołowej względem kątowej liczby falowej określamy:
![]() |
(2.83)
|
Równanie na natężenie prądu elektrycznego płynącego w obwodzie LC w poszczególnym oczku w cewce określamy przez:
![]() |
(2.84)
|
[edytuj] Wahadła matematyczne sprzężone ze sobą za pomocą sprężynek
Równanie ruchu n-tej masy dla układu sprzężonych ze sobą wahadeł matematycznych za pomocą sprężynek opisujemy:
![]() |
(2.85)
|
Równanie (2.85) jest równaniem, którego typ jest już znany, czyli stąd możemy wyliczyć częstotliwość kołową drgań w zależności od liczby falowej:
![]() |
(2.86)
|
Jeśli w prowadzimy częstotliwość drgań wahadła matematycznego, którego kwadrat jest napisany w punkcie (1.6), i je oznaczymy przez
, wtedy nasz związek dyspersyjny (2.86) piszemy w postaci:
![]() |
(2.87)
|
Widzimy, że związek (2.87) zależy od częstotliwości kołowej drgań podstawowych wahadła matematycznego ω0, a także zależy od stałej sprężystości, masy kulek, a także od kątowej liczby falowej "k". Jeśli przyjmować będziemy
, wtedy związek dyspersyjny (2.87), przy założeniu ka<<1, zapisujemy w postaci przybliżonej:
![]() |
(2.88)
|
[edytuj] Związek dyspersyjny dla plazmy
Wyprowadzimy tutaj związek dyspersyjny dotyczącej plazmy, który możemy przestawić jak się przekonamy od częstotliwości podstawowej ωp, a także od długości światła. Wiemy, że kwadrat częstotliwości kołowej przestawiamy w postaci wzoru:
![]() |
(2.89)
|
Widzimy, że pozostało na wyliczyć kwadrat częstotliwości podstawowej, tj.
, który wyprowadzimy w dalszych rozważaniach. Elektryczna obojętna plazma składa się z większości cząsteczek obojętnych i pewnej liczby cząsteczek zjonizowanych. Jonosfera ziemi składa się się z dużej ilości cząsteczek obojętnych N2 i O2. Proces jonizacji tego gazu powodowany jest poprzez pochłoniecie promieniowania pochodzącego od słońca. Największą gęstość jonów zjonizowanych występuje przy odległości 200km do 400km nad powierzchnią ziemi. Przy niższych wysokościach proces jonizacji nie zachodzi ze względu na pochłonięcie promieniowania przez wyższe warstwy atmosfery, i dlatego na tej wysokości gęstość jonów i elektronów maleje do zera. Plazma jest elektrycznie obojętna i nie stanowi źródła pola elektrycznego, ale istnieją w niej warstwy zjonizowane na wskutek działania pola elektrycznego zewnętrznego, że dodatnie jony są przyspieszane w jedną stronę, a elektrony w drugą stronę. W skutek nadmiaru ładunku zostaje zlikwidowane do zera ładunek w wyniku ich przyciągania do siebie, czyli jony i elektrony. A bezwładność odciąga je od stanu równowagi powodując nowy nadmiar ładunku, a także niedobór ładunków. Całkowite natężenie pola elektrycznego wyznaczamy z prawa Gausa, pamiętając przy okazji ze pole na zewnątrz takiego układu jest równe zero, przestawiamy je w zależności od powierzchni rozważanych okładek A:
![]() |
(2.90)
|
Całkowity ładunek znajdujących się w grubości x plazmy przestawiamy w zależności od koncentracji jonów i grubości warstwy przy pomocy wzoru:
![]() |
(2.91)
|
Równanie ruchu plazmy możemy przestawić w zależności od pola działającego ze strony pola w wewnątrz tej plazmy według mechaniki klasycznej:
![]() |
(2.92)
|
Wzór na natężenie pola iksowego elektrycznego Ex (2.90) i wzór na łądunek znajdujących się w plazmie (2.91) podstawiamy do wzoru (2.92) opisującą drugą zasadę dynamiki Newtona:
![]() |
(2.93)
|
Końcowe równanie (2.92) jest równaniem oscylatora harmonicznego, w której kwadrat częstotliwości kołowej drgań, która jest jednocześnie kwadratem częstotliwości podstawowej ωp, który to przestawiamy w zależności od koncentracji jonów i ich masy, jest:
![]() |
(2.94)
|
[edytuj] Drgania wymuszone układów harmonicznych
[edytuj] Oscylator tłumiony w przedstawieniu jednowymiarowym
Wyobraźmy sobie oscylator harmoniczny, który doznaje siły kierującej
pochodzącej od sprężyny o stałej sprężystości
, dla której układ doznaje siły oporu
, którym występuje stała Γ, która jest współczynnikiem tłumienia na jednostkę masy, lub po prostu jest zwana współczynnikiem tłumienia. Na układ też działa siły wymuszająca układ do ruchu F(t). Biorąc wszystkie te uwagi możemy napisać z drugiej zasady dynamiki Newtona równość różniczkową:
![]() |
(3.1)
|
[edytuj] Tłumienie oscylacji harmonicznych bez udziału siły wymuszającej
Równanie (3.1) można napisać przy założeniu, że nie mamy siły wymuszającej F(t), wtedy tak otrzymane równanie możemy podzielić przez masę M, i wszystkie wyrazy przenosząc na jedną stronę, w takim wypadku:
![]() |
(3.2)
|
Rozwiązaniem szczególnym równania jednorodnego liniowego (3.2) jest rozwiązaniem oscylatora tłumionego zależne od stałej τ i częstotliwości kołowej ω:
![]() |
(3.3)
|
Dalszym krokiem jest policzenie pierwszej i drugiej pochodnej zupełnej położenia względem czasu masy M oscylatora tłumionego (3.3), jego pierwsza pochodna jest:
![]() |
(3.4)
|
Wyznaczmy teraz drugą pochodną zupełnej wyrażenia (3.3) względem czasu, czyli pierwszą pochodną pierwszej pochodnej wyrażenia (3.4) względem czasu, zatem do dzieła:
![]() ![]() |
(3.5)
|
Szczególne rozwiązanie (3.3), pierwszą pochodną funkcji x(t) (3.4), a także drugą pochodną funkcji x(t) (3.5) podstawiamy do równania różniczkowego (3.2):
![]() ![]() |
(3.6)
|
Tożsamość (3.6) dzielimy obustronnie przez
, bo ono jest zawsze niezerowe, w ten sposób otrzymujemy tożsamość, w którym grupujemy wyrazy względem sinωt i cosωt:
![]() |
(3.7)
|
Równość (3.7) jest słuszna dla każdego czasu, stąd wynika, że czynniki stojące przy sin(ωt+θ) i cos(ωt+θ) powinny być zawsze równe zero, w ten sposób dostajemy układ dwóch równań:
![]() |
(3.8)
|
Na podstawie układu równań (3.8) otrzymujemy następujące tożsamości na τ i ω2:
|
|
Rozwiązanie (3.3) równania różniczkowego jednorodnego (3.2) możemy przestawić w postaci równoważnej przy pomocy funkcji funkcji trygonometrycznych sinus i kosinus jako ich kombinację liniową pomnożoną przez funkcję eksponencjalną malejącą wraz z czasem do zero dla "t" nieskończonego:
![]() |
(3.11)
|
Możemy wyznaczyć pierwszą pochodną wyrażenia (3.11), w ten sposób mamy wzór na prędkość ciała w chwili "t":
![]() |
(3.12)
|
Wyznaczmy warunki brzegowe dla t=0, tj. dla równania (3.11), który jest położeniem masy w czasie "t", dla równania (3.12), który jest prędkością ciała o masie M w czasie "t", by potem wyznaczyć stałe A i B w zależności od warunków brzegowych:
|
|
Wykorzystując fakty, tzn. (zależność na położenie ciała w czasie t=0) (3.13) i (zależność prędkości ciała w chwili zerowej)(3.14), w ten sposób ogólne rozwiązanie oscylatora harmonicznego tłumionego (3.11) piszemy w postaci:
![]() |
(3.15)
|
Jeśli natomiast
, i wiedząc coś o liczbach zespolonych, które są wprowadzone w algebrze, możemy powiedzieć na podstawie (3.10):
|
|
Zatem równość (3.15) możemy napisać na podstawie tożsamości (3.16) i (3.17), w którym ω jest liczbą zespoloną, a także wykorzystując fakty o związkach pomiędzy funkcjami trygonometrycznymi a funkcjami hiperbolicznymi:
![]() |
(3.18)
|
Wyznaczmy teraz energię całkowitą zmieniającą się wraz z czasem, która jest sumą energii kinetycznej i potencjalnej sprężystości, którego zapis jest przy wykorzystaniu wzoru (3.3), który jest wzorem na położenie ciała w chwili "t" i biorąc przykład słabego tłumienia, wtedy wyraz z Γ jako czynnik w (3.12) możemy pominąć, bo ono jest bardzo małe, zatem w takim wypadku średnia energia całkowita jest:
![]() ![]() ![]() |
(3.19)
|
[edytuj] Drgania stacjonarne podczas działania siły harmonicznej siły wymuszającej
Każdą siłę wymuszającą możemy rozłożyć w szereg Fouriera, których każdy wyraz jest funkcją kosinus o przesunięciu fazowym φ(ω) o amplitudzie F(ω), która ta całkowita siła jest przestawiona:
![]() |
(3.20)
|
W szeregu (3.20) weźmy tylko jeden wyraz, tak by funkcja kosinus nie zawierała wcale przesunięcia fazowego, wtedy tą siłę przestawiamy:
![]() |
(3.21)
|
Całkowite równanie ruchu masy M dla drgań harmonicznych, gdy istnieją w nim straty energii spowodowane niezerowym współczynnikiem tłumienia, na którą to masę działa siła harmoniczna drgająca z częstością ω, przestawiamy:
![]() |
(3.22)
|
Drgania wykonywane przez oscylator opisywany przez równanie (3.22) jest narysowane jako kombinacja funkcji sinus i kosinus:
![]() |
(3.23)
|
Wzór (3.23) podstawiamy do równania różniczkowego (3.22) otrzymując tożsamość:
![]() |
(3.24)
|
Równanie (3.24) możemy pogrupować względem funkcji sinus i kosinus, w ten sposób otrzymać dalsze przekształcenie tej równości:
![]() |
(3.25)
|
Równanie (3.25) jest spełnione dla każdego czasu "t", wtedy współczynniki stojące przy kosinusach i sinusach są równe zero, wtedy mamy układ równań, które będziemy rozwiązywać ze względu na współczynniki A i B:
![]() |
(3.26)
|
Z pierwszego równania (3.26) możemy wyznaczyć stałą B zależną od stałej A, częstotliwości podstawowej ω0, częstotliwości drgań siły wymuszającej i współczynnika tłumienia Γ, którego równanie jest:
![]() |
(3.27)
|
Tożsamość (3.27) możemy podstawić do drugiego równania układu równań (3.26), z którego w dalszych rozumowaniach otrzymamy wzór na stałą A:
![]() |
(3.28)
|
Z tożsamości (3.28) możemy wyznaczyć stałą A, którego wzór jest napisany przy pomocy częstotliwości kołowej drgań podstawowych ω0 i drgań wymuszonych ω:
![]() |
(3.29)
|
Równanie przestawione wzorem (3.29) jest to wzór na amplitudę absorpcyjną . Wzór na A (3.29) podstawiamy do (3.27), w ten sposób otrzymujemy równość:
![]() |
(3.30)
|
Równanie na stałą B (3.30) nazywamy równaniem na amplitudę elastyczną . Wkład energii spowodowanej przez amplitudę elastyczną po uśrednianiu w ciągu całego okresu jest równy zero, bo straty energii są równe sile F0cosωt pomnożonej przez
i dlatego ta amplituda w ogólności nie powoduje strat energii. Można udowodnić, że całkowita strata energii spowodowana przez amplitudę absorpcyjną spowodowaną po uśrednianiu jego w ciągu jednego okresu jest równa:
![]() |
(3.31)
|
A dowód jego polega na uśrednianiu strat energii w czasie t w ciągu jednego okresu, jak widzimy poniżej straty energii są zależne od amplitudy absorpcyjnej, częstotliwości drgań siły wymuszającej i na samym końcu od amplitudy F0 drgań harmonicznej siły zewnętrznej działającej na nasz układ drgający:
![]() ![]() |
(3.32)
|
Całkowita średnia energia tracona uśrednianą w ciągu jednego okresu jest wyrażona przy pomocy amplitudy elastycznej Ael, absorpcyjnej Aab, współczynnika Γ i masy drgającej M:
![]() |
(3.33)
|
Można udowodnić, że wzór (3.33) za pomocą prostych przekształceń, jest równy wzorowi (3.31). Całkowita energia zmagazynowana w oscylatorze tłumionym podtrzymywanym siłą wymuszającą jest wyrażona po uśrednianiu jej w ciągu jednego okresu:
![]() |
(3.34)
|
[edytuj] Rezonans
Wzór (3.31) możemy przepisać w postaci równania na moc traconą podczas oporów ruchów, która jest tak skonstruowana, by dla ω=ω0, by było P(ω0)=P0:
![]() |
(3.35)
|
Stała proporcjonalności P0 została tak dobrana, by dla ω=ω0, całkowita moc promieniowania była równa P0. Wartość P równej połowie maksymalnej rezonansowej jest napisana w takiej formie, by było spełnione:
![]() |
(3.36)
|
Z równości podanej w punkcie (3.36) możemy wyznaczyć częstotliwość ω w zależności od częstotliwości podstawowej ω0 i współczynnika tłumienia, którego postać naszej tożsamości na tą wspomnianą częstotliwość jest:
![]() |
(3.37)
|
Ponieważ częstotliwość kołowa drgań kołowych jest wielkością nieujemną, wtedy wyrażenie (3.37) możemy tak napisać, w którym drugi składnik jest połową współczynnika tłumienia wziętej z minusem lub plusem:
![]() |
(3.38)
|
Wyrażenie z plusem (3.38) możemy odjąć od wyrażenia wziętego z minusem, w ten sposób otrzymujemy tożsamość na szerokość połówkową:
![]() |
(3.39)
|
[edytuj] Amplituda elastyczna i jego zależność od częstości
Amplituda elastyczna jest dana wzorem (3.29), a jego stosunek do amplitudy absorpcyjnej piszemy:
![]() |
(3.40)
|
Gdy częstotliwość siły wymuszającej jest o wiele mniejsza od drgań podstawowych układu, wtedy absorpcja energii jest zbyt mała w porównaniu z absorpcją energii w rezonansie przez układ. Dla drgającego układu daleko od rezonansu masy możemy napisać tożsamość, która jest spełniona w przybliżeniu:
![]() |
(3.41)
|
[edytuj] Krzywe rezonansowe
Dla oscylatora harmonicznego drgania wymuszone możemy opisywać za pomocą kilku wielkości zależne od częstotliwości kołowych, ale napiszmy najpierw wzór na amplitudę absorpcyjną:
![]() |
(3.42)
|
Teraz wyznaczmy kwadrat modułu amplitudy, która jest sumą kwadratów amplitudy elastycznej (3.30) i absorpcyjnej (3.29), wtedy:
![]() ![]() |
(3.43)
|
Energię oscylacji określamy według wzoru (3.31), która jest definicją mocy traconej przez układ, znając definicję amplitudy absorpcyjnej (3.30), piszemy ją:
![]() ![]() |
(3.44)
|
Całkowita energia zmagazynowana przez układ określamy jako sumę energii kinetycznej i potencjalnej posiadanej przez układ, którą piszemy w postaci wzoru (3.34), do której wykorzystamy wzór (3.43), w ten sposób otrzymamy końcową tożsamość opisującą tą naszą wielkość:
![]() |
(3.45)
|
W powyższych wielkościach mianownik jest określony poprzez częstotliwość drgań podstawowych i częstotliwość drgań wymuszonych:
![]() |
(3.46)
|
Definicja D napiszmy, gdy częstotliwość ω jest bliska ω0, dla stanu słabego tłumienia, tzn. dla którego zachodzi Γ<<ω0, przy pomocy tychże dysput sumę ω0+ω możemy przestawić w przybliżeniu w postaci 2ω0, wtedy D opisujemy wzorem przybliżonym na podstawie jej ścisłej definicji (3.46):
![]() |
(3.47)
|
Wtedy możemy wprowadzić współczynnik R(ω), który jest tak dobrany by dla ω0 był równy jeden, czyli zachodziło R(ω0)=1:
![]() |
(3.48)
|
Funkcja (3.48) jest funkcją parzystą ze względu na ω0-ω, następnie udowodnimy, aby wspomniana funkcja była równa połowie jedynki, czyli wtedy powinno zachodzić
, co z tego możemy wywnioskować
, co dalej można łatwo udowodnić, że szerokość połówkowa naszej funkcji jest równa Γ. Funkcja R(ω) w optyce jest nazywana często rozkładem lorentzowskim, a w fizyce jądrowej jest nazywana krzywą rezonansową Breita-Wignera, dla którego zamiast ω i ω0 są zastępowane tam kolejno przez
i
.
[edytuj] Niestacjonarne drgania harmoniczne tłumione
Niestacjonarne drgania harmoniczne, która jest rozwiązaniem równania (3.23), która jest szczególnym rozwiązaniem równania różniczkowego (3.22), a dla równania jednorodnego (3.2), którego rozwiązaniem jest (3.11), to suma tych rozwiązań tutaj rozważanych piszemy:
![]() |
(3.49)
|
Częstotliwość ω1 we wzorze (3.49) jest to wielkość opisana wzorem (3.10).
[edytuj] Układ drgający niestacjonarny spoczywający w stanie początkowym
Aby układ w stanie początkowym t=0 miał położenie zerowe, wtedy musi zachodzić B1=Ael. Dalej będziemy mieli na myśli słabe tłumienie, tzn. czynnik exp(-Γ t) praktycznie się wcale nie zmienia w chwili w ciągu jednego okresu, zatem dla chwili początkowej mamy
, ponieważ interesują nas częstotliwości rezonansowe niezbyt odległe od częstotliwości podstawowej, czyli kładziemy A1=-Aab, czyli powinno chodzić:
![]() |
(3.50)
|
co prędkość drgań początkowych jest równa zero, gdy zachodzi w przybliżeniu ω=ω1, lub gdy Aab=0, co implikuje Γ=0. Na podstawie tychże omawianych wniosków dostajemy:
![]() |
(3.51)
|
Gdy częstość wymuszająca drgania jest równa częstotliwości drgań tłumionych, czyli ω=ω1, wtedy na podstawie (3.49) możemy powiedzieć, że zachodzi własność:
![]() |
(3.52)
|
[edytuj] Oscylator niestacjonarny przy braku tłumienia
Tutaj kładziemy Γ=0, wtedy amplituda absorpcyjna Aab jest równa zero, wtedy wzór na amplitudę elastyczną (3.29) jest wyrażona przy pomocy częstotliwości drgań podstawowych ω0 i częstotliwości drgań wymuszonych ω:
![]() |
(3.53)
|
Wtedy położenie masy M (3.51), przy powyżej kładzionych warunkach i wzoru na amplitudę absorpcyjną (3.53), piszemy:
![]() |
(3.53)
|
Jeśli wprowadzimy związki na częstotliwość średnią
i wolnozmienną amplitudę oscylacji harmonicznych z niską częstością modulacji napisaną
, wtedy położenie ciała w zależności od amplitudy modulacyjnej i amplitudę drgań przestawiamy:
|
|
Energia modulacji jest wprost proporcjonalna do kwadratu amplitudy modulacji określoną w punkcie (3.55) jest określona w zależności od czasu i częstotliwości drgań podstawowych ω0 i częstotliwości drgań wymuszonych:
![]() |
(3.56)
|
Gdy częstotliwość ω jest równa częstotliwości podstawowej ω0, to wtedy położenie ciała (3.53) możemy napisać, korzystając z wiadomości o granicach na funkcjach trygonometrycznych znanych z analizy matematycznej:
![]() |
(3.57)
|
[edytuj] Drgania niestacjonarne nietrwałe i jego dudnienie
Będziemy liczyli energię drgań niestacjonarnych od czasu w przypadku słabego tłumienia, zatem w takim wypadku czynnik eksponencjalny występujący we wzorze (3.51) w ciągu jednego okresu niewiele się zmienia, wtedy możemy rozważmy przypadek, gdy częstotliwość wymuszająca i drgań tłumionych są bardzo bliskie częstotliwości podstawowej ω0, zatem prędkość masy M liczymy jako pochodna wspomnianego wzoru w postaci przybliżonego wzoru:
![]() ![]() |
(3.58)
|
Całkowita energia układu mechanicznego dla słabego tłumienia, piszemy przy pomocy prędkości drgań masy w zależności od czasu (3.58) i położenia drgań tej samej masy w zależności do czasu (3.51):
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
(3.59)
|
Funkcja kosinus mieści się w zakresie wartości od minus jedynki do jedynki, zatem kres dolny rozwiązania (3.59) jest:
![]() |
(3.60)
|
Kres górny rozwiązania (3.56) na podstawie własności funkcji kosinus rozważanej powyżej przyjmuje postać:
![]() |
(3.61)
|
Według rysunku obok, gdy by nie było tłumienia, oscylacje były by draniami oscylatora harmonicznego i takie dudnienia zachodziły by bez końca. Przy istnieniu dudnienia układ swoją fazę oscylacji ustala do siły wymuszającej, by ostatecznie osiągnąć częstotliwość drgań oscylacji równą ω, a względna faza pomiędzy oscylatorem a siłą wymuszającą zostaje ustalona w taki sposób, by ilość energii dostarczonej do układu była równa energii traconej przez sam oscylator przez siłę oporu spowodowanej w równaniu (3.22) przy współczynniku tłumienia Γ.
[edytuj] Rezonans o dwóch stopniach swobody
[edytuj] Drgania wymuszone dwóch wahadeł wymuszonych
Równania ruchu drgań pierwszego i drugiego wahadła matematycznego, które są sprzężone ze sobą za pomocą sprężyny, do drugiej sprężyny jest podłączona siła, która zmienia się kosinusoidalnie z siłą F0cosωt, to równania ruchu takiego układu są przestawione:
![]() |
(3.62)
|
Możemy dodać i odjąć powyższe równanie układu równań, w ten sposób otrzymać układ następny równoważny do poprzedniego układu równań:
![]() |
(3.63)
|
Nasz układ posiada dwie częstotliwości podstawowe, które są przestawione w punkcie (1.125) i (1.126).
Podamy teraz wzory na amplitudę elastyczną i absorpcyjną, wykorzystując przy okazji wzór (3.23), którego amplituda A jest sumą dwóch amplitud absorpcyjnych podanych dla obu częstotliwości podstawowych (3.29):
![]() ![]() |
(3.64)
|
A amplituda elastyczna dla naszego układu zapisujemy podobnie jak przy amplitudzie absorpcyjnej podanej w punkcie (3.64), która ta amplituda B jest sumą amplitud elastycznych dla obu drgań amplitud elastycznych podanych w punkcie (3.30):
![]() |
(3.65)
|
[edytuj] Drgania wymuszone układów izolowanych o wielu stopniach swobody
[edytuj] Wahadła sprzężone ze sobą w przybliżeniu ciągłości
Wahadła matematyczne sprzężone ze sobą za pomocą sprężynek opisujemy przy pomocy równania (2.85), które dla porządku dziennego jeszcze raz powtórzymy.
![]() |
(3.66)
|
Załóżmy, że funkcja ψn jest wolno-zmienną funkcją n. Co wynika z tego, że wszystkie ciężarki znajdujące niemal blisko siebie wykonują w przybliżeniu te same drgania, zatem napiszmy teraz funkcją ψ dla n, n+1 i n-1:
![]() |
(3.67)
|
![]() |
(3.68)
|
![]() |
(3.69)
|
Na podstawie (3.67), (3.68) i (3.69) można otrzymać dwie następujące tożsamości, która ta pierwsza jest różnicą funkcji wychylenia od stanu równowagi dla n+1 i dla n, a ta druga różnica jest różnicą funkcji zależnej od n dla n i n-1:
|
|
Równości (3.70) i (3.71) możemy podstawić do równania różniczkowo-różnicowego (3.66), w ten sposób otrzymać tożsamość różniczkową drugiego stopnia, z której możemy wyznaczyć funkcję ψ, która jest rozwiązaniem równania opisującej ruch drgający:
![]() |
(3.72)
|
[edytuj] Równanie Kliena-Gordona
Jeśli sobie obierzemy funkcję zależną od z i t, którą jest funkcja harmoniczna zależna od czasu, pomnożonej przez amplitudę zależną od z, to jego definicja:
![]() |
(3.73)
|
A jej drugie pochodne cząstkowe liczone względem czasu i położenia, których pierwsza jest wprost proporcjonalna do A(z), a druga do drugiej pochodnej zupełnej względem położenia kulki:
|
|
Wykorzystując tożsamości (3.74) i (3.75) i podstawiając to wszystko do równości różniczkowej (3.72), wtedy możemy napisać końcowe równanie na A(z):
![]() |
(3.76)
|
[edytuj] Fale sinusoidalne
Fale sinusoidalne możemy opisać równaniem (3.76) dla której definicja stałej k2, która jest wprost proporcjonalna do różnicy kwadratów częstotliwości kołowej drgań wymuszonej ω i częstotliwości podstawowej ω0, podajemy w tej samej linijce oba równania, co równanie falowe:
|
|
wtedy A(z) jest kombinacją funkcji sinus i kosinus, których argumentem jest funkcja kz, którą piszemy przy pomocy stałej A i B:
![]() |
(3.79)
|
Stałe A i B występujące w równaniu (3.79) są to stałe opisywane przez warunki brzegowe.
[edytuj] Fale eksponencjalne
Są to fale opisywane na podstawie (3.76), dla którego definicję stałej χ2, która jest wprost proporcjonalna do różnicy kwadratów częstotliwości kołowej podstawowej drgań ω0 i częstotliwości wymuszonej ω, które te związki podajemy w jednej i tej samej linijce co równanie różniczkowe drgań tejże opisywanej fali:
|
|
Rozwiązanie A(z), które jest rozwiązaniem (3.80) piszemy jako kombinacja funkcji eksponencjalnych, których pierwsza jest funkcja malejącą do zera, a druga rosnąca do nieskończoności dla "z" coraz większego:
![]() |
(3.82)
|
Stałe A i B występujące w równaniu (3.82) są stałe opisywane przez warunki brzegowe.
[edytuj] Związki dyspersyjne
Związki (3.78) i (3.81) możemy zapisać na częstotliwość kołową drgań w zależności od częstotliwości kołowej drgań zależnych od czasu i stałej liczbie falowej "k" lub "χ", które te związki piszemy:
|
|
Ośrodkiem dyspersyjnym nazywamy ośrodek w którym rozchodzą się fale sinusoidalne, a ośrodkiem reaktywnym nazywamy taki ośrodek, w którym rozchodzą się fale eksponencjalne. Ten sam ośrodek może być jednocześnie ośrodkiem dyspercyjnym lub reaktywnym dla różnych zakresów częstotliwości kołowej.
[edytuj] Wahadła matematycznych sprzężonych ze sobą przy pomocy jednakowych sprężynek
Równanie rózniczkowe różnicowe (2.85) możemy zapisać wprowadzając częstość kołową drgań podstawowych wahadła matematycznego ω0, którym to w tym równaniu grupujemy wyrazy względem wyrazów z tym samym wskaźnikiem n:
![]() |
(3.85)
|
Wprowadźmy definicję wychylenia od stanu równowagi od czasu, która zależy od czasu, którą przestawiamy przy pomocy funkcji kosinus z argumentu ωt+φ, którą jest suma iloczynu częstotliwości kołowej przez czas "t" i przesunięcia fazowego φ, co ta funkcja trygonometryczna pomnożona jest przez amplitudę An:
![]() |
(3.86)
|
Wzór na wychylenia od stanu równowagi (3.86) podstawiamy do (3.85), w ten sposób mamy równość na związek dyspersyjny na ω, wtedy:
![]() |
(3.87)
|
[edytuj] Zakres dyspersyjny częstości
Amplituda An możemy przestawić jako kombinację funkcji sinus i kosinus względem argumentu kna, której definicja:
![]() |
(3.88)
|
Definicja amplitud An+1 i An-1, na podstawie definicji An (3.88) (tylko we wzorze tym zamiast n podstawiamy kolejno n+1 i n-1) wykorzystując z twierdzenia o rozdzielności mnożenia względem dodawania, by później było można zastosować twierdzenie od kosinusie sumy i kosinusie różnicy dla tych amplitud:
![]() |
(3.89)
|
![]() |
(3.90)
|
Zapiszmy sumę amplitud An+1 i An-1, możemy je zapisać przy wykorzystaniu tożsamości (3.89) i (3.90), by końcowy wzór na tą sumę przestawić w zależności od amplitudy An:
![]() ![]() |
(3.91)
|
Obliczenia (3.91) podstawiamy do (3.88) przy wykorzystaniu tożsamości trygonometrycznych na kątach połówkowych, otrzymujemy:
![]() |
(3.92)
|
[edytuj] Dolny zakresy reaktywny częstości
Amplitudę An dla zakresu reaktywnego możemy przestawić jako kombinację funkcji eksponencjalnych z argumentów różniących się znakiem w postaci:
![]() |
(3.93)
|
Wtedy suma amplitud An+1 i An-1 możemy opisać przy wykorzystaniu tożsamości (3.93), by końcowy wzór na tą sumę przestawić w zależności od amplitudy An:
![]() ![]() ![]() |
(3.94)
|
Obliczenia (3.94) możemy podstawić do równości (3.87), w ten sposób otrzymać tożsamość matematyczną, do którego wykorzystamy definicję funkcji hiperbolicznych:
![]() |
(3.95)
|
Wzór (3.94) możemy zapisać w formie podobnej do (3.92) przy wykorzystaniu tożsamości na funkcjach hiperbolicznych, którego zapis w formie bardziej uproszczonej:
![]() |
(3.96)
|
[edytuj] Górny zakres reaktywny
Przyjmijmy, że amplituda jest opisywana przez zygzakowatą (przemienną) falę eksponencjalną, której definicja:
![]() |
(3.97)
|
Wtedy suma amplitud An+1 i An-1 możemy zapisać przy wykorzystaniu tożsamości (3.97), by końcowy wzór na tą sumę przestawić w zależności od amplitudy An:
![]() ![]() ![]() |
(3.98)
|
Obliczenia (3.98) możemy podstawić do równości (3.87), w ten sposób otrzymać tożsamość, do którego wykorzystamy definicję funkcji hiperbolicznych z znanej analizy matematycznej, doprowadzając końcowe równanie do postaci bardziej uproszczonej:
![]() |
(3.99)
|
[edytuj] Opis fal biegnących
Będziemy rozpatrywali tutaj układy otwarte, tzn. nie mające zewnętrznych granic. Powietrze na dworze dla rozchodzących się fal dźwiękowych zachowuje się jak układ otwarty, pomijamy tutaj efekty powstałe dotyczące odbić fal dźwiękowych od budynków. Rozchodzące się fale dźwiękowe w zamkniętym pomieszczeniu zachowuje się jak układ zamknięty. Jeśli ścianki tego pomieszczenia pokryjemy materiałem całkowicie pochłaniających fale dźwiękowe, to wtedy taki układ zachowuje się jak układ otwarty. Rozpatrzmy tutaj fale, których wychylenia od stanu równowagi piszemy:
![]() |
(4.1)
|
Rozpatrzmy układ, który zachowuje się jak w punkcie z=0 w czasie t', który piszemy w zależności od "t" i położenia, w której rozchodzi się fala z prędkością fazową vf:
![]() |
(4.2)
|
Wzór (4.2) możemy podstawić do wzoru (4.1), w ten sposób otrzymać równość zależną od czasu i położenia zetowego "z", w której znajduje się fala biegnąca:
![]() |
(4.3)
|
Dla fal bięgnących, która jest funkcją kosinus, zapisaną w zależności od liczby falowej i częstotliwości kołowej, przestawiamy wzorem:
![]() |
(4.4)
|
Możemy porównać wzór (4.4) ze wzorem (4.3), które oznaczają to samo, w ten sposób dostajemy wniosek, że liczba falowa jest stosunkiem częstotliwości kołowej drgań harmonicznych i prędkości fazowej drgań:
![]() |
(4.5)
|
co z (4.5) możemy napisać definicję prędkości fazowej w zależności od częstotliwości kołowej i liczby falowej, która jest ich stosunkiem, w tej samej linijce napiszmy dwa następne wzory określającą prędkość fazowa, która za pierwszym razem jest iloczynem długości falowej i częstotliwości, a za drugim razem jest stosunkiem liczby falowej przez okres drgań fal biegnących:
|
|
|
Zdefiniujemy funkcję fazową, która jest zdefiniowana w zależności od częstotliwości kołowej i liczby falowej, a także od zmiennych czasu "t" i położenia "z":
![]() |
(4.9)
|
Podamy teraz inną metodę wyprowadzenia wzoru na prędkość fazową, zatem na podstawie definicji funkcji fazowej (4.9), której różniczka jest napiszemy z definicji różniczki zupełnej, którą wyrazimy w zależności od częstotliwości kołowej ω i liczby falowej "k":
![]() |
(4.10)
|
Rozpatrzmy, gdy różniczka funkcji fazowej jest równa dokładnie zero, w ten sposób dostajemy równość z definicji prędkości znaną z kinematyki klasycznej, którą nazwiemy prędkością fazową:
![]() |
(4.11)
|
[edytuj] Związek dyspersyjny dla wahadeł matematycznych sprzężonych przy pomocy sprężynek
Równanie różnicowo-różniczkowe możemy napisać wzorem (2.85), które dla porządku dziennego tutaj przepisujemy:
![]() |
(4.12)
|
Dla fal biegnących, którego wychylenie od stanu równowagi jest opisane w punkcie (4.4), dla której drugą pochodną cząstkową względem czasu piszemy:
![]() |
(4.13)
|
Równanie (4.12) możemy zapisać, przy pomocy zachodzącej tożsamości dla fal biegnących (4.13), która jest drugą pochodną cząstkową wychylenia podłużnego względem czasu, która jest iloczynem wychylenia od stanu równowagi ψn przez kwadrat częstotliwości kołowej drgań i to wszystko wzięte z minusem:
![]() |
(4.13)
|
[edytuj] Harmoniczne fale biegnące
Dla fal sinusoidalnych fal biegnących wychylenie od stanu równowagi możemy przepisać w zależności od częstotliwości kołowej i liczby falowej w zależności od zmiennej czasowej i położenia z=na:
![]() |
(4.14)
|
Suma wychyleń (4.14) dla n+1 i n-1 piszemy przy pomocy wychylenia od stanu równowagi dla n w sposób:
![]() |
(4.15)
|
Związek (4.15) możemy podstawić do (4.13), w ten sposób dostajemy związek dyspersyjny na kwadrat częstotliwości kołowej w zależności od liczby falowej w postaci:
![]() |
(4.16)
|
[edytuj] Eksponencjalne fale biegnące w układzie otwartym
Dla fal eksponencjalnych biegnących w układzie otwartym wychylenie od stanu równowagi dla z=na zależnego od częstotliwości kołowej i stałej zaniku fal χ, piszemy:
![]() |
(4.17)
|
Suma wychyleń (4.17) dla n+1 i n-1 piszemy przy pomocy wychylenia od stanu równowagi dla n w sposób:
![]() ![]() |
(4.18)
|
Tożsamość napisana w punkcie (4.18) podstawiamy do (4.13), przy wykorzystaniu definicji funkcji hiperbolicznych, w ten sposób mamy związek dyspersyjny na kwadrat częstotliwości kołowej względem wielkości stałej zaniku dla rozchodzących się fal biegnących eksponencjalnych:
![]() |
(4.19)
|
[edytuj] Przemienne eksponencjalne fale biegnące w układzie otwartym
Dla fal eksponencjalnych biegnących w układzie otwartym wychylenie od stanu równowagi, dla z=na zależnego od częstotliwości kołowej i stałej zaniku fal χ, piszemy w postaci:
![]() |
(4.20)
|
Suma wychyleń (4.20) dla n+1 i n-1 piszemy przy pomocy wychylenia od stanu równowagi dla n:
![]() ![]() |
(4.21)
|
Tożsamość napisana w punkcie (4.21) podstawiamy do (4.13), przy wykorzystaniu definicji funkcji hiperbolicznych, w ten sposób mamy związek dyspersyjny na kwadrat częstotliwości kołowej w zależności od stałej zaniku χ dla rozchodzących się fal biegnących:
![]() |
(4.22)
|
[edytuj] Prędkość fazowa fal poprzecznych rozchodzącej się w strunie obciążonymi ciężarkami
Prędkość fazowa rozchodzenia się fal w strunie obciążonej sprężynkami możemy napisać patrząc na definicję kwadratu częstotliwości kołowej (2.71), z którego wyznaczymy samą częstotliwość, co po podzieleniu jego obustronnie przez liczbę falową "k" otrzymamy prędkość fazową:
![]() |
(4.23)
|
Biorąc przypadek dla małej liczby falowej "k", wtedy obliczenia dla wzoru (4.23) przeprowadzamy tak jak w punkcie (2.77), w ten sposób możemy otrzymać wzór na prędkość fazową w przybliżeniu małej liczby falowej wykorzystując przy tym definicję masowej gęstości liniowej, zatem dokonajmy tychże obliczeń:
![]() ![]() |
(4.24)
|
[edytuj] Prędkość fazowa fal podłużnych rozchodzącej się w strunie obciążonymi ciężarkami
Prędkość rozchodzenia się fal w zależności liczby falowej możemy przestawić wykorzystując wzór na kwadrat częstotliwości kołowej (2.79) z którego wyznaczymy samą częstotliwość, co po podzieleniu jego przez liczbę falową "k", w ten sposób dostajemy wzór na prędkość fazową rozchodzenia się fal:
![]() |
(4.25)
|
Dokonując podobnych obliczeń jak w punkcie (4.24) dla k bardzo małego, w ten sposób wykorzystując definicję gęstości liniowej, by otrzymać później tożsamość:
![]() |
(4.26)
|
bo zachodzi KLL=Ka, w którym każda sprężyna ma współczynnik sprężystości równej K, sprężynki wszystkie razem wzięte mają współczynnik sprężystości KL i dlatego tak powinniśmy napisać.
[edytuj] Prędkość rozchodzenia się dźwięku w naczyniu
Będziemy tutaj wyprowadzać równanie rządzącego przy rozchodzącego się dźwięku w danym ośrodku znając pochodną zupełną ciśnienia względem objętości w stanie równowagi danego ośrodka, tutaj będziemy się zajmować powietrzem, ale związek , który tutaj otrzymamy jest zgodny dla innych ośrodków sprężystych a także dla powietrza. Powietrze znajdujące się w naczyniu zachowuje się jak sprężynka o stałej sprężystości KL, którego długość jest L1, która zostaje ściśnięta do długości L, a siła wywierana na ścianki naczynia jest wyrażona wzorem:
![]() |
(4.27)
|
Możemy zróżniczkować tożsamość (4.27), w ten sposób otrzymać równość różniczkową:
![]() |
(4.28)
|
Z drugiej jednak strony powiedźmy, gdy nastąpi zmiana długości sprężyny dL towarzyszy zmiana ciśnienia dp i zmiana siły wywieranej dF, zatem tą siłę piszemy:
![]() |
(4.29)
|
Porównując wzory (4.29) z (4.28) możemy napisać tożsamość na współczynnik sprężystości powietrza:
![]() |
(4.30)
|
Gęstość liniową masy możemy powiązać z gęstością objętościową, która jak udowodnimy gęstość liniowa masy jest równa iloczynowi gęstości objętościowej pomnożonej przez powierzchnię przekroju ośrodka, w której dźwięk rozchodzi się względem linii prostopadłej względem tego przekroju, według związku:
![]() |
(4.31)
|
Wzory (4.30) i (4.31) możemy podstawić do wzoru końcowego (4.24), w ten sposób otrzymujemy końcową równość na kwadrat prędkości rozchodzenia się dźwięku:
![]() |
(4.32)
|
[edytuj] Prędkość rozchodzenia się fal w powietrzu
Biorąc równanie stanu gazu doskonałego dla izotermy, czyli pod stałą temperaturą, wyznaczać z niego będziemy ciśnienie panujące w panujące gazie w zależności od jej objętości:
![]() |
(4.33)
|
Wyznaczmy teraz pierwszą pochodną zupełną ciśnienia (4.33) względem objętości V i kładąc V=V0:
![]() |
(4.34)
|
Wzór (4.34) podstawiamy do (4.32), w ten sposób dostajemy równość, którą wyrażamy w zależności od ciśnienia gazu i jego gęstości objętościowej:
![]() |
(4.35)
|
[edytuj] Poprawka do modelu dotyczącego prędkości rozchodzenia się dźwięku
Model opisujący prędkość dźwięku w gazie opisaną w punkcie (4.35) jest modelem tylko dobrym, tylko jest jedyna różnica, że w tym modelu prędkość dźwięku jest napisana z błędem 15% w porównaniem z prędkością dźwięku policzonego z doświadczenia, zatem ten model zawiera pewne błędy, który świadczy, że zastosowaliśmy model gazu doskonałego dla stałej temperatury układu. Lepszym modelem jest równanie adiabaty. W porównaniu modelu Boyla i adiabaty przy sprężaniu gazu ciśnienie jest nieco większe a przy rozprężaniu jest nieco mniejsze. Podczas rozprężaniu i sprężenia gazu temperatura gazu nie jest w rzeczywiści stała w danym punkcie przestrzeni. Biorąc równanie stanu gazu doskonałego dla adiabaty, czyli układu nie wymieniającego się energii z otoczeniem, wyznaczać z niego będziemy ciśnienie panujące w gazie w zależności od jej objętości:
![]() |
(4.36)
|
Wyznaczmy teraz pierwszą pochodną zupełną ciśnienia (4.37) względem objętości V i kładąc V=V0, by w ten sposób otrzymać równość:
![]() |
(4.37)
|
Wzór (4.37) podstawiamy do (4.32), w ten sposób otrzymać równość na prędkość fazową rozchodzenia się fal dźwiękowych, którą wyrażamy w zależności od ciśnienia gazu i jego gęstości objętościowej, która panuje w stanie równowagi:
![]() |
(4.38)
|
[edytuj] Prędkość fazowa w linii przesyłowym
Częstotliwość fazową linii przesyłowej możemy przestawić poprzez wzór (2.83), który możemy przestawić w granicy małych liczb falowych w postaci:
![]() |
(4.39)
|
Kwadrat prędkości fazowej przestawiamy z jego definicji (4.11) dla małych liczb falowych, która jest przestawiona jako odwrotność iloczynu ilorazu indukcji własnej przez przez długość linii przez iloraz pojemności linii przez jej długość:
![]() |
(4.40)
|
[edytuj] Dolnoprzepustowa linia przesyłowa
Z prawa Gaussa dla elektrostatyki, w zależności od ładunku znajdującego się wewnętrznej powierzchni możemy opisać z tego prawa wniosek:
![]() |
(4.41)
|
Wyznaczmy teraz różnicę potencjału między dolną częścią a jej górną częścią kondensatora, którą możemy przestawić według obliczeń:
![]() |
(4.42)
|
Całkowita pojemność kondensatora (4.40) jest przestawiana z jego definicji jako iloraz ładunku na okładkach kondensatora cylindrycznego przez jego różnicę potencjałów pomiędzy cylindrycznymi okładkami:
![]() |
(4.43)
|
Wyznaczmy indukcję magnetyczną pola magnetycznego znajdującego się w kondensatorze pomiędzy jego okładkami:
![]() |
(4.44)
|
Policzmy strumień indukcji pola magnetycznego z jego formalnej definicji, a później indukcję własną:
![]() |
(4.45)
|
Napiszmy teraz iloczyn pojemności elektrycznej kondensatora C i indukcyjności z jej definicji, tzn. (4.41) i (4.42):
![]() |
(4.46)
|
Jeśli popatrzymy na wzór na kwadrat prędkości fazowej (4.40) i na wzór (4.46) dowiadujemy się, że prędkość fazowa jest równa prędkości światła w dolnoprzepustowej linii przesyłowej.
[edytuj] Linia przesyłowa zrobiona przy pomocy dwóch równoległych pasków przewodzących
Różnica potencjałów między paskami określamy z jego definicji, a także natężenie pola elektrycznego określamy z prawa Gaussa, piszemy je związkami:
|
|
Pojemność elektryczna, która jest ilorazem ładunku znajdującego na kondensatorze przez różnicę potencjału znajdującego się na okładkach, wyraża się:
![]() |
(4.49)
|
Strumień indukcji pola magnetycznego względem powierzchni bocznej możemy określić jako iloczyn igrekowej współrzędnej indukcji pola magnetycznego przez powierzchnię o wartości ga:
![]() |
(4.50)
|
Z prawa Stokesa możemy również powiedzieć, że iloczyn współrzędnej igrekowej pola magnetycznego i długości bocznej kondensatora jest równa iloczynowi przenikalności magnetycznej ośrodka, w którym jest próżnia, przez natężenie prądu elektrycznego płynącego w linii:
![]() |
(4.51)
|
Wzór na indukcję pola magnetycznego a właściwie na jej współrzędną igrekową podstawiamy do (4.50), wykorzystując definicję indukcyjności własnej, możemy pisać:
![]() |
(4.52)
|
Mnożymy wzory (4.49) i (4.52), w ten sposób otrzymujemy odwrotność na kwadrat prędkości światła w próżni:
![]() |
(4.53)
|
Patrząc wzór na kwadrat prędkości fazowej (4.40) i definicję iloczynu ilorazu pojemności linii przesyłowej przez "a" i przez iloraz indukcji pola magnetycznego przez "a". dowiadujemy się, że prędkość fazową linii przesyłowej jest równa prędkości światła w próżni "c".
[edytuj] Zjawisko załamania świata i dyspersja
Z rysunku obok otrzymujemy natychmiast, następujące dwa równania, które wynikają z rysunku obok:
|
|
Dzielimy wzory (4.54) i (4.55) względem siebie, dostajemy:
![]() |
(4.56)
|
Teraz wykorzystamy definicję współczynnika załamania n=c/v dla rozważanego ośrodka pierwszego i drugiego, to na podstawie tego równość (4.56) przepisujemy jako:
![]() |
(4.57)
|
[edytuj] Wyprowadzenie wzoru na współczynnik załamania promieni świetlnych
Względna przenikalność elektryczną możemy przepisać według następującego schematu:
![]() |
(4.58)
|
Dalszym krokiem jest wyprowadzenie wzoru na współczynnik przenikalność elektryczną względną, wiedząc że
, którego przestawienie jest:
![]() |
(4.59)
|
Na podstawie alternatywnej definicji polaryzacji elektrycznej P=Nqx(t) przenikalność względna elektryczna możemy przepisać do:
![]() |
(4.60)
|
Napiszmy teraz równanie ruchu ładunku elektrycznego, który napiszemy za pomocą częstotliwości drgań podstawowych ω0, współczynnika tłumienia Γ, i przy pomocy natężenia pola elektrycznego zmieniającego się według funkcji kosinus E(t)=E0cosωt, wtedy:
![]() |
(4.61)
|
Położenie ciała o masie M możemy zapisać jako rozwiązanie równania (4.61) przy pomocy amplitudy elastycznej Ael i absorpcyjnej Aab według:
![]() |
(4.62)
|
Rozważmy sobie ośrodek przezroczysty, takich jak szkło, woda, są one przezroczyste i nie następuje wtedy absorpcja promieni świetlnych. Dla ośrodków barwnych, nastepuje absorpcja wszystkich częstości oprócz jednej, i interesują nam te częstości, które są daleko od częstości rezonansowych, wtedy amplituda absorpcyjna według (3.29) jest w przybliżeniu równa zero, zatem istnieje tylko amplituda elastyczna, która jest równa dla naszego przypadku według (3.30):
![]() |
(4.63)
|
Powyżej skorzystano z definicji współczynnika załamania, a także z definicji prędkości światła w ośrodku materialnym, dla którego względna przenikalność magnetyczna jest w przybliżeniu równa jeden. Jeśli powiązać kwadrat współczynnika załamania z liczbą falową, z częstotliwością kołową i z prędkością światła w próżni, to wtedy wzór (4.63) jest w postaci:
![]() |
(4.64)
|
Związek (4.64) jest związkiem liczby falowej z częstotliwością kołową, który jest związkiem dyspersyjnych dla światła widzialnego.
[edytuj] Strumień energii i impedancja
Prędkość fazową rozchodzenia się fal w strunie obciążonych punktowymi ciężarkami określamy jako pierwiastek z stosunku siły napięcia działająca na strunę T0 i gęstości liniowej struny ρ0, a impedancję charakterystyczną nazywamy pierwiastek z iloczynu gęstości liniowej struny i siły napięcia, których definicje tychże wielkości określamy poprzez wzory:
|
|
Poniżej podamy również wzór na impedancję w tej samej strunie, ale dla drgań poprzecznych.
[edytuj] Fale poprzeczne biegnące w strunie ciągłej
Całkowita siła iksowa działająca na ciągłą strunę wywołująca jego ruch falowy w postaci fal biegnących określamy przez:
![]() |
(4.67)
|
[edytuj] Wprowadzenie do definicji impedancji charakterystycznej
Każdy element struny ciągłej wykonuje ruch harmoniczny określanych jako funkcja, dla której jest ona wprost proporcjonalna do kosinusa z argument różnicy iloczynu częstotliwości kołowej przez czas i iloczynu liczby falowej przez położenie danego elementu struny, co jego definicja tejże funkcji piszemy przez:
![]() |
(4.68)
|
Pierwsza pochodna wychylenia od stanu równowagi ψ(z,t) względem położenia z i czasu t struny piszemy przez:
|
|
Patrząc na wzory (4.69) i (4.70) możemy napisać wzór łączący pochodną cząstkową wychylenia od stanu równowagi (4.42) względem położenia z i pochodnej cząstkową wychylenia względem czasu w sposób:
![]() |
(4.71)
|
Całkowita siła działająca na strunę ze strony generatora (4.65) możemy napisać wykorzystując wniosek (4.71), wtedy:
![]() |
(4.72)
|
Wprowadźmy teraz impedancję, jest ona stosunkiem siły T0 i prędkości fazowej:
![]() |
(4.73)
|
Przy definicji impedancji (4.73) wzór na siłę iksowa działającą na strunę ciągłą określamy jako iloczyn pochodnej cząstkowej wychylenia od stanu równowagi ψ względem czasu "t" przez impedancję Z, to razem wszystko wzięte z minusem:
![]() |
(4.74)
|
[edytuj] Moc wytwarzana przez generator drgań struny
Całkowita wytwarzana moc w czasie "t" jest iloczynem siły Fx zdefiniowaną w punkcie (4.74) i prędkości wychylenia od stanu równowagi, która jest pochodną cząstkową wychylenia ψ względem czasu:
![]() |
(4.75)
|
Innym sposobem wyprowadzenie mocy wytwarzanej przez generator jest opisany przez wzór w zależności od pierwszej pochodnej cząstkowej wychylenia względem położenia:
![]() |
(4.76)
|
[edytuj] Promieniowanie wytarzane przez układ sprężynek i mas
Niech mamy układ sprężynek, które łączą jednakowe masy, wtedy wielkość odpowiadająca wielkości napięcia w strunie T0 jest wielkość Ka, wtedy impedancję możemy określić jako pierwiastek z wielkości Ka i gęstości masowej liniowej:
![]() |
(4.77)
|
Całkowita moc wytwarzana przez generator podczas wytwarzania fal podłużnych w układzie sprężynek w zależności od czasu "t" lub położenia jest pisana:
![]() |
(4.78)
|
Całkowita siła działająca na układ mas sprężynek przy założeniu, że w stanie równowagi sprężynki pozostają napięte z siłą F0, która nie wpływa na ruch drgający układów kulek, jest określona:
![]() |
(4.79)
|
[edytuj] Fale głosowe
Prędkość fazowa fal dźwiękowej jest określana przez wzór (4.38), porównując ten wzór z (4.26) (w tym drugim mamy gęstość liniową, a w tym pierwszym gęstość objętościową), dochodzimy do wniosku, że rolę napięcia sprężyny spełnia κp0, wtedy impedancję piszemy:
![]() |
(4.80)
|
Natężenia fali dźwiękowej można napisać na podstawie tożsamości (4.76) określoną względem pochodnej wychylenia podłużnego ψ(z,t) względem czasu "t" lub położenia "z":
![]() |
(4.81)
|
Siłą działająca na jednostkę powierzchni określamy względem stałego ciśnienia panującego w ośrodku p0, która to fala dźwiękowa zaburza stałe nasze ciśnienie atmosferyczne przy pomocy członu harmonicznego, który powoduje, że mamy raz ciśnienie wyższe, a za drugim niższe od naszego rozważanego ciśnienia p0.
![]() |
(4.82)
|
[edytuj] Odbicie fal klasycznych
Całkowita siła działająca na falę by je wytworzyć jest wyrażona jako wielkość wprost proporcjonalna do pierwszej pochodnej cząstkowej wychylenia od stanu równowagi i impedancji Z, którą wyrazimy wzorem (4.48). Falę padającą możemy wyrazić przez (4.4), którą tutaj przepiszemy:
![]() |
(5.1)
|
W punkcie z=0 fala padająca jest opisana wzorem (5.1), którą możemy zapisać formułą:
![]() |
(5.2)
|
[edytuj] Współczynnik odbicia
Oznaczmy ostatni punkt struny przez L, a amortyzator przez P, wtedy siła działająca ze strony amortyzatora na falę padającą i odbitą jest zapisana:
|
|
A więc całkowita siła działająca jednocześnie na falę odbitą i padającą jest sumą sił (5.3) i (5.4), i po prostu jest napisana:
![]() |
(5.5)
|
Fala odbita ψodb i padająca ψpad, której superpozycją jest fala ψ, której będziemy liczyli pochodną cząstkową względem czasu, można ją rozłożyć na sumę pochodnych cząstkowych względem czasu fali padającej i odbitej:
![]() |
(5.6)
|
Z drugiej jednak całkowita siła działająca ze strony oporu amortyzatora równa jest sile oporu amortyzatora na strunę, wtedy całkowita siła działająca na "falę" ψ,, dla którego impedancja jest wyrażona przez Z2, jest wyrysowana:
![]() ![]() |
(5.7)
|
Siła (5.5) jest równa sile wyrażona troszeczkę w innej postaci (5.7), wtedy przyrównując do siebie te siły, bo one oznaczają to samo, dostajemy wniosek:
![]() |
(5.8)
|
Wzór (5.8) możemy napisać troszeczkę w innej postaci grupując wyrazy w nim występujące, by po lewej stronie znajdowała się pochodna cząstkowa wychylenie od stanu równowagi fali odbitej względem czasu, a po prawej stronie znajdowała się pochodna cząstkowa wychylenia fali padającej od stanu równowagi względem czasu, co to wszystko możemy wyrazić przez wzór:
![]() |
(5.9)
|
[edytuj] Fale padająca i odbita, superpozycja tychże fal
W równanie ruchu (5.9) możemy przecałkować obie jego strony, w ten sposób otrzymamy wychylenie od stanu równowagi fali odbitej w zależności od amplitudy fali padającej:
![]() |
(5.10)
|
Współczynnik odbicia występującej we wzorze (5.8) nazywamy wielkość określaną:
![]() |
(5.11)
|
Fala odbita jest falą biegnącą w kierunku -z, której jej podstać dla z<0 otrzymujemy na podstawie jej związku dla z=0 zastępując zmienne z=0 i "t" zmiennymi "z", t+z/vf, dla której vf jest prędkością fazową fali. Wzór na falę odbitą piszemy przez:
![]() |
(5.12)
|
Całkowite przemieszenie poszczególnych części fali jest sumą fali padającej i odbitej wyrażony związkiem (5.12):
![]() |
(5.13)
|
[edytuj] Przemieszczenia fali dobitej i padającej, a siła kierująca
Ogólnie siłą kierująca możemy wyrazić przez wzór zdefiniowany w punkcie (4.41), która jak wiemy jest zależna od pierwszej pochodnej cząstkowej przemieszczenia względem położenia "z", która jest równa sile kierującej działającej ze strony sprężyny znajdującej się na lewo od punktu "z" na sprężynę ze strony prawej od tego punktu. Jak można udowodnić współczynnik odbicia pierwszej pochodnej przemieszczenia względem czasu dla tych fal jest taki sam jak dla samych przemieszczeń, a współczynnik odbicia dla pochodnej przemieszczeń względem położenia jest równy minus współczynnika odbicia przemieszczeń, który wyrażamy przy pomocy poprzednich rozważań równej (5.11). Aby to wszystko udowodnić napiszmy poniższe rozważania.
Przemieszczenia fali padającej i odbitej wyrażamy przy pomocy wzoru (5.12):
|
|
Porównując wzory (5.14) i (5.15) dochodzimy do wniosku, że współczynnik proporcjonalności jest R12. Wtedy możemy policzyć pochodne cząstkowe wychyleń względem czasu fali padającej (5.14) i odbitej (5.15):
|
|
Współczynnik proporcjonalności prędkości przemieszczeń między falą odbita i padającą jest to współczynnik napisany jako R12. Pochodne cząstkowe przemieszczeń względem położenia danego elementu układu drgającego możemy napisać:
|
|
Porównując tożsamość (5.18) z (5.20) możemy dojść do wniosku, że współczynnik odbicia dla pochodnej cząstkowej przemieszczenia względem położenia określamy poprzez wzór:
![]() |
(5.20)
|
[edytuj] Ciągłość parametrów fali na granicy dwóch ośrodków
Na granicy między dwoma ośrodkami prędkość, czyli pierwsza pochodna przemieszczenia względem czasu, a także siła kierująca (definicja (4.41)) są wielkościami ciągłymi na granicy dwóch ośrodków. Siła działająca na falę padająca i przechodzącą na tej granicy wziętej razem dwóch ośrodków jest wyrażana przez:
![]() |
(5.21)
|
Skokowi siły napinającej struny według (5.21) towarzyszy kant struny i w rezultacie nachylenie struny nie jest wielkością ciągłą. natomiast wielkością ciągłą jest iloczyn nachylenia struny i siły napinającej.
[edytuj] Współczynnik transmisji
Napiszmy teraz wzór, który określimy przy pomocy współczynnika odbicia R, która jest superpozycją fali padającej i odbitej dla ośrodka pierwszego, którą piszemy poprzez wzór:
![]() |
(5.22)
|
Wzór (5.21) określa wielkości falowe, którymi mogą być przemieszczenie, prędkość, lub nawet siła kierująca, który nasz wspomniany wzór piszemy przy pomocy współczynnika odbicia R, który dla przemieszczeń i prędkości są wyrażane przez R12. Fala rozchodząca się w ośrodku drugim z ośrodka pierwszego jest określana przy pomocy współczynnika transmisji T:
![]() |
(5.23)
|
Współczynnik transmisji możemy wyznaczyć dla z=0, wtedy następująca równość wielkości (5.23) z wielkością (5.22) piszemy według:
![]() |
(5.24)
|
Do wzoru (5.24) możemy podstawić wzory (5.22) i (5.23), by w ten sposób otrzymać tożsamość:
![]() |
(5.25)
|
Z równości (5.25) możemy otrzymać wzór na współczynnik transmisji, który jest sumą jedynki i współczynnika odbicia:
![]() |
(5.26)
|
We wzorze (5.26) R jest równe R12, gdy φ0 jest to przemieszczenie φ lub prędkością ∂ψ/∂ t, natomiast gdy jest równa -R12, co wtedy mamy doczynienia z siła kierującą -T∂ψ/∂ z. (Uwaga: Symbol T jest używany zarówno do oznaczenia napięcia struny, jak i dla współczynnika transmisji).
[edytuj] Impulsy, paczki i modulacje fal
Dotychczas rozważaliśmy drgania, które były drganiami harmonicznymi o jednej tylko częstości ω. Dowiedzieliśmy, że superpozycja fal o jednakowych amplitudach, ale o prawie jednakowych częstościach kołowych prowadzi do zjawiska dudnień. Omówimy tutaj zjawisko dudnień powstałych z nakładania się wielu fal o zbliżonych amplitudach. Jak się okazuje, że fale mogą rozchodzić się z prędkością grupową niosąc ze sobą energię z tą właśnie prędkością.
[edytuj] Prędkość grupowa
Sumę dwóch drgań harmonicznych o jednakowych amplitudach, ale dla różnych liczb falowych i częstotliwości kołowych, piszemy w postaci:
![]() |
(6.1)
|
Definicję częstotliwości modulowanej i liczby falowej modulowanej, a także ich odpowiedniki średnie rysujemy:
|
|
||||
|
|
Wykorzystując powyższe definicję średniej częstotliwości kołowej i średniej liczby dwóch fal o zbliżonych częstościach, możemy określić jako sumę (6.1) w postaci zwartej z wyodrębnionymi amplitudami modulowanymi zależnych od częstotliwości i liczb falowych modulowanych, a także z wyodrębnioną częścią harmoniczną, która zależy od średniej częstotliwości i średniej liczby falowej, co te wszystkie wnioski określimy przez dwa poniższe wzory:
|
|
Funkcję fazową amplitudy (6.7) można napisać w takiej postaci, która po zróżniczkowaniu obu jego stron dla tej funkcji stałej, w ten sposób otrzymujemy wzór, z którego możemy wyznaczyć prędkość rozchodzenia się modulacji:
![]() |
(6.8)
|
Prędkość modulacji możemy rozłożyć w szereg Taylora, w ten sposób dowiadujemy się, że prędkość rozchodzenia się modulacji, która jest ilorazem różnicowym częstotliwości kołowej i liczby falowej k, jest ona w przybliżeniu równa pochodnej zupełnej częstotliwości kołowej względem liczby falowej:
![]() |
(6.9)
|
Na podstawie tożsamości (6.9) możemy określić prędkość grupową, którą nazywamy pochodną zupełną częstotliwości kołowej względem liczby falowej:
![]() |
(6.10)
|
[edytuj] Składanie się fal o prostokątnym rozkładzie widmowym
Superpozycja N fal, których amplitudy są jednakowe, ale o częstościach kołowych równomiernie rozłożonych pomiędzy częstościami ω1, a ω2, których ta całkowita fala ψ(t) jest zdefiniowana:
![]() |
(6.11)
|
Wielkość δω jest tak zdefiniowana w taki sposób jako iloraz różnicy częstotliwości kołowych ω2 i ω1 podzielonej przez N-1, co matematycznie:
![]() |
(6.12)
|
Najlepszym sposobem jest przestawienie sumy drgań harmonicznych (6.11) przy pomocy funkcji eksponencjalnych, których częścią rzeczywistą jest ta nasza rozważana suma drgań harmonicznych (6.11), którego przestawienie jest:
![]() |
(6.13)
|
Funkcja (6.13) jest szeregiem geometrycznym o ilorazie exp(-δωt), wtedy zapisujemy go w postaci zwartej:
![]() ![]() |
(6.14)
|
Całkowita funkcja falowa zapisana w postaci zwartej (6.11) jest częścią rzeczywistą obliczeń napisaną w punkcie (6.14), którą rysujemy w postaci funkcji zależnej od czasu "t", częstotliwości średniej drgań ωśr:
![]() |
(6.15)
|
Amplituda fali ψ(t), występująca we wzorze (6.15) jest zależna od czasu i jest napisana przy pomocy funkcji sinus, która występuje w mianowniku i liczniku, której licznik jest zależny od liczby przegród, jest wykreślona:
![]() |
(6.16)
|
Gdy N=2, wtedy równość na funkcję ψ(t) (6.15), którą definiujemy zamienieniając 1/2ωt przez zmienną x, wtedy mamy:
![]() |
(6.17)
|
Wzór na funkcję ψ(t) dla N=2 otrzymaliśmy jak w punkcie (1.116), co jest zgodnie z naszymi oczekiwaniami. Przy liczeniu funkcji (6.16) dla t=0 należy się dokładnie przyjrzeć, bo licznik i mianownik tego wyrażenia jest równy zero, wtedy należy dokonać zamiany zmiennych według
, i obliczyć granicę wykorzystując wiadomości o granicach z analizy matematycznej.
![]() |
(6.18)
|
Amplitudę (6.16) możemy przestawić jako funkcję A(t) przy A=A(0)/N, dla którego tą wielkość piszemy:
![]() |
(6.19)
|
Przejdźmy teraz do granicy, dla którego N jest bardzo duże, wtedy odstępy pomiędzy składowymi drgań harmonicznych δω są na tyle małe, że ich doświadczalnie nie możemy rozróżnić, wtedy mamy doczynienia z fizyką, a nie z matematyką. Amplitudę (6.16) możemy napisać dla warunku δω bardzo małego i przy zachodzącej tożsamości (6.12):
![]() |
(6.20)
|
Wychylenie od stanu równowagi (6.11) jako superpozycja N fal harmonicznych, przy wykorzystaniu (6.12), piszemy jako:
![]() |
(6.21)
|
Równość (6.21) możemy zapisać w postaci całkowej, przy założeniu, że δω jest bardzo małe, praktycznie zerowe, jako całkę od ω=ω1, aż do ω=ω2:
![]() |
(6.22)
|
[edytuj] Szereg Fouriera liczona względem czasu
Szeregiem Fouriera nazywamy takie wyrażenie, której każdą funkcję F(t) możemy rozłożyć ten nasz szereg, którego postać jest zależna od współczynników B0 i współczynników zależnych od n, czyli od An i od Bn:
![]() |
(6.23)
|
Wówczas jak możemy udowodnić, że współczynniki An i Bn możemy napisać przy pomocy funkcji F(t), jak to robiliśmy też dla szeregu Fouriera zdefiniowanego w punkcie (2.35), tylko tam była definicja względem liczby falowej, a tutaj względem częstotliwości kołowej:
|
|
![]() |
(6.26)
|
[edytuj] Przejście od szeregu Fouriera w jej postać całkową
Rozważmy kilka pierwszych wyrazów w szeregu Fouriera, tzn. A1sinω1t+B1cosω1t, A1sinω1t+B1cosω1t. Pierwsze wyrazy są tak małe w podanych wyrazach, że można je zaniedbać. Widzimy, że T1, która jest okresem. Sztucznie utworzona funkcja zależy od naszego okresu, więc ten okres możemy zwiększać o dowolną wielokrotność, zatem okres T1 jest dowolny, to częstotliwość kołową ω1=2π/T1 możemy dowolnie zmniejszać. W gruncie rzeczy T1 możemy wziąć na tyle duże by wyrazy An i Bn można było zaniedbać. Rozważmy n na tyle duże by nie można było zaniedbać współczynników An , Bn, wtedy szereg (6.23) możemy zapisać biorąc jak na razie dla uproszczenia wyrazy z sinusami:
![]() |
(6.27)
|
Obierzmy sobie teraz funkcję ω zmiennej n i ω1, która jest iloczynem liczby n i częstotliwości ω:
![]() |
(6.28)
|
Niech przyrost δω będzie taki, że różniczka z n jest ilorazem z delty częstotliwości kołowej ω przez częstotliwość ω1:
![]() |
(6.29)
|
Załóżmy, że w pasmie n do n+δ n wszystkie współczynniki An są sobie równe, wtedy wszystkie te wyrazy możemy przegrupować względem tychże współczynników, by później było można napisać:
![]() ![]() |
(6.30)
|
W końcu na podstawie obliczeń (6.30) możemy napisać całkę Fouriera, którego postać piszemy względem współczynników A(ω) i B(ω) (powyżej było przyjęte, że B(ω)=0, a poniżej jest dla różnego od zera dla ogólności), której ogólna postać naszej całki:
![]() |
(6.31)
|
Na podstawie definicji współczynnika An (6.26) i Bn(6.25), a także z własności ω1T1=2π, a także będziemy przyjmować t0=-∞, a T0=∞, wtedy możemy napisać definicję współczynników A(ω) i B(ω) przy wcześniejszych poczynionych uwagach:
|
|
[edytuj] Całki Fouriera dla tłumionego oscylatora tłumionego
Wychylenie od stanu równowagi dla tłumionego oscylatora od stanu równowagi przy zerowym przesunięciu fazowym przyjmuje postać:
![]() |
(6.34)
|
Wtedy patrząc na tożsamość (3.10) możemy napisać tożsamość wiążąca częstotliwość drgań tłumionych ω1 względem drgań harmonicznych podstawowych i współczynnika tłumionego, z którego skorzystamy w dalszym kroku obliczeń. Policzmy tutaj dwie całki poniższe, które będą nam bardzo potrzebne w dalszych rozważaniach:
![]() ![]() |
(6.35)
|
Dalej policzmy drugą z kolei całkę:
![]() ![]() |
(6.36)
|
Policzmy teraz współczynnik A(ω) według wzoru (6.32) występującego jako współczynnik we wzorze (6.31) korzystując przy tym z tożsamości fizycznej (6.36):
![]() ![]() ![]() ![]() ![]() |
(6.37)
|
Dalszym krokiem jest policzenie B(ω) według wzoru (6.33) występujące jako współczynnik we (6.31) wykorzystując przy tym z tożsamości fizycznej (6.34):
![]() ![]() ![]() ![]() ![]() |
(6.38)
|
Wyznaczmy wielkość I(ω), która jest sumą kwadratów wielkości 2πA(ω)(6.38) i 2πB(ω)(6.38):
![]() ![]() ![]() ![]() |
(6.39)
|
[edytuj] Całka Fouriera dla fal biegnących w ośrodku dyspersyjnym jednorodnym
Całka Fouriera dla fal biegnących w ośrodku dyspersyjnym jednorodnym dla fal dyspersyjnych dla z=0 piszemy w postaci zależnej od częstotliwości kołowej:
![]() |
(6.40)
|
Dla k≠0 całkę Fouriera możemy przestawić w zależności od liczby falowej, która z kolei zależy od liczby falowej dla fal biegnących, co można uzyskać zamieniając zmienną "t" na primowane, i podstawiając dalej za t' w tożsamości (6.40) wyrażenie t-z/v i pamiętając, że ψ(0,t')=ψ(z,t):
![]() |
(6.41)
|
Wykorzystajmy definicję prędkości fazowej, którego definicja jest podana w punkcie (4.11), wtedy równość (6.41) dla liczby falowej zależnej od częstotliwości kołowej możemy przepisać do:
![]() |
(6.42)
|
[edytuj] Klasyczne równanie falowe a fale ulegające dyspersji
Załóżmy, że mamy wychylenie od stanu równowagi, którego równanie piszemy przy pomocy częstotliwości drgań ω, liczby falowej zależnej od częstotliwości kołowej k(ω):
![]() |
(6.43)
|
Równanie falowe różniczkowe, którego rozwiązaniem jest (6.43), jest w postaci zależnej od prędkości fazowej, którego z kolei zależy od częstotliwości kołowej, możemy otrzymać licząc drugie pochodne wychylenia naszego przemieszczenia względem czasu i położenia:
|
|
Równości (6.44) i (6.45) porównujemy je względem siebie, w ten sposób otrzymujemy wzór na równanie falowe zależne od drugich pochodnych przemieszczenia ψ względem czasu i położenia "z":
![]() |
(6.46)
|
Do równości (6.46) wykorzystamy definicję prędkości fazowej (4.11), w ten sposób otrzymujemy równanie różniczkowe, które jest zależne od prędkości fazowej, a która z kolei zależy od częstotliwości kołowej lub mówiąc równoważnie od liczby falowej, bo we związku dyspersyjnym częstotliwość jest zależna od liczby falowej:
![]() |
(6.47)
|
[edytuj] Fale dwu- i trójwymiarowe rozchodzące się w przestrzeni
Dotychczas rozważaliśmy fale opisując je w jednym kierunku. Wprowadzimy fale biegnące trójwymiarowe opisując je jako fale płaskie odpowiednio obracając układ współrzędnych. Przekonamy się, że przejście do trzech wymiarów daje nam coś więcej niż zwykła zamianę zmiennych dla fal płaskich rozchodzącej się względem starego układu wzdłuż osi "z", na podstawie tego otrzymujemy dodatkowe stopnie swobody. Dla trzech wymiarów możemy mieć falę elektromagnetyczną, którą dla jednego kierunku mamy czystą falę biegnącą.
[edytuj] Wektor propagacji dla płaskich fal harmonicznych
Fala płaska rozchodząca się w jednym wymiarze, która jest falą harmoniczną dla z=0, jest funkcją o stałej amplitudzie A i częstotliwości kołowej ω, którego równanie fali dla tego "z" piszemy:
![]() |
(7.1)
|
Gdy we wzorze (7.1) podstawimy t'=t-z/v, wtedy otrzymamy fale płaską określoną w płaszczyźnie z:
![]() |
(7.2)
|
Iloczyn liczby falowej i położenia z, czyli kz, możemy napisać dla układu obróconego względem tego układu o pewien kąt:
![]() |
(7.3)
|
W ogólnym przypadku iloczyn wektora propagacji i wektora położenia możemy zapisać względem jego współrzędnych w postaci:
![]() |
(7.4)
|
Równanie fali płaskiej rozchodzącej się w trójwymiarowej przestrzeni, nazywamy wychylenie od stanu równowagi, wykorzystując przy tym fakt (7.3), jest określane przez:
![]() |
(7.5)
|
[edytuj] Prędkość fazowa
Funkcję fazową napisaną dla fali płaskiej rozchodzącej się w kierunku
nazywamy funkcję zależną od częstotliwości kołowej i liczby falowej:
![]() |
(7.6)
|
Będziemy się tutaj zajmować przypadkiem, gdy φ jest funkcją stałą:
![]() |
(7.7)
|
Interesuje nasz przypadek, gdy
, bo fala rozchodzi się zgodnie z kierunkiem wektora propagacji, zatem w takim przypadku prędkość fazową definiujemy jako stosunek częstotliwości kołowej ω i liczby falowej k, która jest długością wektora propagacji.
![]() |
(7.8)
|
[edytuj] Przykłady fal biegnących w trójwymiarze
Fala elektromagnetyczna rozchodząca się w przestrzeni, której częstość kołowa jest zależna od składowych wektora propagacji, posiada związek dyspersyjny:
![]() |
(7.9)
|
Dla fal rozchodzących się w ośrodku dyspersyjnym z prędkością "v", w którym możemy zapisać jako stosunek prędkości światła i stałej załamania, w którym to możemy zastąpić w (7.8) "c" przez prędkość "v", wtedy mamy związek dyspersyjny dla tej rozważanej prędkości:
![]() |
(7.10)
|
Rozchodząca się fala w jonosferze, która ma częstość podstawową równą ωpl, dla którego częstość kołowa jest zależna od tej częstości i od wektora propagacji, ma związek dyspersyjny:
![]() |
(7.11)
|
[edytuj] Fale stojące
Fale biegnące w dwóch osobnych kierunkach i mające takie same amplitudy i częstotliwości kołowe, w wyniku nakładania się takich fal, otrzymamy związek:
![]() |
(7.12)
|
Pamiętając, że
, to stosując odpowiednie przekształcenia trygonometryczne możemy napisać na podstawie (7.12) związek:
![]() |
(7.13)
|
[edytuj] Trójwymiarowe równanie falowe według klasycznego równania fali
Wyznaczmy teraz drugą pochodną cząstkową równania fali (7.5) względem czasu:
![]() |
(7.14)
|
A teraz policzmy drugie pochodne wspomnianego równania fali względem położeń w trzech możliwych sposobach, tzn. względem trzech możliwych współrzędnych:
|
|
|
Całkowite równanie falowe na podstawie związku (7.14), (7,16) i (7.17), przy wykorzystaniu definicji długości wektora propagacji, którego długość jest liczbą falową, przy definicji prędkości falowej (7.8), i definicji laplasjanu, czyli operatora Δ, jest napisane według:
![]() |
(7.18)
|
[edytuj] Przykłady równań fali
Równanie falowe dla fali elektromagnetycznej rozchodzącej się w próżni, przy wniosku, że prędkość fazowa fali jest równa prędkości światła w próżni równej "c", jest:
![]() |
(7.19)
|
Dla fal elektromagnetycznych rozchodzącego się w ośrodku dyspersyjnym należy wprowadzić równania falowe, dla prędkości rozchodzącej się fali równej c/n, dla której równanie różniczkowe falowe zapisujemy w formie:
![]() |
(7.20)
|
Dla fal rozchodzących się w jonosferze równanie falowe jest natomiast zależne dodatkowo od częstotliwości podstawowej ωpl, która jest przedstawiona w postaci:
![]() |
(7.21)
|
[edytuj] Fale elektromagnetyczne rozchodzące się w falowodzie
Równanie falowe dla przestrzeni dwuwymiarowej, w której rozchodzą się fale elektromagnetyczne (7.19) określanych względem współrzędnych igrekowej i zetowej, jest napisane według:
![]() ![]() |
(7.22)
|
Całkowite równanie fali elektromagnetycznej rozchodzącej się w falowodzie, określamy jako funkcję y i z, jest:
![]() |
(7.23)
|
Związek (7.23) na podstawie (7.22) zachodzi, gdy mamy związek dyspersyjny:
![]() |
(7.24)
|
Funkcja fali (7.23) na obu końcach naszego falowodu jest tak skonstruowane by w tych punktach była równa zero, co dla y=0 zachodzi automatycznie, ale znów dla y=b już nie, aby to zachodziło musi być spełniony warunek:
![]() |
(7.25)
|
Będziemy teraz rozpatrywać częstość progową, czyli tzn. dla której zachodzi m=1, wtedy warunek dyspersyjny (7.24) dla której zachodzi wniosek na ky, zapisujemy na w sposób:
![]() |
(7.26)
|
Widzimy, że związek dyspersyjny (7.26) ma kształt podobny do związku dyspersyjnego w jonosferze (7.11) albo do związku dyspersyjnego dla wahadeł sprzężonych ze sobą (2.86) dla małych k:
![]() |
(7.27)
|
Dla częstości kołowej poniżej częstości πc/b równanie fali przestawiamy w zależności od stałej zaniku χ:
![]() |
(7.28)
|
Związek dyspersyjny dla równania fali (7.28), który można otrzymać podstawiając go do (7.19), otrzymujemy wtedy inny wzór do (7.26), który jest w postaci:
![]() |
(7.29)
|
[edytuj] Fale biegnące krzyżujące się nawzajem
Równanie fali (7.23) możemy zapisać jako superpozycja fal biegnących o takich samych częstotliwościach, ale różnych liczbach falowych, którego schemat krzyżowania się tychże fal przestawiana jest jako:
![]() |
(7.30)
|
dla której we wzorze (7.30) mamy związek dla wektorów propagacji:
![]() |
(7.31)
|
Krzyżowanie się fal wynika z tego, że
i
dla związków (7.31) mają przeciwne skierowane składowe igrekowe tychże omawianych wektorów propagacji.
[edytuj] Prędkość fazowa i grupowa fal rozchodzących się fal falowodzie
Ilustracja dla fal biegnących w falowodzie pozwala na ustalić związek pomiędzy prędkością grupową a fazową w wyniku nakładania się dwóch fal krzyżujących się ze sobą. Odpowiednik
znosi składową igrekową z
z tą sama składową, ale oba te wektory propagacji mają tą sama składową zetową. Podczas gdy promień przebywa odległość ct, w tym czasie czoło fali przebyło odległość równą przy wykorzystaniu na własnościach o trójkącie prostokątnym:
![]() |
(7.32)
|
Podobnie mówimy dla prędkości grupowej uzmysławiając sobie, że zachodzi:
![]() |
(7.33)
|
Wyznaczmy teraz iloczyn prędkości fazowej i grupowej, którego przestawienie mówi, że ta wielkość jest równa kwadratowi prędkości światła:
![]() |
(7.34)
|
Drugi bardzo ważny związek wynika skorzystania z udowodnionego związku (7.34):
![]() |
(7.35)
|
jeśli do związku (7.24) podstawimy związek (7.25), w ten sposób otrzymujemy związek dyspersyjny w zależności do liczby dyskretnej dodatniej całkowitej m i współrzędnej zetowej wektora propagacji:
![]() |
(7.36)
|
Co potwierdza zgodność (7.35) z (7.36).
[edytuj] Odbijanie światła i jego przepuszczanie padającego ze szkła do powietrza
Określmy teraz dwa ośrodki, którym pierwszym jest powietrze a drugim szkło, szkło rozciąga się z minus nieskończoności do zera, a druga część przestrzeni stanowi próżnia. Wartości wektorów propagacji określamy z definicji prędkości fazowej dla próżni dla światła w sposób:
|
|
Współrzędne igrekowe wektorów propagacji są sobie równe w próżni i w szkle, i je zapisujemy na w sposób:
![]() |
(7.39)
|
Związek otrzymany w punkcie (7.39) podstawiamy do związku na całkowitą wartość wektora propagacji dla przestrzeni dwuwymiarowej, z którego możemy otrzymać wzór na kwadrat współrzędnej zetowej wektora propagacji:
![]() |
(7.40)
|
Gdy promienia załamanego nie ma, wtedy kwadrat współrzędnej zetowej wektora falowego dla fali w próżni jest równy zero. Ten kąt przy którym to pojawia się jest zwany kątem granicznym padania całkowitego odbicia wewnętrznego i przejawia się on według (7.40):
![]() |
(7.41)
|
[edytuj] Wyprowadzenie prawa Snelliusa
Weźmy sobie dwa dowolne ośrodki, w którym pierwszy jest ośrodek o współczynniku załamania n2, a drugi o współczynniku załamania n1. Mając na myśli, że igrekowe współrzędne promienia przed i po załamaniu w zależności od całkowitej długości liczby falowej przestawiamy jako:
|
|
Porównując te dwie wielkości na współrzędne igrekowe dla ośrodka pierwszego jak i drugiego, które się nie zmieniają, to możemy określić pewien związek pomiędzy kątami padania i załamania:
![]() |
(7.44)
|
[edytuj] Całkowite wewnętrzne odbicie
Rozważmy przypadek, gdy we wzorze (7.40) wyrażenie to przyjmuje wartość ujemną dla pewnych katów większych od kąta granicznego, zatem weźmy sobie na tyle duży kąt by stała
miała wartość ujemną, wtedy należy napisać:
![]() |
(7.45)
|
Wtedy związek (7.40), który jest napisany w zależności od kąta padania na granicy dwóch ośrodków θ1, na postawie (7.45), jest:
![]() |
(7.46)
|
Fala elektromagnetyczna rozchodząca się w próżni jest fala zanikającą, która opisywana jest przy pomocy funkcji proporcjonalnej do funkcji wykładniczej:
![]() |
(7.47)
|
Widzimy, że wtedy funkcja falowa jest falą biegnącą wzdłuż osi igrekowej, a falą wykładniczą wzdłuż osi zetowej. Patrząc na wzór (7.47), to średnia gęstość energii jest wprost proporcjonalna do jego średniej czasowej kwadratu funkcji falowej.
[edytuj] Teoretyczny opis fal na wodzie
Zawsze nas interesowały fale na wodzie. Pominiemy niektóre właściwości rzeczywistej wody, miedzy innymi jej lepkość, dla której Prof. Richard Feynman określa je mianem suchej wody. Przemieszczenie na wodzie poszczególnych punktów, które będziemy opisywać jej jako położenie poszczególnych elementów wody zależne od położenia i czasu, określamy je przez:
![]() |
(7.48)
|
Chwilowa prędkość wody o współrzędnych x, y w stanie równowagi, która jest określana jako pochodna cząstkowa przemieszczenia (7.48), jest pisana:
![]() |
(7.49)
|
Lokalne prawo zachowania masy mówi, że suma dywergencji iloczynu prędkości przez gęstość i pochodnej cząstkowej gęstości cząstkowej względem czasu jest równa zero, a my będziemy rozpatrywać, gdy gęstość wody nie zmienia się w czasie i też w przestrzeni, wtedy dywergencja wektora prędkości powinna być równa zero, co do tak uzyskanego prawa podstawiamy (7.49):
![]() |
(7.50)
|
Warunek braku pęcherzyków we wzorze (7.50) indukuje, że dywergencja wektora przemieszczenia poszczególnych elementów wody indukuje natomiast fakt:
![]() |
(7.51)
|
Brak wirów w wodzie indukuje wskazuje, że rotacja wektora prędkości (7.49) pokazuje, a co za tym idzie, że rotacja przemieszczenia jest wielkością stałą i równą zero, tzn. że powinno zachodzić przy braku pęcherzyków :
![]() |
(7.52)
|
Przyjmijmy, że funkcje ψx i ψy, które są przemieszczeniami w zależnymi od czasu i położenia iksowego, która względem tych argumentów zmieniają się w sposób sinusoidalny lub kosinusoidalny, przedstawiają się:
|
|
[edytuj] Zależności pomiędzy ruchem poziomym a pionowym
Wzory (7.53) i (7.54) możemy podstawić do (7.51) i (7.52), w ten sposób otrzymujemy warunki w postaci dwóch równań różniczkowych względem dwóch nieznanych funkcji, które chcemy wyznaczyć:
|
|
Równanie (7.55) różniczkujemy obustronnie względem "y" i wykorzystujemy do niego równość (7.56), otrzymujemy równanie drugiego rzędu, którego rozwiązanie podamy w jednej linijce:
![]() |
(7.57)
|
Końcowe rozwiązanie f(y) podane w punkcie (7.57) podstawiamy do wzoru różniczkowego za f(y) (7.55), w ten sposób otrzymujemy równość na rozwiązanie g(y):
![]() |
(7.58)
|
[edytuj] Warunki brzegowe na dnie zbiornika wodnego
Na dnie zbiornika wodnego nie ma ruchu pionowego wody, zatem powinno zachodzić dla y=-h warunek ψy(-h)=0, co indukuje od razu warunek brzegowy zależności stałej B od stałej A, który piszemy wzorem
, co daje nam równania na przemieszczenia względem osi iksowej i igrekowej, co je piszemy wzorami:
|
|
Gdy mamy do czynienia z głęboką wodą, wtedy przemieszczenia (7.59) i (7.60) możemy napisać po pominięciu wyrazu z e-2kh, w ten sposób otrzymujemy wzory na ψy i na ψx:
|
|
Gdy mamy doczynienia z wodą płytką, wtedy należy pominąć wyrazy z dokładnością do wyrazów liniowych względem przestawienia
, w ten sposób dostajemy wniosek:
|
|
[edytuj] Związek dyspersyjny dla fal grawitacyjnych rozchodzącej się w wodzie
Całkowita siła działając na nieskończenie mały element objętości wody o masie Δm, która działa na powierzchni piszemy z definicji siły pochodzącej od ciśnienia, jest:
![]() ![]() ![]() |
(7.65)
|
Druga zasada dynamiki Newtona dla naszego elementu wody o masie ΔM dla niezrównoważonej siły Fx wprost proporcjonalnej do drugiej pochodnej igrekowego przemieszczenia względem czasu jest określana według:
![]() |
(7.66)
|
Jeśli połączymy związki (7.65) i (7.66) przy istniejących funkcjach jako przemieszczenia iksowego (7.60) i przemieszczenia igrekowego (7.59), w ten sposób dostajemy związek dyspersyjny kwadratu częstotliwości kołowej względem liczby falowej i głębokości wody:
![]() |
(7.67)
|
Patrząc na wzór (7.67) możemy napisać wzory na związki dyspersyjne na kwadrat częstotliwości kołowych dla wody głębokiej i płytkiej:
|
|
|
||||
|
|
[edytuj] Fale biegnące na wodzie
Rozważmy teraz fale na wodzie, które przestawiamy względem funkcji sinus i kosinus przy argumencie ωt-kx, których definicja przemieszczenia igrekowego i iksowego przestawiamy koleino:
|
|
Związki na f(y) i g(y) dla ψy (7.70) i ψx (7.71) są (7.55) i (7.56), zatem dalsze obliczenia co do fal rozchodzących się po wodzie są takie same jak dla fal stojących, więc nie należy ich powtarzać.
[edytuj] Fale elektromagnetyczne w próżni
Zestaw równań obowiązujących w próżni, dla którego gęstość objętościowa ładunku i gęstość objętościowa prądu elektrycznego jest równa zero, określamy wedle:
|
|
|||||||||||
|
|
|
|
|
||||||||
|
|
|
|
||||||||
Rozpatrzmy teraz równanie fali elektromagnetycznej rozchodzącej się w próżni, dostajemy gdy równość (7.75) zróżniczkujemy obustronnie cząstkowo względem czasu, i do niego podstawiając wzór (7.74), otrzymujemy:
![]() |
(7.76)
|
Tutaj skorzystamy z tożsamości różniczkowej, które tutaj nie będziemy udowodniać i ją przestawiamy w postaci:
![]() |
(7.77)
|
Równość (7.76) możemy napisać na podstawie zachodzącej tożsamości różniczkowej (7.77) i pierwszego równania różniczkowego Maxwella (7.72) w postaci:
![]() |
(7.78)
|
Równanie (7.78) jest równoważne trzem równaniom różniczkowym, które przepisujemy trzema związkami, które są równaniami falowymi:
|
|
|
Równania (7.79), (7.80) i (7.81) są to klasyczne równania falowe dla każdej współrzędnej z osobna. Podobne związki możemy otrzymać, gdy mamy doczynienia z falami magnetycznymi, które też otrzymujemy w postaci trzech związków. Całkowity wektor fali elektrycznej i magnetycznej, gdy mamy poszczególne współrzędne, które są równaniami falowymi, przestawiamy:
|
|
Z warunku (7.72) wynika, że składowa zetowa natężenia pola elektrycznego nie zależy od współrzędnej zetowej, czyli stąd wynika, że mamy doczynienia z falami poprzecznymi, napiszemy teraz dowód, która mówi, że fala (7.82) nie zależy od współrzędnej iksowej i igrekowej, co matematycznie zapisujemy:
![]() |
(7.84)
|
Z (7.84) wynika, co chcieliśmy wykazać. Następnie wykażmy czy pole magnetyczne jest polem poprzecznym, w tym celu należy wykorzystać drugie prawo Maxwella (7.73) i z niezależności współrzędnych wektora wektora indukcji magnetycznej (7.82), z której wynika od razu poprzeczność naszej fali magnetycznej. Wykażemy teraz, że zetowa współrzędna fali elektromagnetycznej nie zależy od czasu dla pola elektrycznego, w tym celu należy wykorzystać tożsamość (7.75) dla fali magnetycznej zapisaną w punkcie (7.82) w postaci pewnego wektora, stąd wnioskujemy, że współrzędna zetowa fali elektrycznej nie zależy od czasu, a matematyczny zapis tego dowodu jest:
![]() |
(7.85)
|
Co z czwartego różniczkowego prawa Maxwella wynika nasz upragniony dowód o niezależności współrzędnej zetowej od czasu. Podobnie wykazujemy, że współrzędna zetowa pola magnetycznego też nie zależy od czasu, w tym celu należy skorzystać ze wzoru (7.74), co tutaj dowodzimy jak w punkcie (7.85).
[edytuj] Polaryzacja liniowa i eliptyczna pola elektromagnetycznego
Współrzędne pola elektrycznego nie są sprzężone równaniami różniczkowymi Maxwella Ex i Ey, one są niezależne od siebie, tzn. możemy wytworzyć je w taki sposób by Ex było zawsze różne od zera, a składowa Ey było zawsze równe zero, mówimy wtedy, że fala elektryczna jest spolaryzowana w kierunku wektora
. Układ jest spolaryzowany kołowo, gdy natężenia współrzędnych fali elektrycznej spełniają zależności:
|
|
Widzimy, że według przestawienia współrzędnych iksowej (7.86) i igrekowej (7.87) pola elektrycznego zmieniają się z częstością kołową ω, co takie przestawienie składowych pola elektrycznego nazywamy polaryzacją kołową, to te współrzędne pola możemy powiązać równaniem:
![]() |
(7.88)
|
Natomiast, gdy pole jest związane zależnościami dla drgających składowych fali pola elektrycznego z częstością kołową ω:
|
|
wtedy na podstawie zależności współrzędnej iksowej (7.89) i igrekowej (7.90) pola elektrycznego, to takie przestawienie składowych pola elektrycznego nazywamy polaryzacją eliptyczną, możemy je powiązać równaniem:
![]() |
(7.91)
|
Widzimy, że polaryzacja kołowa przestawiona wzorem (7.88) jest szczególnym przypadkiem polaryzacji eliptycznej (7.91) dla E0x=E0y=E0.
[edytuj] Biegnące Fale elektromagnetycznego harmoniczne
Przyjmijmy, że mamy falę biegnącą harmoniczną przestawiającą drgania pola elektrycznego w czasie drgającą z częstością kołową ω i posiadającej liczbę falową równą "k" przestawianą przy założeniu, że ta fala jest fala liniową.
![]() |
(7.92)
|
Biorąc warunki koleino (7.75) i (7.76), wtedy przy pomocy równania fali (7.92) przestawiających falę elektryczną biegnącą harmoniczną, które możemy wykorzystać do przestawienia poniższych równości:
![]() |
(7.93)
|
![]() |
(7.94)
|
Widzimy na podstawie wzorów (7.93) i (7.94), że pochodna By względem czasu "t" i położenia "z" jest wprost proporcjonalna do odpowiednich pochodnych fali elektrycznej Ex, zatem wnioski stąd wynikające:
|
|
|
W powyższych rozważaniach pola elektryczne i magnetyczne biegły w zgodnych kierunkach zgodnie z osią "z", którego równania dla składowej pola magnetycznego igrekowego powstaje po podzieleniu Ex przez prędkość światła, co (7.92) dla pola elektrycznego razem z równania fali pola indukcji magnetycznej napisaną poniżej jest zgodne z równaniami Maxwella:
![]() |
(7.98)
|
Ale również możemy rozpatrywać falę By, dla której ta fala biegnie w kierunku przeciwnym do osi "z", którego równanie fali jest ze znakiem przeciwnym niż u (7.92) i to równanie fali wyrażamy według:
![]() |
(7.99)
|
to równanie razem z (7.92) również jest zgodne z równaniami Maxwella. Oba te kierunki są uwzględnione od razu w równaniach (7.95), (7.96) i (7.97).
[edytuj] Stojąca fala elektromagnetyczna
Zajmijmy się teraz falami elektromagnetycznymi, zatem w nim występująca fala elektryczna może się rozchodzić w dwóch przeciwstawnym kierunkach, tzn. z kierunkiem +z i kierunkiem -z, wtedy suma amplitud tejże fali, dla której wykorzystamy wiadomości z trygonometrii, jest opisywana:
![]() |
(7.100)
|
Wektor indukcji magnetycznej jest napisany w zależności od natężenia pola elektrycznego, którą piszemy w postaci iloczynu wektorowego wektora jednostkowego w kierunku propagacji fali przez natężenie pola elektrycznego i ta całość podzielona przez prędkość światła w próżni "c".
![]() |
(7.101)
|
Dla fali magnetycznej w zależnej od fali elektrycznej biegnącej zgodnie lub niezgodnie z osią zetową, to ich superpozycja jest wyrażona:
![]() |
(7.102)
|
Jeśli patrząc na wzór (7.102), to można powiedzieć, że współrzędne igrekowe fali magnetycznej lecącą w kierunku +z odpowiadające fali elektrycznej iksowej lecącą w kierunku +z, a także odpowiadające współrzędne fali magnetycznej igrekowej lecącą do -z, która odpowiadają iksowej fali elektrycznej lecącą zgodnie -z, jeśli te dwie składowej pola magnetycznego lecącego zgodnie i odwrotnie do zwrotu osi zetowej dodamy je do siebie, w ten sposób otrzymujemy stojącą fale magnetyczną w postaci:
![]() |
(7.103)
|
[edytuj] Strumień rozchodzenia się energii fali płaskiej przez pewną powierzchnię
Gęstość energii występujących w danym punkcie nieskończenie małym w próżni jest wyrażona w zależności od natężenia pola elektrycznego i indukcji pola magnetycznego poprzez wzór:
![]() |
(7.104)
|
Spróbujmy wyrazić energię występującym w nieskończenie małym elemencie powierzchni prostopadłej do osi z i grubości nieskończenie małej Δz, wtedy ta wielkość jest wyrażona przez:
![]() |
(7.105)
|
Wielkość (7.105) zróżniczkujmy obustronnie względem czasu wykorzystując definicję o pochodnej iloczynu dwóch funkcji, w ten sposób otrzymujemy wyrażenie mówiące ile jest energii dostarczanej do objetości nieskończenie małej w ΔzA w jednostce czasu.
![]() |
(7.106)
|
Teraz wykorzystajmy tożsamości (7.93) i (7.94), by wyliczyć wartość energii wpływającej do powierzchni A o grubości Δz.
![]() ![]() ![]() |
(7.107)
|
Współrzędna zetowa ilości przekazywanej energii przez powierzchnię A w jednostce czasu i jednostce powierzchni jest to iloczyn wektorowy wektora natężenia fali elektrycznej przez wektor natężenia fali magnetycznej, jest wyrażona na podstawie (7.107) jako:
![]() ![]() |
(7.108)
|
Wektorem Poyntinga nazywamy wektor, którego wartość jest ilością energii przekazywanej w jednostce czasu przez jednostkę powierzchni, a kierunek jest prostopadły do tej powierzchni, a zwrot jest na zewnątrz jej.
![]() |
(7.109)
|
Jeśli wyrazić współrzędną iksową natężenia pola elektrycznego i igrekową współrzędną indukcji pola magnetycznego przy pomocy funkcji kosinus, w ten sposób mamy dwa równania:
|
|
Gęstość energii (7.104), u której występuje wektor natężenia pola elektrycznego, który ma tylko składową niezerową iksową, a dla pola magnetycznego posiada tylko składową igrekową, jest:
![]() |
(7.112)
|
Ilość energii całkowitej wypływającej w jednostce czasu przez jednostkę powierzchni jest narysowana przez:
![]() |
(7.113)
|
[edytuj] Gęstość i strumień energii dla elektromagnetycznej stojącej
Mając wzory dla fali stojącej pola elektrycznego (7.100) i pola magnetycznego (7.103), to wtedy gęstość energii pola elektrycznego i gęstość energii pola magnetycznego mają maksima w odstępach 1/4 okresy i w odstępach położenia równej 1/4 długości fali. Energia pola elektrycznego oscyluje z częstością kołową równej 2ω od zera do podwojonej wartości średniej, podobnie dzieje się z energią pola magnetycznego, a wiec energia oscyluje z tematu, w którym mamy czystą energią pola elektrycznego do tematu w którym mamy czystą energię magnetyczną, z maksimami odległymi od siebie równej 1/4λ. Przypomina to energie oscylatora harmonicznego, w której w jednej chwili mamy czystą energię potencjalną do chwili, gdy mamy czystą energię kinetyczną.
[edytuj] Ciśnienie promieniowania
Foton, który zostaje pochłonięty przez daną jednostkę powierzchni przekazuje jej pewien pęd równą pędowi cząstki, a jeśli foton całkowicie zostanie odbity od powierzchni przekazuje jej podwójny pęd równy zmianie pędu fotonu po odbiciu i przed odbiciem. Aby wyznaczyć zależność ciśnienia promieniowania promieniowania w zależności od energii danego fotonu musimy napisać siłę Lorentza w postaci:
![]() |
(7.114)
|
Weźmy sobie średnią czasową wyrażenia (7.114). Wyraz
ma średnią wartość równą zero, także ostatnio wyraz rozważanego wyrażenia jest równy zero, dzieje się to tak dlatego, że przyrost prędkości w kierunku osi z można zaniedbać. Pozostaje nam jeden wyraz, dla którego będziemy liczyli średnią powyższego wyrażenia, co możemy przepisać go w postaci:
![]() |
(7.115)
|
Z drugiej jednak strony moc wykonywana nad przesunięciem ładunku jest iloczynem siły (7.114) i prędkości:
![]() |
(7.116)
|
Ponieważ cBy=Ex, wtedy możemy połączyć równości (7.115) z (7.116), w ten sposób otrzymać związek między przekazywanym pędem a mocą wykonywaną nad ładunkiem:
![]() |
(7.117)
|
Ze wzoru (7.117) możemy wyznaczyć pęd fotonu w zależności od jego energii, który to pęd jest ilorazem energii fotonu W przez prędkość światła "c", lub inaczej energii fali W odpowiada pęd fali
:
![]() |
(7.118)
|
Energia W możemy również wyznaczyć z poniższego wzoru, która jest to energia policzona ze związku między energią a pędem ze szczególnej teorii względności dla masy jego spoczynkowej równy zero (jeśli tam podstawić M=0), co w rezultacie otrzymujemy w wyniku końcowych perypetiach równość (7.118).
![]() |
(7.119)
|
[edytuj] Moment pędu przekazywanej przez fale biegnącą fali płaskiej
Fala elektromagnetyczna posiada nie tylko energię i pęd liniowy, ale również moment pędu. Będziemy rozpatrywać fale elektryczne, które wprawiają w ruch obrotowy ładunek q, co jest niemożliwe przy polu spolaryzowanym liniowo. Rozważmy falę elektryczną o stałym natężeniu obracającego się z pewną prędkością kołową, w którym pole magnetyczne w zależności od pola magnetycznego piszemy poprzez wzór
. Pole elektryczne rozpędza ładunek q, a pole magnetyczne zakrzywia jego tor ruchu naszej cząstki. Moment skręcający
oddziałujący na ładunek q z siłą, która wynosi
, której ta wspomnianą wielkość piszemy pomnożoną przez częstotliwość kołową ω, wtedy:
![]() |
(7.120)
|
iloczyn wektorowy
ma kierunek zgodny z
, a
ma kierunek
. A ponieważ średnia czasowa prędkości w ciągu jednego okresu jest równe zero, to momentu całkowity pochodzący od pola magnetycznego nic nie wnosi do całkowitego momentu skręcającego, ale do niego wnosi tylko pole elektryczne:
![]() |
(7.121)
|
wtedy średnia czasowa momentu skręcającego w ciągu jednego okresu jest wyrażona przez wyrażenie:
![]() |
(7.122)
|
Moment pędu fali elektromagnetycznej spolaryzowanej kołowo możemy wyrazić od energii fali W na podstawie (7.122):
![]() |
(7.123)
|
[edytuj] Fale elektromagnetyczne rozchodzące się w ośrodku jednorodnym liniowym
Prędkość światła w zależności od prędkości światła w próżni i względem stałych przenikalności elektrycznej i magnetycznej ośrodka liniowego wyrażamy wzorem:
![]() |
(7.124)
|
Wtedy kwadrat liczby falowej jest równy stosunkowi kwadratów częstotliwości kołowych drgań przez kwadrat prędkości światła pomnożonej przez iloczyn względnych przenikalności elektrycznej εr i magnetycznej μr ośrodka liniowego:
![]() |
(7.125)
|
[edytuj] Promieniowanie pochodzące od przyspieszającego ładunku
Załóżmy, że ładunek q znajduje się w chwili początkowej w spoczynku od czasu t=-∞ do t=0, wtedy pole elektryczne tego ładunku jest to pole opisywane przez prawo Coulumba. W chwili t=0 cząstka doznaje chwilowego przyspieszenia przez króŧki czas Δt, a później porusza się ze stała prędkością v=aΔt. W chwili t=0 pole elektryczne dane jest przez prawo Coulumba:
![]() |
(7.126)
|
Określmy teraz czas t długi w porównaniu z czasem Δt. Gdy odległość r jest mniejsza niż c(t-Δt) od ładunku q załamanie już mineło i pole elektryczne jest wtedy wywoływane przez ładunek poruszających się z prędkością v, jest określone jest przez wzór w zależności od kąta θ, w którym jest to kąt pomiędzy wektorem wodzącym punktu, w którym następuje obserwacja, którego początek ma w punkcie położenia chwilowego ładunku q, względem prędkości źródła pola elektrycznego:
![]() |
(7.127)
|
Gdy prędkość ładunku jest o wiele mniejsza niż prędkość ładunku, czyli v<<c, wtedy pole elektryczne jest opisywane przez (7.126). Oznaczmy przez
i
, które są to wartości bezwzględne składowych
koleino prostopadłej i równoległej do kierunku wektora
, które jest polem zajętym przez załamanie, pochodzących od ładunku, który poruszał się z przyspieszeniem. Prawo zachowania strumienia linii sił pola elektrycznego pociąga ze sobą oczywiście ciągłość linii sił i dlatego stosunek poprzecznej składowej pola elektrycznego
do składowej równoległej (podłużnych) da się przestawić jako stosunek odległości
przez
, który piszemy poprzez:
![]() |
(7.128)
|
Ponieważ pole elektryczne jest opisywane przez równanie (7.126), również przez ten sam wzór jest opisujemy składową równoległą (radialna) pola elektrycznego, tzn.:
![]() |
(7.129)
|
Możemy połączyć równanie (7.128) z równaniem na składową radialną pola elektrycznego (7.129), w ten sposób otrzymujemy równość na składową prostopadłą pola elektrycznego:
![]() |
(7.130)
|
Całkowity wektor natężenia pola elektrycznego promieniowania określamy w zależności od wektora
jako:
![]() |
(7.131)
|
[edytuj] Energia promieniowania wysyłana przez ładunek punktowy
Wiemy, że natężenia pola elektrycznego jest prostopadłe do wektora jednostkowego propagacji fali elektrycznej
, i wykorzystując definicję wektora Poyntinga, który jest natężeniem wysyłanym przez ładunek punktowy w danym punkcie na jednostkę powierzchni i jednostkę czasu, to wtedy możemy napisać wzór na ten wektor:
![]() |
(7.132)
|
Kwadrat przyspieszenia prostopadłego jest wyrażona od całkowitego przyspieszenia z definicji trójkąta prostokątnego w postaci:
![]() |
(7.133)
|
Tożsamość (7.133) możemy podstawić do (7.132), w ten sposób mamy wzór na różniczkę mocy przekazywaną w jednostce czasu przez powierzchnię dS poprzez promieniowanie:
![]() |
(7.133)
|
Całkowita moc przekazywana w jednostce czasu przez powierzchnię zamkniętą wyrażamy poprzez całkę po całym przebiegu zmienności θ i φ, otrzymujemy:
![]() |
(7.134)
|
[edytuj] Całkowite promieniowanie dipola elektrycznego
Załóżmy, że ładunek q porusza się ruchem harmonicznym, którego funkcję położenia i przyspieszenia podamy w jednej linijce:
|
|
Całkowita moc promieniowania dipola elektrycznego P, która jest zależna od częstotliwości kołowej drgań, która jest wysyłana do nieskończoności, wyliczamy go podstawiając (7.136) do (7.134), w ten sposób:
![]() |
(7.137)
|
[edytuj] Szerokość naturalna promieniowania linii atomu dla światła
Zależność całkowitej średniej energii w zależności do czasu wyrażamy przy pomocy wzoru (3.19), zatem możemy udowodnić tożsamość, że pochodna zupełna energii układu względem czasu wziętej z minusem i po podzieleniu przez energie układu równa jest odwrotności współczynnika zaniku drgań τ, czyli współczynnikowi tłumienia Γ:
![]() |
(7.138)
|
Całkowita energia średnia oscylatora harmonicznego tłumionego dla słabego tłumienia zależnego od czasu (wtedy ω≈ω0) możemy napisać poprzez:
![]() |
(7.139)
|
Całkowita moc promieniowania wysyłana przez pewną powierzchnię zamkniętą określamy przez wzór (7.137) i na podstawie tego możemy powiedzieć, że pochodna zupełna energii układu względem czasu wziętej z minusem wyrażamy:
![]() |
(7.140)
|
Do wzoru (7.138) postawiamy wzory (7.140) i (7.139), w ten sposób otrzymujemy wzór na naturalną szerokość połówkową linii światła:
![]() |
(7.141)
|
[edytuj] Prawo Rayleigha
Weźmy sobie teraz elektron występujący w cząsteczce molekuły, który drga pod wpływem harmonicznego natężenia pola fali elektrycznej, którego współrzędna iksowa jest:
![]() |
(7.142)
|
Niech ładunek q będzie związany z w cząsteczką za pomocą sprężyny o współczynniku sprężystości
, który porusza się ruchem drgającym, wtedy równanie ruchu wynikające z drugiej zasady dynamiki jest:
![]() |
(7.143)
|
Ale ponieważ zachodzi
, co jest słuszne dla ruchu drgającego harmonicznego, wtedy równość (7.143) możemy przepisać, po zastąpieniu iksowego przyspieszenia przez iloczyn kwadratu częstotliwości kołowej drgań i położenia tej cząstki względem osi iksowej, wtedy równanie ruchu elektronu przestawiamy:
![]() |
(7.144)
|
Do wzoru (7.137) na całkowitą moc wypromieniowanej do nieskończoności podstawiamy do tożsamości (7.144):
![]() |
(7.145)
|
W omawianym przypadku, gdy częstość kołową drgań podstawowych ω0 jest o wiele większa niż częstość drgań harmonicznego pola elektrycznego, czyli ω0>>ω, co można przy pomocy tychże dysput powiedzieć:
![]() |
(7.146)
|
Wzór (7.146) jest tzw. prawem Rayleigha, która jest odwrotnie proporcjonalna do czwartej potęgi długości fali promieniowania wysyłany przez ładunek.
[edytuj] Czemu niebo jest niebieskawe
Światło w zależności od długości fali jest rozpraszane w taki sposób, że światło o większej długości fali promieniowania elektromagnetycznego jest rozpraszane słabiej niż światło o mniejszej długości, to rozpraszanie jest dane przez prawo Rayleigha, którego to prawo obowiązuje dla równości (7.146), widzimy że światło niebieskie o długości fali 440-490nm jest lepiej rozpraszane niż światło czerwone, którego długość fali jest 630-780nm, i dlatego w dzień widzimy niebo niebieskawe, a wieczorem niebo czerwonawe.
[edytuj] Polaryzacja liniowa, kołowa, eliptyczna fal
[edytuj] Wstęp do opisu polaryzacji
Dla wszystkich ruchów falowych badaliśmy przemieszczenie, które określamy przy pomocy pewnej funkcji względem położenia równowagi. W ogólnym przypadku przemieszczenie określamy jako funkcje z "x","y","z", czyli
. Zwykle badamy fale płaskie, które określamy jako funkcje ψ(z,t), gdzie z jest liczone wzdłuż rozchodzenia się fali płaskiej, często interesującymi wielkościami są
i
. Dla fal płaskich przemieszczenie będziemy określać przez wektor przemieszczenia w postaci:
![]() |
(8.1)
|
W przypadku fal poprzecznych wektor przemieszczenia
ma tylko składowe x i y. O takich falach mówimy, że one mają polaryzacje poprzeczną. Fala dźwiękowa jest to fala mająca przemieszczenie wzdłuż osi zetowej. Nie ma poprzecznych fal dźwiękowych, chociaż możemy je sobie wyobrazić jako fale odbijające, które są jako fale podłużne odbijające się od ścianek wewnętrznych rury, i w rezultacie dla fali rozchodzącej się wzdłuż rury fala dźwiękowa oprócz składowych podłużnych mają składowe poprzeczne. Fala elektromagnetyczna, czyli suma fal elektrycznych i magnetycznych, ma tylko składowe poprzeczne. Fala elektromagnetyczna ma też składowe podłużne, jeśli ona rozchodzi się w falowodzie. Falą poprzeczną nazywamy falę, dla której składowa zetowa przemieszczenia nigdy nie występuje, zatem drgania poprzeczne w analogii do (8.1) są napisane:
![]() |
(8.2)
|
[edytuj] Promieniowanie ładunku punktowego
Innym przykładem fali płaskiej są drgania wyemitowane w początku układu współrzędnych, które daleko od punktu emisji obserwujemy jako fale płaskie, których przemieszczenie ładunku q jako źródła fal opisujemy przez:
![]() |
(8.3)
|
Natężenie pola elektrycznego fali (7.134) dla odległości r=z od początku układu współrzędnych wytworzone przez ładunek poruszający się z pewnym przyspieszeniem na podstawie (8.3) piszemy:
![]() |
(8.4)
|
Biorąc współczynnik proporcjonalności
, wtedy przemieszczenie fali może być uważane za natężenie pola elektrycznego w czasie przedwczesnym
.
[edytuj] Polaryzacja liniowa
W przypadku drgań poprzecznych, którego drgania są tam i z powrotem, to te drgania można opisywać w dwóch niezależnych kierunkach:
|
|
Superpozycją drgań (8.5) i (8.6) w dwóch osobnych kierunkach nazywamy falę:
![]() |
(8.7)
|
Całkowita długość drgań amplitudy opisywanej wzorem (8.7) piszemy jako pierwiastek sumy kwadratów amplitud drgających w dwóch prostopadłych kierunkach:
![]() |
(8.8)
|
kierunkiem polaryzacji, patrząc na wzory (8.7) i (8.8), nazywamy kierunek, którego wektor jest napisany:
![]() |
(8.9)
|
łatwo sprawdzić, że wektor
jest wektorem jednostkowym, co wynika z prostopadłości wektorów
i
, które z kolei są wektorami jednostkowymi.
[edytuj] Liniowo spolaryzowaną falę biegnącą
Liniowo spolaryzowaną falę biegnącą nazywamy falę przy ustalonych płaszczyźnie drgań biegnącego w kierunku +z:
![]() |
(8.10)
|
[edytuj] Liniowo spolaryzowana fala stojąca
Załóżmy, że mamy dwie fale biegnące w kierunku +z, a druga w kierunku -z, to superpozycja tychże fal możemy opisać:
![]() ![]() |
(8.101)
|
[edytuj] Polaryzacja kołowa
Tutaj przemieszczenie fali poprzecznej ulega ruchowi kołowemu. Równanie tejże fali opisujemy przy pomocy funkcji sinus i kosinus z ωt, których drganie w kierunku
wyprzedza drganie w kierunku
o kąt 90o;
![]() |
(8.12)
|
Podobnie, gdy fala drgań w kierunku
od
wyprzedza o kąt 90o, to fala spolaryzowana kołowo jest opisywana:
![]() |
(8.13)
|
[edytuj] Spolaryzowana kołowo fala biegnąca
Spolaryzowana kołowo fala biegnącą nazywamy falę biegnącą w kierunku osi zetowej (+z), którego przepis jest zależny od częstotliwości kołowej ω i liczby falowej "k" :
![]() |
(8.14)
|
[edytuj] Spolaryzowana kołowo fala stojąca
Kołowo spolaryzowaną falę stojącą nazywamy dwie fale biegnące zdefiniowane w punkcie (8.14), którego biegną kierunkach +z(pierwsza fala) i w kierunku -z (druga fala), to wyniku superpozycji tychże fal otrzymujemy falę tzn. stojącą, który wyjdzie nam z obliczeń:
![]() ![]() ![]() |
(8.15)
|
[edytuj] Polaryzacja eliptyczna
Falę spolaryzowaną eliptycznie nazywamy falę, których amplitudy drgań w kierunkach
i
mają różne wartości amplitudy i wzory na fale stojącą dla której superpozycja dwóch fali biegnących w tych samych kierunkach, ale w innych płaszczyznach, piszemy:
![]() |
(8.16)
|
Eliptycznie spolaryzowana fale stojącą nazywamy dwie fale biegnące zdefiniowane w punkcie (8.16), którego biegną kierunkach +z(pierwsza fala) i w kierunku -z (druga fala), to wyniku superpozycji tychże fal otrzymujemy falę tzn. stojącą, który wyjdzie nam z obliczeń w sposób:
![]() ![]() ![]() |
(8.15)
|
[edytuj] Superpozycja dwóch fal o niejednakowych amplitudach
Weźmy sobie dwie fale o jednakowej częstotliwości kołowej ω i o niejednakowych amplitudach, których superpozycja jest opisywana:
![]() |
(8.18)
|
Opiszmy falę biegnącą opisywaną według (8.18), którego będziemy rysowali na wykresach poniżej, których współrzędną iksową jest jedna z tych fali, a osią igrekową jest druga z tychże fali, dla różnych ich przesunięć fazowych:
[edytuj] Liczby zespolone w opisie fal
Zastosowanie liczb zespolonych z ilustrujmy na przykładzie fali elektromagnetycznej, a właściwie dla fali pola elektrycznego, którego superpozycja fal biegnących w dwóch prostopadłych płaszczyznach jest w postaci wektora natężenia fali pola elektrycznego, gdy te dwie fale są przesunięte względem siebie o pewne przesunięcie fazowe φ1-φ2, to ich superpozycja jest:
![]() |
(8.19)
|
Równanie (8.19) możemy przestawić w postaci zespolonej, którego częścią rzeczywistą jest równanie (8.19):
![]() |
(8.20)
|
Zespolone natężenie funkcji falowej możemy napisać jak superpozycja dwóch fal o funkcjach falowych
i
:
![]() |
(8.21)
|
Funkcje falowe
i
są to funkcje, których definicje są wyrysowane w zależności od częstotliwości kołowej i liczby falowej zapisane są przy pomocy funkcji eksponencjalnej:
|
|
Amplitudy występujące w równaniu są to amplitudy zależne od przesunięć fazowych φ1 i φ2 zapisane pod funkcją eksponencjalną, jak widzimy, te amplitudy są wielkościami zespolonymi:
|
|
[edytuj] Superpozycja funkcji ortonormalnych
Funkcję, które opisywaliśmy w punkcie (8.22) i (8.23) są to funkcje ortonormalne i tworzą zupełną bazę ortonormalną zupełną bazowych funkcji falowych, które przestawiamy:
|
|
Aby sprawdzić, czy funkcję
i
zdefiniowane w poprzednim rozdziale tworzą zupełny układ bazowy należy na tych funkcjach wykonać następujące działania, które są opisywane wzorami (8.26) i (8.27), zatem do dzieła:
![]() |
(8.28)
|
![]() |
(8.29)
|
Ponieważ nasze rozważane funkcje bazowe są funkcjami ortonormalnymi, zatem zdefiniujmy moduł funkcji natężenia pola elektrycznego
zdefiniowanego w punkcie (8.21):
![]() ![]() |
(8.30)
|
[edytuj] Inne postacie ortonormalnych funkcji bazowych
Aby napisać wektor natężenia fali pola elektrycznego (8.21) w innych bazach niż tutaj obranego pierwotnie, których z oczywistych powodów jest ich nieskończenie wiele, które można tak obrać przez obrót pierwotnego układu współrzędnych o pewien kąt φ, których nowe wersory naszej bazy są napisane względem odpowiedników starej bazy w sposób:
|
|
Na podstawie przestawienia wersorów w (8.31) (8.32) dowiadujemy się, że one są do siebie ortonormalne, wtedy ich funkcje falowe są ortonormalne:
|
|
[edytuj] Przestawienie polaryzacji fali płaskiej
Falę biegnącą fali płaskiej, która jest spolaryzowana kołowo, którą możemy przedstawić w zależności od częstości kołowej i liczby falowej w postaci poniżej:
![]() ![]() |
(8.35)
|
Końcowy wynik (8.35) możemy przestawić w postaci zespolonej za pomocą funkcji eksponencjalnych:
![]() |
(8.36)
|
Wykorzystamy tutaj pewne tożsamości na liczbach zespolonych, która wynika z definicji eix jako liczby zespolonej zapisaną w postaci eksponencjalnej, które te tożsamości piszemy wzorami przy rozwinięciu jego do postaci trygonometrycznej:
|
|
Na podstawie obliczeń (8.37) i (8.38) wzór (8.36) możemy przepisać po odpowiednich jego przekształceniach do wzoru, z którego będzie jasne jak wybierzemy funkcję ortogonalne:
![]() |
(8.39)
|
Wybierzmy sobie teraz funkcję ortonormalne spełniające warunki (8.26) i (8.27):
|
|
Jeśli dodatkowo weźmiemy, że liczby A+ i A-, których definicje są takie by przy wektorach ortonormalnych (8.40) i (8.41) była spełniona tożsamość (8.39):
![]() |
(8.42)
|
Na podstawie wzoru (8.36) i przedstawienia bazy ortogonalnej, tzn. (8.40) i (8.41) przy definicji jego współczynników (8.42) równanie na natężenie fali elektromagnetycznej piszemy w postaci formalnej:
![]() |
(8.43)
|
[edytuj] Dyfrakcja jako superpozycja bardzo dużej liczby fal
Będziemy się tutaj zajmować zjawiskiem interferencji lub dyfrakcji. Zjawiskiem dyfrakcji nazywamy takie nakładanie się fal w wyniku superpozycji, która powstaje na detektorze w przypadku dojścia tam fali wytworzonych przez dwa źródła lub dwie szczeliny, które natomiast są pobudzone do wytwarzania tychże fal przez pewne prawdziwe źródło fal, jakie może być na przykład żarówka.
[edytuj] Superpozycja nakładania się dwóch lub więcej źródeł fal spójnych
Źródłami spójnymi w przypadku dwóch źródeł wytwarzających dwie fale, nazywamy falę, gdy różnica faz jest wielkością niezależną od czasu, co jest jedyne możliwe, gdy mamy doczynienia z falami o jednakowych częstotliwościach kołowych. Interferencją konstruktywną nazywamy falę wyniku nakładania się maksimów lub minimum tychże fal. A wynik interferencji dekonstruktywnej nazywamy takie nakładanie się dwóch fal, które powstaje wyniku nakładania się doliny pierwszej fali z grzbietem drugiej fali. Obraz interferencyjny nazywamy obraz, który powstaje w wyniku obszarów maksimum i minimum, który powstaje na detektorze. Jeśli będziemy rozpatrywać nakładanie się dwóch fal pochodzących od dalekich źródeł, wtedy mówimy, że detektor jest w dalekim polu źródeł.
Z rysunku obok wynika z definicji trójkąta prostokątnego:
![]() |
(9.1)
|
Korzystając ze wzoru skróconego mnożenia znanego z algebry dla różnicy kwadratów, co można napisać wzór wynikający z (9.1):
![]() |
(9.2)
|
Interesuje nas przypadek, gdy L2P i L1P są prawie równe, ale ich różnica nie przekracza
, wtedy obliczenia (9.2) można dokończyć przy naszych dysputach według schematu:
![]() |
(9.3)
|
Patrząc na wzór (9.3) można napisać, że przybliżenie dalekiego pola mamy, gdy iloczyn odległości do detektora od źródła fal przez długość fali światła jest równa w przybliżeniu kwadratowi odległości między dwoma źródłami, co piszemy:
![]() |
(9.4)
|
[edytuj] Źródło pola elektromagnetycznego fali podłużnej
Gdy fala jest wytwarzana przez oscylujący ładunek, która jak wiadomo porusza się z przyspieszeniem zależnych od czasu w postaci drgań harmonicznych tego ładunku, to fala elektromagnetyczna pola elektrycznego wykorzystując wzór (7.131) jest napisana w zależności od czasu:
![]() |
(9.5)
|
Widzimy na podstawie (9.5) dla fali pola elektrycznego, która jest falą z pewną amplitudą, jest:
|
|
Jeśli mamy N źródeł, to koleino oznaczamy natężenie fali pola elektrycznego przy pomocy wskaźnika jeden lub dwa, aż do N.
[edytuj] Różnica faz spowodowana różnicą dróg optycznych
Widzimy, że według rysunku obok różnica między promieniem r2, a r1 jest równa Δs, i jest wyrażona względem odległości dwóch szczelin i kąta pomiędzy poziomem, a promieniem pierwszym lub drugim, bo ten kąt dla promienia pierwszego i drugiego jest jednakowy:
![]() |
(9.8)
|
Różnica faz między promieniem drugim a pierwszym, bo promień drugi przebywa dłuższą drogę optyczną niż pierwszy, zapisujemy w postaci:
![]() |
(9.9)
|
[edytuj] Równanie interferencyjne dla dwóch szczelin
Weźmy sobie teraz sobie falę płaską wytwarzaną przez źródło pierwsze jak i drugie, która jest falą pola elektrycznego, przy założeniu, że te dwie fale wytwarzane są przez dwa źródła o takiej samej częstotliwości i przesunięciu fazowym, dodajmy te dwie fale do siebie w wyniku superpozycji tychże fal daleko od źródła fal, otrzymujemy stąd wniosek:
![]() ![]() |
(9.10)
|
Aby wystąpiły maksima fali w wyniku interferencji dwóch fal i w wyniku tego drugi czynnik musi być równy plus jeden lub minus jeden, co jest możliwe gdy:
![]() |
(9.11)
|
Fala napisana wzorem (9.10), aby przyjmowała wartość zero, czyli wtedy natężenie pola elektrycznego jest równe zero, to musi być spełnione:
![]() |
(9.12)
|
[edytuj] Warunek spójności dwóch źródeł
Wykażemy tutaj, że jeśli mamy źródła fal a,b,c, to wtedy działa ono jako efektywne punktowe źródło światła, jeśli jest spełniony pewien warunek, który mamy zamiar wyprowadzić. Według rysunku obok możemy powiedzieć:
![]() |
(9.13)
|
Wzór (9.13) podstawiamy do (9.11) dla n=1, w ten sposób otrzymujemy tożsamość:
![]() |
(9.14)
|
W rzeczywistości d możemy być on wiele mniejsze niż to wynika z równania końcowego (9.14), zatem nasz warunek na spójność wielu źródeł piszemy:
![]() |
(9.15)
|
[edytuj] Obraz pojedynczej szczeliny w obrazie dyfrakcyjnym
Posłużmy się teraz konstrukcją Huygensa, którą zastępujemy falę od odległego źródła falą płaską, która pada na N szczelin, w których mamy oscylujące ładunki, które wytwarzają fale kuliste, którego natężenie fali jest opisane w przestrzeni rzeczywistej, mamy:
![]() ![]() |
(9.16)
|
Superpozycję fali zapisaną w punkcie (9.16) możemy również zapisać w postaci zespolonej w przestrzeni zespolonej przy pomocy eksponentów, którego odpowiednik wspomnianego równania jest:
![]() |
(9.17)
|
Poszczególne odległości r1, r2,..,rN piszemy w zależności od ilości N szczelin, kąta padania θ i odległości między szczelinami:
![]() |
(9.18)
|
Wzory (9.18) wsadzamy do tożsamości (9.17), stąd w ostatecznych rozrachunkach otrzymujemy;
![]() ![]() |
(9.19)
|
Dalszym krokiem jest wykorzystanie wzoru na zwartą sumę szeregu geometrycznego, który ma iloraz
, w ten sposób (9.19) przyjmuje postać:
![]() ![]() |
(9.20)
|
[edytuj] Przykład ciągły nieskończenie wielu szczelin na odległości D
Weźmy sobie N szczelin, których liczba jest nieskończona przy stałym D przy d dążących do zero, zatem w takim przypadku mamy prawo powiedzieć na podstawie wzoru (9.17), że wypadkową amplitudę piszemy jako:
![]() |
(9.21)
|
Całkowite przemieszczenie na podstawie obliczeń przeprowadzonych w punkcie (9.21) i na podstawie definicji naszego przemieszczenia (9.17) jest określone:
![]() |
(9.22)
|
Według obliczeń (9.22) możemy obliczyć średni strumień energii, która tutaj po uśrednieniu po czasie mamy:
![]() |
(9.23)
|
[edytuj] Nieoznaczoność liczby falowej i położenia i dyfrakcyjna kątowa szerokość wiązki
Patrząc na rysunek powyżej szerokość połówkowa wiązki nazywamy przedział, który rozciąga się od
do
i jego wartość określona jest przez:
![]() |
(9.24)
|
Całkowita wartość liczby falowej jest tak określona, by jego kwadrat był
, zdefiniujmy w nim współrzędne iksowe kx=k cosθ, a kz=kcosθ i na samym końcu ky, która jest równa zero, zatem stosując przybliżenie małych kątów, wtedy kosinus zastąpimy przez jedynkę, a sinus kąta przez sam kąt, wtedy współrzędną iksową liczby falowej możemy zapisać w przybliżeniu wedle sposobu:
![]() |
(9.25)
|
Zatem nieoznaczoność liczby falowej i położenia w ostatecznych rozrachunkach, przy wykorzystaniu wzorów (9.23) i (9.24) zastępując D przez Δx, jest:
![]() |
(9.26)
|
Można również udowodnić podobnie jak powyżej, że nieoznaczoności zachodzą również dla współrzędnej igrekowej i zetowej, a także dla częstotliwości kołowej, które możemy napisać jako związki:
|
|
|
|
[edytuj] Optyka geometryczna
[edytuj] Odbicie od płaszczyzny fali płaskiej
Odległość dróg fazowych dwóch promieni padającej pierwszej względem promienia drugiego według rysunku obok możemy zapisać według:
![]() |
(9.31)
|
Różnica dróg fazowych między promieniem drugim a pierwszym, bo promień drugi przebywa dłuższą drogę optyczną niż pierwszy, możemy zapisać:
![]() |
(9.32)
|
Równanie fal materii pierwszego i drugiego promienia, uwzględniają przesunięcie fazowe tychże fal przed ugięciem fal, a także przesunięcie drugiej fali względem pierwszej po ugięciu, pisząc je wedle:
|
|
Według zasady Huygensa musimy dodać fale (9.33) do (9.34), które ulegają superpozycji w bardzo dużej odległości od kryształu, co można zapisać w przybliżeniu, że te dwie fale poruszają się po liniach prostych:
![]() ![]() |
(9.35)
|
W wyrażeniu (9.58) kosinus przyjmuje wartość jeden lub minus jeden, gdy wyrażeniu pod kosinusem jest
, wtedy moduł wspomnianego wyrażenia przyjmuje wartość maksymalną, gdy zachodzi zależność:
![]() |
(9.36)
|
Przyjmujemy, że λ>d, wtedy
, ale ponieważ lewa strona równania jest zawsze mniejsza niż jeden, a prawa strona jest większa niż jeden dla m>0, co stąd wynika, że jedyną możliwością jest m=0, co na podstawie (9.36) dostajemy α=β. Otrzymujemy, ze promień odbicia jest równy promieniowi padania.
[edytuj] Załamanie fali płaskiej
Odległość AB i DC możemy policzyć przy pomocy właściwości trójkąta równobocznego, które tutaj piszemy wzorami:
|
|
Różnica dróg między promieniem drugim a pierwszym, bo promień drugi przebywa dłuższą drogę optyczną niż pierwszy, możemy zapisać:
![]() ![]() ![]() |
(9.39)
|
Funkcje falowe występujące w zjawisku załamania mają taki sam wygląd jak w zjawisku odbicia, czyli funkcje (9.33) i (9.34), zatem złożenie tychże funkcji piszemy wedle:
![]() ![]() |
(9.40)
|
W wyrażeniu (9.37) kosinus przyjmuje wartość jeden lub minus jeden, gdy wyrażeniu pod kosinusem jest
, wtedy moduł wspomnianego wyrażenia przyjmuje wartość maksymalną, gdy zachodzi zależność:
![]() |
(9.41)
|
Dla dowolnie małego d występujące w tożsamości (9.41) prawa jego strona jest dowolnie duża dla m>0, a lewa skończona jest skończona, zatem jedyną możliwością jest m=0, zatem co wynika ze wspomnianego równania wniosek:
![]() |
(9.42)
|
Prawo opisane wzorem (9.40) nazywamy prawem załamania Snelliusa.
[edytuj] Zasada Fermata a prawo odbicia
Całkowity czas jaki promień leci na powierzchnię graniczną z punktu A, który pada na powierzchnię i odbija się od punktu B i leci do punktu C, jest wyrażona przy stałych parametrach y1, y2:
![]() |
(9.43)
|
Zróżniczkujmy obie strony równości (9.43) i przyrównajmy różniczkę czasu do zera, by otrzymać później prawo odbicia, wiedząc jednocześnie, że x1=L-x2.
![]() |
(9.44)
|
Ale ponieważ z definicji sinusów znanej z trygonometrii dla trójkąta prostokątnego możemy otrzymać wnioski:
|
|
Wykorzystując tożsamości trygonometryczne (9.45) i (9.46), które podstawiamy do (9.44), otrzymujemy:
![]() |
(9.47)
|
W końcowym wniosku w zasadzie Fermata (9.47) wnioskujemy równość kątów α i β, co jest treścią ostateczną zasady Fermata dla zjawiska odbicia.
[edytuj] Prawo załamania fali a zasada Fermata
Całkowity czas jaki promień leci w pierwszym ośrodku o współczynniku załamania n1 padający na powierzchnię graniczną pomiędzy dwoma ośrodkami z punktu A, który załamuje się w punkcie punkcie B i leci do punktu C w drugim ośrodku o współczynniku załamania n2, jest wyrażony przy pomocy stałych y1, y2:
![]() |
(9.48)
|
Zróżniczkujmy obie strony równości (9.48) i przyrównajmy różniczkę czasu do zera, wtedy możemy otrzymać prawo załamania w dalszych niż poniżej dysputach, wiedząc jednocześnie, że x1=L-x2:
![]() |
(9.49)
|
Ale ponieważ z definicji sinusów znanej z trygonometrii dla trójkąta prostokątnego mamy wnioski:
|
|
Wykorzystując tożsamości trygonometryczne (9.50) i (9.51), które podstawiamy do (9.49), w ten sposób otrzymujemy równość:
![]() |
(9.52)
|
Wzór końcowy wynikający z obliczeń (9.52) jest tzw. prawem Snelliusa dla zjawiska załamania.
[edytuj] Zwierciadło wklęsłe
Rozważmy sobie teraz zwierciadło wklęsłe, zbadajmy jaki jest obraz danego przedmiotu występującego przed tym obiektem. Wysokość obrazu jest h i znajduje się w odległości od zwierciadła o x. Od przedmiotu promień świetlny leci do zwierciadła względem poziomu pod kątem α, i do zwierciadła dociera ono na wysokość równą wysokości hx przedmiotu minus długość względem pionu na jaką promień się obniżył, która jest iloczynem odległości przedmiotu od zwierciadła pomnożonej przez kąt z jakim leci nasz promień względem poziomu:
![]() |
(9.53)
|
Kąt γ według rysunku obok określamy jako iloraz wysokości h1 napisaną wzorem (9.15) przez promień krzywizny zwierciadła R.
![]() |
(9.54)
|
Całkowity kąt β jest równy sumie kata α i γ, który jest kątem pomiędzy odcinkiem łączący R z promieniem świetlnym w punkcie granicznym zwierciadła, który się odbija od zwierciadła w tymże punkcie:
![]() |
(9.55)
|
Nachylenie promienia odbitego od zwierciadła względem poziomu określamy jako sumę katów β i γ, który jest jednocześnie ilorazem sumy wysokości powstałego obrazu i wysokości nad osią główną zwierciadła, od której to wysokości odbija się promień świetlny wychodzący od przedmiotu:
![]() |
(9.56)
|
Załóżmy, że mamy dwa wzory (9.56), dla różnych α, tzn. dla α1 i α2 i odejmijmy od siebie te wzory i tak otrzymane równanie podzielmy przez α1-α2, otrzymując w ten sposób tożsamość:
![]() |
(9.57)
|
Na podstawie obliczeń przeprowadzonych powyżej wzór na odległość ogniska od zwierciadła jest równa połowie promienia krzywizny zwierciadła, a także wzory na zwierciadło są w takim razie opisywane przez:
|
|
Wzór (9.58) podstawiamy do (9.56), w ten sposób otrzymujemy następną tożsamość, z której napiszemy jak zmienia się wysokość obrazów względem ich odległości od zwierciadła, czyli wzór na powiększenie obrazu względem przedmiotu, który jest stosunkiem odległości obrazu i przedmiotu:
![]() ![]() |
(9.60)
|
[edytuj] Soczewka dwuwypukła
Narysujmy soczewkę dwuwypukłą o promieniach krzywizny R1 i R2, i przedmiot będący po jego lewej stronie, a przedmiot po jego prawej strony, ale ten obiekt jest po przeciwnej stronie niż przedmiot. Na podstawie rysunku obok wyprowadzimy wzór soczewki i soczewkowy, a także wzór na powiększenie obrazu względem wielkości przedmiotu.
W poniższych obliczeniach skorzystaliśmy z małości kątów, tzn. dla których zachodzi tgα≈sinα≈α, a także z warunku ze soczewka jest cienka ze wzoru h1≈h3. Wysokość na jaką pada promień pochodzący od przedmiotu z wysokości hx. który leci na soczewkę jest sumą wysokości przedmiotu i długości jaką światło leciało pod kątem do góry.
![]() |
(9.61)
|
Wyznaczmy teraz kąt z jaką promień pochodzący z przedmiotu pada na soczewkę z jego lewej strony względem R1, który przedstawiamy jako sumę kąta promienia lecącego z przedmiotu względem poziomu i nachylenia promienia krzywizny względem poziomu:
![]() |
(9.62)
|
Kąt pod jakim promień świetlny wszedł do wnętrza w soczewce względem promienia krzywizny R1 jest określany za pomocą prawa Snelliusa dla małych kątów padania i odbicia:
![]() |
(9.63)
|
Kąt pod jakim pada promień z wnętrza na granicę soczewki względem R2, która chce wylecieć na zewnątrz jej jest określana z definicji dowolnego trójkąta, że suma kątów jest równa 180 stopni:
![]() ![]() |
(9.64)
|
Kąt pod jakim wyszedł promień świetlny z soczewki względem R2 jest opisany przez wzór wynikający z prawa załamania (9.49) dla małych katów padania i załamania:
![]() |
(9.65)
|
Kąt pod jakim promień wyszedł z soczewki względem poziomu z oczywistych powodów piszemy przez wzór:
![]() ![]() ![]() |
(9.66)
|
Kąt pomiędzy promienień wchodzący w obraz a poziomem jest wyrażony poprzez iloraz wysokości jakiej promień musi przebyć do obrazu pod kątem przez odległość obrazu od soczewki:
![]() |
(9.67)
|
Weźmy sobie dwa promienie, które wylatują z przedmiotu na wysokości hx pod kątami: α1 i α2, tworzymy w ten sposób dwa równania (9.67), które odejmiemy je od siebie i dzieląc tak powstały wzór przez α1-α2, otrzymujemy:
![]() |
(9.68)
|
Z obliczeń (9.68) otrzymujemy równanie soczewki i równanie soczewkowe, czyli które przepisujemy przy pomocy dwóch wzorów, oznaczając przez "f" długość ogniskowej soczewki:
|
|
Tożsamość (9.69) podstawiamy w miejsce (9.70) we wzorze (9.67), bo one oznaczają to samo, w ten sposób możemy otrzymać równość, z którego wynika, że powiększenie obrazu jest stosunkiem położenia obrazu "y" i położenia przedmiotu "x" względem soczewki:
![]() ![]() |
(9.70)
|
[edytuj] Fale w mechanice kwantowej i elektrodynamice klasycznej
[edytuj] Przykład cząsteczki amoniaki jako układ drgający w mechanice kwantowej
Rozważmy sobie cząsteczkę, w którym atomy wodoru H tworzą trójkąt równoboczny. Dla azotu N mamy dwa możliwe drgania, w których może on drgać na dwa sposoby odpowiadające dwom wahadłom a i b sprzężonym ze sobą. Jedno drganie jest po pierwszej stronie płaszczyzny H3, nazwijmy je przez "a", a drugie po drugiej stronie i nazwijmy je przez b. W mechanice klasycznej te położenia są położeniami równowagi trwałej, więc nie możliwe jest przejście z jednego stanu do drugiego. W mechanice kwantowej wprowadza się sprzężenie pomiędzy tymi stanami pozwalający na przenikanie płaszczyzny bariery potencjału. Załóżmy, że cząstka początkowo znajdowała się w stanie "a" w chwili t=0, wtedy zachodziło by |ψa|2=1, |ψb|2=0, znaczy to, że prawdopodobieństwo, że cząstka drga w stanie "a" jest równe pewności, ale N nie drga w stanie "b" w tymże czasie. Prawdopodobieństwa znalezienia cząstki w stanie "a" lub "b" możemy przestawić przez dwa równania poprzez czas "t":
|
|
- gdzie ω1 i ω2 są częstościami postaci drgań normalnych.
Równania (10.1) i (10.2) są kolejno bardzo podobne do (1.133) i do (1.134). Całkowite prawdopodobieństwo znalezienia cząstki w stanie "a" lub "b" jest oczywiście równe jedności. Załóżmy, że cząstka w stanie 2 drga z nieco większą częstością ω2>ω1 niż w stanie jeden, atom N znajdujący się w tym stanie, który jest stanem nietrwałym i wyniku promieniowania elektromagnetycznego przechodzi ze stanu 2 zwany stanem wzbudzonym do stany 1 zwany stanem stanem podstawowym. To promieniowanie można wykryć w wyniku jego częstości dudnień ν2-ν1, która ma wartość
, co odpowiada to długości fali równą około 1,5 cm, która jest w stanie radarowym lub mikrofalowym. Jeśli przez cząsteczki amoniaku będziemy przepuszczać o tej częstości promieniowanie, to okazuje się, że niektóre z atomów azotu N przejdzie ze stanu podstawowego do wzbudzonego.
[edytuj] Przykład układu sprzężonych ze sobą obojętnych mezonów
Bardzo interesującym przykładem jest układ sprzężonych ze sobą układu dwóch mezonów
i
, które zachowują się jak układ sprzężonych ze sobą wahadeł, które jako sprężynki spełniają rolę piony π, które każdy drogą słabego oddziaływania może oddziaływać z dwoma mezonami. Mamy tutaj dwa drgania proste zwanych mezonami
,
. W przypadku, gdy te dwie postacie drgań prostych nie rozpadają się na cząstki, to prawdopodobieństwo, że układ znajduje się w jednych tych stanów z dwóch jest określone przez wzory (10.1) i (10.2). Jeśli weźmiemy przykład tłumienia drgań dla tych dwóch cząstek, co jest równoważne rozkładowi tych mezonów na inne cząstki, to prawdopodobieństwo, że układ znajduje się w pierwszym lub drugim stanie, jest określane przez:
![]() |
(10.3)
|
![]() |
(10.4)
|
[edytuj] Wyprowadzenie związku dyspersyjnego dla fal de Broglie'a
Funkcja falowa drgań harmonicznych o częstości ω przestawiamy jako funkcję eksponencjalną, której amplituda jest zależna od położenia z:
![]() |
(10.5)
|
Jeśli ten związek podstawimy do równania mechaniki kwantowej zależnego od czasu, w ten sposób uzyskamy zależność różniczkową na funkcję f(z):
![]() |
(10.6)
|
Ze związku (10.6) otrzymujemy natychmiast funkcję f(z), która jest kombinacją liniową funkcji sinus i kosinus, którym argumentem jest iloczyn liczby falowej i położenia cząstki kwantowej, wtedy całą funkcję falową ψ(z,t) (10.5) możemy określić:
![]() |
(10.7)
|
Jeśli jeszcze raz podstawimy (10.7) do równania zależnego od czasu, wtedy otrzymujemy związek dyspersyjny częstotliwości kołowej ω w zależności od liczby falowej k:
![]() |
(10.8)
|
[edytuj] Fale stojące na kwantowej strunie skrzypcowej
Mając związek w postaci funkcji falowej (10.7), która na jej dwóch końcach ta funkcja nie przedstawia żadnych drgań, wtedy dla jednego końca z=0, otrzymujemy natychmiast B=0, mamy:
![]() |
(10.9)
|
Ponieważ długość sprężynki jest skończona i wynosi L, to aby dla prawego końca funkcja (10.9) była równa zero, to musi być na pewno spełniony związek:
![]() |
(10.10)
|
Unormujmy funkcję (10.9) do jedynki licząc całkę z kwadratu modułu (10.9) wykorzystując warunki brzegowe (10.10):
![]() ![]() |
(10.11)
|
Ostateczną funkcję falową (10.9) na podstawie obliczeń (10.11) jest funkcja zależna od "z" przy warunkach brzegowych (10.10):
![]() |
(10.12)
|
Związek dyspersyjny jaki rządzi równaniem (10.12), który wynika z (10.8), jest:
![]() |
(10.13)
|
[edytuj] Kwantowy i klasyczny ośrodek niejednorodny
Rozważmy kwantowy ośrodek, w którym panuje niejednorodny potencjał zależny od z, czyli V(z), wtedy możemy napisać w analogii do (10.6) wzór:
![]() |
(10.14)
|
W mechanice klasycznej dla niejednorodnej struny (2.57) równanie, w której w strunie wszędzie panuje takie same naprężenie równe T0(z)=T0, powstałe po podstawieniu funkcji falowej (10.5), jest:
![]() |
(10.15)
|
[edytuj] Obszar zakazany w sensie klasycznym i przenikanie do niego cząstki kwantowej i klasycznej
Energia klasyczna całkowita cząstki krążącej pomiędzy barierami potencjału, która ta cząstka przez nie nie może przeniknąć jest określona przez:
![]() |
(10.16)
|
Jeśli ona porusza się w przedziale o zerowym potencjale klasycznym i chce się przedostać do nieskończonego potencjału klasycznego, to cząstka w nim przebywająca ma energię kinetyczną ujemną, co jest pozbawione sensu fizycznego w mechanice fizycznej, ale nie w mechanice kwantowej. W fizyce klasycznej cząstka uderzająca w barierę potencjału zmienia swój pęd na przeciwny. Związek dyspersyjny rządzący pomiędzy częstością kołową drgań cząstki kwantowej, a liczbą falową uzyskujemy podstawiając do (10.16) zależności
i
:
![]() |
(10.17)
|
[edytuj] Przechodzenie przez barierę potencjału a analogia z sprzężonymi wahadłami w mechanice klasycznej
W mechanice klasycznej istnieje związek dyspersyjny dla małych liczb falowych dla nieciągłej struny (2.86):
![]() |
(10.18)
|
Inny związek dla częstości reaktywnych, który jest przestawiony związkiem (3.96) dla małych współczynników wygaszania, jest:
![]() |
(10.19)
|
W mechanice kwantowej związek dyspersyjny dla fal eksponencjalnych możemy przedstawić przy pomocy związku zapisanego przy pomocy współczynnika wygaszania χ:
![]() |
(10.20)
|
[edytuj] Prędkość fazowa i grupowa według mechaniki klasycznej i szczególnej teorii względności
Biorąc związek słuszny w mechanice kwantowej klasycznej (10.17), to policzmy jego prędkość fazową i grupową poruszającej się cząstki w mechanice kwantowej falowej:
|
|
Biorąc związek znany ze szczególnej teorii względności możemy zapisać związek pomiędzy energią a pędem cząstki, to zależność dyspersyjna pomiędzy częstotliwością a liczbą falową przestawiamy według związku:
![]() |
(10.23)
|
Prędkość grupową cząstki według mechaniki kwantowej relatywistycznej możemy przedstawić jako pochodna częstotliwości kołowej wynikającej z równania (10.23) względem liczby falowej, to wszystko określamy poprzez:
![]() |
(10.24)
|
Widzimy, że według mechaniki kwantowej klasycznej (10.22) i według mechaniki relatywistycznej (10.24) prędkość grupowa jest równa prędkości cząstki.
[edytuj] Wyprowadzenie równań falowych fal de Broglie'a według hamiltonianu klasycznego i relatywistycznego
Funkcję falową określmy dla fali o częstości kołowej ω i liczbie falowej "k" jako fale biegnące w kierunku zgodnym z osią "z" i niezgodnym z tą osią przestawiamy:
![]() |
(10.25)
|
Policzmy teraz drugą pochodną wyrażenia (10.25) względem czasu najpierw licząc jego pierwszą pochodną, dla której dowiemy się, że ona jest wprost proporcjonalna do kwadratu częstotliwości kołowej wziętej z minusem i do funkcji falowej ψ(z,t):
![]() |
(10.26)
|
A teraz policzmy drugą pochodną funkcji falowej (10.26), ale najpierw licząc jego pierwszą pochodną, dla której dowiemy się, że ona jest wprost proporcjonalna do kwadratu liczby falowej wziętej z minusem i do funkcji falowej ψ(z,t):
![]() |
(10.27)
|
Będziemy tutaj korzystać ze wzoru na całkowitą energie układu składającą się z jednej cząstki (10.16) i podstawimy do niego wzór na energię równą
i na pęd
, wtedy po podstawieniu do niego wzoru na częstotliwość kołową wynikającego z (10.26) i na kwadrat liczby falowej wynikającego z (10.27), w ten sposób otrzymujemy równanie falowe mechaniki kwantowej klasycznej wychodząc z zasady zachowania energii dla mechaniki Newtona:
![]() |
(10.28)
|
Powyższe równanie mechaniki falowej klasycznej jest spełnione dla V stałego, ale można go uogólnić dla pewnej niestałej funkcji potencjału zależną od z. To równanie tak uogólnione zwane jest jednowymiarowym równaniem Schrödingera. Mając wzór (10.23), który można go tak rozpisać w troszeczkę w innej postaci wykonując potęgowanie w jego drugim wyrazie, w ten sposób otrzymujemy równość mechaniki relatywistycznej kwantowej Kliena-Gordona:
![]() |
(10.29)
|
Równanie falowe (10.28) dla m=0 jest równaniem falowym opisujących fotony dla fal nieulegających dyspersji pędzących z prędkością równą "c".
[edytuj] Promieniowanie elektromagnetyczne wysyłane przez jednowymiarowy atom w mechanice falowej
Weźmy sobie teraz elektron krążący w studni potencjału, którego funkcja falowa jest superpozycją dwóch stanów pierwszego podstawowego i pierwszego wzbudzonego:
![]() |
(10.30)
|
Wyznaczmy teraz kwadrat modułu wyrażenia (10.30), którego definicja znana jest z algebry, przedstawiony w postaci:
![]() |
(10.31)
|
Policzmy teraz całkę kwadratu modułu funkcji falowej (10.29), czyli funkcji (10.30), wykorzystując wzór na sinus i kosinusa podwojonego kąta, co według powyższych rozważań tą naszą całkę przedstawiamy:
![]() |
(10.32)
|
Wyznaczmy teraz całkę z iloczynu funkcji "z" i kwadratu modułu funkcji falowej (10.31), którą możemy przecałkować względem "z" i wyłączając przed nawias funkcję niezależną od "z" i wykorzystując wzór na sumę sinusów znanego ze szkoły średniej:
![]() ![]() ![]() ![]() ![]() |
(10.33)
|
Całkowita średnia prędkość cząstki określamy przy pomocy wzorów (10.30) i (10.31), która jak się przekonamy zależy od czasu "t" i różnicy częstotliwości kołowych ω1 i ω2 i amplitud A1 i A2:
![]() |
(10.34)
|
Położenie elektronu krążącego, który ma ładunek q=-e oscyluje według wzoru (10.34) z częstością dudnień ω2-ω1, który to wysyła promieniowanie w postaci lecących fotonów ze wspomnianą częstością kołową.
[edytuj] Źródła spójne w czasie przy częstości dudnień optycznych
Zapiszmy teraz w postaci zespolonej całkowite zespolone natężenie fali, w której pierwsza fala jest napisana przy częstości kołowej ω1 i o przesunięciu fazowym φ1, a druga fala jest o częstości ω2 i o przesunięciu fazowym φ2, które razem to przepisujemy dokonując superpozycję tychże fal:
![]() |
(10.35)
|
Każdą liczbę zespoloną możemy rozłożyć jako wektor w płaszczyźnie zespolonej, zatem długość wektora jako liczby zespolonej (10.35) możemy przedstawić jako moduł z tej wielkości:
![]() |
(10.36)
|
Średnią energię przedstawiamy jako wielkość uśrednianą w ciągu jednego okresu wielkości (10.36) dla zespolonego natężenia fali (10.35), jak widzimy w tak otrzymanym wzorze na strumień energii występuje częstość dudnień równą ω1-ω2, która powinna przyjmować stosunkowo niską wartość. Czas koherencji, który jest odwrotnością szerokości pierwszej fali dla zmieniającej się stałej fazowej w ciągu jednego okresu, której również amplitudy tychże drgań też się zmieniają w sposób przypadkowy, przedstawiamy:
|
|
Aby można było zaobserwować dudnienia zapisane według ν1-ν2, to czasy (okresy) koherencji powinny być duże w porównaniu w porównaniu ze wspomnianą częstością dudnień, zatem aby eksperyment obserwacji tychże fal udał się, to czasy dudnień powinny spełniać następujące dwa warunki:
|
|
[edytuj] Czemu niebo jest jasne
W rozdziale "Czemu niebo jest niebieskawe" rozpatrywaliśmy czemu niebo w dzień jest niebieskie, a wieczorem czerwone. A teraz rozpatrzmy czemu niebo jest jasne. Rozpatrzymy teraz rozumowanie, które mówi dlaczego niebo powinno być niewidzialne. Rozpatrzmy elektron znajdujący w cząsteczce, który wysyła promieniowanie we wszystkich kierunkach znajdujący się pierwszym atomie, a teraz weźmy sobie drugi atom, który jest w pół długości fali dalej od obserwatora, oba te cząsteczki są pobudzane tymi samymi przesunięciami fazowymi, i tymi samymi amplitudami, co w wyniku nakładania się tych dwóch fali daje nam zero w położeniu obserwatora. Dla rozpraszania pod kątem 900 możemy spełnić te żądania dla fazy i amplitudy tak by cząsteczce pierwszej odpowiadała cząsteczka druga. Dla rozpraszania bliskiego mimo, że cząsteczki znajdujące od obserwatora są pobudzane o pół okresu wcześniej i są oddalone o półokresu dalej od obserwatora, ale one nie powodują one interferencji dekonstruktywnej. Biorąc pod uwagę, że sześcian zawiera
cząsteczek, gęstość powietrza maleje wykładniczo wraz z wysokością, co powinno wystarczyć do zaistnienia interferencji dekonstruktywnej, ale to nie wystarcza, dlaczego? Powyższe przewidywania przeczą doświadczeniu, bo natężenie jest niemal takie jakie są wyniki rozpraszania na N cząsteczkach i napisane jest w postaci sumy natężeniem rozpraszania przez te cząsteczki, które są przyczynkami do tego natężenia. I nie wiadomo nadal dlaczego interferencja dekonstruktywna nie następuje. Jeśli obierzemy wodę zamiast powietrza interferencja dekonstruktywna następuje dla rozpraszania 900. Wiązka z latarki przechodzi przez czystą wodę przy nieznacznej stracie natężenia, jeśli pominąć nieznaczne rozszerzenia wiązki spowodowanej przez dyfrakcję. Rozproszenie przy 900 przy przejściu 10 m jest małe w porównaniu z powietrzem przy takiej samej wysokości. Przy tym rozproszeniu w wodzie dodaje się do siebie amplitudy fali rozproszeniowej, a w przypadku powietrza ich natężenia, dlatego w wodzie występuje interferencja dekonstruktywna. Tajemnicza kryje się w tym, że cząsteczki wody są rozłożone równomierne i znajdzie się zawsze cząsteczka wody numer dwa, wyniku której przewidziana interferencja jest dekonstruktywna, a w powietrzu na pierwszą cząsteczkę powietrza przypada średnio jedna ale druga cząsteczka powietrza, czasem w cale, zatem to wyjaśnia dlaczego w powietrzu sumują się natężenia, a wodzie amplitudy.
Weźmy sobie teraz obszar pierwszy i drugi znajdujący w tej samej odległości od słońca w przybliżeniu, ale obszar drugi znajduje się pół okresu dalej, a każdy z tych obszarów jest mały w porównaniu z długością fali światła monochromatycznego tutaj rozpatrywanego. Weźmy wszystkie cząsteczki znajdujące się w obszarze pierwszym, w którym każda cząsteczka wnosi przyczynek
, a w obszarze dwa przyczynek
, czyli całkowite natężenie pola w pierwszym obszarze jest
, a w obszarze drugim jest
, jeśli przyjmować będziemy falę biegnącą, który drugi obszar jest pół okresu dalej niż pierwszy, czyli powinno wtedy zachodzić
, zatem całkowite natężenie fali docierające do obserwatora jest:
![]() |
(10.41)
|
Pole wytwarzane przez jedną pobudzoną cząsteczkę jest opisywane przez funkcję:
![]() |
(10.42)
|
A całkowite natężenie fali wytarzanej przez cząsteczki w obszarze pierwszym i drugim jest opisywane poprzez wzór:
![]() |
(10.43)
|
Biorąc natężenie fali wytwarzanej prze pojedynczą cząsteczkę (10.41) i przez cząsteczki znajdujące się obszarze drugim (10.43), wstawiając to wszystko do wzoru na całkowite natężenie fali w punkcie obserwatora, po tak otrzymanej tożsamości i podzieleniu przez kosinus, otrzymujemy zależność pomiędzy amplitudami:
![]() |
(10.44)
|
We wzorze n1 jest raz mniejsze, a za drugim większe od n2, a średnio rzecz mówiąc oba te wielkości średnio są sobie równe, wtedy średnio powinna występować interferencja dekonstruktywna pod kątem 900. Całkowite natężenie fali obserwowanego przez naszego obserwatora jest wyrażona poprzez wzór:
![]() |
(10.45)
|
Weźmy sobie, że liczba cząsteczek znajdujących się w pierwszym obszarze jest
, a w drugim obszarze jest
, wtedy pierwszy czynnik występujący po prawej stronie pod średnią możemy przepisać wiedząc, że średnia wartość wartości n1 jest
, a także średnia liczba cząsteczek n2 jest równa
, wtedy można napisać:
![]() |
(10.46)
|
Policzmy teraz średnią wartość trzeciego wyrażenia występującego w punkcie (10.46), która jak się przekonamy jest równa zero, a oto jego dowód:
![]() |
(10.47)
|
Biorąc rozkład Poissona, w którym wartość oczekiwana jest równa wariancji, zatem na podstawie tego możemy napisać dwa wnioski wynikające z tego rozkładu:
|
|
Biorąc wnioski (10.48) i (10.49), a także wniosek (10.46) i (10.47), wtedy (10.45) możemy zapisać do postaci:
![]() |
(10.50)
|
Dla wody warunki (10.48) i (10.49), które są dla powietrza, odpowiadają warunkom:
|
|
Patrząc na warunek (10.50), który odpowiada ośrodkowi, który jest powietrzem opisanym na podstawie wzorów (10.51) i (10.52), natomiast dla wody odpowiada warunek:
![]() |
(10.53)
|
[edytuj] Fale elektromagnetyczne rozchodzące się w ośrodkach materialnych
[edytuj] Równania Maxwella
|
|
|||||||||||
|
|
|
|
|
||||||||
|
|
|
|
||||||||
|
|
|||||||||||
|
|
|
||||||||||
|
|
||||||||||
[edytuj] Prowadzenie do teorii o ośrodkach liniowo izotropowych
Siła działający na ładunek q określamy poprzez wzór napisany przy pomocy natężenia pola elektrycznego i indukcji pola magnetycznego w postaci:
![]() |
(10.60)
|
Siły działająca na ładunki znajdujące się w ośrodku pochodzą od natężenia pola elektrycznego
i indukcji pola magnetycznego
. Ośrodek nazywamy izotropowym, jeśli polaryzacja
ma kierunek wzdłuż wektora
, a magnetyzacja ma kierunek zgodny z wektorem indukcji magnetycznej
. Z powyższych warunków możemy wywnioskować wnioski, że jeśli polaryzacja jest równa zero, to natężenie pola też jest równe zero, a także jeśli magnetyzacja jest równa zero, to wektor indukcji jest równy zero. Ośrodek izotropowy jest taki, że jeśli iksowa składowa polaryzacji zależy tylko od iksowej składowej iksowej pola elektrycznego, a nie od igrekowego, ten związek tej zależności jest:
![]() |
(10.61)
|
Przy dostatecznie małych natężeniach pola elektrycznego w wyrażeniu (10.61) możemy zaniedbać wyrazy wyższe niż liniowe. Polaryzację i magnetyzację liniową możemy przedstawić jako:
|
|
Przenikalności elektryczne i magnetyczne zdefiniowane poprzez wzory (10.58) i (10.59) są napisane:
|
|
Względne przenikalności elektryczne i magnetyczne zdefiniowane są na podstawie wzorów (10.64) i (10.65):
|
|
Uogólnienie wniosku (10.62) dla pól zależnych od czasu, w którym polaryzacja jest równoległa do natężenia pola elektrycznego dla pól rozprzestrzeniających się w czasie prowadzi do fałszywych wniosków, że powinno zachodzić dla częstości drgań fali elektrycznej ω:
![]() |
(10.68)
|
W tym celu należy wprowadzić podatność elektryczną elastyczną i absorpcyjną, dzięki któremu możemy wprowadzić polaryzację iksową, którego definicja jest sumą iloczynu przenikalności elektrycznej w próżni przez podatność elastyczną i przez Ex i iloczynu przenikalności elektrycznej w próżni przez podatność elektryczną absorpcyjną i przez natężenia fali pola elektrycznego przesuniętego w fazie o 900:
![]() |
(10.69)
|
Natężenie fali elektromagnetycznej zmienia się według funkcji proporcjonalnej do kosinusa, tzn.:
![]() |
(10.70)
|
wtedy polaryzację wyrażoną wzorem (10.69) przy definicji natężenia iksowego fali pola elektrycznego (10.70) możemy przepisać:
![]() |
(10.71)
|
[edytuj] Prosty model ośrodka izotropowego liniowego i wyznaczanie podatności elastycznej i magnetycznej
Rozważmy prosty model ośrodka liniowego i izotropowego, w której jądra są znacznie cięższe od rozważanego ładunku, który jest połączony z jądrem sprężyście o stałej sprężystości równą
. Ruch jądra jest pomijalny, wiec nie mamy magnetyzacji. Na ten ładunek działa siła tłumiąca jego ruch o stałej tłumienia równą
, zatem równanie ruchu ładunku q, na którą działa pole elektryczne o przebiegu kosinusoidalnej, jest:
![]() |
(10.72)
|
Natężenie fali pola elektrycznego, w zależności od częstotliwości występujące w równaniu (10.72), jest:
![]() |
(10.73)
|
Równanie ruchu położenia cząstki, który określamy poprzez amplitudę elastyczną Eel i amplitudę absorpcyjną Aab, piszemy:
![]() |
(10.74)
|
Amplitudę absorpcyjną Aab i elastyczną Ael, które już określaliśmy dla F0=qE0 w punktach (3.29) i (3.30), które przepiszemy tutaj dla przejrzystości wykładu:
|
|
Ponieważ polaryzację definiujemy jako iloczyn koncentracji N, ładunku q i położenia tegoż ładunku (10.74) przy definicji natężenia fali pola elektrycznego (10.73):
![]() ![]() |
(10.77)
|
Możemy porównać wzory (10.71) z (10.77) na polaryzację iksową, wtedy dostajemy wzory na podatność elektryczną elastyczną i absorpcyjną:
|
|
[edytuj] Liczby zespolone w równaniach elektrodynamiki klasycznej i prowadzenie definicji podatności zespolonej
Równania elektrodynamiki klasycznej Maxwella nie zawierają w sobie liczb zespolonej, ale też jego rozwiązania, ale możemy prowadzić podatność elektryczną jako wielkość zespoloną jako:
![]() |
(10.80)
|
Polaryzację zespoloną możemy przedstawić w zależności od natężenia pola elektrycznego zespolonego dla fali elektrycznej o częstotliwości kołowej ω wedle:
![]() |
(10.81)
|
Jeśli przedstawimy natężenie zespolone pola elektromagnetycznego w postaci wzoru zależnego od częstotliwości kołowej zapisywaną w konwencji używanej w mechanice kwantowej:
![]() |
(10.82)
|
wtedy zapis polaryzacji zespolonej (10.81) przy zastosowaniu definicji podatności elektrycznej (10.80) i wzoru na natężenie fali pola elektrycznego (10.82) prowadzi do:
![]() ![]() |
(10.83)
|
Część rzeczywista równania na polaryzację zespoloną (10.83) jest w naszym przypadku polaryzacją elektryczną ośrodka, która jest taka sama jak wzór (10.71), wtedy definicja podatności zespolonej elektrycznej (10.80) jest definicją poprawną i zgadza się z naszymi rozważaniami.
[edytuj] Prosty model ośrodka liniowego izotropowego i jego zespolona stała dielektryczna
Względną zespoloną stałą dielektryczną możemy liczyć z bardzo podobnej definicji do (10.66), przy definicji podatności elektrycznej równą (10.80):
![]() |
(10.84)
|
[edytuj] Rozwiązywanie ruchu ładunku q w ośrodku tłumionym o wymuszonym pod działaniem konsinosidalnej siły elektrycznej przy wprowadzeniu liczb zespolonych
Równanie ruchu danego ładunku w ośrodku przestawiamy analogicznie do (10.73), który jest przedstawiony w przestrzeni rzeczywistej, a tutaj przedstawimy je w przestrzeni zespolonej:
![]() |
(10.85)
|
Będziemy tutaj rozpatrywać ruch ładunków ośrodka, którego rozwiązanie przepuszczalne jest dla x=x0e-iωt, wtedy podstawiając to do (10.85) i pamiętając przy tym
, i
, wtedy powyższe równanie ruchu cząstki obdarzonej ładunkiem q w ośrodku jest:
![]() |
(10.86)
|
Równanie (10.86), a właściwie jej część rzeczywista jest równoważna z rozwiązaniem (10.74) przy definicji amplitudy elastycznej Ael (10.75) i absorpcyjnej Aab (10.86). Stałą podatność piszemy w postaci wzoru wynikającą z (10.81), w której wykorzystamy definicję polaryzacji elektrycznej (10.77) do zapisania tej wielkości:
![]() |
(10.87)
|
Zespoloną względna przenikalność elektryczną liczmy ze wzoru (10.87), w ten sposób otrzymujemy analogiczny do (10.84) wzór:
![]() |
(10.88)
|
Jak można wykazać, wzór (10.88) jest równoważny ze wzorem (10.84), co tutaj nie będziemy wykazywali, co jest trywialne.
[edytuj] Równania elektrodynamiki klasycznej dla ośrodka liniowego izotropowego
|
|
|||||||||||
|
|
|
|
|
||||||||
|
|
|
|
||||||||
Będziemy tutaj rozpatrywali, gdy gęstość objętościowa i prąd objętościowy ładunku swobodnych są te wielkości równe zero, jeśli będziemy rozpatrywali fale elektromagnetyczne.
[edytuj] Równanie falowe
Mnożymy obustronnie lewostronnie przez
wzór (10.92), dalej po skorzystaniu ze wzoru (10.91), w ostateczności otrzymujemy:
![]() |
(10.93)
|
Podobne mnożymy obustronnie lewostronnie przez
wzór (10.91), dalej po skorzystaniu ze wzoru (10.92), otrzymujemy:
![]() |
(10.94)
|
Będziemy tutaj wykorzystywali tożsamość, którą przepiszemy jako bez dowodu, którego przeprowadzenie jest bardzo łatwe.
![]() |
(10.95)
|
Jeśli skorzystamy z tożsamości (10.95) do (10.93) i (10.94), to w ten sposób otrzymujemy dwa równania falowe na natężenia pola elektrycznego i indukcję pola magnetycznego:
|
|
Równania (10.96) i (10.97) są to równania fali, których jest trzy dla fali pola elektrycznego i trzy dla fali pola magnetycznego, których razem mamy sześć równań dla fali pola elektromagnetycznego, które w ogólnym przypadku zapisujemy:
![]() |
(10.98)
|
dla którego kwadrat liczby falowej przepisujemy do postaci zależnej od względnej przenikalności dielektrycznej i względnej przenikalności magnetycznej:
![]() |
(10.99)
|
Dla przenikalności elektrycznej względnej i przenikalność magnetycznej, które te dwie wielkości są zespolone, wtedy (10.99) możemy zapisać w postaci kwadratu zespolonej liczby falowej dla rzeczywistej częstotliwości kołowej jako:
![]() |
(10.100)
|
Każdą falę możemy przedstawić w postaci zespolonej według wzoru:
![]() |
(10.101)
|
Dla fal elektromagnetycznych zachodzi bardzo podobny związek do (10.1101), który tutaj dla naszego przypadku nie będziemy przepisywać dla zwięzłości wykładu.
[edytuj] Związek pomiędzy natężeniem pola elektrycznego i magnetycznego dla fali płaskiej elektromagnetycznej
Związek pomiędzy falami magnetycznymi i elektrycznymi na podstawie (10.83) i (10.84) dla iksowego pola elektrycznego i igrekowego pola magnetycznego przepisujemy w postaci dwóch równań:
|
|
Napiszmy teraz natężenie iksowe pola elektrycznego i magnetycznego w postaci równań fali biegnącej w kierunku osi z, jest ona superpozycją fali biegnącej w kierunku +z i w kierunku -z:
|
|
Równania (10.104) i (10.105) podstawiamy do (10.102), w ten sposób mamy tożsamość, którego postać po obu stronach jest jednakowa:
![]() |
(10.106)
|
Powyższą równość jest tożsamością, jak można było udowodnić przy definicji współczynnika załamania n=c/v, następnie przejdźmy do dowodu równości (10.103):
![]() |
(10.107)
|
Powyższa równość staje się tożsamością przy definicji współczynnika załamania n=c/v.
[edytuj] Współczynnik transmisji i odbicia
Natężenie pola elektrycznego padająca na granicę dwóch ośrodków, na której tej granicy fala padająca jest odbijana przy współczynniku odbicia R12, a także część tej fali przechodzi do drugiego ośrodka, którego to równania fali dla natężenia fali elektrycznej padającej i odbitej przedstawiamy:
|
|
Dla pola magnetycznego zachodzą bardzo podobne związki do (10.108) i (10.109) przy definicji współczynnika odbicia i transmisji:
|
|
Ciągłość w punkcie z=0, daje nam, że suma współczynnika odbicia R12 i jedynki jest równa współczynnikowi transmisji dla równań (10.108) i (10.109). A także w wyniku tożsamości
, której natężenie pola magnetycznego jest wielkością ciągłą dla równań (10.101) i (10.111) dla tego samego "z", możemy zapisać:
|
|
Zdefiniujemy teraz impedancję charakterystyczną, która jest w zależności od przenikalności elektrycznej i magnetycznej ośrodka, którego początkowa definicja jest ilorazem przenikalności magnetycznej ośrodka przez współczynnik załamania:
![]() |
(10.114)
|
Jeśli wykorzystamy definicje impedancji (10.114), a także wykorzystamy do tego celu tożsamości (10.112) i (10.113), w ten sposób możemy wyznaczyć wzór na współczynnik odbicia, dalej przepiszemy jego postać dla przenikalności magnetycznej równą jeden:
![]() |
(10.115)
|
Jeśli współczynnik załamania ośrodka pierwszego jest równy jedności, tzn. n1=1, wtedy współczynnik odbicia (10.115), przy definicji impedancji (10.114), jest równy:
![]() |
(10.116)
|
Jeśli w współczynniku odbicia R12podstawimy współczynnik załamania przez n=nR+inl dla współczynnika załamania n1 równy jeden i dla współczynnika załamania n2 równą po prostu n. Wtedy można przedstawić R12 i jego kwadrat modułu w jednej linijce:
|
|
[edytuj] Uproszony prosty model przewodnika, przenikalność i przewodność elektryczna
Weźmy sobie uproszczony model przewodnika, w którym każdy elektron już nie oddziaływuje z jądrem na w sposób sprężysty, na te ładunki w przewodniku działa siła pochodząca od pola elektrycznego, jego równanie ruchu jest:
![]() |
(10.119)
|
Weźmy teraz stacjonarne pola elektryczne, to (10.119) możemy rozwiązać dokonując w nim podstawienia
:
![]() ![]() |
(10.120)
|
Jeśli przyjmować będziemy przykład dużego tłumienia w przewodniku, wtedy wyrażenie (10.120) możemy przepisać w sposób przybliżony:
![]() |
(10.121)
|
Wyobraźmy sobie elektron, który oddziałuje z polem o charakterze fali, dla Γ>>ω, wtedy na podstawie tego współczynnik elastyczny Ael jest równy zero, wtedy współczynnik absorpcyjny przybiera postać:
![]() |
(10.122)
|
Pochodna cząstkowa położenia ładunku q określimy jako pochodna zupełna położenia elektronu (10.73):
![]() |
(10.123)
|
Gęstość prądu elektrycznego, z którego wyprowadzimy różniczkowe prawo Ohma, definiujemy:
![]() |
(10.124)
|
Do tego problemu podejdźmy teraz z innej strony rozpatrując, że pole elektryczne drga i zapiszemy go w postaci eksponencjalnej, wtedy gęstość prądu elektrycznego możemy przepisać w formie:
![]() |
(10.125)
|
Z posługiwania się tożsamością (10.125) możemy wywnioskować, że przewodność elektryczną definiujemy:
![]() |
(10.126)
|
Zamiast posługiwania się powyższym wzorem na przewodność elektryczną, tzn. jej częścią rzeczywistą i urojoną, można się posługiwać zespolonymi wielkościami, dla których zespoloną podatność elektryczną i zespoloną przewodność elektryczną dla zerowania się częstotliwości drgań podstawowych piszemy przy wykorzystaniu (10.87):
|
|
Gdy wybierzemy przykład ω<<Γ, wtedy zespoloną podatność elektryczną (10.127) i zespoloną przewodność elektryczną (10.128) piszemy:
|
|
Wyznaczmy teraz współczynnik załamania światła, który możemy przedstawić z zależności od podatności elektrycznej, gdy współczynnik przenikalności magnetycznej względny wynosi jeden:
![]() |
(10.131)
|
[edytuj] Rozrzedzony ośrodek oporowy
Dla rozrzedzonego ośrodka oporowego, zachodzą warunki, z których będziemy skorzystać:
![]() |
(10.132)
|
wtedy współczynnik załamania światła (10.131) jest napisany:
![]() |
(10.133)
|
Przy liczeniu współczynnika załamania pomijaliśmy wyrazy wyższego rzędu. Policzmy zespoloną liczbę falową wykorzystując (10.133) przy użyciu definicji kwadratu ωpl i definicji przewodności elektrycznej σ (10.124):
![]() |
(10.134)
|
Część urojona liczby falowej (10.134) powoduje, że fala zanika wraz zwiększającą odległością zetową "z", ta część funkcji falowej dana jest przez e-Im(k)z, natężanie opisywanej fali jest wprost proporcjonalne do kwadratu tej części wspomnianej, czyli do e-2Im(k)z. Odległość zaniku natężenia fali elektromagnetycznej, określamy przez d=(2Im(k))-1, piszemy przy pomocy współczynnika przewodności:
![]() |
(10.135)
|
Mając współczynnik załamania zapisanego wzorem (10.133), wtedy kwadrat modułu współczynnika odbicia (10.118) określamy poprzez:
![]() |
(10.136)
|
Do końcowego wzoru na współczynnik odbicia (10.136) podstawiamy cześć zespoloną wyrażenia (10.133), otrzymujemy:
![]() |
(10.137)
|
[edytuj] Gęsty ośrodek oporowy
W przypadku gęstego ośrodka oporowego zachodzą warunki, z których będziemy korzystali:
![]() |
(10.138)
|
Biorąc wzór na współczynnik załamania n (10.133) skorzystamy tutaj z tożsamości, którą można udowodnić korzystając z wiadomości z algebry, tj.
, co odpowiada
, by potem było można oczywiście powiedzieć:
![]() |
(10.139)
|
Zespoloną liczbę falową, na podstawie obliczeń (10.139) i przy skorzystaniu (10.131), możemy przedstawić w postaci zależnej od współczynnika przewodności elektrycznej:
![]() |
(10.140)
|
Mając zdefiniowany zespolony współczynnik załamania (10.139) możemy policzyć kwadrat modułu współczynnika odbicia (10.118):
![]() ![]() |
(10.141)
|
Ponieważ energia fali padającej do ośrodka jest wprost proporcjonalna do kwadratu funkcji falowej, zatem współczynnik "d" przedstawiający zanik natężenia fali wkraczający w ośrodek jest opisywany poprzez wzór:
![]() |
(10.142)
|
Jak widzimy z powyższych rozważań środek oporowy rozrzedzony jest w zasadzie czarnym, co daje niemal zupełną absorpcję docierającego do niego światła. A przeciwnie ośrodek oporowy gesty zachowuje się jak dyskretny skupiony opór falowy, który daje niemal stukrotne odbicie.
[edytuj] Całkowicie elastyczny obszar częstości
Weźmy sobie równanie (10.114) przy definicji położenia danego przez funkcję eksponencjalną e-iωt:
![]() |
(10.143)
|
Rozwiązanie równania (10.136) w postaci funkcji prędkości dla ω>>Γ jest funkcją zależną od Ex, częstotliwości kołowej ω, masy ładunku i ładunku o wartości "q":
![]() |
(10.144)
|
Jak widzimy cząstki, dla której prędkość względem siły jest przesunięta o ±900, wtedy praca wykonywana nad ładunkiem jest równa zero. Gęstość prądu elektrycznego określamy przy pomocy definicji prędkości zdefiniowaną w punkcie (10.144):
![]() |
(10.145)
|
Przestawimy teraz w jednej linijce jaka jest definicja przewodności elektrycznej według (10.145) (która jest liczbą zespoloną) i kwadrat współczynnika załamania, patrząc na równanie (10.126):
|
|
[edytuj] Przedział dla częstości dyspersyjnych
Obszar częstości dyspersyjnej możemy zdefiniować warunkami:
![]() |
(10.148)
|
Na podstawie warunku dyspersyjnego dla kwadratu współczynnika załamania (10.147) można powiedzieć, że ona mieści się w przedziale od zera do jedynki, stąd wynika, że sam współczynnik załamania też leży w tym samym przedziale:
![]() |
(10.149)
|
Widzimy na podstawie obliczeń (10.149), że współczynnik n jest rzeczywisty. Ten ośrodek spełniającej (10.149) jest ośrodkiem przezroczystym, który nie ma wcale absorpcji, a względne natężenie światła jest opisywane poprzez współczynnik odbicia (10.118).
[edytuj] Przedział dla częstości reaktywnych
Obszar częstości reaktywnej możemy zdefiniować warunkami:
![]() |
(10.150)
|
wtedy na podstawie (10.147) i (1.150) zachodzi warunek:
![]() |
(10.151)
|
Współczynnik załamania dla częstości kołowych reaktywnych na podstawie (10.151) możemy przepisać, ale nie tym razem dla kwadratu, ale dla samego współczynnika załamania, w zależności od częstotliwości kołowej fali:
![]() |
(10.152)
|
Policzmy czemu jest równa zespolona liczba falowa, której przedstawienie w zależności od współczynnika załamania (10.152) podamy w postaci urojonej w zależności od częstotliwości kołowej, modułu współczynnika załamania i prędkości światła:
![]() |
(10.153)
|
Patrząc na (10.153), to całkowite natężenie fali w ośrodku reaktywnym jest przestawiane przy pomocy funkcji eksponencjalnych oraz względem sinusa i kosinus (bo eksponens z liczby urojonej możemy rozłożyć na jej część rzeczywistą przy pomocy kosinusa i na część urojoną przy pomocy sinusa) przy pomocy czasu:
![]() |
(10.154)
|
Policzmy teraz kwadrat modułu współczynnik odbicia przy pomocy współczynnika załamania (10.152), która jest ona wyrażona przy pomocy liczby stojącej przy jednostce urojonej i bez części rzeczywistej:
![]() |
(10.155)
|























![\simeq-{{2Kx}\over{M}}\left[1-{{a_0}\over{a}}\left(1-{{1}\over{2}}{{x^2}\over{a^2}}+...\right)\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/8/e/4/8e4b5583c8921b835e979580c890c81f.png)





















































































































![E_a={{1}\over{2}}E\left[1+\cos(\omega_1-\omega_2)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/8/e/b/8eb23d0ef4803e5d752b05aee4d8f59a.png)
![E_b={{1}\over{2}}E\left[1-\cos(\omega_1-\omega_2)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/1/b/d/1bdfae70abe3d2be84545b357ff4ea74.png)







, który jest równy tangensowi nachylenia struny względem osi z dla chwili t dla wycinka struny o położeniu z. Siłę 













![\psi(x,t)=\cos\left(\omega t+\varphi\right)\left[A\sin\left(2\pi {{z}\over{\lambda}}\right)+B\cos\left(2\pi{{z}\over{\lambda}}\right)\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/c/d/a/cda812d429ce35abed3e40ead1903250.png)
![\psi(0,t)=\cos(\omega t+\varphi)[0+B]=0\;](http://upload.wikimedia.org/wikibooks/pl/math/0/7/b/07bd840a2d994bf99c2d3fc85554950a.png)
















![F(z)=\sum_{n=0}^{\infty}\left[A_n\sin n{{2\pi}\over{\lambda_1}}z+B_n\cos n{{2\pi}\over{\lambda_1}}z\right]=B_0+\sum_{n=1}^{\infty}A_n\sin n{{2\pi}\over{\lambda_1}}z+\sum_{n=0}^{\infty}B_n\cos n{{2\pi}\over{\lambda_1}}z=\;](http://upload.wikimedia.org/wikibooks/pl/math/c/8/9/c8921aeb526a014d547b0fa5608d8f9b.png)














![={{2}\over{\lambda_1}}(2)\int_0^L\left[\sin \left({{1}\over{2}}mk_1(2z+L)\right)\sin\left({{1}\over{2}}mk_1L\right)\right]=](http://upload.wikimedia.org/wikibooks/pl/math/2/5/7/2577ca9e853429744d4af6adbf7a8e21.png)

![=-{{4}\over{\lambda_1}}\sin\left({{1}\over{2}}m{{2\pi}\over{2L}}L\right){{1}\over{mk_1}}\left[\cos \left({{1}\over{2}}mk_1(2z+L)\right)\right]_0^L=\;](http://upload.wikimedia.org/wikibooks/pl/math/c/5/a/c5a80cfea13880cd4655be825a261857.png)
![=-{{2}\over{\pi m}}\sin \left(m{{\pi}\over{2}}\right)\left[\cos(m{{3}\over{2}}\pi)-\cos(m{{1}\over{2}}\pi)\right]={{4}\over{\pi m}}\sin \left(m{{\pi}\over{2}}\right)\sin \left(m\pi\right)\sin\left(m{{\pi}\over{2}}\right)=0\;](http://upload.wikimedia.org/wikibooks/pl/math/b/1/1/b11af8b4d2a9fd2d7751f9ac5521913b.png)
![A_m={{2}\over{\lambda}}\int_{z_1}^{z_1+\lambda_1}F(z)\sin mk_1z dz={{2}\over{\lambda_1}}\left[\int_{0}^L\left(\sin mk_1z-\sin mk_1(z+L)\right)dz\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/a/2/8/a28989f182972508cc21e8053e71df31.png)
![=
-{{4}\over{\lambda_1}}\left[\int_0^L\cos \left({{1}\over{2}}mk_1(2z+L)\right)\sin\left({{1}\over{2}}mk_1L\right)\right]=\;,](http://upload.wikimedia.org/wikibooks/pl/math/1/e/1/1e1abfc76cc8f970f6e4d3ed558822a9.png)
![=-{{4}\over{\lambda_1}}\sin ({{1}\over{2}}mk_1L){{1}\over{mk_1}}\left[\sin {{1}\over{2}}mk_1(2z+L)\right]_0^L=](http://upload.wikimedia.org/wikibooks/pl/math/0/0/8/0081d6e4d8be89f370d460b3eee4509b.png)
![=-{{2}\over{m\pi}}\sin m{{\pi}\over{2}}\left(\sin 3m{{\pi}\over{2}}-\sin m{{\pi}\over{2}}\right)=-{{4}\over{m\pi}}\sin\left( m{{\pi}\over{2}}\right)\left[\cos\left(m\pi\right)\sin \left(m{{\pi}\over{2}}\right)\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/c/2/d/c2d2aaffd394f59e0c2e03f31fb22370.png)


![={{4}\over{\pi}}\left[\sin k_1z+{{1}\over{3}}\sin 3k_1z+{{1}\over{5}}\sin 5k_1z+...\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/a/e/5/ae553364d0888c57e35a6b3c75a2a00d.png)








![F_x=\Delta z{{\partial}\over{\partial z}}\left[T_0(z){{\partial\psi(z,t)}\over{\partial z}}\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/6/9/7/697731afa7f27d32ecbc4dc0251bbd04.png)
![\rho(z)\Delta z{{\partial^2\psi(z,t)}\over{\partial t}}=\Delta z{{\partial}\over{\partial z}}\left[T_0(z){{\partial\psi(z,t)}\over{\partial z}}\right]\Rightarrow {{\partial^2\psi(z,t)}\over{\partial t^2}}={{1}\over{\rho(z)}}\left[T_0(z){{\partial\psi(z,t)}\over{\partial z}}\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/b/e/d/bed5c238bd78d9996919495947dd9ab4.png)


![{{1}\over{\rho(z)}}{{d}\over{dz}}\left[T_0(z){{dA(z)}\over{dz}}\right]=-\omega^2A(z)\;](http://upload.wikimedia.org/wikibooks/pl/math/b/9/4/b946ffabc8367b22bc5c9e3a4b6e176d.png)















![\omega(k)=2\sqrt{{{T_0}\over{Ma}}}\left[{{1}\over{2}}ka-{{1}\over{48}}(ka)^3+....\right]=\sqrt{{{T_0a}\over{M}}}k\left[1-{{1}\over{48}}(ka)^2+....\right]\Rightarrow \omega(k)\simeq\sqrt{{{T_0a}\over{M}}}k\;](http://upload.wikimedia.org/wikibooks/pl/math/0/e/3/0e348b48fb3d19b4e5f41923c4fc94c5.png)



























![\cos(\omega_1 t+\theta)\left[{{1}\over{4\tau^2}}-\omega_1^2-{{\Gamma}\over{2\tau}}+\omega_0^2\right]+\sin(\omega_1 t+\theta)\left[
{{\omega_1}\over{\tau}}-\Gamma\omega_1\right]=0](http://upload.wikimedia.org/wikibooks/pl/math/3/e/4/3e4f2dfeafe7292e85566aa65bc8bc0d.png)







![x(t)=e^{-{{1}\over{2}}t\Gamma}\left[x(0)\cos\omega_1 t+\left(\dot{x}(0)+{{1}\over{2}}\Gamma x(0)\right){{\sin\omega_1 t}\over{\omega_1}}\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/3/e/e/3ee37ca9cfe3b8f6d2110c80bbf3e49f.png)


![x(t)=e^{-{{1}\over{2}}t\Gamma}\left[x(0)\cosh|\omega_1| t+\left(\dot{x}(0)+{{1}\over{2}}\Gamma x(0)\right){{\sinh|\omega_1| t}\over{|\omega_1|}}\right]](http://upload.wikimedia.org/wikibooks/pl/math/e/b/f/ebf49591ab569ee361c1cd99d416191a.png)

![={{1}\over{2}}Me^{-\Gamma t}\left[\omega_1^2\langle(A\cos\omega_1 t-B\sin\omega_1 t)^2\rangle+\omega_0^2\langle(A\sin\omega_1 t+\cos\omega_1 t)^2\rangle\right]=](http://upload.wikimedia.org/wikibooks/pl/math/f/5/6/f56a09bb0e50afce9aeb4aa784b68f2a.png)

![F(t)=\sum_{\omega}F(\omega)\cos[\omega t+\varphi(\omega)]\;](http://upload.wikimedia.org/wikibooks/pl/math/c/3/1/c3110f9dc601ce760ac35964639bb9bf.png)












![=
F_0\omega\left[A_{abs}\langle\cos^2\omega t\rangle+A_{el}\langle\cos\omega t\sin\omega t\rangle\right]={{1}\over{2}}F_0\omega A_{ab}\;](http://upload.wikimedia.org/wikibooks/pl/math/c/a/0/ca0949b0b34e64480534f51feaa5a2b8.png)
![P_t=M\Gamma\langle \dot{x}_s^2\rangle={{1}\over{2}}M\Gamma\omega^2\left[A_{ab}^2+A_{el}^2\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/b/5/e/b5e2261046ce381e0e034969f8c35749.png)
![E={{1}\over{2}}M\langle\dot{x}^2\rangle+{{1}\over{2}}M\omega_0^2\langle x_s^2\rangle+{{1}\over{2}}M\omega_0^2={{1}\over{2}}M\left(\omega^2+\omega_0^2\right)\left[{{1}\over{2}}A_{ab}^2+
{{1}\over{2}}A_{el}^2\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/9/7/6/9765aab3dfc75092449c651dc3f92c93.png)








![|A(\omega)|^2={{\Gamma^2\omega^2}\over{[(\omega_0^2-\omega^2)^2+\Gamma^2\omega^2]^2}}+
{{(\omega_0^2-\omega^2)^2}\over{[(\omega_0^2-\omega^2)^2+\Gamma^2\omega^2]^2}}=\;](http://upload.wikimedia.org/wikibooks/pl/math/7/1/1/71195b24af103a459b6f1352c1ca3f26.png)





![D\simeq4\omega_0^2\left[(\omega_0-\omega)^2+{{\Gamma^2}\over{4}}\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/3/b/2/3b25d0b2844ab8877e514ce3fbaf43c4.png)




![x(t)=\left[1-e^{-{{1}\over{2}}\Gamma t}\right]\left(A_{ab}\sin\omega t+A_{el}\cos\omega t\right)=\left[1-e^{-{{1}\over{2}}\Gamma t}\right]x_s(t)\;](http://upload.wikimedia.org/wikibooks/pl/math/7/d/1/7d1fb150486effff23fde7ef1504638c.png)




![E(t)=E_0\sin^2{{1}\over{2}}(\omega_0-\omega)t={{1}\over{2}}E_0\left[1-\cos(\omega_0-\omega)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/2/5/e/25ed27d32739c35adeab24ea0a97648b.png)


![\simeq\omega_0\left[A_{ab}\left(\cos\omega t-e^{-{{1}\over{2}}\Gamma t}\cos\omega_1 t\right)-A_{el}\left(\sin\omega t-e^{-{{1}\over{2}}\Gamma t}\sin\omega_1 t\right)\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/a/2/4/a249dcc0a03504ccb8884c5898015169.png)

![={{1}\over{2}}M\omega_0^2\Bigg\{\left[A_{ab}\left(\cos\omega t-e^{-{{1}\over{2}}\Gamma t}\cos\omega_1 t\right)-A_{el}\left(\sin\omega t-e^{-{{1}\over{2}}\Gamma t}\sin\omega_1 t\right)\right]^2+\;](http://upload.wikimedia.org/wikibooks/pl/math/7/0/3/70394009f9633c17ee0a283b4ca1da72.png)
![+\left[A_{ab}\left(\sin\omega t-e^{-{{1}\over{2}}\Gamma t}\sin\omega_1 t\right)+A_{el}\left(\cos\omega t-e^{-{{1}\over{2}}\Gamma t}\cos\omega_1 t\right)\right]^2\Bigg\}=\;](http://upload.wikimedia.org/wikibooks/pl/math/d/a/9/da91a30280bad9f777a6b4db550613ab.png)






![={{1}\over{2}}M\omega_0^2\left[A_{ab}^2+A_{el}^2+(A_{ab}^2+A_{el}^2)e^{-\Gamma t}-
2(A_{ab}^2+A_{el}^2)e^{-{{1}\over{2}}\Gamma t}\cos(\omega-\omega_1)t\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/1/5/3/153d9b45125a04d059aace7f76ef85b3.png)
![=
E\left[1+e^{-\Gamma t}-2e^{-{{1}\over{2}}\Gamma t}\cos(\omega-\omega_1)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/b/0/8/b084914f4d3f3e95c9585ebb006457bf.png)
![E(t)=E\left[1+e^{-\Gamma t}-2e^{-{{1}\over{2}}\Gamma t}\right]=E\left(1-e^{-{{1}\over{2}}\Gamma t}\right)^2\;](http://upload.wikimedia.org/wikibooks/pl/math/d/f/d/dfd213f16d5870eb07fc9222fc810c26.png)
![E(t)=E\left[1+e^{-\Gamma t}+2e^{-{{1}\over{2}}\Gamma t}\right]=E\left(1+e^{-{{1}\over{2}}\Gamma t}\right)^2\;](http://upload.wikimedia.org/wikibooks/pl/math/1/9/5/195f836598ff13454bce8ddf9170c3b6.png)




































![\omega^2=\omega_0^2+{{2K}\over{M}}\left[1-{{1}\over{2}}\left(e^{\chi a}+e^{-\chi a}\right)\right]=\omega_0^2+{{2K}\over{M}}\left(1-\cosh \chi a\right)\;](http://upload.wikimedia.org/wikibooks/pl/math/a/3/e/a3efe0be3b5bb9441221e4ed01c9f24c.png)





![\omega^2=\omega_0^2+{{2K}\over{M}}\left[1+{{1}\over{2}}\left(e^{\chi a}+e^{-\chi a}\right)\right]=\omega_0^2+{{2K}\over{M}}\left(1+\cosh ka\right)=\omega_0^2+{{4K}\over{M}}\cosh^2{{\chi a}\over{2}}\;](http://upload.wikimedia.org/wikibooks/pl/math/c/b/3/cb311c590490f970ff4aae5d3b828a9c.png)














































































![P(t)=F_x{{\partial\psi}\over{\partial t}}=\left[-T_0{{\partial\psi}\over{\partial z}}\right]{{\partial\psi}\over{\partial t}}=
\left[-T_0{{\partial\psi}\over{\partial z}}\right]\left[-v_f{{\partial\psi}\over{\partial z}}\right]={{v_f}\over{T_0}}\left[-T_0{{\partial \psi}\over{\partial z}}\right]^2={{1}\over{Z}}\left[-T_0{{\partial\psi}\over{\partial z}}\right]^2\;](http://upload.wikimedia.org/wikibooks/pl/math/a/e/2/ae298d73095e14d326507a4860a2c370.png)

![P(z,t)=Z\left[{{\partial\psi(z,t)}\over{\partial t}}\right]^2={{1}\over{Z}}\left[-Ka{{\partial\psi(z,t)}\over{\partial z}}\right]^2\;](http://upload.wikimedia.org/wikibooks/pl/math/f/c/2/fc2006817117c0d840ee309103f0dce5.png)


![I(z,t)=Z\left[{{\partial\psi(z,t)}\over{\partial t}}\right]^2={{1}\over{Z}}\left[-\kappa p_0{{\partial \psi(z,t)}\over{\partial z}}\right]^2\;](http://upload.wikimedia.org/wikibooks/pl/math/f/f/e/ffe0482d0fe8f96fc2c63634a8577091.png)










![{{\partial\psi_{odb}(0,t)}\over{\partial t}}=\left[{{Z_1-Z_2}\over{Z_1+Z_2}}\right]{{\partial\psi_{pad}(0,t)}\over{\partial t}}\;](http://upload.wikimedia.org/wikibooks/pl/math/d/9/9/d99ae54f0c6e2938a03b5d830c46c6e1.png)


![\psi_{odb}(z,t)=R_{12}A\cos\left[\omega\left(t+{{z}\over{v_f}}\right)\right]=R_{12}A\cos(\omega t+kz)\;](http://upload.wikimedia.org/wikibooks/pl/math/7/c/1/7c182849950054081c95fec7f821df15.png)

































![\psi(t)={{A(0)}\over{\Delta\omega}}\left[\delta\omega\cos\omega_1 t+\delta\omega\cos(\omega_1+\delta\omega)t+...+\delta\omega\cos\omega_2 t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/9/b/0/9b0384697c27a2e6cefe32faebb859db.png)



















![=\int_0^{\infty}e^{-{{1}\over{2}}\Gamma t}[\sin(\omega+\omega_1)t+\sin(\omega-\omega_1)t]dt=\;](http://upload.wikimedia.org/wikibooks/pl/math/8/b/2/8b2329b184b502ccdce2e33ab561aea2.png)




![=\int_0^{\infty}e^{-{{1}\over{2}}\Gamma}[\cos(\omega+\omega_1)t+\cos(\omega-\omega_1)t]dt=\;](http://upload.wikimedia.org/wikibooks/pl/math/c/c/4/cc4b039b0cb271da7b30b206e9749857.png)



![I(\omega)=[2\pi A(\omega)]^2+[2\pi B(\omega)]^2=\left({{2\omega(\omega^2-\omega_0^2)+\omega\Gamma^2}\over{(\omega_0^2-\omega^2)^2+\Gamma^2\omega^2}}\right)^2+\left({{\Gamma(\omega^2+\omega_0^2)}\over{(\omega_0^2-\omega)^2+\Gamma^2\omega^2}}\right)^2=\;](http://upload.wikimedia.org/wikibooks/pl/math/2/0/a/20a74c6ae4f9989b9b7f177f458d2ad1.png)


![=
{{4\omega^2\left[(\omega^2-\omega^2_0)^2+\Gamma^2\omega^2\right]+\Gamma^2\left[\omega^2\Gamma^2+(\omega^2-\omega_0^2)^2\right]}\over{\left((\omega_0^2-\omega)^2+\Gamma^2\omega^2\right)^2}}={{4\omega^2+\Gamma^2}\over{(\omega^2-\omega^2_0)^2+\Gamma^2\omega^2}}](http://upload.wikimedia.org/wikibooks/pl/math/4/b/6/4b6afa36f50d2d6741e04b1b16e72b51.png)
![\psi(0,t)=\int_{0}^{\infty}\left[A(\omega)\sin\omega t+B(\omega)\cos\omega t\right]d\omega\;](http://upload.wikimedia.org/wikibooks/pl/math/e/8/5/e859f014678568ecd882500139db10a4.png)
![\psi(z,t)=\int_{0}^{\infty}\left[A(\omega)\sin\omega \left(t-{{z}\over{v}}\right)+B(\omega)\cos\omega\left(t-{{z}\over{v}}\right)\right]d\omega\;](http://upload.wikimedia.org/wikibooks/pl/math/a/5/a/a5ae478b33c95bfcb32b7566e2d0bb82.png)
![\psi(z,t)=\int_{0}^{\infty}\left[A(\omega)\sin(\omega t-k(\omega)z)+B(\omega)\cos(\omega t-k(\omega)z)\right]d\omega\;](http://upload.wikimedia.org/wikibooks/pl/math/0/c/8/0c8345f9b50a115da891a18f9b72388c.png)
![\psi(z,t)=A\cos[\omega t-k(\omega)z]\;](http://upload.wikimedia.org/wikibooks/pl/math/c/7/6/c766b9fb834631715a4e714b55b95022.png)
![{{\partial^2\psi(z,t)}\over{\partial t^2}}=-\omega^2A\cos[\omega t-k(\omega)z]\;](http://upload.wikimedia.org/wikibooks/pl/math/8/c/7/8c7a17ff5e5f95832cd11b2b0d20ea51.png)
![{{\partial^2\psi(z,t)}\over{\partial z^2}}=-k^2A\cos[\omega t-k(\omega)z]\;](http://upload.wikimedia.org/wikibooks/pl/math/7/6/4/764c0cb41639b087284cf6dca5b7e186.png)































































![\psi_y=2A\cos\omega t\sin kx[k(y+h)]\;](http://upload.wikimedia.org/wikibooks/pl/math/e/a/9/ea9d1ab65e77e7532d44d10ea1bd76f2.png)

![F_x=-L\Delta y[p(x+\Delta x)-p(x)]=\;](http://upload.wikimedia.org/wikibooks/pl/math/3/1/0/310373bbb961e479c6e3f5de5387de29.png)
![=-L\Delta y\rho g\left[\psi_y(x+\delta x)-\psi_y(x)\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/1/b/8/1b855ff9f03c132cb365ab659217eb28.png)
![=L\delta y\Delta x\rho g{{\partial\psi_y}\over{\partial y}}=-(\Delta M)g\left[{{\partial\psi_y}\over{\partial y}}\right]_{y=0}\;](http://upload.wikimedia.org/wikibooks/pl/math/4/b/d/4bd0c6f1909e602002304ae19af7968a.png)
![F_x=(\Delta M)\left[{{\partial^2\psi_y}\over{\partial t^2}}\right]_{y=0}\;,](http://upload.wikimedia.org/wikibooks/pl/math/d/7/8/d78692c087ddfd76b6c157a411cbee14.png)




































![B_y(z,t)={{1}\over{2}}{{A}\over{c}}\cos(\omega t-kz)-{{1}\over{2}}{{A}\over{c}}\cos(\omega t+kz)={{-A}\over{c}}\sin\omega t\sin kz={{-1}\over{c}}E_x\left[z-{{1}\over{4}}\lambda,t-{{1}\over{4}}T\right]
\;](http://upload.wikimedia.org/wikibooks/pl/math/2/7/4/2745796fbe5a9968ec33a61315ef32c3.png)































![\hat{S}(\vec{r},t)={{1}\over{\mu_0}} \vec{E}\times\vec{B}={{1}\over{\mu_0}} \vec{E}\times\left(\hat{r}\times\vec{E}\right)={{1}\over{\mu_0c}}E^2\hat{r}={{1}\over{\mu_0}}\left[-{{q\vec{a}_{\bot}}\over{4\pi\epsilon_0 rc^2}} \right]^2\hat{r}={{\mu_0q^2\vec{a}_{\bot}^2}\over{16\pi^2cr^2}}\hat{r}](http://upload.wikimedia.org/wikibooks/pl/math/3/e/2/3e23166346349fe145aed623b2673540.png)






![{{1}\over{E}}\left[-{{dE}\over{dt}}\right]={{1}\over{\tau}}\;](http://upload.wikimedia.org/wikibooks/pl/math/b/7/6/b76c42686d77ae5d663f20d8b4ecb005.png)





















![\vec{\psi}(t)=\left[\hat{x}A\cos(\omega t-kz)+\hat{y}A\sin(\omega t-kz)\right]+\left[\hat{x}A\cos(\omega t+kz)+A\sin\left(\omega t+kz\right)\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/9/b/c/9bc6467befd199b37f491d47a99fb941.png)
![=\left[\hat{x}A(\cos(\omega t-kz)+\cos(\omega t+kz))+\hat{y}A(\sin(\omega t-kz)+\sin(\omega t+kz)\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/3/7/6/37636edce225a0ff0691937a096a2d65.png)
![=\left[\hat{x}A\cos\omega t+\hat{y}A\sin\omega t\right]\cos kz\;](http://upload.wikimedia.org/wikibooks/pl/math/3/e/3/3e3a800e8639e710321dc897c24c99d0.png)

![\vec{\psi}(t)=\left[\hat{x}A_1\cos(\omega t-kz)+\hat{y}A_2\sin(\omega t-kz)\right]+\left[\hat{x}A_1\cos(\omega t+kz)+A_2\sin\left(\omega t+kz\right)\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/c/0/7/c07471e37ff44b4fa356a34875cbcb55.png)
![=\left[\hat{x}A_1(\cos(\omega t-kz)+\cos(\omega t+kz))+\hat{y}A_2(\sin(\omega t-kz)+\sin(\omega t+kz)\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/8/5/6/856ccc30b16293517d13668a741a7112.png)
![=\left[\hat{x}A_1\cos\omega t+\hat{y}A_2\sin\omega t\right]\cos kz\;](http://upload.wikimedia.org/wikibooks/pl/math/d/a/d/dadfb8e5a55c775092abfdbd488fb45e.png)


















![\vec{E}=\hat{x}A\cos(\omega t-kz)={{A}\over{2}}\left\{\hat{x}\cos(\omega t-kz)+\hat{y}\cos\left[\left(\omega t-{{\pi}\over{2}}\right)-kz\right]\right\}+\;](http://upload.wikimedia.org/wikibooks/pl/math/4/5/6/4569ac4f7f1e9eee4aac550190e8a224.png)
![+{{A}\over{2}}\left\{\hat{x}\cos(\omega t-kz)-\hat{y}\cos\left[\left(\omega t-{{\pi}\over{2}}\right)\right]\right\}\;](http://upload.wikimedia.org/wikibooks/pl/math/a/e/b/aebab174ea65164e77ba878bceb11b45.png)
![\vec{E}=A\hat{x}e^{-i(\omega t-kz)}={{A}\over{2}}\left[\hat{x}e^{-i(\omega t-kz)}+\hat{y}e^{-i(\omega t-kz-{{\pi}\over{2}})}\right]+
{{A}\over{2}}\left[\hat{x}e^{-i(\omega t-kz)}-\hat{y}e^{-i(\omega t-kz-{{\pi}\over{2}})}\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/c/1/9/c194cb1e3a2d982f8804aedd1a4d89e2.png)














































![=
2Ae^{-i(\omega t-ks+\delta_0-{{\Delta\phi}\over{2}})}\cos\left({{\Delta\phi}\over{2}}\right)=
2Ae^{-i(\omega t-ks+\delta_0-{{\Delta\phi}\over{2}})}\cos \left[{{\pi}\over{\lambda}}d(\sin\alpha-\sin\beta)\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/3/4/f/34fa70aeddb41691d4a40a85fba727e7.png)






![=
2Ae^{-i(\omega t-ks+\delta_0-{{\Delta\phi}\over{2}})}\cos\left({{\Delta\phi}\over{2}}\right)=
2Ae^{-i(\omega t-ks+\delta_0-{{\Delta\phi}\over{2}})}\cos\left[{{\omega d}\over{c}}\left(n_1\sin\alpha-n_2\sin\beta\right)\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/9/7/a/97a90e4e9679210a5f11a43dcc619a4c.png)

































![|\psi_a|^2={{1}\over{2}}\left[1+\cos(\omega_1-\omega_2)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/0/b/3/0b321466ea870c9fb7e74c87f59ee1c6.png)
![|\psi_b|^2={{1}\over{2}}\left[1-\cos(\omega_1-\omega_2)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/b/8/2/b827fc4c4a73f0b152ca7460c7664b5c.png)
![|\psi(K^0_1)|^2={{1}\over{4}}\left[e^{-{{t}\over{\tau_1}}}+e^{-{{t}\over{\tau_2}}}+2e^{-{{1}\over{2}}\left({{t}\over{\tau_1}}+{{t}\over{\tau_2}}\right)}\cos(
\omega_1-\omega_2)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/5/c/e/5ce3f3ec3b3ca810e521851b0ed91b6d.png)
![|\psi(K^0_1)|^2={{1}\over{4}}\left[e^{-{{t}\over{\tau_1}}}+e^{-{{t}\over{\tau_2}}}-2e^{-{{1}\over{2}}\left({{t}\over{\tau_1}}+{{t}\over{\tau_2}}\right)}\cos(
\omega_1-\omega_2)t\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/3/f/e/3fe04003ba130500794e3d982da47c6d.png)

![{{d^2f(z)}\over{dz^2}}={{2m}\over{\hbar^2}}\left[V-\hbar\omega\right]f(z)\;](http://upload.wikimedia.org/wikibooks/pl/math/c/1/a/c1a2626c6ed32e2568863b9c3d0d30b3.png)
![\psi(z,t)=\left[A\sin kz+B\cos kz\right]e^{i\omega t}\;](http://upload.wikimedia.org/wikibooks/pl/math/3/0/8/30894b5a3788a41a3557c0456c29db37.png)



![1=\int_0^L|\psi|^2dz=A^2\int_0^L\sin^2kz dz={{1}\over{2}}A^2\int_0^L\left[1-\cos 2kz\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/c/0/b/c0b34934ac919531c5a0eb0cb4a3d977.png)
![{{1}\over{2}}A^2\left[1-{{1}\over{2k}}\cos2kz\Bigg|_0^{n{{\pi}\over{L}}}\right]={{1}\over{2}}A^2L\;](http://upload.wikimedia.org/wikibooks/pl/math/f/9/1/f911f5e14e7126317c6e1daefd058e8e.png)


![{{d^2f(z)}\over{dz^2}}={{2m}\over{\hbar^2}}\left[V(z)-\hbar\omega\right]f(z)\;](http://upload.wikimedia.org/wikibooks/pl/math/9/1/0/91073750ea11e5d56542335487c495e5.png)


















![\int_{-{{L}\over{2}}}^{{{L}\over{2}}}z|\psi|^2dz=A_1A_2\cos(\omega_2-\omega_1)t\int_{-{{L}\over{2}}}^{{{L}\over{2}}}z\left[\sin(k_2+k_1)z+\sin(k_2-k_1)z\right]dz=\;](http://upload.wikimedia.org/wikibooks/pl/math/4/0/e/40efd18164e2fca4fb7eb43d4b05163c.png)
![=-A_1A_2\cos(\omega_2-\omega_1)t\Bigg\{\left[z{{1}\over{k_1+k_2}}\cos(k_2+k_1)z+z{{1}\over{k_1-k_2}}\cos(k_2-k_1)z\right]_{-{{L}\over{2}}}^{{{L}\over{2}}}-\;,](http://upload.wikimedia.org/wikibooks/pl/math/f/d/4/fd48547fb2f32de99eff7f2c7cc829f1.png)
![-\left[{{1}\over{(k_1+k_2)^2}}\sin(k_2+k_1)z+{{1}\over{(k_2-k_1)^2}}\sin(k_2-k_1)z\right]_{-{{L}\over{2}}}^{{{L}\over{2}}}\Bigg\}=\;](http://upload.wikimedia.org/wikibooks/pl/math/7/1/6/7162b5ffe1b306c2570de1718e4400d2.png)
![=A_1A_2\cos(\omega_2-\omega_1)t\left[-{{2L^2}\over{9\pi^2}}+{{2L^2}\over{\pi^2}}\right]=A_1A_2\cos(\omega_2-\omega_1)t\left[{{-2L^2+18L^2}\over{9\pi^2}}\right]=\;](http://upload.wikimedia.org/wikibooks/pl/math/e/c/f/ecf64e9d84e5439f049eb6b812dbd6e3.png)



![|E_c(t)|^2=E_1^2+E_2^2+2E_1E_2\cos\left[(\omega_1-\omega_2)t+(\varphi_1-\varphi_2)\right]\;](http://upload.wikimedia.org/wikibooks/pl/math/1/a/6/1a6de28ae745ad927c63029545f1bbdd.png)









![(n_1-n_2)^2=[(n_1-\overline{n}_1)-(n_2-\overline{n}_2)]^2=(n_1-\overline{n}_2)^2+(n_2-\overline{n}_2)^2-2(n_1-\overline{n}_1)(n_2-\overline{n}_2)\;](http://upload.wikimedia.org/wikibooks/pl/math/4/b/b/4bbba3237538c24aeaf15321c66a1067.png)




















































![\psi(z,t)=\left[Ae^{ikz}+Be^{-ikz}\right]e^{-i\omega t}\;](http://upload.wikimedia.org/wikibooks/pl/math/7/2/6/726798ede9f666517a6451a44612660f.png)























































