function dy=vdp1000(t,y) % 建立函数，函数名字是vdp，t，y表示输入值，dy是输出值
dy=zeros(2,1); %建立空矩阵 2*1的0矩阵
dy(1)=y(2); %第一个微分方程
dy(2)=1000*(1-y(1)^2)*y(2)-y(1);%第二个微分方程
建立m文件的作用就是在[T,Y]=ode15s('vdp1000',[0 3000],[2 0]); 调用这个vdp.m函数。
如不明白请留言...
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024