Question

matlab微分方程组的 数值解怎么求啊= =

具体 如下 [x y]=dsolve('D2x=-50*1.225*exp((1000-y)/10.7)*(Dx)^2',... 'D2y=9.8-100*1.225*exp((1000-y)/10.7)*(Dy)^2','y(0)=10000,Dx(0)=222,Dy(0)=0','t')解析解 求不出来。。不知道数值解怎么求呃= = 急啊

Answer 1

>> syms t

exp(-2*t)*(sin(2*t+sym(pi)/3)+cos(3*t))
ans =
exp(-2*t)*(sin(2*t+1/3*pi)+cos(3*t))
>> y=dsolve(['D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=',...
'exp(-2*t)*(sin(2*t+1/3*pi)+cos(3*t))'],'y(0)=1','y(1)=3','y(pi)=2',...
'Dy(0)=1','Dy(1)=2')

追问

你这是什么东西。。对不上号= =