刚性微分方程与非刚性的大概区分(自看,求指点)
(这里的内容是作为自己的一个粗略的总结,不确定是否正确,希望有大佬能够更明确的指出其中的错误,作出指导)
对于刚性和非刚性微分方程的区分,可以简单的转变为在将原方程转换为常微分方程组后,进行一个简单的系数判断:
例如:
y''' - 3y'' - y'y = 0, y(0) = 0, y'(0) = 0, y''(0) = -1;
在这里可以设 y1 = y, y2 = y', y3 = y'', 有
y1' = y2, y1(0) = 0,
y2' = y3, y2(0) = 1,
y3' = 3y3 + y2y1, y3(0) = -1,
这里可以简单判断出方程组的右侧系数矩阵值差异不大,得到的特征值差异随之也不大,可以简单判断为非刚性微分方程。
MATLAB中解非刚性微分方程常用 ode45 ;
又例如:
y'' - 1000(1 - y^2)y' + y = 0, y(0) = 2, y'(0) = 0;
在这里可以设y1 = y, y2 = y', 有
y1' = y2, y1(0) = 2,
y2' = 1000(1 - y1^2)y2 - y1, y2(0) = 0,
这里可以简单判断出方程组的右侧系数矩阵值差异较大,得到的特征值差异随之较大,可以简单判断为刚性微分方程。
MATLAB中解刚性微分方程常用 ode15s , ode23s , ode23t , ode23tb ;
对于所有的初值问题解方程方法均有 solver( 'f( t, y )' , [t0 tend], y0 ) :
--solver : 指代所有的方法调用;
--f( t, y ) : 为M文件定义的微分方程 y' = f(x, y) 右端的函数, 即使用不到两个参数,也要保留两个参数,且返回的结果应是列向量;
--[t0 tend] : 为参数中变量 x / t 的取值范围;
--y0 : 为参数中 y 的初始值,要求为列向量
转载于:https://www.cnblogs.com/yanmingjiang/p/9426810.html
刚性微分方程与非刚性的大概区分(自看,求指点)相关推荐
- 像“今日头条”这样的模块滑动特效,大概咋写,求指点 ?---酷课堂iOS交流群问答整理(201805期)
酷课堂iOS交流群 我们是一个什么样的组织: 酷课堂iOS交流群,聚集了一群热爱技术.有趣.有料,平均Q龄在10年以上的"老司机",他们遍布在全国/球各地,有知名企业iOS工程师. ...
- MATLAB 数学应用 微分方程 常微分方程 求解非刚性ODE
本文介绍两个使用 ode45 来求解非刚性常微分方程的示例.MATLAB拥有三个非刚性 ODE 求解器. ode45 ode23 ode113 对于大多数非刚性问题,ode45 的性能最佳.但对于允许 ...
- 几类常微分方程的matlab求解方法 | 刚性微分方程、隐式微分方程、微分代数方程
目录 微分方程的转换 一.单个高阶常微分方程 二.高阶常微分方程组 刚性微分方程求解 隐式微分方程求解 微分代数方程求解 微分方程的转换 根据微分方程求解的标准型,要得到微分方程的数值解,应该先将该方 ...
- 非刚性配准(Non-rigid ICP )
原文https://blog.csdn.net/linmingan/article/details/79270874?utm_medium=distribute.pc_relevant.none-ta ...
- 《Master Opencv...读书笔记》非刚性人脸跟踪 II
上一篇博文中,我们了解了系统的功能和模块,明确了需要采集哪些类型的样本点及利用类的序列化的保存方式.这次将介绍几何约束模块,通过统计形态分析法(Statistical Shape Analysis, ...
- ASTM D 4804非刚性固体塑料可燃性试验
ASTM D 4804测定非刚性固体塑料可燃性的标准试验方法– 标准名称 ASTM D 4804 Standard Test Method for Determining the Flammabili ...
- 信捷伺服刚性调整_伺服电机刚性怎么调整
伺服电机是指在伺服系统中控制机械元件运转的发动机.是一种补助马达间接变速装置.可使控制速度,位置精度非常准确.将电压信号转化为转矩和转速以驱动控制对象. 调整伺服电机刚性 1.伺服主要靠脉冲来定位,基 ...
- 烽火HG680-L_HG680-LV_hg680-lc_hg680-ly高安及非高安版本区分及刷机说明
烽火HG680-L_HG680-LV_hg680-lc_hg680-ly高安及非高安版本区分及刷机说明 峰火hg680-L分高安和非高安版本: 河北的100%是高安,其他地区基本是非高安: 高安又分S ...
- 想学非标机械设计的朋友看过来!不然后悔终生!
想学非标机械设计的朋友看过来!不然后悔终生! 一.非标机械设计具体是什么? 简单的说就是非标准设备以及非标准件. 非标准工装的设计. 举例说: C616 C620 C630 C6125 M7475 B ...
最新文章
- C# 打好基础图文详解
- any() missing 1 required positional arguments: dim
- jenkins+Docker+springcloud微服务持续集成
- wxWidgets:wxLogFormatter类用法
- tl494组成的超声波发射电路_透彻分析利用TL494组成的逆变器电路
- Hibernate中对增删改查的小结
- 眼图观测实验报告_眼图观察测量实验
- CTF工具-gdb插件peda,pwndbg安装
- ACdream区域赛指导赛之手速赛系列(2)
- html标签的引号嵌套,使用YQL多查询& XPath解析HTML,如何转义嵌套引号?
- msyql数据类型对照java类型
- 汇编语言上机考试三星题——加密的key和明文字符串
- python settings_python settings如何配置sqlite?
- java用dda算法绘制直线_使用DDA算法画出给定两点的直线
- 2020switch电信最快的dns_《2020switch电信最快的dns》电影_2020switch电信最快的dns正片免费观看-深圳市创业发科技有限公司...
- python实现C4.5
- 理解物理 CPU 与逻辑 CPU
- 人工智能期末考试复习(贲可荣 张彦铎)
- 使用scrapy图片管道下载图片
- VMware安装Win11+WSA子系统和使用教程
热门文章
- 功能测试非常重要的43个功能测试点...
- MacBook 安装 Docker 与 docker-compose 套装
- 半双工通信和全双工通信的区别
- 苹果服务器维护不能刷机,iphone刷机失败不开机报错维修方法分享
- python支持的编程方式包括,python的两种编程方式是什么
- Qt之调用FFTW3实现音频频谱(实现)
- poj2976 01分数规划
- 使用字符映射完成简单的对称加密
- python贪吃蛇控制台_通过 PRA 的 python-pygame 玩转贪吃蛇
- java毕业设计二次元信息分享平台(附源码、数据库)