Question

ode45求解微分方程组时，出现如下问题

%建立dyfun.m文件
function dy=dyfun(t,y)
dy=zeros(6,1);
dy(1)=y(2);
dy(2)=0.5*(-3*y(2)+1.5*y(4)-100*y(1)+50*y(3)=2*sin(3.754*t));
dy(3)=y(4);
dy(4)=0.5*(1.5*y(2)-3*y(4)+1.5*y(6)+50*y(1)-100*y(3)+50*y(5)-2*cos(2.2*t));
dy(5)=y(6);
dy(6)=0.5*(-1.5*y(4)-3*y(6)+50*y(3)-100*y(5)+1*sin(2.8*t));
end

%在windows窗口运行
y0=[1,1,1,1,1,1];
tspan=[0 30];
[t,y]=ode45(@dyfun,tspan,y0);
plot(t,y(:,1))

错误使用 feval
未定义与 'double' 类型的输入参数相对应的函数 'dyfun'。

出错 odearguments (line 87)
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

出错 ode45 (line 113)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

Answer 1

这个程序的问题在于：dy(2)=0.5*(-3*y(2)+1.5*y(4)-100*y(1)+50*y(3)=2*sin(3.754*t));其中有两个＝号，把后面的那个改为减号就行，可改为：
dy(2)=0.5*(-3*y(2)+1.5*y(4)-100*y(1)+50*y(3)-2*sin(3.754*t));

追问

这里是我写错了，后来我自己找到问题了，谢谢！