Matlab中如何求解隐式常微分方程

你可真搞笑

于 2022-04-01 22:28:13 发布

阅读量3.8k

点赞数

分类专栏： Matlab 文章标签： matlab 线性代数

本文链接：https://blog.csdn.net/weixin_43568289/article/details/123885727

版权

Matlab 专栏收录该内容

55 篇文章 105 订阅 ¥39.90 ¥99.00

订阅专栏

超级会员免费看

本文介绍了在Matlab中如何使用ode45和ode15i求解隐式常微分方程。尽管ode45可以通过与其他函数配合使用来解决隐式方程，但其对初始值敏感，容易导致计算错误。相比之下，ode15i是专为隐式方程设计，能提供更准确的结果。建议在求解隐式式子时选择ode15i。

摘要由CSDN通过智能技术生成

在matlab中求解隐式常微分方程，通常有两种方法，一种是使用matlab中的内嵌的ode15i等求解隐式微分方程的函数，另一种是仍使用matlab中求解常微分方程的显式函数，如ode45等等。但是小编想要表达的绝不是这两种方法，而是这两种方法之间的差异：

一、ode45

使用这种方法求解隐式微分方程的思路就是结合其他函数求解隐式方程，例如fzero等等。

以求解如下隐式方程为例：
$\dot{y}^2-5/(0.5*\dot{y}-2/y)=0$

可以发现这个式子很复杂，没有显式式子，因此只能使用隐式函数求解。代码如下所示：

% Matlab

%% ode45
tspan = [0 1];                                                             % 求解时间域，此处假定x=t，为时间变量
y0 = 1.924;                                                                % 初始值

[t,y] = ode45(@(t,y)odefcn(t,y),tspan,y0);