2021-01-13 Matlab求解微分代数方程 (DAE)

Matlab求解微分代数方程 (DAE)

什么是微分代数方程？

微分代数方程是一类微分方程，其中一个或多个因变量导数未出现在方程中。方程中出现的未包含其导数的变量称为代数变量，代数变量的存在意味着不能将这些方程记为显式形式 y′=f(t,y)。

ode15s 和 ode23t 求解器可以使用奇异质量矩阵 M(t,y)y′=f(t,y) 来解算微分指数为1的线性隐式问题，包括以下形式的半显式 DAE

y′0=f(t,y,z)

0 =g(t,y,z)

在此形式中，由于主对角线存在一个或多个零值，因此代数变量的存在会产生奇异质量矩阵。

默认情况下，求解器会自动检验质量矩阵的奇异性，以检测 DAE 方程组。如果提前知道奇异性，则可将 odeset 的 MassSingular 选项设为 'yes'。对于 DAE，还可以使用 odeset 的 InitialSlope 属性为求解器提供 y′(0) 的初始条件估计值。

举个例子

其中x1(0)=0.8;x2(0)=x3(0)=0.1;

1）方程写成DAE形式

2）编程求解

%% clcclearclose allodefun = @(t,x)[-0.2*x(1)+x(2)*x(3)+0.3*x(1)*x(2);                2*x(1)*x(2)-5*x(2)*x(3)-2*x(2)^2;                x(1)+x(2)+x(3)-1]; %微分方程M = [1 0 0;0 1 0;0 0 0]; % 质量矩阵options=odeset('mass',M); % 定义mass属性x0=[0.8;0.1;0.1];[t,x]=ode15s(odefun,[0 10],x0,options);figureplot(t,x(:,1),t,x(:,2),t,x(:,3))grid onlegend('x1','x2','x3')

2021-01-13 Matlab求解微分代数方程 (DAE)相关推荐

matlab求解微分代数方程组,微分代数方程(DAE)的Matlab 解法.PDF
微分代数方程(DAE)的Matlab 解法微分代数方程(DAE)的Matlab解法所谓微分代数方程,是指在微分方程中,某些变量满足某些代数方程的约束.假设微分方程的更一般形式可以写成前面所介绍 ...
用matlab求解线性代数方程
用matlab求解线性代数方程 % 功能:线性代数方程求解:AX=B % 解的判定矩阵:C=[A B] % 求解方法分三种情况讨论: % 唯一解:A为非奇异方阵,x=inv(A)*B % 无穷解:ra ...
龙格库塔法matlab求解微分方程组,微分方程组的龙格库塔公式求解matlab版.pdf
微分方程组的龙格库塔公式求解matlab版微分方程组的龙格-库塔公式求解matlab版南京大学王寻 1. 一阶常微分方程组考虑方程组     y'f x,y,z , y x y ...
基于MATLAB的微分代数方程解法（附完整代码）
目录一. 微分代数方程求解例题1 二. 全隐式微分方程三. 延迟微分方程求解例题2 一. 微分代数方程求解例题1 初始条件: 求数值解: 解: ①方法1求解矩阵形式表示该微分代数方程: ( ...
Matlab求解微分方程组
我们采用ode方法: (1)求解普通微分方程组:使用ode45方法 1. 创建一个函数文件eq2.m,在函数文件中描述这个解的微分方程组: %eq2.m文件 %描述微分方程组function dy=e ...
数学建模学习（29）：matlab求解微分方程组详细讲解，代码+案例讲解，学不会找我！
文章目录前言求解微分方程组求解矩阵微分方程组总结前言上一篇我已经详细讲过求微分方程,这一篇是对上一篇的补充,也就是变得稍微复杂一点,就是要求方程组了,如果你学会了上一篇,那么求解方程组其实 ...
python解隐式方程_求解隐式ODE（微分代数方程DAE）
如果代数操作失败,可以对约束进行数值求解,例如在每个时间步运行fsolve:import sys from numpy import linspace from scipy.integrate imp ...
MATLAB怎么解方程解,怎么用MATLAB求解微分方程组并画出解函数图？
!using["XSLSF"]; //使用命名空间XSLSF //数组xArray存放x的值:ti为当前有效值的个数:tmax为ti对应的时间:tmi ...
（三十七：2021.01.13）Pre-MICCAI 2019学习（二）《前列腺近距离放射治疗中，检测粒子在CT中的3D位置和方向》
<Automatic 3D seed location and orientation in CT images for prostate brachytherapy> 讲在前面摘要论 ...

2021-01-13 Matlab求解微分代数方程 (DAE)

Matlab求解微分代数方程 (DAE)

什么是微分代数方程？

2021-01-13 Matlab求解微分代数方程 (DAE)相关推荐

最新文章

热门文章