MATLAB ode函数使用过程中微分方程创建方法

MATLAB ode函数使用过程中微分方程创建方法

例题一范德波尔方程

以上给出了方程以及方程的初试条件

function dy = vdp1000(t,y)  
dy = zeros(2,1);
dy(1)=y(2);     % 初试设定是 dy(1)=y(1)';dy(2)=y(2)';     而此句规定y(1)' = y(2)  从而有  dy(2) = y(2)'=y(1)''
dy(2)=1000*(1-y(1)^2)*y(2)-y(1); %注意到方恒的左边都是关于dy(n)的表达式，右边都是关于y(n)的表达式
end

方程设置完毕，求解见博客：https://blog.csdn.net/weixin_41855010/article/details/104194829

例题二 多个变量

求解过程见https://wenku.baidu.com/view/77d72d4733687e21af45a99c.html

本文主要想引出一个观点（不知道是不是对）：

y是解向量，y(1),y(2)…表示的是解向量中的变量，而系统默认dy(1)=y(1)’;dy(2)=y(2)’；以这个思路去建立微分方程。