Question

Matlab 解微分方程

这是我的程序，药解这个微分方程，运行总出错Not enough input arguments.求大神帮忙改正！谢谢！
function ydot=vdpol(t,y)
% Van der Pol微分方程
mu=2;
ydot=[y(2);mu*(1-y(1)^2)*y(2)-y(1)];
%solve_vdpol.m
%用常微分方程方法求解Van del Pol
%积分的时间区域
tspan=[0 20];
y0=[2;0];
disp('ode45 Van der pol')
[t,y]=ode45(@vdpol,tspan,y0);
%绘制结果
plot(t,y(:,1),t,y(:,2),'--');

Answer 1

function ydot=vdpol(t,y)
% Van der Pol微分方程
mu=2;
ydot=[y(2);mu*(1-y(1)^2)*y(2)-y(1)];

这段要单独到保存到 vdpol.m文件中

追问

我已经单独保存到m文件，可是我运行以下程序时需要调用此前定义的函数
%function main vdpol
%solve_vdpol.m
。。。。。
依然出错，请问大神怎么调用保存的函数？

追答

function ydot=vdpol(t,y)
% Van der Pol微分方程
mu=2;
ydot=[y(2);mu*(1-y(1)^2)*y(2)-y(1)];

这几句独立保存到vdpol.m文件中，保存在matlab当前工作的文件夹

tspan=[0 20];
y0=[2;0];
disp('ode45 Van der pol')
[t,y]=ode45(@vdpol,tspan,y0);
%绘制结果
plot(t,y(:,1),t,y(:,2),'--');
这几句内容写在另一个m文件中，或者直接粘贴在命令行运行