二维空间的抛物型方程-含时间项的PDE差分法

已于 2024-04-16 10:30:19 修改
阅读量1.6k 收藏 5
点赞数 1
分类专栏： 传统数值方法 文章标签： 矩阵 线性代数
于 2020-09-23 20:35:57 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forrestguang/article/details/108760841
版权
传统数值方法 专栏收录该内容
35 篇文章 18 订阅 ¥199.90 ¥299.90
订阅专栏
本文探讨了二维空间抛物型方程的解法，重点关注使用迭代法来解决含时间项的偏微分方程（PDE）。文章提到了四种不同的差分格式：不稳定且不可用的向前欧拉格式、向后欧拉格式、BDF方法以及CN（Crank-Nicolson）方法。特别地，矩阵迭代在向后欧拉和CN格式中被详细阐述。
摘要由CSDN通过智能技术生成

使用迭代法求解

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

在这里插入图片描述

向前欧拉格式（注意，不稳定，不可用）

import numpy as np
import matplotlib.pyplot as plt
了解本专栏 订阅专栏 解锁全文
Galerkin码农选手
关注
  • 1
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
专栏目录
订阅专栏
微分方程的热传导数值解法及MATLAB实现：向前向后欧拉、C-N、二阶BDF，并比较欧拉、梯形公式和改进欧拉的差异
xOYUDrxns的博客
04-12 754
综上所述，通过对微分方程一维抛物热传导方程的向前向后欧拉法、C-N格式和二阶BDF格式的分析，以及五点差分、九点差分和紧差分格式的比较，我们可以得出关于热传导过程数值计算的详细方法说明。源码可以是MATLAB编写的求解微分方程的代码，读者可以根据提供的源码进行实际计算和分析。在数值计算中，常用的求解一维抛物型热传导方程的方法有向前欧拉法、向后欧拉法、Crank-Nicolson格式和二阶BDF格式。在文章中，可以手动编写微分方程的公式，以及求解过程中所使用的差分格式和数值方法的示意图。
二维抛物方程的初边值问题的数值解法（附Matlab程序，六点对称格式，Du Fort-Frankel 格式，ADI格式，LOD格式）
Sophiayinqianbao的博客
04-21 1861
二维抛物方程的初边值问题的数值解法（附Matlab程序，六点对称格式，Du Fort-Frankel 格式，ADI格式，LOD格式）
ADI(交替方向隐格式)求解二维抛物方程(含matlab程序)
05-17
ADI(交替方向隐格式)求解二维抛物方程(含matlab程序)其中包括了ADI算法的解析步骤和计算例子，最后附上MATLAB程序以供参考
抛物型差分（二维—ADI格式）
whuawell的博客
12-15 3441
抛物型差分（二维—ADI格式） example : 模拟四方边界 ， 上：100°、下：0°、左：75°、右：50°   //getStart.m // simulate function getStart() clc X_INTERVAL = [0 20]; Y_INTERVAL = [0 30]; T = [ 0 10]; deltax = 0.5; delta
偏微分方程数值解—ADI格式求解二维抛物型方程
徐思孟的博客
11-11 8213
二维热传导的ADI格式求解
解二维抛物型方程初边值问题的交替方向隐式方法_周维奎.pdf
01-19
解二维抛物型方程初边值问题的交替方向隐式方法.详细地研究了各种情形的求解二维抛物型方程初边值问题的交替方向隐式差分法。这 种方法可以将二维隐式方法归结为求解三对角线性方程组。与一维情况类似 ,可以继续利用追 赶法求解。它具有运算速度快 ,存储量小 ,无条件稳定等优点。该方法是求解二维抛物型方程的 有效方法 ,必将得到更广泛的应用。
抛物型偏微分方程的有限差分法（3）
01-26
抛物型偏微分方程（PDEs...总的来说，理解并熟练掌握抛物型PDE的有限差分法是数值计算的重要组成部分。选择合适的离散方法和边界处理策略，可以有效求解复杂的物理或工程问题，并为理解和预测现实世界的现象提供基础。
抛物型偏微分方程的有限差分法.zip_偏微分 差分_向后差分格式_向后差分法_抛物型_抛物型 隐
07-15
综上所述，"抛物型偏微分方程的有限差分法.zip" 文件很可能包含了一系列关于如何使用向前差分显格式和向后差分隐格式求解抛物型PDE的教程、代码示例或者分析文档，帮助读者理解和应用这两种方法。这些资源对于学习和...
Matlab实现解抛物型方程求解
07-01
抛物型方程是偏微分方程（PDE）的一种，广泛应用于热传导、波动传播以及流体动力学等领域。掌握在Matlab中解决这类问题的方法对于本科和硕士阶段的科研学习至关重要。 首先，我们来看`adi.m`文件，这是一个实现了 ...
LAB12_EDP:使用 Crank-Nicolson 方法求解抛物线方程-matlab开发
05-31
这种方法属于有限差分法的一种，具有良好的稳定性和精度，尤其适用于时间依赖的问题。下面将详细介绍Crank-Nicolson方法及其在MATLAB中的实现。 一、Crank-Nicolson方法基础 Crank-Nicolson方法是由John Crank和...
解二维抛物型方程的恒稳高精度格式 (1999年)
06-18
建立了解二维抛物型方程的一族含参数的绝对稳定的高精度的差分格式,进而,在特殊情况(θ=0,r=16)下,得到显式差分格式ωn+1=(1+136
二维抛物线方程交替方向隐格式 matlab程序
04-30
简单的二维抛物线方程例子 采用二维交替方向隐格式求解 并且附有matlab程序 适合借鉴
matlab的欧拉方法代码-efficient-PIDE-2D:为二维空间上非结构化网格上的抛物线积分微分方程的数值解提供了简短的MATLAB
05-26
matlab的欧拉方法代码高效PIDE-2D 作者：Gujji Reddy，Alan Seitenfuss，Debora Medeiros，Luca Meacci，MiltonAssunção和Michael Vynnycky 为二维空间上非结构化网格上的抛物线积分微分方程的数值解提供了简短的MATLAB实现。 三角形上的分段线性有限元空间用于空间离散化，而时间离散化则基于后向欧拉法和Crank-Nicolson方法。 选择正交规则以离散化Volterra积分项，以便与时间步长方案一致。 此外，在组装过程中使用矢量化技术介绍了该代码的有效版本，并进行了比较研究。 数值例子证明了该代码的灵活性。 有六个zip文件，其中包含必要的MATLAB文件： BE-LRR-unvectorized.zip （后向Euler，左矩形规则，未向量化）； BE-LRR-vectorized.zip （向后Euler，左矩形规则，矢量化）； BE-RRR-unvectorized.zip （后向Euler，直角矩形规则，未向量化）； BE-RRR-vectorized.zip （向后Euler，直角矩形规则
抛物型方程,抛物型方程的标准形式,matlab源码.zip
10-15
在MATLAB中，解决抛物型方程通常涉及数值方法，如有限差分法、有限元法或谱方法。这些方法将连续问题转化为离散问题，通过代数系统求解近似解。例如，有限差分法会将空间和时间域离散化，用网格点上的函数值来近似原...
【有限差分法】（三）一维和二维抛物方程CN格式以及长时间稳定性分析（附算例与Python代码）
qq_39482771的博客
11-28 8006
前言 我们知道差分里的CN格式是无条件稳定的。但是最近在学习有限元结合CN格式算长时间抛物问题的时候（在时间方向用差分空间方向用有限元）发现稳定性却不能保证，其数值解和真解误差会随着时间越来越大。并且实际能算的时间只有几秒，这几秒钟几乎是没有实际意义的。所以就想到一个问题，纯CN差分格式在长时间计算的时候，是否也是个理论看还行，实践臭弟弟的“花瓶”呢？ 在以前做差分课本上的算例大多也是只有几秒，还真没注意到长时间稳定性的问题(¬_¬ )。所以本文找了最简单的一维抛物问题，把时间层增加至1000S，来看看C.
python pde adi(抛物型差分（二维—ADI格式）)
whuawell的博客
05-02 3140
#coding:utf-8 from mpl_toolkits.mplot3d import axes3d import matplotlib.pyplot as plt import numpy as np import time def createMatrix( m, n): A = np.zeros( (n + 2,m + 2)) Up = np.ones( (m+2,1)) * 1
拋物型方程向前差分matlab,不规则区域上Laplace方程基本解方法的MATLAB实现
weixin_36048512的博客
03-19 841
Advances in Applied MathematicsVol.06No.04(2017), Article ID:21282,6pages10.12677/AAM.2017.64055The Method of Fundamental Solutions for Laplace Equation on Irregular Domain by Using MATLABKang Shan, ...
抛物型差分（二维-1）
whuawell的博客
03-10 1520
// ImpMx.m function ImpMx = getImpMx() clear;clc;   M = 4; N = 3; Up = zeros(1,M); Down = Up; Lf = zeros(1,N); Rt = Lf; rx = 0.1; ry = 0.3; r = rx + ry; ImpMatrix = diag( ones(1, M * N)) * ( 1 + 2 * r
Matlab偏微分方程快速上手：使用pde有限元工具箱求解二维偏微分方程
热门推荐
weixin_47006934的博客
02-01 2万+
Matlab偏微分方程快速上手：使用pdetool工具箱求解二维偏微分方程，适用于数学建模、数学实验，简单的偏微分方程数值计算与工程问题。
二维抛物型方程预-校格式
最新发布
06-18
二维抛物型方程预-校格式（Petrov-Galerkin预处理）通常用于求解偏微分方程（PDEs）中的线性或非线性问题，特别是那些涉及时间依赖和空间复杂性的方程。这种方法结合了有限元素法（FEM）和 galerkin 概念，常用于数值计算中。 在二维抛物型方程的预-校格式中，过程大致分为以下几个步骤： 1. **离散化**：首先，我们将连续的二维空间划分为多个子区域（元素），每个元素内用一个简单的函数（如多项式或正弦函数）近似。这一步通常使用有限元基函数。 2. **弱形式**：将偏微分方程转化为弱形式，这意味着方程在每个元素上乘以一个测试函数并积分，这样可以将问题转换为矩阵形式，更便于计算。 3. **预处理**（Petrov-Galerkin）：这里的关键在于选择合适的测试函数（也称为权重函数）和试函数。预处理包括： - 选择适当的正交性条件，使得测试函数和有限元基函数在每个元素上是正交的。这有助于简化计算，因为这会导致系数矩阵是对称的。 - 应用边界条件，可能需要对有限元的边界进行特定的处理，以保证解决方案的正确性。 4. **列主元方法（Row-reduction）**：为了提高矩阵的稀疏性和计算效率，预处理常常涉及到列主元方法，这可以简化求解过程，特别是在大规模计算中。 5. **迭代求解**：使用预处理后的矩阵和初始条件，应用线性代数算法（如直接求解器或迭代方法）求解系统。 6. **后处理**：得到解之后，可能还需要进行后处理，如插值、平滑或其他可视化操作。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

Galerkin码农选手

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

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

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

打赏作者

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

抵扣说明：

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

余额充值