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/05 14:52:45
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','没有解析解,请问数值解怎 么求
xTNA~Shkٝ=֠NwMcvЮ.zE4*1$<DDcWpmbQz1q/v증W/6޾ɇtK텪@ 00@ݵK@U#IK Jf15㉊+ٖ V):'x9Kӵ`?QkF'D껓ƽ4cְCcL_Z2/|&l>^9e>\hˍʁs!ͺZδICɘeͳrNJ(ZƼd%F"pT{~,>oov;eݳ^XƌvX֦9ThR80L7&WirO4Gkݜlcq}"R;b[7x!--#a?'s0NNJx,d[TIVa|ukOv@TikJ*}֭Qr}wLS l)ZNl!3=-ZjwQ.+EdX*Ē C$&C!L$QflJI'颤jh"h@C=a&(F

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

 

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

收起