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','
没有解析解，请问数值解怎 么求
或ODE23
即龙格库塔法
话说没有分么？追问
你能解出来才行啊
追答试着解了一下：
>> dy=inline('[y(2);2783.83*(2*0.059*cos(105.52)-1)*y(1)]','t','y');
>> ode45(dy,1,[1 0])
y-t的图像为：
多谢，还可以再问你个问题吗，矩阵形式怎么解：
clear all
n=2;
EI=4.2063e6;
L=7;
M=[m11,m12;m21,m22];
K=[k11,k12;k21,k22];
S=[s11,s12;s21,s22];
DL=inline('[x(n+1:end,1); inv(M)*(-K*x(1:n,1)-(98264*cos(79.865*t))*S*x(1:n,1))]',...
't','x','flag','n','M','K','S');
请看私信。。。，多谢了
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')
%%%%%%%%%%%%%%%%%%
不是这个曲线，function dy=chbq(t,y) 方程右边不对，不能含有y1 y2 之外的变量
追答不懂
追问不是这个曲线，function dy=chbq(t,y) 方程右边不对，不能含有y1 y2 之外的变量
我这是用simulink算出来的,但是怎么编程解不会
什么变量啊，不能说清楚吗，搞不懂你
'y(0)=1' 你的图上是=1 吗？
多谢~
