rr struny zagadnienie Dirichleta
: 11 cze 2024, o 03:04
Dzień dobry proszę o sprawdzenie i pomoc w zrozumieniu rozwiązania zagadnienia Dirichleta struny. Jest ono niejednorodne. To są moje notatki, więc z pewnym pp są momentami błędne.
\(\displaystyle{ \begin{cases} u''_{tt}(x,t)=a^{2}u''_{xx}(x,t)\;\;\; (1) \\ u'_{x}(0,t)=u'_{x}(\pi ,t)=0 \;\;\; (2)\; brzegowe\\u(x,0)=\cos 2x \;\; u'_{t}(x,0)=7\cos 5x \;\; początkowe \end{cases} }\)
\(\displaystyle{ t}\) jest nieujemnym czasem, \(\displaystyle{ x}\) jest położeniem na strunie \(\displaystyle{ x\in [0,\pi ] }\)co znaczy, że struna ma długość \(\displaystyle{ \pi}\) .
Najpierw sprawdzamy czy warunki zgodności pasują, to ogarniam, że jak są niedobre to całe zagadnienie do kosza.
Brzegowe są oba zerowe pasują. Jeszcze z tego pierwszego początkowego liczymy pochodną cząstkową po iksach w punktach \(\displaystyle{ 0}\) i \(\displaystyle{ \pi}\). \(\displaystyle{ u'_{x}(x,0)=-2\sin 2x}\) i to jest oczywiście równe \(\displaystyle{ 0}\) w punktach \(\displaystyle{ 0}\) i\(\displaystyle{ \pi}\).
Pasują, idziemy dalej.
\(\displaystyle{ h(t,x)=\cos 3x}\) jest niejednorodnością. Jak zawsze najpierw rozwiązujemy rr jednorodne a potem dopiero niejednorodne.
Zakładamy, że rozwiązanie rr jednorodnego to \(\displaystyle{ u(x,t)=g(t)f(x)}\) czyli iloczyn dwóch funkcji i to mają być rozwiązania nietrywialne, czyli jak jakiś element iloczynu jest \(\displaystyle{ 0}\) to takie odrzucamy. Dalej \(\displaystyle{ A}\),\(\displaystyle{ B}\) są stałymi.
\(\displaystyle{ u_{tt}=g''(t)f(x)}\) i \(\displaystyle{ u_{xx}=a^{2}g(t)f''(x)}\)
\(\displaystyle{ g''(t)f(x)=a^{2}g(t)f''(x)}\) dzielimy stronami, bo interesują nas rzeczy (w jednym punkcie co najmniej) niezerowe.
\(\displaystyle{ \frac{g''(t)}{a^{2}g(t)} = \frac{f''(x)}{f(x)} }\)Tutaj zauważamy, że skoro lewa strona nie zależy od \(\displaystyle{ x}\) a prawa od \(\displaystyle{ t}\), to oba ułamki są stałe.
Więc między\(\displaystyle{ f''(x)}\) a \(\displaystyle{ f(x)}\) jest zależność liniowa \(\displaystyle{ f''(x)=\lambda f(x)}\), lambda jest stałą. Jest to rr liniowe autonomiczne 2 rzędu. Korzystamy z warunków brzegowych
\(\displaystyle{ f'(0)=0=f'(\pi )}\) są dwa, a więc można to rozwiązać. Ale lambda nie wiadomo jaka jest, więc rozważamy trzy przypadki.
Wielomian charakterystyczny jest taki sam \(\displaystyle{ r^{2}=\lambda }\). Gdy lambda jest:
1) dodatnia \(\displaystyle{ r= \pm sqrt{\lambda }}\)
ogólne: \(\displaystyle{ f(x)=Ae^{\sqrt { \lambda x}} +Be^{-\sqrt { \lambda x}}}\)
\(\displaystyle{ f'(0)=f'(\pi )=0}\)
\(\displaystyle{ f'(x)=A\sqrt { \lambda }e^{\sqrt { \lambda x}}-B\sqrt { \lambda }e^{-\sqrt { \lambda x}}}\)
\(\displaystyle{ f'(0)=0=A\sqrt{\lambda }-B\sqrt{\lambda }}\)
\(\displaystyle{ A=B}\)
\(\displaystyle{ 0=f'(\pi )=A\sqrt { \lambda }e^{\sqrt { \lambda \pi }}-A\sqrt { \lambda }e^{-\sqrt { \lambda \pi }}}\) funkcja wykładnicza i lambda są obie dodatnie więc
\(\displaystyle{ A=0}\) rozwiązanie trywialne odrzucam.
2) zerowa \(\displaystyle{ r=0}\) dwukrotnie -> \(\displaystyle{ f(x)=A+Bx}\)
\(\displaystyle{ f'(0)=f'(\pi )=0=B}\)
\(\displaystyle{ B}\) jest zerowe, \(\displaystyle{ A}\) dowolne niezerowe, więc \(\displaystyle{ f(x)}\) jest funkcją stałą niezerową.
Niech \(\displaystyle{ f(x)=f_{0}(x)=A\cos 0x}\) <- można tak napisać?
3) ujemna \(\displaystyle{ r= \pm sqrt{-\lambda}i}\) <- tu się dzieją sztuczki, chociaż jako tako ogarniam
\(\displaystyle{ f(x)=A\cos (\sqrt (-\lambda )x)+B\sin (\sqrt (-\lambda )x)}\) tak się rozwiązuje gdy wartość własna jest zespolona.
\(\displaystyle{ f'(0)=f'(\pi )=0}\)
\(\displaystyle{ f'(x)=-A\cos(\sqrt {\lambda }x)+B\cos (\sqrt{-\lambda }\pi )}\)
\(\displaystyle{ B=0}\)
\(\displaystyle{ A}\) dowolne (niezerowe?) gdy \(\displaystyle{ \sin (\sqrt{-\lambda }\pi )=0}\)
\(\displaystyle{ \sqrt{-\lambda }=k}\), \(\displaystyle{ k\in \mathbb{N}}\)
\(\displaystyle{ \lambda =-k^{2}}\)
\(\displaystyle{ f_{k}(x)= A\cos (\sqrt{ -k^{2}} x) =f_{n}=A\cos (nx)}\)
I tutaj się okazało, że nie ma sinusa w funkcji \(\displaystyle{ f}\). Czy zawsze tak będzie w rr struny?
Okeeeej no i tutaj dzieje się magia. Bo my chcemy znaleźć \(\displaystyle{ u}\), a wiadomo, że rozwiązanie ogólne to kombinacja liniowa rozwiązań pojedynczych, ale co to w ogóle jest. Co się stało z numeracją to skomplikowane, ale tu jest (6). Rozumiem, że skoro \(\displaystyle{ f_{0}}\) było funkcją stałą, to można ją zespolić z \(\displaystyle{ g_{0}}\) i pominąć to oznaczenie.
\(\displaystyle{ u(x,t)=g_{0}(t)+ \sum_{n=1}^{\infty } g_{n}\cos (nx)}\) \(\displaystyle{ (6)}\)
\(\displaystyle{ g}\) z indeksem \(\displaystyle{ 0}\) rozdzieliliśmy, bo będziemy je oddzielnie liczyć. Idea jest taka, że \(\displaystyle{ f_{n}}\) ma być najprostszą możliwą funkcją, a wszelkie dodatkowe rzeczy mają być w \(\displaystyle{ g_{n}}\)???
Wracamy do rr \(\displaystyle{ (1)}\) niejednorodnego. Rozwijamy \(\displaystyle{ h(t,x)}\) w szereg fouriera względem układu \(\displaystyle{ (cosnx)^{\infty }_{n=0}}\), ale \(\displaystyle{ h}\) jest domyślnie rozwinięta jak jest cosinusem \(\displaystyle{ h(t,x)=f_{3}(x)}\). To akurat widać.
Liczymy pochodne cząstkowe drugiego rzędu z \(\displaystyle{ u}\) żeby wstawić do \(\displaystyle{ (1)}\). To ma sens no bo jak inaczej.
\(\displaystyle{ u''_{tt}(x,t)=g''_{0}(t)+ \sum_{n=1}^{\infty } g''_{n}(t)\cos (nx) }\)
\(\displaystyle{ u''_{xx}(x,t)= \sum_{n=1}^{\infty }g_{n}(t)(-n^{2}\cos (nx)) }\)
Oba razem z \(\displaystyle{ (1)}\)
\(\displaystyle{ g''_{0}(t)+ \sum_{n=1}^{\infty } g''_{n}(t)\cos (nx) =a^{2}u''_{xx}(x,t)= \sum_{n=1}^{\infty }g_{n}(t)(-n^{2}\cos (nx)) +7\cos 3x}\)
Teraz korzystamy z własności, że \(\displaystyle{ \forall x\in D \;\; af(x)=bf(x) \Rightarrow a=b}\). Mogę tak napisać? Muszę zakładać coś poza tych, że f\(\displaystyle{ }\) nie jest stale równa \(\displaystyle{ 0}\)? Oddzielnie traktujemy indeksy \(\displaystyle{ n}\), które wystąpiły wcześniej a całą resztę grupowo. W tym momencie wystąpiło prawdopodobnie nieprzeliczalnie dużo błędów rachunkowych, także proszę mnie sprawdzić.
Właśnie to całe porównywanie współczynników jest najgorsze i ja nie ogarniam.
\(\displaystyle{ n=0}\): \(\displaystyle{ g''_{0}(t)=0}\) funkcje liniowe czyli\(\displaystyle{ g_{0}=A_{0}+B_{0}t}\)
\(\displaystyle{ n=3}\): \(\displaystyle{ g''_{3}(t)=-a^{2}\cdot 3^{3}g_{3}(t)+1\;\;(7)}\) liniowe rzędu 2 niejednorodne
\(\displaystyle{ n}\) ogólne:\(\displaystyle{ g''_{n}(t)=-a^{2}n^{2}g_{n}(t)\; (8)}\) rr liniowe autonomiczne drugiego rzędu i już były analogiczne przykłady, więc analogicznie liczymy.
\(\displaystyle{ (8)\; -a^{2}n^{2}<0}\) no bo zarówno \(\displaystyle{ a= \neq 0}\) i \(\displaystyle{ n \neq 0}\)
\(\displaystyle{ g_{n}(t)=A_{n}\cos (ant) + B_{n}\sin (ant)}\)
co do \(\displaystyle{ (7)}\) to jak wyżej
\(\displaystyle{ g_{3oj}=A_{3}\cos (3at) +B_{3}\sin (3at)}\)
Potrzebujemy dodać rozwiązanie szczególne. Jako że niejednorodność to funkcja stała to przewiduję, że rozwiązanie szczególne to też jakaś stała \(\displaystyle{ c}\).
\(\displaystyle{ c''=-a^{2}\cdot 9 \cdot c +1}\)
\(\displaystyle{ c= \frac{1}{9a^{2}} }\)
więc \(\displaystyle{ g_{3}=A_{3}\cos 3at++B_{3}\sin (3at)+\frac{1}{9a^{2}}
}\)
I to wszystko podstawiamy do \(\displaystyle{ (6)}\). Dlaczego \(\displaystyle{ g_{3}}\) jest mnożone przez \(\displaystyle{ \cos 3x}\) ??
\(\displaystyle{ (9)}\)
\(\displaystyle{ u(x,t)=A_{0}+B_{0}t+ \sum_{n=1}^{\infty } (A_{n}\cos (ant) +B_{n}\sin (ant))\cos (nx) + \frac{1}{9a^{2}}\cos 3x }\)
Wracamy do warunków początkowych
\(\displaystyle{ u(x,0)=\cos 2x}\), \(\displaystyle{ u'_{t}(x,t)=7\cos (5x)}\)
dokończę potem
\(\displaystyle{ \begin{cases} u''_{tt}(x,t)=a^{2}u''_{xx}(x,t)\;\;\; (1) \\ u'_{x}(0,t)=u'_{x}(\pi ,t)=0 \;\;\; (2)\; brzegowe\\u(x,0)=\cos 2x \;\; u'_{t}(x,0)=7\cos 5x \;\; początkowe \end{cases} }\)
\(\displaystyle{ t}\) jest nieujemnym czasem, \(\displaystyle{ x}\) jest położeniem na strunie \(\displaystyle{ x\in [0,\pi ] }\)co znaczy, że struna ma długość \(\displaystyle{ \pi}\) .
Najpierw sprawdzamy czy warunki zgodności pasują, to ogarniam, że jak są niedobre to całe zagadnienie do kosza.
Brzegowe są oba zerowe pasują. Jeszcze z tego pierwszego początkowego liczymy pochodną cząstkową po iksach w punktach \(\displaystyle{ 0}\) i \(\displaystyle{ \pi}\). \(\displaystyle{ u'_{x}(x,0)=-2\sin 2x}\) i to jest oczywiście równe \(\displaystyle{ 0}\) w punktach \(\displaystyle{ 0}\) i\(\displaystyle{ \pi}\).
Pasują, idziemy dalej.
\(\displaystyle{ h(t,x)=\cos 3x}\) jest niejednorodnością. Jak zawsze najpierw rozwiązujemy rr jednorodne a potem dopiero niejednorodne.
Zakładamy, że rozwiązanie rr jednorodnego to \(\displaystyle{ u(x,t)=g(t)f(x)}\) czyli iloczyn dwóch funkcji i to mają być rozwiązania nietrywialne, czyli jak jakiś element iloczynu jest \(\displaystyle{ 0}\) to takie odrzucamy. Dalej \(\displaystyle{ A}\),\(\displaystyle{ B}\) są stałymi.
\(\displaystyle{ u_{tt}=g''(t)f(x)}\) i \(\displaystyle{ u_{xx}=a^{2}g(t)f''(x)}\)
\(\displaystyle{ g''(t)f(x)=a^{2}g(t)f''(x)}\) dzielimy stronami, bo interesują nas rzeczy (w jednym punkcie co najmniej) niezerowe.
\(\displaystyle{ \frac{g''(t)}{a^{2}g(t)} = \frac{f''(x)}{f(x)} }\)Tutaj zauważamy, że skoro lewa strona nie zależy od \(\displaystyle{ x}\) a prawa od \(\displaystyle{ t}\), to oba ułamki są stałe.
Więc między\(\displaystyle{ f''(x)}\) a \(\displaystyle{ f(x)}\) jest zależność liniowa \(\displaystyle{ f''(x)=\lambda f(x)}\), lambda jest stałą. Jest to rr liniowe autonomiczne 2 rzędu. Korzystamy z warunków brzegowych
\(\displaystyle{ f'(0)=0=f'(\pi )}\) są dwa, a więc można to rozwiązać. Ale lambda nie wiadomo jaka jest, więc rozważamy trzy przypadki.
Wielomian charakterystyczny jest taki sam \(\displaystyle{ r^{2}=\lambda }\). Gdy lambda jest:
1) dodatnia \(\displaystyle{ r= \pm sqrt{\lambda }}\)
ogólne: \(\displaystyle{ f(x)=Ae^{\sqrt { \lambda x}} +Be^{-\sqrt { \lambda x}}}\)
\(\displaystyle{ f'(0)=f'(\pi )=0}\)
\(\displaystyle{ f'(x)=A\sqrt { \lambda }e^{\sqrt { \lambda x}}-B\sqrt { \lambda }e^{-\sqrt { \lambda x}}}\)
\(\displaystyle{ f'(0)=0=A\sqrt{\lambda }-B\sqrt{\lambda }}\)
\(\displaystyle{ A=B}\)
\(\displaystyle{ 0=f'(\pi )=A\sqrt { \lambda }e^{\sqrt { \lambda \pi }}-A\sqrt { \lambda }e^{-\sqrt { \lambda \pi }}}\) funkcja wykładnicza i lambda są obie dodatnie więc
\(\displaystyle{ A=0}\) rozwiązanie trywialne odrzucam.
2) zerowa \(\displaystyle{ r=0}\) dwukrotnie -> \(\displaystyle{ f(x)=A+Bx}\)
\(\displaystyle{ f'(0)=f'(\pi )=0=B}\)
\(\displaystyle{ B}\) jest zerowe, \(\displaystyle{ A}\) dowolne niezerowe, więc \(\displaystyle{ f(x)}\) jest funkcją stałą niezerową.
Niech \(\displaystyle{ f(x)=f_{0}(x)=A\cos 0x}\) <- można tak napisać?
3) ujemna \(\displaystyle{ r= \pm sqrt{-\lambda}i}\) <- tu się dzieją sztuczki, chociaż jako tako ogarniam
\(\displaystyle{ f(x)=A\cos (\sqrt (-\lambda )x)+B\sin (\sqrt (-\lambda )x)}\) tak się rozwiązuje gdy wartość własna jest zespolona.
\(\displaystyle{ f'(0)=f'(\pi )=0}\)
\(\displaystyle{ f'(x)=-A\cos(\sqrt {\lambda }x)+B\cos (\sqrt{-\lambda }\pi )}\)
\(\displaystyle{ B=0}\)
\(\displaystyle{ A}\) dowolne (niezerowe?) gdy \(\displaystyle{ \sin (\sqrt{-\lambda }\pi )=0}\)
\(\displaystyle{ \sqrt{-\lambda }=k}\), \(\displaystyle{ k\in \mathbb{N}}\)
\(\displaystyle{ \lambda =-k^{2}}\)
\(\displaystyle{ f_{k}(x)= A\cos (\sqrt{ -k^{2}} x) =f_{n}=A\cos (nx)}\)
I tutaj się okazało, że nie ma sinusa w funkcji \(\displaystyle{ f}\). Czy zawsze tak będzie w rr struny?
Okeeeej no i tutaj dzieje się magia. Bo my chcemy znaleźć \(\displaystyle{ u}\), a wiadomo, że rozwiązanie ogólne to kombinacja liniowa rozwiązań pojedynczych, ale co to w ogóle jest. Co się stało z numeracją to skomplikowane, ale tu jest (6). Rozumiem, że skoro \(\displaystyle{ f_{0}}\) było funkcją stałą, to można ją zespolić z \(\displaystyle{ g_{0}}\) i pominąć to oznaczenie.
\(\displaystyle{ u(x,t)=g_{0}(t)+ \sum_{n=1}^{\infty } g_{n}\cos (nx)}\) \(\displaystyle{ (6)}\)
\(\displaystyle{ g}\) z indeksem \(\displaystyle{ 0}\) rozdzieliliśmy, bo będziemy je oddzielnie liczyć. Idea jest taka, że \(\displaystyle{ f_{n}}\) ma być najprostszą możliwą funkcją, a wszelkie dodatkowe rzeczy mają być w \(\displaystyle{ g_{n}}\)???
Wracamy do rr \(\displaystyle{ (1)}\) niejednorodnego. Rozwijamy \(\displaystyle{ h(t,x)}\) w szereg fouriera względem układu \(\displaystyle{ (cosnx)^{\infty }_{n=0}}\), ale \(\displaystyle{ h}\) jest domyślnie rozwinięta jak jest cosinusem \(\displaystyle{ h(t,x)=f_{3}(x)}\). To akurat widać.
Liczymy pochodne cząstkowe drugiego rzędu z \(\displaystyle{ u}\) żeby wstawić do \(\displaystyle{ (1)}\). To ma sens no bo jak inaczej.
\(\displaystyle{ u''_{tt}(x,t)=g''_{0}(t)+ \sum_{n=1}^{\infty } g''_{n}(t)\cos (nx) }\)
\(\displaystyle{ u''_{xx}(x,t)= \sum_{n=1}^{\infty }g_{n}(t)(-n^{2}\cos (nx)) }\)
Oba razem z \(\displaystyle{ (1)}\)
\(\displaystyle{ g''_{0}(t)+ \sum_{n=1}^{\infty } g''_{n}(t)\cos (nx) =a^{2}u''_{xx}(x,t)= \sum_{n=1}^{\infty }g_{n}(t)(-n^{2}\cos (nx)) +7\cos 3x}\)
Teraz korzystamy z własności, że \(\displaystyle{ \forall x\in D \;\; af(x)=bf(x) \Rightarrow a=b}\). Mogę tak napisać? Muszę zakładać coś poza tych, że f\(\displaystyle{ }\) nie jest stale równa \(\displaystyle{ 0}\)? Oddzielnie traktujemy indeksy \(\displaystyle{ n}\), które wystąpiły wcześniej a całą resztę grupowo. W tym momencie wystąpiło prawdopodobnie nieprzeliczalnie dużo błędów rachunkowych, także proszę mnie sprawdzić.
Właśnie to całe porównywanie współczynników jest najgorsze i ja nie ogarniam.
\(\displaystyle{ n=0}\): \(\displaystyle{ g''_{0}(t)=0}\) funkcje liniowe czyli\(\displaystyle{ g_{0}=A_{0}+B_{0}t}\)
\(\displaystyle{ n=3}\): \(\displaystyle{ g''_{3}(t)=-a^{2}\cdot 3^{3}g_{3}(t)+1\;\;(7)}\) liniowe rzędu 2 niejednorodne
\(\displaystyle{ n}\) ogólne:\(\displaystyle{ g''_{n}(t)=-a^{2}n^{2}g_{n}(t)\; (8)}\) rr liniowe autonomiczne drugiego rzędu i już były analogiczne przykłady, więc analogicznie liczymy.
\(\displaystyle{ (8)\; -a^{2}n^{2}<0}\) no bo zarówno \(\displaystyle{ a= \neq 0}\) i \(\displaystyle{ n \neq 0}\)
\(\displaystyle{ g_{n}(t)=A_{n}\cos (ant) + B_{n}\sin (ant)}\)
co do \(\displaystyle{ (7)}\) to jak wyżej
\(\displaystyle{ g_{3oj}=A_{3}\cos (3at) +B_{3}\sin (3at)}\)
Potrzebujemy dodać rozwiązanie szczególne. Jako że niejednorodność to funkcja stała to przewiduję, że rozwiązanie szczególne to też jakaś stała \(\displaystyle{ c}\).
\(\displaystyle{ c''=-a^{2}\cdot 9 \cdot c +1}\)
\(\displaystyle{ c= \frac{1}{9a^{2}} }\)
więc \(\displaystyle{ g_{3}=A_{3}\cos 3at++B_{3}\sin (3at)+\frac{1}{9a^{2}}
}\)
I to wszystko podstawiamy do \(\displaystyle{ (6)}\). Dlaczego \(\displaystyle{ g_{3}}\) jest mnożone przez \(\displaystyle{ \cos 3x}\) ??
\(\displaystyle{ (9)}\)
\(\displaystyle{ u(x,t)=A_{0}+B_{0}t+ \sum_{n=1}^{\infty } (A_{n}\cos (ant) +B_{n}\sin (ant))\cos (nx) + \frac{1}{9a^{2}}\cos 3x }\)
Wracamy do warunków początkowych
\(\displaystyle{ u(x,0)=\cos 2x}\), \(\displaystyle{ u'_{t}(x,t)=7\cos (5x)}\)
dokończę potem