微分方程数值解
文章平均质量分 82
用python求解常微分方程、偏微分方程
浪子膏丷
这个作者很懒,什么都没留下…
展开
-
【有限差分法】(三)一维和二维抛物方程CN格式以及长时间稳定性分析(附算例与Python代码)
前言我们知道差分里的CN格式是无条件稳定的。但是最近在学习有限元结合CN格式算长时间抛物问题的时候(在时间方向用差分空间方向用有限元)发现稳定性却不能保证,其数值解和真解误差会随着时间越来越大。并且实际能算的时间只有几秒,这几秒钟几乎是没有实际意义的。所以就想到一个问题,纯CN差分格式在长时间计算的时候,是否也是个理论看还行,实践臭弟弟的“花瓶”呢?在以前做差分课本上的算例大多也是只有几秒,还真没注意到长时间稳定性的问题(¬_¬ )。所以本文找了最简单的一维抛物问题,把时间层增加至1000S,来看看C.原创 2021-11-28 17:53:46 · 6316 阅读 · 0 评论 -
【有限差分法】(一)有限差分法的基本流程与常用格式
2.有限差分法的引入与介绍有限差分法:通过将求解区域离散化得到离散的微分算子,将微分方程定解问题转化为代数方程组求解问题。有限差分法的原理:我们在用数值方法求解微分方程时,由于计算机不能处理连续的数学方程,所以在用计算机处理问题时,必须先将一个无限的连续问题通过一些特殊的方法离散化,使之成有限形式的线性代数方程组,再将这个线性代数方程组交给计算机求解以达到我们的目的。在偏微分方程中,通过离散化区域获得导数(或偏导数)的手段就构成了我们的有限差分法。3.有限差分法的流程以一维情况为例来说明一个有限原创 2021-11-25 22:12:45 · 19448 阅读 · 1 评论 -
【微分方程数值解】常微分方程(一)欧拉方法和改进欧拉方法(附python算例,封装类)
aa原创 2021-04-20 22:19:51 · 4685 阅读 · 0 评论 -
【微分方程数值解】有限差分法(二)两点边值问题数值算例(附python代码)
什么是一维椭圆问题形如u′′(x)=af(x)的形式。标准形式一边是关于u关于x的二阶导数,一边是关于x的函数。求解思路我们现在用有限差分法解下面一个简单的一维弦震动问题−u′′(x)=16π2sin(4πx)其中边值条件为u(0)=0,u(1)=0.根据我们之前的内容,可以将解此题的内容分为四步:对求解域进行区间剖分用差分的格式代替对方程中对各阶导数进行边界的处理解代数方程...原创 2019-05-25 10:30:29 · 10925 阅读 · 0 评论 -
【微分方程数值解】常\偏微分方程及其常用数值解法概述
博主在这段时间将不断更新有限差分法在偏微分方程中的应用,并以python为工具讲解三类偏微分方程的有限差分法实现。 直通链接: ... 未完待续有限差分法的引入与介绍==========微分方程(model)是描述现实问题的重要模型,很多重要的现象都是用微分方程描述刻画的,例如记忆棉受力分析、声波传播、弦振动问题等等。将一个问题转化为微分方程后,如何得到尽可能逼近真解的数值解,是我们...原创 2019-03-31 10:29:24 · 10787 阅读 · 0 评论