Question

MATLAB求解微分方程组 调用 ode45解法的编程问题

各位大神，你们好：我是MATLAB 的初学者，想编个小程序，求解一阶微分方程组，并输出图像。编程如下：创建了M-文件：function dy=rigid(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2);[T,Y]=ode45(@rigid,[0,12],[0,1,1]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')但是在command windows窗口中输入rigid时显示如下：??? Input argument "y" is undefined.Error in ==> rigid at 3dy(1)=y(2)*y(3);请问是什么原因？请大神们多多帮忙，谢谢

Answer 1

，你这里为什么会有两个输入？而且t根本就没用。你这样定义引用的时候肯定会出错了，你改成function
dy=rigid(y)试试。