偏微分方程算法之克兰克-尼科尔森法（Crank-Nicolson法）

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

一、研究对象

       从偏微分方程的向前欧拉法和向后欧拉法的计算结果分析可知，向前欧拉法虽然计算简单但稳定性差，收敛要求高，只有选择合理的时间、空间步长时才能确保数值计算收敛；向后欧拉法则不需要考虑时间、空间步长的选取，但需要在每个时间层上求解一个线性方程组，计算量大。同时，两种方法在时间、空间步长均减半的情况下，计算误差无法确定性的减少1/4，而是减少1/2到1/4。为了实现误差与步长同比例变化，我们可以通过对时间的一阶偏导数采用二阶中心差商的方式，提升收敛速度。

       我们继续以与向前、向后欧拉法相同的抛物型方程初边值问题为研究对象：

了解本专栏 订阅专栏 解锁全文 超级会员免费看
猿何试Bug个踌
关注
  • 10
    点赞
  • 踩
  • 14
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
专栏目录
订阅专栏
Crank-Nicolson格式求解二维抛物线方程
SDU MATH
11-25 2141
本文给出了CN格式求解二维抛物线方程的MATLAB代码以及验证了收敛阶
一维薛定鄂Schrodinger方程解法Crank-Nicholson 法
09-09
一维薛定鄂Schrodinger方程解法Crank-Nicholson 法 复数三对角矩阵追赶法
热传导偏微分方程Crank-Nicloson格式附MATLAB
03-02
热传导偏微分方程差分，使用Crank-Nicloson格式，附MATLAB实例
Numerical-Solution-PDE-.rar_PDE算法_crank method_elliptic equation
09-14
偏微分方程的数值解法的MATLAB程序介绍了应用最为广泛的椭圆型、双曲型、抛物型偏微分方程的数值解法，而且还详细编程实现了每种方程的多种常见数值解法。 附件使用MATLAB编程来实现这些算法
Crank-Nicolson隐式方法解线性双曲型方程
最新发布
程序员的博客
04-28 356
线性双曲型方程的Cranck-Nicolson隐式解法
隐式克兰克—尼科尔森方法
08-28
安德森《计算流体力学及其应用》第九章不可压库埃特流的数值解 9.3节 数值方法：隐式克兰克—尼科尔森方法。
【转载】crank-nicolson求解一维扩散方程
cyberck的博客
08-25 946
【Matlab】抛物方程Crank-Nicolson格式
热门推荐
See_Star的博客
12-26 1万+
1 题目 (1) 编制用Crank-Nicolson格式求抛物方程数值解的通用程序。 {∂u∂t−a∂2u∂x2=f(x,t)(0<x<1,0≤t≤ T)u(x,0)=φ(x)(0≤x≤1)u(0,t)=α(t),u(1,t)=β(t)(0<t≤T) \left\{\begin{array}{ll} \frac{\partial u}{\partial t}-a \frac{\partial^{2} u}{\partial x^{2}}=f(x, t) & (0<
crank-Nicolson方法
weixin_42245449的博客
07-06 9703
http://blog.sina.com.cn/s/blog_61c5b0020100zkw7.html 追赶法 https://wenku.baidu.com/view/885b27a40029bd64783e2c8a?fr=step_zhidao
持续不定期更新：CFD&C++之不可压库埃特流（Coutte）的数值解（隐式克兰克-尼科尔森方法）
陈十七的博客
01-13 1707
之前一篇https://blog.csdn.net/weixin_39124457/article/details/101038218用到的是显式方法，这次尝试用隐式方法解决不可压库埃特流的问题。 资料参考《计算流体力学基础及其应用》第9章第一节，电子书可自行搜索下载。 代码自己打的，可下载：https://download.csdn.net/download/weixin_39124457/...
Crank-Nicolson
renbaifen的博客
11-08 593
function [error,h,k] = heateq_CN(a,b,ti,tf,m,kappa) % a = 0; % b = 1; % m = 100; % ti = 0;%initial time % tf = 1.1;%final_time % kappa = 0.02; % Discretization set up h = abs(b-a)/(m+1); % h = delta x k = 4*h; .
抛物型偏微分方程的Crank-Nicolson 方法； Richardson 外推法；紧差分法
liuqihang11的博客
01-07 4577
抛物型偏微分方程的Crank-Nicolson 方法； Richardson 外推法；紧差分法
para-crank Nicolson scheme_抛物型方程_抛物方程_crank-nicolson_
09-29
Crank-Nicolson格式求解抛物型方程
偏微分方程数值解法的MATLAB源码--古典显式格式求解抛物型偏微分方程等
11-05
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)
使用CN方法获得一维热传导方程的稳态解。：使用Crank-Nicolson方法求解一维热方程并绘制等高线图-matlab开发
06-01
我们有兴趣使用 CN 方法获得一维热传导方程的稳态解。 边界条件是：在 x=0 和 0.3 m 处 T=300 K，在所有其他内部点处 T=100 K。 α = 〖3*10〗^(-6) m-2s-1 . 这里，t=30 分钟，Δx=0.015m 和 Δt=20 秒
偏微分方程数值解CH2PPT教案学习.pptx
10-02
在讲解的隐式差分格式中，比如尼科尔森（Nicolson）和克兰克-尼科尔森（Crank-Nicolson）格式，是用于求解热传导等扩散问题的常用方法。克兰克-尼科尔森格式是一种时间-空间半离散的隐式格式，它可以提供比显式格式...
热传导方程有限差分法的MATLAB实现
07-02
有限差分法是一种数值分析方法，广泛应用于求解偏微分方程，特别是在物理、工程和科学计算领域，如热传导问题。MATLAB作为强大的数值计算工具，提供了一个理想的平台来实现这种算法。 首先，热传导方程描述了热量...
crank-nicolson差分格式原理
06-28
### 回答1： Crank-Nicolson差分格式是一种常用的数值方法，用于求解偏微分方程的数值解。实际应用中，很多偏微分方程都难以求出解析解，因此需要通过数值计算的方式来求得近似解。 Crank-Nicolson差分格式的原理是将时间轴上的偏微分方程离散化，然后通过差分运算得到一个方程组，进而求解该方程组得到数值解。在离散化过程中，将时间和空间都离散化为若干个区间，其间隔大小分别为Δt和Δx，以此来逼近连续的偏微分方程。对于每个时间步长，我们分别在时间和空间方向上进行更新，进而得到一个更新后的解。 具体来说，Crank-Nicolson差分格式采用的是一个二阶中心差分格式，该格式是一个隐式的数值差分方式，对数值解的数值稳定性和收敛性都有极好的保证。其基本思想是将一个解的时间和空间方向上的差分分别计算出来，然后取两者的平均值作为解的近似值。通过这种方式，能够减少误差，提高计算精度。 总之，Crank-Nicolson差分格式是一种基于二阶中心差分的隐式数值方法，它能够很好地解决偏微分方程计算的误差和收敛性问题，因此在实际应用中受到广泛的欢迎。 ### 回答2： Crank-Nicolson差分格式是一种二阶精度的有限差分格式，被广泛应用于偏微分方程的数值求解，尤其适用于热传导方程、扩散方程等具有稳定解的情况。 其基本思想是将时间步长分为两个部分，用前一半的时间步长采用向前差分(displacement)的方法，用后一半时间步长采用向后差分(Bacward)的方式，即将每个方程都变为两个方程的平均值，这样可以将误差降到二次项。由于这个方法同时使用向前差分和向后差分，所以它同时考虑了未知量的前后信息，从而更加准确。 在时间上离散之后，可以得到离散的差分方程组，然后通过求解线性方程组得到数值解。由于Crank-Nicolson差分格式所计算的每一步都是一个算术平均值，因此结果更加平滑，具有良好的稳定性和收敛性。 Crank-Nicolson差分格式的特点在于其数值解更加准确，因为它同时考虑了未知量的前后信息。其缺点是相对于其它差分格式耗时较大，因为必须计算两个半步长的方程。此外，该方法存在计算的矩阵会是个隐式矩阵，所以需要通过专门的技巧求解线性方程组。 总之，Crank-Nicolson差分格式是一种准确稳定的有限差分格式，适用于求解偏微分方程，特别是求解热传导方程、扩散方程等具有稳定解的情况。 ### 回答3： Crank-Nicolson差分格式是一种常见的数值求解偏微分方程的方法，特别是在热传导等问题中。其原理是使用离散化的空间间隔和时间步长，将偏微分方程转化为线性代数方程组的形式，从而可以通过迭代方法求解数值解。 首先，我们将偏微分方程中的空间变量离散化成网格点，时间变量则离散化成时间步长。这样我们就得到了一个网格，其中每个网格点的值都代表了这个时刻的偏微分方程解。我们通过离散化之后，可以得到一个类似于有限差分的方程组。我们将网格点的值表示为$u_{i,j}$，其中$i$表示空间网格点的编号，$j$表示时间步长的编号。因此，我们将偏微分方程的解离散化为了一个二维矩阵。 接下来，我们将偏微分方程中的时间偏导数用向前差分和向后差分的形式离散化，同时用中心差分的形式离散化空间偏导数。由于Crank-Nicolson方法同时采用了前向差分和后向差分，在时间上具有更高的精度。因为它的中心时间步长和二阶时间导数的中心值相等，从而得到的离散方程形式为： $\frac {1}{2 k}(u_{i,j+1}-u_{i,j-1})=\frac{1}{2h^2}(u_{i+1,j}-2u_{i,j}+u_{i-1,j})+f_{i,j}$ 其中$k$表示时间步长，$h$表示空间步长，$f_{i,j}$是源项。 将上式简化可得： $-r u_{i-1,j}+(2+2r)u_{i,j}-r u_{i+1,j} =s_{i,j-1}+s_{i,j}$ 其中： $r=\frac{k}{2 h^{2}}$ $s_{i,j-1}=\frac{u_{i+1,j-1}-2 u_{i,j-1}+u_{i-1,j-1}}{h^{2}}$ $s_{i,j}=f_{i,j}+\frac{u_{i+1,j-1}-2 u_{i,j-1}+u_{i-1,j-1}}{h^{2}}$ 我们使用上述离散方程组来解决偏微分方程的数值解。初始条件可以通过直接给定或者通过求解另一个初始条件的常微分方程得到。对于边界条件，可以使用指定的边界条件，或者通过用空间点表示数值解的导数来自然地得到。此外，由于离散化的时间步长和空间步长的选择对数值结果有很大影响，因此确定合适的时间和空间步长是使用Crank-Nicolson方法的一个重要考虑因素。 总之，Crank-Nicolson差分格式是一种非常有效的数值求解方法，特别是在处理非线性、高阶或多维偏微分方程时。通过将偏微分方程转化为线性方程组的形式，我们可以应用广泛的数值方法来求解。
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

猿何试Bug个踌

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

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

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

打赏作者

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

抵扣说明：

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

余额充值