偏微分方程算法之高精度算法（Richardson法+紧差分法）

已于 2024-04-13 10:46:25 修改
阅读量1.3k 收藏 25
点赞数 30
分类专栏： 偏微分方程（抛物型） 文章标签： 算法
于 2024-04-13 10:41:06 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/L_peanut/article/details/137628912
版权

目录

一、研究对象

二、理查德森外推法（Richardson法）

2.1 理论推导

2.2 算例实现

三、紧差分法

2.1 理论推导

2.2 算例实现

一、研究对象

        在常微分方程算法中求解关于两点狄利克雷（Dirichlet）边值问题时，我们通过利用理查德森外推法、紧差分法在几乎不增加计算成本的情况下，将二阶精度的数值解通过线性组合构造了四阶精度的数值解。同样，我们也可以将这些高精度算法应用于偏微分方程计算中。

        研究对象与向前、向后欧拉法相同，为抛物型方程初边值问题：

了解本专栏 订阅专栏 解锁全文 超级会员免费看
猿核试Bug愁
关注
  • 30
    点赞
  • 踩
  • 25
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
专栏目录
订阅专栏
常微分方程算法之高精度算法（Richardson法+紧差分法）
欢迎各位大佬光临“寒舍”！请多多指点，共同进步！
04-08 1805
常微分方程算法-提高求解精度的方法及C++实现
matlab开发-RichardsonExtrapolation
08-23
matlab开发-RichardsonExtrapolation。一个有许多钟和哨声的广义理查森外推程序。
Matlab实现解抛物型方程求解
07-01
1.版本：matlab2019a，不会运行可私信 2.领域：基础教程 3.内容：Matlab实现解抛物型方程求解 4.适合人群：本科，硕士等教研学习使用
偏微分方程数值解法的MATLAB源码--古典显式格式求解抛物型偏微分方程等
06-23
1、古典显式格式求解抛物型偏微分方程（一维热传导方程） 2、古典隐式格式求解抛物型偏微分方程（一维热传导方程） 3、Crank-Nicolson隐式格式求解抛物型偏微分方程 4、正方形区域Laplace方程Diriclet问题的求解 如： function [U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) %古典显式格式求解抛物型偏微分方程 %[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) % %方程：u_t=C*u_xx 0 <= x <= uX,0 <= t <= uT %初值条件：u(x,0)=phi(x) %边值条件：u(0,t)=psi1(t), u(uX,t)=psi2(t)
抛物型偏微分方程的Crank-Nicolson 方法； Richardson 外推法；紧差分法
liuqihang11的博客
01-07 4594
抛物型偏微分方程的Crank-Nicolson 方法； Richardson 外推法；紧差分法
理查森外推法 matlab,数值代数–理查森外推法.doc
weixin_39559277的博客
03-25 1492
数值代数–理查森外推法实验四一、实验名称理查森外推算法二、实验目的与要求：实验目的：掌握理查森外推算法。实验要求：1. 给出理查森外推算法思路，2. 用C语言实现算法，运行环境为Microsoft Visual C++。三、算法思路：1. 假设函数泰勒展开式可表示为 和，将两式相减，消去偶数项，则，整理得到下式，记L表示，表示微分形式，则有 (1)用h/2代替h，有 (2)，由(1)(2)两式子...
理查生外推（Richardson extrapolation）
Angelaboy的博客
01-05 883
理查生外推（Richardson extrapolation）
理查森导数外推法：点中单值实函数的一阶和二阶导数的理查森外推法。-matlab开发
05-29
函数 RICHARDSONDER 在给定中心差分公式的初始步长 H_IN 和外推的阶数 N 的情况下，实现了理查森外推算法，用于逼近 X0 点单值实函数 F 的一阶和二阶导数。 输出是一阶和二阶导数 F_PRIME_X0 和 F_SECOND_X0 的两个外推表。 每个表的第一列是使用的不同步长值的数组，而其他列包含导数外推； 从第 3 列开始的每一列的固定 h 的错误顺序比其左侧的低（表中的 0 个条目是未计算的值）。 因此，最佳近似值是每个表的最后一个元素（最小误差和最小步长）。 有关该函数的用法示例，请参阅 RICHARDSONSCRIPT。
Richardson外推法计算给定点处的一阶和二阶导数
u011378809的专栏
08-13 3908
//计算给定点处的一阶和二阶导数 #include #include using namespace std; class deriv { private:  int k;  double d1, d1_1, d1_2, d2, d2_1, d2_2, d1_new, d2_new;  double f_x, f_xhm, f_x2hm, f_xhp, f_x2hp, h,
数值分析中理查德森外推求导数
05-31
利用理查德森外推求函数的导数,避免人工求导的过程，程序稳定，调试成功。
理查森外推：具有许多花里胡哨的广义理查森外推例程。-matlab开发
06-01
一个通用的理查森外推例程，可以执行多个外推步骤，可以容纳任意领先的误差项，可以处理任何统一的细化方案（例如 h/2、h/3、h/5 等），并且可以报告完整的外推表，包括与已知的良好解决方案进行比较。 能够处理矢量输入。 与 MATLAB 和 Octave 兼容。 包括 Octave 风格的单元测试。 一个人为的用法示例： 假设您的近似序列是 [1, 2, 3, 4]，其中每个连续值都将网格细化了 3 倍。假设近似中的三个主要误差项是 O(h^2)、O(h^4 ) 和 O(h^5) 并且您知道确切值是 4.14。 然后你会调用 Octave:78> [result, normtable] = richardson([1,2,3,4],[2,4,5], 3, 4.14) 结果 = 4.14163223140496 常态 = 3.13999999999999968 NaN NaN NaN
广义理查德森迭代法matlab编程
07-04
广义理查德森迭代法引入了n个迭代参数，其收敛条件为A是对称正定矩阵，可用于求解AX=b的解！
理查森外推法：它根据理查森外推法计算函数的导数。-matlab开发
05-30
它将理查森外推法应用于泰勒级数，以使用“n”次迭代来逼近任何函数 f(x) 在 x_0 处的导数。 这是一个 O(n^2) 算法，可以在“数值数学和计算，Ward Cheney 和 David Kincaid，第 6 版”第 4.3 节中找到。
非线性偏微分方程几种经典算法
12-11
对非线性偏微分方程几种经典算法进行了综述，给出了推导过程
抛物型偏微分方程的有限差分法（2）
01-26
【抛物型偏微分方程的有限差分法（2）】 抛物型偏微分方程在数值计算中通常使用有限差分法来求解。本篇内容主要介绍了Crank-Nicolson方法，这是一种结合了向前欧拉显格式和向后欧拉隐格式优点的数值解法。向前欧拉...
MATLAB常用算法程序集 偏微分方程的数值解法.zip
最新发布
08-19
这个压缩包包含了多个MATLAB程序，专门针对不同类型的偏微分方程进行数值求解，对于学习和应用偏微分方程的数值方法具有很高的参考价值。 1. **椭圆型偏微分方程（Elliptic PDEs）**： - `peEllip5.m` 和 `peEllip...
常微分方程的数值解法（5）
01-26
Richardson 外推法的主要优点是可以提高低阶差分法的计算精度，从而提高计算的速度和精度。此外，Richardson 外推法也可以应用于各种常微分方程的数值解法。 5. Richardson 外推法的缺点 Richardson 外推法的主要...
Matlab偏微分方程的数值解法常用程序-偏微分方程的数值解法_程序.rar
08-13
3. peParabImp.m：这个名字暗示了该程序可能实现了隐式方法，如有限差分法或有限元法，用于求解抛物型偏微分方程。隐式方法的优点在于其稳定性，即使在大的时间步长下也能保持数值解的稳定。 4. EllIni2Uyr.m：这个...
Matlab数值计算方法程序源代码 算法设计及其MATLAB实现 共93页.pdf
05-29
有限差分法是数值分析中的基础工具，用于离散化偏微分方程。 第四章讨论了方程求根的各种算法。二分法是最基础的连续函数零点寻找方法，开方法则是基于牛顿-拉弗森迭代的一种改进。Newton下山法，也称为牛顿法，...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

猿核试Bug愁

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值