MATLAB解微分方程 :a=105.52 b=0.059 dsolve('D2y+2783.83*(1-2*b*cos(a*t))*y=0','y(0)=1','Dy(0)=0','没有解析解,请问数值解怎 么求
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/26 06:39:06
MATLAB解微分方程 :a=105.52 b=0.059 dsolve('D2y+2783.83*(1-2*b*cos(a*t))*y=0','y(0)=1','Dy(0)=0','没有解析解,请问数值解怎 么求
MATLAB解微分方程 :a=105.52 b=0.059 dsolve('D2y+2783.83*(1-2*b*cos(a*t))*y=0','y(0)=1','Dy(0)=0','
没有解析解,请问数值解怎 么求
MATLAB解微分方程 :a=105.52 b=0.059 dsolve('D2y+2783.83*(1-2*b*cos(a*t))*y=0','y(0)=1','Dy(0)=0','没有解析解,请问数值解怎 么求
ODE45
或ODE23
即龙格库塔法
话说没有分么?
function dy=chbq(t,y) % 在当前工作文件夹下保存为 chbq.m 文件 a=105.52;b=0.059; dy=...
全部展开
function dy=chbq(t,y) % 在当前工作文件夹下保存为 chbq.m 文件 a=105.52;b=0.059; dy=[y(1);-2783.83*(1-2*b*cos(a*t))*y(1)]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%% 运行以下: [t,y]=ode45('chbq',[0 2],[1;0]); plot(t,y(:,1)) title('y-t') figure plot(t,y(:,2)) title('dy-t') %%%%%%%%%%%%%%%%%%
收起