matlab里的ode15s函数怎么使用

`ode15s` 是 Matlab 中的一个常用函数，用于求解刚体的常微分方程组。以下是 `ode15s` 函数的基本使用方法：

1. 定义微分方程：将微分方程写成一组形如 dy/dx = f(x,y) 的一阶微分方程组形式。

2. 定义初值：确定初值 y0。

3. 定义求解区间：确定求解区间 tspan，例如 [t0, tf]。

4. 调用 `ode15s` 函数：使用以下语法调用函数：

   [t,y] = ode15s(odefun,tspan,y0,options)

其中，odefun 是定义微分方程的函数句柄；t 和 y 分别是求解区间 tspan 上的时间向量和解向量；options 是选项结构体，可用来设定求解精度、终止条件等。

下面是一个简单的例子，用 `ode15s` 函数求解常微分方程：

% 定义微分方程
odefun = @(t,y) [-2*y(1)+y(2); 3*y(1)-4*y(2)];

% 定义初值和求解区间
y0 = [1;0];
tspan = [0,2];

% 调用 ode15s 函数
[t,y] = ode15s(odefun,tspan,y0);

% 绘制解的图像
plot(t,y(:,1),'r-',t,y(:,2),'b--')
legend('y_1','y_2')
xlabel('t')

这个例子中，我们定义了一个二阶微分方程组，使用 `ode15s` 函数求解后，得到了在求解区间 [0,2] 上的解向量 y，然后绘制了 y1 和 y2 随时间的变化曲线。