3.3
, < < , t > 0, (3.38)
.
(3.38) , t.
:
, < < , (3.39)
.
, , t 0, (3.38) (3.39).
(3.38) , , , t, .
,
. (3.40)
ϳ (3.38), :
.
ϳ (3.38), :
,
.
˳ t, , , , , (). :
. (3.41)
T(t) ().
, (3.42)
. (3.43)
(3.42) :
, (3.44)
1 .
(4.43) :
, (3.45)
2, 3 , , , .
ϳ (3.44) (3.45) (3.22), (3.20)
, (3.46)
.
:
, .
. (3.47)
, , , , (3.20) (3.47), :
. (3.48)
(3.48) (3.20).
, , (3.48) (3.21). (3.21), :
,
(3.48), :
. (3.49)
, :
. (3.50)
(3.49) (3.50),
,
(3.51)
.
ϳ (3.51) (3.48), :
,
. (3.52)
(3.52), (3.20), (3.21).
|
|
(3.52) , . , , - , (3.52) :
. (3.53)
(3.53), .
) -
. (3.54)
:
. (3.54/)
, , , . (3.54/) :
. (3.55)
(3.55) ,
,
.
0
.
.
,
.
, :
. (3.56)
(3.56)
.
ϳ (3.56), :
.
. (3.57)
) :
,
. (3.58)
) :
. (3.59)
(3.59) , :
. (3.59/)
(3.59/) , :
.
, , :
.
,
.
:
.
, :
,
. (3.60)
. , . ,
.
,
.
ϳ (3.60), :
. (3.61)
(3.61) (3.53), :
.
,
,
,
.
, (3.20), (3.21), :
,
. (3.62)
(3.62) - , .
. , . , .
3.4
3.1 .
, < < , t>0 (3.63)
. (3.64)
(3.62), (3.63), (3.64) :
.
:
.
(4.64), :
,
. (3.65)
.
ϳ (3.64),
,
, (3.66)
.
(3.66). :
|
|
.
.
:
= .
,
=
= .
,
,
, (3.67)
.
, , :
, (3.68)
.
, (3.69)
z , , , . . >5, >5) = 0,5. - , . , .
(3.67) :
. (3.70)
, = 0 , .
3.2
, < < , t>0 (3.71)
. (3.72)
, , , (3.62):
,
, (3.73)
- , .
3.3 , , , , .
.
> Eq:=diff(u(x,t),t$2)=a^2*diff(u(x,t),x$2);
.
> pdsolve(Eq,HINT=X(x)*T(t));
, , where , () T(t) .
(), .
> Eq1:=diff(X(x),`$`(x,2))=-lambda^2*X(x);
(), , : U(0,t)=X(0)T(t)=0, X(0)=0. .
> dsolve({Eq1,X(0)=0},X(x));
, X(l)=0. ( ), , , .
> _EnvAllSolutions:=true;
> solve(sin(lambda*l)=0,lambda);
" " .
, .
> nu:=n->Pi*n/l;
, .
> X:=(x,n)->sin(x*nu(n));
T(t), .
> Eq2:=diff(T(t),`$`(t,2))=-a^2*lambda^2*T(t);
> dsolve({Eq2,D(T)(0)=0},T(t));
> T:=(t,n)->cos(nu(n)*a*t);
.
> U:=(x,t)->Sum(A[n]*X(x,n)*T(t,n),n=1..infinity);
:
> U:=(x,t)->Sum((C[n]*cos(a*Pi*n*t/l)+ +D[n]*sin(a*Pi*n*t/l))*sin(Pi*n*x/l),n=1..infinity);
> a:=1;l:=Pi;
> U(x,t);