使用变量的ode23 ode求解器

Question

我想对系统使用ode23求解器。我的系统有三个方程，第一个方程依赖于上面代码中声明的变量'z‘。下面是示例代码：

clc; clear;
z=1;

function Fv=funsys(t,Y,z); 
Fv(1,1)=2*Y(1)+Y(2)+5*Y(3)+exp(-2*t) + z; 
Fv(2,1)=-3*Y(1)-2*Y(2)-8*Y(3)+2*exp(-2*t)-cos(3*t);
Fv(3,1)=3*Y(1)+3*Y(2)+2*Y(3)+cos(3*t); 
end

[tv,Yv]=ode23('funsys',[0 pi/2],[1;-1;0]); 

变量'z‘未声明的错误。那么，有可能解决这个问题吗？(函数'funsys‘需要依赖于'z')

Answer 1

是的，看看文档中的将额外的参数传递给ODE函数：

您可以通过在函数外部定义参数并在指定函数句柄时传入它们，从而传入额外的参数。

在你的例子中，它应该是这样的：

[tv,Yv]=ode23(@(t,y) funsys(t,y,z),[0 pi/2],[1;-1;0]);