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)’;以这个思路去建立微分方程。