• Nie Znaleziono Wyników

Spistreści Równaniafalowe

N/A
N/A
Protected

Academic year: 2021

Share "Spistreści Równaniafalowe"

Copied!
15
0
0

Pełen tekst

(1)

Równania falowe

Spis treści

1 Przykładowe rozwiązania 2

2 Zestaw zadań do ćwiczenia samodzielnego 12

(2)

1 Przykładowe rozwiązania

Zadanie 1.1.

Rozważmy problem drgania struny jednorodnej o zamocowanych końcach:

uxx = 1 a2utt,

gdzie a jest stała różn a od 0. Jest to rodzaj równania falowego, przy czym funkcja u(x, t) jest funkcj a wychylenia struny z położenia równowagi (w położeniu równowagi zajmuje ona przedział 0, π na osi x). Określamy warunki brzegowe

u(0, t) = u(π, t) = 0, (1)

co oznacza, że końce struny sa nieruchome, oraz warunki pocz atkowe

u(x, 0) = φ(x), ut(x, 0) = ψ(x) (2)

dla t  0 i x ∈ 0, π , przy czym równości

φ(0) = φ(π) = 0, φ(0) = φ(π) = 0, ψ(0) = ψ(π) = 0

sa warunkami zgodności mi edzy warunkami pocz atkowymi (2) a brzegowymi (1). Funkcja φ(x) określa poczatkowy kształt struny, a ψ(x) określa rozkład pr edkości pocz atkowej.

Znaleźć rozwiazania równania. Co należy założyć o funkcjach φ i ψ, aby uzyskać rozwi azanie w sensie klasycznym?

Metoda szukania nietrywialnego rozwiazania polega na zostosowaniu techniki rozdzielenia zmiennych i poszukiwania rozwiazań u(x, t) w postaci:

u(x, t) =

 k=1

vk(x) · wk(t).

Jeśli każdy składnik tej sumy spełnia równanie różniczkowe, to vk(x)wk(t) = 1

a2vk(x)wk(t), vk(x)

vk(x) = 1 a2

wk(t) wk(t).

Stosunek ten jest zatem stały. Niech wynosi λk. Wtedy dostajemy dwa równania:

vk(x)

vk(x) = λk i 1 a2

wk(t)

wk(t) = λk. (3)

(3)

Jeśli zażadamy, by każdy składnik sumy spełniał warunki brzegowe Dirichleta (1), to ponieważ 0 = vk(0)wk(t), to vk(0) = 0, a ponieważ 0 = vk(π)wk(t), to vk(π) = 0. Stad wniosek, że funkcja v k

jest rozwiazaniem równania różniczkowego zwyczajnego z warunkami brzegowymi typu Dirichleta: vk(x) − λkvk(x) = 0, vk(0) = vk(π) = 0. (4) Oczywiście takie rozwiazanie istnieje dla dowolnej stałej λ k : vk(x) ≡ 0. Interesuja nas takie warto- ściλk, dla których mamy rozwiazania nietrywialne. Niech v k(x) = eαx. Wtedy z równania (4) mamy α2eαx − λkeαx = 0, czyli α2 = λk.

Jeśli λk > 0,to rozwiazaniem s a funkcje v k(x) = C1eλkx + C2eλkx. Z warunków brzegowych (4) otrzymamy C1 = C2 = 0 i stad v k(x) ≡ 0. Podobnie dla λk = 0. Jeśli natomiast λk < 0, to vk(x) = C1cos√

−λkx + C2sin√

−λkx. Z warunku vk(0) = 0 mamy C1 = 0, a z warunku vk(π) = 0 mamy

C2sin

−λkπ = 0.

Jeśli wiec chcemy, by C 2 = 0, to musimy wziać λ k = −k2, n = 1, 2, ... , by funkcje vk(x) = sin kx były nietrywialnymi rozwiazaniami (4). Zatem szukana postać rozwi azania jest teraz nast epuj aca:

u(x, t) =



k=1sin kxwk(t).

Zajmijmy sie teraz drugim równaniem (3) z λ k = −k2. Ma ono postać wk(t) = −a2k2wk(t),

czyli jest to znowu równanie różniczkowe liniowe drugiego rzedu o stałych współczynnikach. Mamy, jak poprzednio:

wk(t) = Cksin akt + Dkcos akt, (5)

gdzie Ck, Dk ∈ R sa dowolnymi stałymi. Trzeba tak dobrać te stałe, by funkcja u (ostatnio wyliczona) spełniała warunki poczatkowe (2). Otrzymujemy:

φ(x) = u(x, 0) =

k=1

sin kxwk(0),

ψ(x) = ut(x, 0) =

 k=1

sin kxwk(0).

Widać wiec, że w k(0) i wk(0) sa współczynnikami rozwini ecia w szereg Fouriera funkcji φ i ψ, odpowiednio, czyli

wk(0) = 2 π

 π

0 φ(t) sin kt dt, wk(0) = 2 π

 π

0 ψ(t) sin kt dt. (6)

(4)

Wykorzystujac teraz (5) i (6), mamy

wk(0) = Cksin 0 + Dkcos 0 = Dk, wk(0) = akCkcos 0 − akDksin 0 = akCk, wiec D k = wk(0) i Ck = wkak(0). Stad otrzymujemy postać funkcji u

u(x, t) =

 k=1

sin kx

wk(0)

ak sin akt + wk(0) cos akt



,

gdzie wk(0) i wk(0) sa dane wzorami (6). Jest to oczywiście tylko kandydat na rozwi azanie, bo nie wiemy jeszcze, czy ta funkcja jest ciagła oraz czy szereg (wyst epuj acy po prawej stronie) dopusz- cza dwukrotne różniczkowanie wyraz po wyrazie z zachowaniem niemal jednostajnej zbieżności dla t  0, x ∈ 0, π .

Załóżmy, że φ i ψ daja si e rozszerzyć na cał a oś rzeczywist a do funkcji okresowych o okresie 2π, nieparzystych i φ ∈ C4, ψ ∈ C3. Wtedy całkujac cztery razy przez cz eści, dostajemy

wk(0) = 2 πk4

 π

0 φ(4)(t) cos kt dt, czyli dla pewnych A, B > 0

|wk(0)|  2

πk4π sup |φ(4)(t)| = A k4 i analogicznie, całkujac trzy razy przez cz eści, możemy oszacować:

|wk(0)|  2π

πk3π sup |ψ(3)(t)| = B k3. W takim razie możemy oszacować wyrazy szeregu



sin kx

wk(0)

ak sin akt + wk(0) cos akt B ak4 + A

k4  B + aA k4 ,

a szeregk=1 k14 jest zbieżny jako harmoniczny rzedu 4. Zatem na podstawie kryterium Weierstrassa szereg wyjściowy jest jednostajnie zbieżny. Analogicznie można oszacować pochodne (wzgledem x) wyrazów szeregu:



k cos kx

wk(0)

ak sin akt + wk(0) cos akt B + aA ak3

i otrzymamy jednostajna zbieżność szeregu pochodnych, a wi ec możliwość różniczkowania szeregu wyraz po wyrazie. Dla kolejnego różniczkowania mamy



−k2sin kx

wk(0)

ak sin akt + wk(0) cos akt B + aA ak2 .

Analogicznie postepujemy, by wykazać możliwość dwukrotnego różniczkowania szeregu wyraz po wyrazie wzgledem t:



sin kx

wk(0)

ak ak cos akt − wk(0)ak sin akt  B + aA k3 ,

(5)

sin kx−akwk(0) sin akt(ak)2wk(0) cos akt  B + aA k2 . Zatem

u(x, t) =

 k=1sin kx

wk(0)

ak sin akt + wk(0) cos akt



jest rozwiazaniem równania w sensie klasycznym, u ∈ C 2 i jest to rozwiazanie jedyne.

Warunek o możliwości odpowiedniego rozszerzenia funkcji φ i ψ jest spełniony, gdy φ ∈ C4((0, π)) i ψ ∈ C3((0, π)) oraz np. istnieja granice:

x→0lim+

djφ(x)

dxj , lim

x→π

djφ(x)

dxj , dla 0  j  4,

x→0lim+

djψ(x)

dxj , lim

x→π

djψ(x)

dxj , dla 0  j  3, przy czym zachodza równości:

x→0lim+φ(x) = lim

x→0+

d2φ(x)

dx2 = lim

x→0+

d4φ(x) dx4 = 0,

x→0lim+ψ(x) = lim

x→0+

d2ψ(x) dx2 = 0.

Założenia na funkcje φ i ψ sa dosyć restrykcyjne i nie stanowi a one warunku koniecznego, by rozwiazanie u było klasyczne. Można je osłabić, inaczej wykonuj ac szacowania (por.[17]).

Zadanie 1.2.

Zastosować metode rozdzielania zmiennych do równania niejednorodnego: utt = a2uxx + Ae−tcos π

2lx (7)

przy warunku brzegowym

ux(0, t) = u(l, t) = 0 (8)

i poczatkowym

u(x, 0) = ut(x, 0) = 0 (9)

w półpasie 0  x  l, t  0.

(6)

Zauważmy od razu, że warunki poczatkowe i brzegowe s a zgodne. Jeśli zaniedbamy funkcj e Ae −tcos 2lπx, to rozwiazań możemy poszukiwać, jak poprzednio w postaci v (x)·w (t), czyli dojdziemy do tożsamości v (x)w(t) = a2v(x)w (t), która oznacza

v(x)

v (x) = λ = 1

a2 · w(t) w (t) .

Otrzymujemy wiec równanie v (x)λv (x) = 0, dla którego wyznaczymy warunki brzegowe z (8):

0 = u(l, t) = v (l)w (t) =⇒ v (l) = 0, 0 = ux(0, t) = v(0)w (t) =⇒ v(0) = 0.

Rozwiażemy wi ec nast epuj ace równanie:

v(x)λv (x) = 0, v (l) = 0 = v(0). (10) Jak poprzednio, dla λ  0 dostajemy tylko rozwiazania zerowe, a dla λ < 0 mamy:

v (x) = C1cos√

−λx + C2sin√

−λx, zatem v(0) = √

−λC2 = 0, gdy C2 = 0, a v (l) = C1cos√

−λl = 0, gdy

−λl = π2 + 2kπ, k jest liczba naturaln a. St ad mamy wartości własne: λ k = −(2k+1)4l22π2 i funkcje własne:

vk(x) = cos2k+12l πx , k = 0, 1, 2, ... . Bedziemy wi ec poszukiwać rozwi azań problemu brzegowo- poczatkowego w postaci

u(x, t) =



k=0wk(t) cos

2k + 1 2l πx



.

Zakładajac, że dozwolone jest różniczkowanie szeregu wyraz po wyrazie, to z (7)



k=0vk(x)wk(t) =



k=0a2vk(x)wk(t) + Ae−tcos π 2lx, czyli



k=0wk(t) cos

2k + 1 2l πx



= −a2

 k=0

2k + 1 2l π

2

wk(t) cos

2k + 1 2l πx



+ Ae−t cos π

2lx, (11) Zauważmy, że ostatni składnik po prawej stronie tego równania można zapisać jako:

Ae−t cos π

2lx = Ae−tcos2k + 1

2l πx, dla k = 0.

Zatem równanie (11) można zapisać w postaci układu, porównujac kolejno wyrazy szeregów, dla k = 0 :

w0(t) cos

2 · 0 + 1

2l πx = −a2

2 · 0 + 1

2l π 2w0(t) cos

2 · 0 + 1

2l πx + Ae−tcos2 · 0 + 1 2l πx

(7)

i dla k > 0 :

wk(t) cos

2k + 1 2l πx



= −a2

2k + 1 2l π

2

wk(t) cos

2k + 1 2l πx



. Po skóceniu otrzymujey:

w0(t) = −a2

1

2lπ 2w0(t) + Ae−t, wk(t) = −a2

2k + 1 2l π

2

wk(t).

Z warunków poczatkowych (9) mamy teraz 0 = u(x, 0) =



k=0wk(0) cos

2k + 1 2l πx



,

0 = ut(x, 0) =



k=0wk(0) cos

2k + 1 2l πx



,

czyli wk(0) = wk(0) = 0 jako rozwiniecia w szereg Fouriera funkcji zerowych. Otrzymujemy wi ec do rozwiazania dwa rówania:

w0(t) + a2

1

2lπ 2w0(t) = Ae−t, w0(0) = w0(0) = 0 (12) i

wk(t) + a2

2k + 1 2l π

2

wk(t) = 0, wk(0) = wk(0) = 0. (13) Zajmijmy sie najpierw (13). Jest to równanie jednorodne, a szukan a funkcj a jest

wk(t) = B1cos

2k + 1 2l πat



+ B2sin

2k + 1 2l πat



.

Z warunków poczatkowych mamy wtedy 0 = w k(0) = B1 oraz 0 = wk(0) = B22k+12l , czyli B2 = 0.

Zatem rozwiazaniem tego równania jest w k(t) ≡ 0 dla k > 0. Rozwiażemy teraz (12). Zauważmy, że rozwiazaniem równania jednorodnego jest

˜

w0(t) = B1cos

2l t

+ B2sin

2l t

.

Mo zemy teraz zastosować technike przewidywań lub uzmienniania stałych. Przewidujemy rozwi azanie szczególne w postaci:

¯

w0(t) = Me−t

dla pewnej stałej M. Wtedy ¯w0(t) = ¯w0(t) i po wstawieniu do równania (12) dostajemy:

Me−t +

2l

2

Me−t = Ae−t,

(8)

czyli stała M musi być równa (2l)A(2l)2+(aπ)2 2. Całka szczególn a równania jest wi ec

¯

w0(t) = A(2l)2

(2l)2+ (aπ)2e−t.

Ostateczne rozwiazanie jest sum a rozw azania ogólnego i szczególnego, czyli w0(t) = ˜w0(t) + ¯w0(t) = B1cos

2l t

+ B2sin

2lt

+ A(2l)2

(2l)2+ (aπ)2e−t. Z warunków poczatkowych, dostajemy teraz

B1= − A(2l)2

(2l)2+ (aπ)2, B2 = A(2l)3 ((2l)2+ (aπ)2) aπ.

Ostatecznie postać szeregu w rozwiazaniu u redukuje si e tylko do składnika z k = 0, wi ec u(x, t) = w0(t) cos

1

2lπx = (2l)2A

(2l)2+ (aπ)2 cos

1

2lπx − cos 2l t

+ 2l sin

2l t

+ e−t

. Jak widać jest to rozwiazanie klasyczne (nie ma szeregu, a wszystkie funkcje wyst epuj ace po prawej stronie sa klasy C .)

–300 –200 –100 100 200 300

–4 –2

2 4

t

–40 –50 –20 –30

–10 x

Wykres rozwiązania

Na koniec przypomnijmy tylko, że wspomniana wcześniej matoda uzmienniania stłych dla równania (12) polega na rozwiazaniu układu:

0 = B1(t) cos2lt + B2(t) sin2lt ,

Ae−t = B1(t)2l sin2lt + B2(t)2l cos2lt i scałkowaniu znalezionych B1 i B2.

(9)

Zadanie 1.3.

Rozwiazać zagadnienie

uxx = utt (14)

z warunkami poczatkowymi

u(x, 0) = sin x, ut(x, 0) = 0 (15)

i niejednorodnymi warunkami brzegowymi

u(0, t) = t2, u(π, t) = t3 (16)

dla 0  x  π i t  0.

Należy poszukiwać rozwiazań w postaci

u(x, t) = W (x, t) + V (x, t),

gdzie funkcja W (x, t) = r0(x)t2 + r1(x)t3 jest tak dobrana, by spełniała warunki brzegowe niejed- norodne (czyli wtedy warunki na V bed a jednorodne i b edzie można j a znaleźć metod a rozdzielania zmiennych). Zwykle r0 i r1 sa funkcjami liniowymi zmiennej x. Możemy wyznaczyć warunki na te funkcje, skoro W ma spełniać warunki brzegowe (16):

t2 = W (0, t) = r0(0)t2+ r1(0)t3, czyli r0(0) = 1 i r1(0) = 0, a z

t3 = W (π, t) = r0(π)t2+ r1(π)t3,

otrzymamy r0(π) = 0 i r1(π) = 1. Skoro wiec s a to funkcje liniowe zmiennej x, to z postaci r 0(x) = ax + b i z warunków na te funkcj e, otrzymujemy, że a = − π1, b = 1, a z postaci r1(x) = cx + d i z warunków na te funkcj e: c = π1, d = 0. Mamy wiec

W (x, t) =

−x π + 1

t2+ x πt3. Łatwo sprawdzić teraz, że V spełnia jednorodne warunki brzegowe.

Przekształcimy teraz równanie (14), używajac znalezionej postaci funkcji W . Ponieważ Wxx+ Vxx = Wtt + Vtt,

(10)

to dostajemy równanie niejednorodne:

Vxx = Vtt + 2

1 − x π

+ 6x

πt (17)

z jednorodnymi warunkami brzegowymi

V (0, t) = V (0, π) = 0 (18)

i z warunkami poczatkowymi

V (x, 0) = u(x, 0) − W (x, 0) = sin x, Vt(x, 0) = ut(x, 0) − Wt(x, 0) = 0. (19) Zatem rozwiazania równania (17) szukamy metod a rozdzielania zmiennych; V (x, t) = v (x)w (t) dla równania jednorodnego i otrzymujemy, jak poprzednio

v(x)

v (x) = λ = w(t) w (t),

z warunkami v (x) = v (π) = 0 dla równania z funkcja v : v (x) − λv (x) = 0. Rozwiazuj ac ten ostatni problem, dostajemy ciag wartości własnych λ k = −k2 i odpowiadajacy mu ci ag funkcji vk(x) = sin kx, k = 1, 2, ... . Zatem mamy rozwiazanie w postaci

V (x, t) =



k=1wk(t) sin kx.

Po wstawieniu do równania różniczkowego (17) dostaniemy

 k=1

(−k2)wk(t) sin kx =

 k=1

wk(t) sin kx + 2

1 − x π

+ 6x

πt. (20)

Ostatnie dwa składniki wystepuj ace po prawej stronie należy teraz rozwin ac w szereg wzgl edem sinusów, aby móc wykonać porównanie wyraz po wyrazie. Mamy

2

1 − x π

+ 6x πt =

 k=1

aksin kx, czyli

ak = 2 π

 π

0



2

1 − x π

+ 6x πt



sin kx dx.

Całkujac przez cz eści, mamy: ak = 2

π



2 −2x π + 6x

π t

−1

k cos kx 

π 0 + 1

k

 π

0

−2 π + 6

πt

cos kx dx



, czyli

ak = 2

kπ(−6t cos kπ + 2) . Wracajac do równania (20), otrzymujemy teraz:



k=1(−k2)wk(t) sin kx =



k=1wk(t) sin kx +

 k=1

2

kπ(−6t cos kπ + 2) sin kx.

(11)

Porównujac teraz wyrazy szeregu, dostajemy równanie różniczkowe zwyczajne niejednorodne wk(t) + k2wk(t) = 2

kπ(−6t cos kπ + 2) z warunkami wyznaczonymi z (19):

sin x = V (x, 0) =



k=1wk(0) sin kx, 0 = Vt(x, 0) =



k=1wk(0) sin kx.

Wynika stad, że w k(0) = 0, k = 1, 2, ... , w1(0) = 1 i wk(0) = 0 dla k = 1. Mamy wiec do rozwiazania dwa równania:

w1(t) + 12w1(t) = 2

1 · π (6t cos(1 · π) − 2) , w1(0) = 1, w1(0) = 0, (21)

wk(t) + k2wk(t) = 2

kπ(6t cos(kπ) − 2) , wk(0) = wk(0) = 0, dla k = 1, 2, ... . (22) Zauważmy, że równanie jednorodne w obu przypadkach jest takie samo:

wk(t) + k2wk(t) = 0, wiec jego rozwi azaniem jest

˜

wk(t) = A sin kt + B cos kt.

Stosujemy metode przewidywań (po prawej stronie, jednakowej dla obu obu równań niejednorodnych, wystepuje funkcja liniowa zmiennej t):

¯

wk(t) = Mt + N.

Wtedy mamy:

k2(Mt + N) = 12 cos kπ t − 4

, wiec N = 12 cos kπk3π i N = −k43π. Ostatecznie mamy

wk(t) = ˜wk(t) + ¯wk(t) = A sin kt + B cos kt + 12 cos kπ

k3π t − N = − 4 k3π.

Z (22) 0 = wk(0) = B − k43π, wiec B + k43π, a 0 = wk(0) = kA + 12 cos kπk3π , wiec A = − 12 cos kπk4π . Natomiast z (21) mamy 1 = w1(0) = B −π4, wiec B = 1 + π4 i 0 = w1(0) = A −12π implikuje A = 12π.

Mamy wiec teraz postać wszystkich w k i możemy zapisać kolejno rozwiazania: wk(t) = 4

k3π



−3

kcos kπ sin kt + cos kt + 3t cos kπ − 1



, k = 2, 3, ... ,

w1(t) = 1

π(12 sin t + (4 + π) cos t − 12t − 4),

(12)

V (x, t) = w1(t) sin x +k=2wk(t) sin kx =

= π1sin x(12 sin t + (4 + π) cos t − 12t − 4)+

+ k=2 k43πsin kx3k cos kπ sin kt + cos kt + 3t cos kπ − 1 i ostatecznie

u(x, t) = W (x, t) + V (x, t) =

=

x π + 1

t2+ xπt3+ 1πsin x(12 sin t + (4 + π) cos t − 12t − 4)+

+ k=2 k43πsin kx3k cos kπ sin kt + cos kt + 3t cos kπ − 1.

2 Zestaw zadań do ćwiczenia samodzielnego

1. Zastosować metode¸ Fouriera do wyznaczenia drgań podłużnych pre¸ta o długości jednostkowej przy naste¸puja¸cych warunkach pocza¸tkowych:

u(x, 0) = ψ(x), ut(x, 0) = 0, gdzie

ψ(x) =

0 dla x ∈ 0, a,

h(x−a)

b−a dla x ∈ (a, b), h > 0, 0 dla x ∈ b, 1.

Udowodnić, że otrzymany szereg:

(i) jest zbieżny na całej płaszczyźnie xt, przy czym zbieżność ta jest niemal jednostajna na każdym prostoka¸cie otwartym ograniczonym prostymi: x ± t = 2k − b, −x ± t = 2k − b k ∈ Z. (Wsk.

Można udowodnić, że jeśli (an) da¸ży monotonicznie do 0, ton=1ansin nx jest jednostajnie zbieżny w każdym przedziale domknie¸tym nie zawieraja¸cym punktów x = 2kπ, k ∈ Z);

(ii) jego suma spełnia warunki pocza¸tkowe w całym przedziale 0, 1 z wyjątkiem punktu x = b;

(iii) sprawdzić, czy otrzymany szereg jest dwukrotnie różniczkowalny wyraz po wyrazie w każdym z prostoka¸tów z i).

2. Wyznaczyć drgania podłużne pre¸ta o długości jednostkowej przy warunkach pocza¸tkowych u(x, 0) = 0, ut(x, 0) = ψ(x), gdzie ψ jest określona, jak w zadaniu poprzednim. Udowodnić, że szereg otrzy- many metoda¸ Fouriera:

(i) jest zbieżny w całej płaszczyźnie xt,

(ii) jego suma spełnia warunki pocza¸tkowe w całym przedziale 0, 1 z wyja¸tkiem punktu x = b,

(13)

(iii) stanowi rozwia¸zanie zagadnienia, czyli szereg jest zbieżny jednostajnie i u(x, t) jest klasy C2.

3. Rozwia¸zać naste¸puja¸ce równania:

(i) utt = uxx + bx(x − l) przy jednorodnych warunkach pocza¸tkowych i brzegowych u(0, t) = 0, u(l, t) = 0;

(ii) ∂t2u2 = ∂x2u2 + x(x − l)t2 przy jednorodnych warunkach pocza¸tkowych i brzegowych u(0, t) = 0, u(l, t) = 0;

(iii) ∂t2u2 = tα ∂∂x2u2, (α > −1) spełniaja¸ce warunki pocza¸tkowe u|t=0 = ϕ0(x), ∂u∂t |t=0 = ϕ1(x) oraz jednorodne warunki brzegowe u|x=0= 0, u|x=a = 0.

4. W półpasie a < x < b, t > 0 rozwia¸zać zagadnienie brzegowe uxx = utt, u(a, t) = u(b, t) = 0.

Czy jest to rozwia¸zanie jedyne?

5. W półpasie a < x < l, t > 0 rozwia¸zać zagadnienie mieszane dla równania utt = a2uxx przy warunkach:

(i) u(0, t) = u(l, t) = 0, u(x, 0) = 0, ut(x, 0) = sinl x;

(ii) u(0, t) = ux(l, t) = 0, u(x, 0) = x, ut(x, 0) = sin2lπx + sin2lx.

6. W półpasie a < x < l, t > 0 rozwia¸zać zagadnienie mieszane dla równania utt = a2uxx + f (x) przy warunkach:

(i) u(0, t) = α, u(l, t) = β, u(x, 0) = ϕ(x), ut(x, 0) = ψ(x),

(ii) ux(0, t) = α, ux(l, t) + hu(l, t) = β, u(x, 0) = ut(x, 0) = 0, h > 0.

7. Posługuja¸c sie¸ rozkładem u(x, t) = V (x, t) + W (x, t), rozwia¸zać:

(i) uxx = utt u(0, t) = µ(t), u(l, t) = ν(t), u(x, 0) = ϕ(x), ut(x, 0) = ψ(x);

(ii) uxx = utt+ f (x, t) u(0, t) = µ(t), ux(l, t) + hu(l, t) = ν(t), h > 0, u(x, 0) = 0, ut(x, 0) = ψ(x).

8. W półpasie a < x < l, t > 0 rozwia¸zać zagadnienie mieszane dla równania utt = a2uxx+ f (x, t) przy warunku pocza¸tkowym u(x, 0) = 0, ut(x, 0) = 0 i naste¸puja¸cym warunku brzegowym i funkcji f :

(i) u(0, t) = u(l, t) = 0, f (x, t) = Ae−tsinπlx, (ii) u(0, t) = u(l, t) = 0, f (x, t) = Axe−t, (iii) u(0, t) = ux(l, t) = 0, f (x, t) = A sin t, (iv) u(0, t) = u(l, t) = 0, f (x, t) = Ae−t cos2lπx.

(14)

9. Rozwia¸zać naste¸puja¸ce zagadnienia mieszane:

(i) uxx = utt z warunkami: u(0, t) = t2, u(π, t) = t3, u(x, 0) = sin x, ut(x, 0) = 0, gdzie 0 < x < π, t > 0,

(ii) uxx = utt z warunkami: u(0, t) = e−t, u(π, t) = t, u(x, 0) = sin x cos x, ut(x, 0) = 1, gdzie 0 < x < π, t > 0,

(iii) uxx = utt z warunkami: u(0, t) = t, ux(π, t) = 1, u(x, 0) = sin 12x, ut(x, 0) = 1, gdzie 0 < x < π, t > 0,

(iv) utt = a2uxx + sin 2t z warunkami: ux(0, t) = t2, ux(l, t) = 2asin2la sin 2t, u(x, 0) = 0, ut(x, 0) = −2 cos2xa, gdzie 0 < x < l, t > 0.

10. Napisać zagadnienie, do którego w wyniku rozdzielenia zmiennych u(x, y , t) = v (x, y )w (t) sprowadza sie¸ mieszane zagadnienie brzegowe:

uxx+ uyy − utt = 0,

u(x, y , t) = 0, t > 0, (x, y ) ∈ C ,

u(x, y , 0) = ϕ(x, y ), ut(x, y , 0) = ψ(x, y ), (x, y ) ∈ G ,

gdzie G jest obszarem płaszczyzny zmiennych x i y o brzegu C , a ϕ(x, y ) i ψ(x, y ) sa¸ danymi funk- cjami cia¸głymi.

11. Wykazać, że jeżeli u(x, t) jest rozwiązaniem równania uxx − utt = 0, to rozwiązaniem tego równania jest także funkcja

v (x, t) = u

x

x2− t2, t x2− t2

wszędzie tam, gdzie ona jest określona.

12. Wykazać, że jeżeli u(x, t) o ciągłych pochodnych cząstkowych trzeciego rzędu jest rozwiąza- niem równania uxx− utt = 0, to rozwiązaniem tego równania jest także funkcja

v (x, t) = u

x

x2− t2, t x2− t2

.

Bibliografia

[1] W. I. Arnold, Metody matematyczne mechaniki klasycznej, PWN, Warszawa 1981.

(15)

[2] W. I. Arnold, Równania różniczkowe zwyczajne, PWN, Warszawa 1975.

[3] W. I. Arnold, Teoria równań różniczkowych, PWN, Warszawa 1983.

[4] A. W. Bicadze, Równania fizyki matematycznej, PWN, Warszawa 1984.

[5] A. W. Bicadze, D. F. Kaliniczenko, Zbiór zadań z równań fizyki matematycznej, PWN, Warszawa 1984.

[6] P. Biler Prof. dr hab.- redakcja naukowa, Warsztaty z równań różniczkowych czastkowych, Toruń 2003. [7] Birkholc A. Analiza matematyczna. Funkcje wielu zmiennych, Wydawnictwo Naukowe PWN, Warszawa 2002.

[8] D. Bleecker, G. Csordas, Basic Partial Differential Equations, Chapman & Hall, Oxford 1995.

[9] L. Evans, Równania różniczkowe czastkowe, PWN, Warszawa 2002. [10] Fichtenholz G.M. Rachunek różniczkowy i całkowy, PWN, Warszawa 1980.

[11] J. Jost, Postmodern Analysis, Springer-Verlag,Berlin-Heidelberg-New York 2002.

[12] W. Kołodziej, Wybrane rozdziały analizy matematycznej, PWN, Warszawa 1982.

[13] H. Marcinkowska, Wstep do teorii równań różniczkowych cz astkowych, PWN, Warszawa 1972. [14] J. Musielak, Wstep do analizy funkcjonalnej, PWN, Warszawa 1976.

[15] Ockendon J., Howison S., Lacey A., Movxhan A., Applied Partial Differential Equations, Oxford University Press, 2003.

[16] J. Ombach, Wykłady z równań różniczkowych wspomagane komputerowo -Maple, Wydawnictwo Uniwersytetu Jagiellońskiego, Kraków 1999.

[17] B. Przeradzki, Równania różniczkowe czastkowe. Wybrane zagadnienia, Wydawnictwo Uniwersytetu Łódz- kiego, Łódź 2000.

[18] B. Przeradzki, Teoria i praktyka równań różniczkowych zwyczajnych, Wydawnictwo Uniwersytetu Łódzkiego, Łódź 2003.

[19] M. M. Smirnow, Zadania z równań różniczkowych czastkowych, PWN, Warszawa 1970.

[20] P. Strzelecki, Krótkie wprowadzenie do równań różniczkowych czastkowych, Wydawnictwo Uniwersytetu War- szawskiego, Warszawa 2006.

[21] B. W. Szabat, Wstęp do analizy zespolonej, PWN, Warszawa 1974.

[22] Whitham G.B., Lecture notes on wave propagation , Springer-Verlag, Berlin-Heidelberg-New York 1979.

[23] Zauderer, Partial Differential Equations of Applied Mathemathics, John Wiley & Sons, Singapore-New York- Chichester-Brisbane-Toronto 1989.

Cytaty

Powiązane dokumenty

Nawet jeżeli dla pewnej funkcji f rozwiązanie istnieje to nie zależy w sposób ciągły od parametrów zadania (czyli funkcji f ).. 4.4

Zmniejszenie kroku h istotnie polepsza dokładność metody łamanych, przy czym należy pamiętać, że nadmierne zmniejszenie kroku daje efekt odwrotny do spodziewanego.

Jak widać, wyniki otrzymane metodą Eulera i metodą Rungego-Kutty są do siebie podobne, aby jednak podobieństwo to stało się wyraźniejsze (i aby dokończyć rozwiązywania

Jeśli maksymalny rząd pochodnych funkcji u jest 2 (czyli w równaniu pojawia się przynajmniej jedna pochodna cząstkowa drugiego rzędu i nie ma pochodnych wyższego rzędu), to

Strzelecki, Krótkie wprowadzenie do równań różniczkowych cz astkowych, Wydawnictwo Uniwersytetu War-  szawskiego, Warszawa 2006..

[r]

• Na ocenę z przedmiotu składa się wynik kolokwium (warte 40 punktów), projekt labo- ratoryjny (warty 40 punktów) oraz aktywność na ćwiczeniach (10 punktów) i

Należy w każdym zadaniu wykonać jedynie podpunkt zgodny z numerem na liście obecności na zajęciach (osoby nieobecne proszone są o kontakt mailowy w celu ustalenia numeru)..