Lutik
Сообщение
#46655 4.12.2009, 19:43
Сделал систему уравнений, метод исключения:
y'=2y-4z
z'=y-3z+e^x
далее я выразил из 2-го уравнения у:
-у=-z'-3z+e^x
y=z'+3z-e^x
тогда нашёл производную и подставил её в 1-ое уравнение:
y'=z''+3z'-e^x
z''+3z'-e^x=2y-4z
z''+3z'+4z=2y-e^x
нашёл характеристическое уравнение,
если
z=e^kx
z'=k*(e^(kx))
z''=k^2*(e^(kx))
k^2*(e^(kx))+3*k*(e^(kx))+4*e^(kx)=0
e^(kx)(k^2+3k+4)=0
k^2+3k+4=0
k1=(-3+корень7i)/2
k2=(-3-корень7i)/2
комплексные корни, значит z1=e^(-3x)cos(7x), z2=e^(-3x)sin(7x)
общее решение z=С1*e^(-3x)cos(7x)+С2*e^(-3x)sin(7x)
Значит потом нужно найти производную z' и подставить в z'=y-3z+e^x и из него найти у?
граф Монте-Кристо
Сообщение
#46660 4.12.2009, 19:55
Нет.
Вам нужно в Ваше уравнение:
z''+3z'+4z=2y-e^x
z''+3z'+4z=2y+e^x
подставить y,который Вы выражали раньше:
y=z'+3z-e^x,
чтобы справа была функция,которая зависит только от x.
Lutik
Сообщение
#46661 4.12.2009, 20:09
y=z'+3z-e^x
y=-3*e^(-3x)*c1*cos(7x)-7*c1*e^(-3x)*sin(7x)-3*c2*e^(-3x)*sin(7x)+7*c2*e^(-3x)*cos(7x)+3*c1*e^(-3x)cos(7x)+3*c2*e^(-3x)sin(7x)-e^x
y=-7*c1*e^(-3x)*sin(7x)+7*c2*e^(-3x)*cos(7x)+3*c1*e^(-3x)cos(7x)+-e^x
Тогда z и y известны значит получается ответ
y=-7*c1*e^(-3x)*sin(7x)+7*c2*e^(-3x)*cos(7x)+3*c1*e^(-3x)cos(7x)+-e^x
z=С1*e^(-3x)cos(7x)+С2*e^(-3x)sin(7x)
граф Монте-Кристо
Сообщение
#46662 4.12.2009, 20:21
Вы не поняли. Перед тем,как решать, нужно подставить в уравнение
z''+3z'+4z=2y+e^x
функцию y, выраженную через z,z' и x, т.е.
z''+3z'+4z=2(z'+3z-e^x)+e^x.
Dimka
Сообщение
#46663 4.12.2009, 20:24

Вы не исправимый.
Lutik
Сообщение
#46664 4.12.2009, 20:29
z''+3z'+4z=2(z'+3z-e^x)+e^x
z''+3z'+4z-2z'-3z=0
z''+z'+z=0
z=e^kx
z'=k*e^kx
z''=(k^2)*e^(kx)
(k^2)*e^(kx)+k*e^kx+e^kx=0
e^(kx)*(k^2+k+1)=0
k1=(-1+корень15i)/2
k2=(-1-корень15i)/2
комплексные числа:
z1=e^(-x)cos(15x)
z2=e^(-x)sin(15x)
z=e^(-x)cos(15x)+e^(-x)sin(15x)
Dimka
Сообщение
#46665 4.12.2009, 20:33
ошибка в решении квадратного уравнения
Lutik
Сообщение
#46666 4.12.2009, 20:36
Спасибо, вместо 4*1 умножил 4*4, глючит меня

z''+3z'+4z=2(z'+3z-e^x)+e^x
z''+3z'+4z-2z'-3z=0
z''+z'+z=0
z=e^kx
z'=k*e^kx
z''=(k^2)*e^(kx)
(k^2)*e^(kx)+k*e^kx+e^kx=0
e^(kx)*(k^2+k+1)=0
k1=(-1+корень
5i)/2
k2=(-1-корень
5i)/2
комплексные числа:
z1=e^(-x)cos
(5x)z2=e^(-x)sin
(5x)z=e^(-x)cos(5x)+e^(-x)sin
(5x)
Dimka
Сообщение
#46667 4.12.2009, 20:37
не угадали.
Lutik
Сообщение
#46669 4.12.2009, 20:47
Извините крыша едет
z''+3z'+4z=2(z'+3z-e^x)+e^x
z''+3z'+4z-2z'-3z=0
z''+z'+z=0
z=e^kx
z'=k*e^kx
z''=(k^2)*e^(kx)
(k^2)*e^(kx)+k*e^kx+e^kx=0
e^(kx)*(k^2+k+1)=0
k1=(-1+корень3i)/2
k2=(-1-корень3i)/2
комплексные числа:
z1=e^(-x)cos(3x)
z2=e^(-x)sin(3x)
z=e^(-x)cos(3x)+e^(-x)sin(3x)
Dimka
Сообщение
#46670 4.12.2009, 20:50
z=e^(-x)(С1cos(3x)+С2sin(3x))
дальше ищите y
Проверку не забудьте сделать.
Lutik
Сообщение
#46672 4.12.2009, 21:01
В это уравнение подставляю производную z': y=z'+3z-e^x
z'=-e^(-x)*С1*cos(3x)-3*e^(-x)*С1sin(3x)-e^(-x)*С2*sin(3x)+3*e^(-x)*С1*cos(3x)
y=-e^(-x)*С1*cos(3x)-3*e^(-x)*С1*sin(3x)-e^(-x)*С2*sin(3x)+3*e^(-x)*С2*cos(3x)+3*e^(-x)*С1*cos(3x)+3*e^(-x)*С2*sin(3x)-e^x
y=2*e^(-x)*С1*cos(3x)-3*e^(-x)*С1*sin(3x)+2*e^(-x)*С2*sin(3x)+3*e^(-x)*С2*cos(3x)-e^x
вроде так
Dimka
Сообщение
#46674 4.12.2009, 21:19
Цитата(Lutik @ 4.12.2009, 23:47)

Извините крыша едет
z''+3z'+4z=2(z'+3z-e^x)+e^x
z''+3z'+4z-2z'-3z=0
z''+z'+z=0
преобразование с ошибкой сделаны. Соответственно дальнейшее решение неверное. Переделывайте заново.
Lutik
Сообщение
#46676 4.12.2009, 21:28
Это моя проблема - невнимательность
z''+3z'+4z=2(z'+3z-e^x)+e^x
z''+3z'+4z-2z'-6z=-e^x
z''+z'-2z=-e^x
z=e^kx
z'=k*e^kx
z''=(k^2)*e^(kx)
(k^2)*e^(kx)+k*e^(kx)-2*e^(kx)=0
e^(kx)*(k^2+k-2)=0
k^2+k-2=0
D=9
k1=(-1+3)/2=1
k2=(-1-3)/2=-2
z1=e^0x=1
z2=e^-2x
z=c1+c2*e^-2x
граф Монте-Кристо
Сообщение
#46678 4.12.2009, 22:05
Цитата(Lutik @ 5.12.2009, 0:28)

k1=(-1+3)/2=1
k2=(-1-3)/2=-2
z1=e^0x=1
z2=e^-2x
Сколько же можно в таких мелочах косячить?
Lutik
Сообщение
#46695 5.12.2009, 10:33
Извините, исправил
k1=(-1+3)/2=1
k2=(-1-3)/2=-2
z1=e^x
z2=e^-2х
z=c1*e^x+c2*e^-2x
Далее подставляем в y=z'+3z-e^x производную z' и получаем y, что требовалось нашли, т.е z и y
граф Монте-Кристо
Сообщение
#46701 5.12.2009, 10:51
Нет, это Вы только общее решение нашли, теперь нужно ещё найти частное.
Lutik
Сообщение
#46709 5.12.2009, 11:46
Да, забыл, что там ещё есть -e^x
z''+z'-2z=-e^x
частное решение имеет вид z*=(Ax^2+Bx)*e^x, тогда нахожу
z*'=(e^x)*A*(x^2)+2*(e^x)*Ax+Bx*e^x+B*e^x
z*''=(e^x)*A*(x^2)+4*(e^x)*A*x+2*(e^x)*A+(e^x)*B*x+2*(e^x)*B
(e^x)*A*(x^2)+4*(e^x)*A*x+2*(e^x)*A+(e^x)*B*x+2*(e^x)*B+(e^x)*A*(x^2)+2*(e^x)*Ax
+Bx*e^x+B*e^x-2*(A*(x^2)*e^x+B*x*e^x)=-e^x
6*(e^x)*A*x+2*(e^x)*A+3*(e^x)*B=-e^x
6A=-1
A=-1/6
тогда В=-2/9
наверно не так, потому что там если находим А, то 6А=-1, а при нахождении В наверно ошибка
Вроде нашёл ошибку:
6A=0 => A=0
3*B=-1
b=-1/3
z*=(-1/3)*x*e^x
граф Монте-Кристо
Сообщение
#46710 5.12.2009, 11:48
Цитата(Lutik @ 5.12.2009, 14:46)

частное решение имеет вид z*=(Ax^2+Bx)*e^x, тогда нахожу...
Почему так?
Lutik
Сообщение
#46711 5.12.2009, 12:01
Если сравнивать правую часть -e^x с f(x)=P(x)e^Lx, L совпадает с корнем k1=1, то тогда (Ax^2+Bx)*e^x,
Lutik
Сообщение
#46712 5.12.2009, 12:01
Если сравнивать правую часть -e^x с f(x)=P(x)e^Lx, L совпадает с корнем k1=1, то тогда (Ax^2+Bx)*e^x,
граф Монте-Кристо
Сообщение
#46714 5.12.2009, 12:04
k=1 - корень кратности один,значит, если справа стоит P(x)*e^(kx), общее решение следует искать в виде z=x*Q(x)*e^(kx).
Lutik
Сообщение
#46715 5.12.2009, 12:12
Q(x) тогда чему будет равен не пойму в уравнении же только -e^x
граф Монте-Кристо
Сообщение
#46717 5.12.2009, 12:19
P(x)=-1 - многочлен нулевой степени, значит и Q(x) - тоже нулевой,т.е. Q(x)=C.
Lutik
Сообщение
#46720 5.12.2009, 12:30
z=x*A*e^x
z'=A*x*e^x+A*e^x
z''=2*A*e^x+x*A*e^x
z''+z'-2z=-e^x
2*A*e^x+x*A*e^x+A*x*e^x+A*e^x-2*x*A*e^x=-e^x
3*A*e^x=-e^x
A=-1/3
z=(-1/3)*x*e^x - частное решение
граф Монте-Кристо
Сообщение
#46721 5.12.2009, 12:40
Да.
Lutik
Сообщение
#46724 5.12.2009, 12:47
потом подставляем в y=z'+3z-e^x , производную и всё система решена
z'=(-1/3)*x*(e^x)-(1/3)*(e^x)
y=(-1/3)*x*(e^x)-(1/3)*(e^x)+3*(-1/3)*x*(e^x)-(e^x)
y=(-4/3)*x*(e^x)-(4/3)*(e^x)
вроде всё.
Спасибо за помощь!!!
Это текстовая версия — только основной контент. Для просмотра полной версии этой страницы, пожалуйста,
нажмите сюда.