matlab中ode23的基本使用方法(秒懂)

  • Post author:
  • Post category:其他



题目如下图


解题过程以及注释如下

function [y1,y2]= myfun3
tspan = 0:0.1:20;%计算结果的数值区间
y0 =[20;12];%y1和y2的初值,即y1(0)和y2(0)
[t,r]= ode23(@dy, tspan,y0);%调用ode23函数
%注意,这里使用的调用函数的方法是@函数名,而非'函数名',因为该函数并非.m文件,不能用引号调用。
y1 = r(:,1);%结果赋值
y2 = r(:,2);%结果赋值

function z =dy(t,y)%因为是方程组,必须用矩阵哦
z = [(1 - 0.013*y(2))*y(1);%存放y1't
    (-1 + 0.011*y(1))*y(2)];%存放y2't



版权声明:本文为weixin_52739840原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。