特殊微分方程的数值解——隐式微分方程（ode15i）

最新推荐文章于 2023-08-28 22:34:14 发布
最新推荐文章于 2023-08-28 22:34:14 发布
阅读量5.5k 收藏 33
点赞数 2
分类专栏： Matlab解微分方程 文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45024585/article/details/107601153
版权

隐式微分方程的数学描述为：
在这里插入图片描述
这是一个关于时间变量t，状态变量，状态变量一阶导数的函数。
已知状态变量的初值，状态变量一阶导数的初值。

相容的初始条件的计算：
在这里插入图片描述
得到相容的初始条件之后，就可以调用ode15i来求解：
在这里插入图片描述
例题：
在这里插入图片描述
选完状态变量之后，按照常规的方法，我们就可以把方程描述为这种形式：

最低0.47元/天 解锁文章
进击的阿尔法猿
关注
  • 2
    点赞
  • 踩
  • 33
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
专栏目录
Matlab中如何求解隐式常微分方程
warnerchang的博客
04-01 3780
在matlab中求解隐式常微分方程，通常有两种方法，一种是使用matlab中的内嵌的ode15i等求解隐式微分方程的函数，另一种是仍使用matlab中求解常微分方程的显式函数，如ode45等等。但是小编想要表达的绝不是这两种方法，而是这两种方法之间的差异： 一、ode45 使用这种方法求解隐式微分方程的思路就是结合其他函数求解隐式方程，例如fzero等等。 以求解如下隐式方程为例： 可以发现这个式子很复杂，没有显式式子，因此只能使用隐式函数求解。代码如下所示： % Matlab %% ode4
特殊微分方程的数值解——刚性微分方程（ode15s）
醉后不知天在水，满船清梦压星河
07-26 4815
我们在微分方程的标准型转换中已经提到的Van der Pol方程就是一种刚性微分方程。 刚性微分方程可能有几个状态变量，其中某个状态变量变化的特别快，远远的快于其他的状态变量，这样就会使得常规算法中，选择步长出现问题。 再谈Van der Pol方程： 计算mu=1000时的Van der Pol方程的数值解。 h_opt=odeset; h_opt.RelTol=1e-6; x0=[2;0]; t_final=3000; f=@(t,x,mu)[x(2);-mu*(x(1)^2-1)*x(2)-x(1)
matlab使用教程(27)—微分代数方程(DAE)求解
最新发布
配电网和matlab的博客
08-28 1370
如果您的方程微分指数为 2 或更高，则需要将方程重写为微分指数为 1 的等效 DAE 方程组。通过对方程进行微分，可以消除代数变量，并且如果执行此操作的次数足够多，这些方程将呈现为显式 ODE 方程组。DAE 方程组的微分指数是为了将方程组表示为等效的显式 ODE 方程组必须执行的求导次数。此方程组的微分指数为 1，因为只需y3的一个导数就能使其成为 ODE 方程组。0 指定了不一致的初始条件，则求解器会将这些值作为估计值进行处理，尝试计算接近估计值的一致值，并继续解算该问题。然而，一个明显的例外是使用。
2.4 一阶隐式微分方程与参数表示
universe_1207的博客
01-16 4019
文章目录类型一：y=f(x,dydx)y=f(x,\frac{dy}{dx})y=f(x,dxdy​)类型二：x=f(y,dydx)x=f(y,\frac{dy}{dx})x=f(y,dxdy​)一个题解类型三：不显含y的方程F(x,y′)=0F(x,y')=0F(x,y′)=0类型四：不显含x的方程F(y,y′)=0F(y,y')=0F(y,y′)=0 一阶隐式微分方程通常表示：F(x,y,y′...
ode45、ode15s、ode15i怎么选？
qq_40380771的博客
06-27 3524
%主程序-求解器：ode45 clear all close all clc y0 = 0; tspan = [0:0.01:10]; [t,y] = ode45(@subfun,tspan,y0); figure(1) plot(t,y);% 画出每次积分结果与时间的图形 %主程序-求解器：ode15s clear all close all clc y0 = 0; yp0 = 1; tspan = [0:0.01:10]; options = odeset('RelTol',1e-12); [t
特殊函数的常微分方程
weixin_38060850的博客
12-04 400
隐式微分方程求解Matlab
村头
08-24 3716
在matlab求解微分方程有很多中，ode45是最常见的，当然，ode45求解有一定条件，就是需要把高阶微分方程转化为1阶微分方程。即使用降阶法。降阶法原理就是引入中间变量，然后实现降阶。F(t,y,y′,y′′)=0令x1​=y,x2​=y′,那么就可以将原来的二阶微分方程降阶为两个一阶微分方程。dtdx1​​=x2​dtdx2​​=G(t,x1​,x2​)​高阶微分方程组同理！！！
python解隐式方程_求解隐式ODE（微分代数方程DAE）
weixin_29955195的博客
01-30 1854
如果代数操作失败，可以对约束进行数值求解，例如在每个时间步运行fsolve：import sysfrom numpy import linspacefrom scipy.integrate import odeintfrom scipy.optimize import fsolvey0 = [0, 5]time = linspace(0., 10., 1000)F_lon = 10.mass = ...
特殊微分方程的数值解——微分代数方程
醉后不知天在水，满船清梦压星河
08-04 2660
所谓的微分代数方程（DAE），是指在微分方程中，某些变量间满足某些代数方程的约束。 微分代数方程的一般形式为： 其中，M（t，x）矩阵的为奇异矩阵。在求解微分方程时，两边不能同时乘以M矩阵的逆矩阵M‘，所以这类方程并不能变成普通的微分方程取求解。 例题1： 假设我们这里有三个微分方程，第一个和第二个是常规形式的微分方程，第三个方程退化成了一个代数方程，所以这类方程又称为微分代数方程。 给定初值条件，x1(0)=0.8, x2(0)=x3(0)=0.1 用矩阵形式表示该微分代数方程： f=@(t,x)[
使用MATLAB求解常微分方程数值解
10-15
"使用MATLAB求解常微分方程数值解" 本文将从基本概念开始，逐步深入探讨使用MATLAB求解常微分方程数值解的方法和原理。常微分方程是数学中一种重要的数学模型，广泛应用于物理、生物、经济等领域。MATLAB作为一个...
ODE15.rar_ode15s_tank4gv_刚性_刚性微分_向后微分
07-15
ODE15s可以处理刚性微分方程。采用向后微分公式，效率不是很高，但是对于一些ODE45会失效或计算效果很差的情况下，可以使用。可能是由于方程刚性引起的。
matlab 常微分方程数值解法 源程序代码.zip_adams多步法_bvp微分方程组_matlab 解方程组_ode_刚性方
07-14
　11.4 微分方程组的数值解 397 　11.4.1 Euler方法 397 　11.4.2 经典四阶Runge-Kutta方法 398 　11.4.3 高阶方程组的求解 399 　11.5 刚性方程组的数值解 401 　11.5.1 梯形公式 401 　11.5.2 隐式Runge-Kutta方法...
MATLAB求解常微分方程数值解.docx
10-09
MATLAB求解常微分方程数值解 本文主要介绍了使用 MATLAB 求解常微分方程的数值解法，包括 Euler 法和 Runge-Kutta 法两种方法，并对其进行了详细的讲解和对比。 一、常微分方程的数值解法 常微分方程是数学中的一...
MATLAB求解常微分方程数值解.doc
10-09
本文主要介绍了使用 MATLAB 求解常微分方程数值解的方法，包括 Euler 法、Runge-Kutta 法和 MATLAB 自带的 ode45 函数。通过对这些方法的详细介绍和比较，可以帮助读者更好地理解常微分方程的数值解法，并掌握使用 ...
隐函数求微分
ResumeProject的博客
12-24 7412
链式求导法则(隐函数求导法则)： 将方程两边都对x求导，有y的地方，先当成y的函数，对y求导，然后再将y对x求导。最后解出dy/dx，dy=y’dx。 直接微分(因为具有形式的不变性，所以可以不考虑变量之间的内在关系)：将x、y看成等同地位，谁也不是谁的函数，方程两边微分，解出dy即可。 x2+y2=1直接微分:2xdx+2ydy=0⇒dy=−xydx链式法则：2x+yy′=0⇒y′=−xy,dy=y′dx x^2+y^2=1\\ 直接微分:2xdx+2ydy=0 \Rightarrow dy=\frac{
c#解微分方程_使用C＃开发符号表达库。 微分，简化，方程求解等
cullen2012的博客
09-07 2759
c#解微分方程 Why does programming a calculator seem to be a task, which every beginner undertakes? History might have the answer — computers were created for this exact purpose. Unlike the beginners, we w...
几类常微分方程的matlab求解方法 | 刚性微分方程、隐式微分方程、微分代数方程
weixin_42019459的博客
03-07 9527
参考：《高等应用数学问题的MATLAB求解（第四版）》 本文涵盖微分方程的转换、刚性微分方程、隐式微分方程、微分代数方程求解方法。
matlab中ode15s函数的用法,求解n个常微分方程构成的方程组，如何使用ode15s
weixin_34221892的博客
03-18 3226
我的方程组涉及的方程重复性较高，而且数量较多，故采用i和for循环来编写方程，不知道哪里有问题，下面是我的代码function dy = ad_mass (t,y)N=100;%间隔数e=0.44;%床孔隙率Dl=11;%轴向弥散系数，cm^2/sus=40.7;%表观速度,cm/su=us/e;%间隙速度,cm/sleng=25.5;%吸附剂填充长度,cmz=leng/N;%间隔长度,cmR=8...
matlab ode15,matlab用ode15数值计算微分代数方程（DAE）的问题
weixin_30082413的博客
03-21 2584
代码如下：前四行都是会用到的参数，odefun中有3个微分方程与1个代数方程lamdap=975e-9;lamdas=1064e-9;t=1e-3;ap=2100e-27;ep=2200e-27;A=2.826e-7;R1=0.6;as=5e-27;es=331e-27;alfap=8e-3;alfal=6e-3;tp=0.8e-3;N=2.1e28;R2=1;h=6.626e-34;c=3e8;...
用MATLAB解线性微分方程组的刚性比
05-16
用MATLAB解线性微分方程组的刚性问题主要体现在数值方法上。在数值方法中，常规的显式方法（如欧拉方法）对于刚性问题的求解效果较差，即使取很小的步长也不能保证数值解的稳定性和精度。而隐式方法（如Crank-Nicolson方法、后向欧拉法、龙格-库塔法等）在求解刚性问题时表现更为优秀。 MATLAB提供了许多数值解法用于解决刚性问题，比如ode15s、ode23s等，它们采用了隐式方法来求解微分方程组，能够有效地解决刚性问题。此外，MATLAB还提供了一些工具箱用于求解刚性问题，如ode15i、ode15t等。 总之，用MATLAB解线性微分方程组的刚性问题需要选用适当的数值方法，通过调整步长和求解参数等方式来保证数值解的稳定性和精度。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值