Здравствуйте, vera-nika!
Предлагаю решение 4 задачи:
u
t=4u
xxu(x,0)=x+3
u(0,t)=u(1,t)=0
Решение.
Будем искать (не равное нулю) решение уравнения в виде произведения двух функций, одна из которых зависит только от х, а другая – только от t, т.е.
u(x,t)=X(x)T(t).
Подставляя это выражение в уравнение
u
t=a
2u
xxимеем
X(x)T'(t)=a
2X''(x)T(t).
Здесь a
2=4. После деления на X(x)T(t) получим:
1/a
2*T'(t)/T(t)=X''(x)/X(x).
Это равенство двух отношений, зависящих только от х и только от t, возможно только в случае, если оба отношения равны постоянному числу – [$955$]. ([$955$]>0):
1/a
2*T'(t)/T(t)=X''(x)/X(x)= - [$955$]
т.е.
X''(x)+[$955$]X(x)=0
T'(t)+[$955$]*a
2*T(t)=0
Первое уравнение системы с граничными условиями
X''(x)+[$955$]X(x)=0
X(0)=X(1)=0
представляет собой задачу Штурма-Лиувиля на отыскание собственных функций и собственных значений дифференциального оператора. Из общего решения уравнения
X(x)=A*cos([$8730$]([$955$]*x))+B*sin([$8730$]([$955$]*x))
использованием краевых условий
X(0)=A*cos(0)+B*sin(0)=A=0
X(1)=A*cos([$8730$]([$955$]))+B*sin([$8730$]([$955$]))=B*sin([$8730$]([$955$]))=0
находим собственные значения
[$955$]=[$955$]
n=(п*n)
2и собственные функции
X
n(x)=B
n*sin(n*п*x).
Второе уравнение системы
T'(t)+[$955$]*a
2*T(t)=0
имеем решение
T(t)=C*e
-[$955$]*a^2*t.
Таким образом, общее решение исходного уравнения принимает вид
u
n(x,t)=X(x)T(t)=C
n*e
-(п*n*a)^2*t*sin(п*n*x).
Поскольку при любых n полученная функция является решением нашего дифференциального уравнения, то и сумма этих решений так же будет являться решением исходного дифференциального уравнения:
u(x,t)=[$8721$]C
n*e
-(п*n*a)^2*t*sin(п*n*x).
Подставляя сюда начальное условие
u(x,0)=[$8721$]C
n*sin(п*n*x)=x+3
Последняя формула показывает, что величины C
n являются коэффициентами разложения функции x+3 в ряд Фурье по синусам в интервале (0,1):
C
n=2/l*
01[$8747$]f(x)sin(пnx/l)dx,
C
n=2*
01[$8747$](x+3)sin(пnx)dx
Интегрируя два раза по частям
C
n=2/(nп)*(4*( -1)
n+1+3)
Получим окончательный ответ
u(x,t)=[$8721$]2/(nп)*(4*( -1)
n+1+3)*e
-4(п*n)^2*t*sin(п*n*x).
Текс ответа в прикрепленном файле.
Будут вопросы обращайтесь в мини-форум.
Удачи