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/05/05 08:16:31
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','
没有解析解,请问数值解怎 么求

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')

 

%%%%%%%%%%%%%%%%%%

收起