函数 > 求解和最优化 > 微分方程求解器 > 示例:Stiffb 和 Stiffr 的增广雅可比矩阵
示例:Stiffb 和 Stiffr 的增广雅可比矩阵
使用 Jacob 函数计算常微分方程 (ODE) 的增广雅可比矩阵,然后将它输入求解器 StiffbStiffr
1. 定义一个具有 4 个未知数的组:
单击可复制此表达式
单击可复制此表达式
* 
yi 变量为 x 的函数。
2. 定义该组中的常数。
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
3. 定义初始值。
单击可复制此表达式
4. 定义与该组右侧对应的矢量函数 D(x,y)
单击可复制此表达式
5. 使用 augment 函数创建增广雅可比矩阵:
单击可复制此表达式
6. 调用 StiffbStiffr 函数:
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
* 
返回的矩阵包含 5 列,分别与点数和 4 个未知数的解对应。
7. 从返回的 StiffbStiffr 矩阵提取 4 个未知数的解:
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
8. 绘制和比较从这两个函数返回的每个未知数的解:
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
单击可复制此表达式
绘图显示两个 ODE 求解器返回相同的解。
这对您有帮助吗?