[原创]偏微分方程数值解法的MATLAB源码【更新完毕】
说明：由于偏微分的程序都比较长，比其他的算法稍复杂一些，所以另开一贴，专门上传偏微分的程序
谢谢大家的支持！
其他的数值算法见：
..//Announce/Announce.asp?BoardID=209&id=8245004
1、古典显式格式求解抛物型偏微分方程（一维热传导方程）
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)
%
%输出参数：U -解矩阵，第一行表示初值，第一列和最后一列表示边值，第二行表示第2层……
% x -空间变量
% t -时间变量
%输入参数：uX -空间变量x的取值上限
% uT -时间变量t的取值上限
% phi -初值条件，定义为内联函数
% psi1 -边值条件，定义为内联函数
% psi2 -边值条件，定义为内联函数
% M -沿x轴的等分区间数
% N -沿t轴的等分区间数
% C -系数，默认情况下C=1
%
%应用举例：
%uX=1;uT=0.2;M=15;N=100;C=1;
%phi=inline('sin(pi*x)');psi1=inline('0');psi2=inline('0');
第十章 偏微分方程数值解法
偏微分方程问题，其求解十分困难。除少数特殊情况外，绝
大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。
§1 差分方法的基本概念
1.1 几类偏微分方程的定解问题
椭圆型方程：其最典型、最简单的形式是泊松（Poisson）方程
),(2222yxfyuxuu
特别地，当0),(yxf时，即为拉普拉斯（Laplace）方程，又称
为调和方程
02222yuxuu
Poisson方程的第一边值问题为
),(),(),(),(),(2222yxyxuyxyxfyuxuyx
其中为以为边界的有界区域，为分段光滑曲线，
称为定解区域，),(yxf，),(yx分别为，上的已知连续函数。
第二类和第三类边界条件可统一表示为
),(),(yxuuyxn
其中n为边界的外法线方向。当0时为第二类边界条件，
0时为第三类边界条件。
抛物型方程：其最简单的形式为一维热传导方程
220(0)uuaatx
方程可以有两种不同类型的定解问题：
初值问题
xxxuxtxuatu)()0,(,0022
初边值问题
221200,0(,0)()0(0,)(),(,)()0uuatTxltxuxxxlutgtultgttT
其中)(x，)(1tg，)(2tg为已知函数，且满足连接条件 )0()(),0()0(21glg
一、 实验目的
1、 了解并掌握不同差分格式的稳定性；
2、 能够掌握比较不同差分格式数值效应的能力。
二、 实验问题
取a=1,2,4, h=0.1, 0.08,用以下几种差分格式求解对流方程
01,0(0,)()0,0txuauxuxfxx
得t=4时数值结果。用图示说明算法的稳定性和间断点附近的计算效果，并进行相应的数值分析。
迎风格式(upwind):
111()nnnnjjjjuuauu
Lax-friedrichs 格式：
1111111()()22nnnnnjjjjjuuuauu
Lax-wendroff 格式：
122111111()(2)22nnnnnnnjjjjjjjuuauuauuu
修正迎风格式:
11(1)nnnjjpjpududu
这里p=[ a],d= a-[ a],/0.8h.为网格比，记号[x]表示不超过x的最大整数。
三、 实验原理
首先取[10,10]x，[0,4]t。按照h=0.1, 0.08划分网格。
再由各种差分格式通过已知的第1层网格点数值可以求出第2层网格点的数值。以此类推，通过逐层的信息最终求得在第4层网格点数值结果。
注意：这里x，t的取值范围应当包含间断点。同时，在所需求的第4层也应当包含间断点。这点要求可以通过初始估算得出。
四、 实验过程
根据要求将网格中x划分为200格，t划分为50格。首先通过matlab构造一个大的零矩阵。然后分别利用差分格式代入并作图。
1、 取a=1
迎风格式(upwind):
Lax-friedrichs 格式：
Lax-wendroff 格式：
修正迎风格式：
2、 取a=2
简析《偏微分方程数值解》教学改革
偏微分方程数值解是本校理学院数学系一门专业课,主要分为有限差分法和有限元法两部分。它为工科专业课程学习和解决实际问题提供必要的数学基础知识及常用的数值解法。随着科技的发展,偏微分方程数值解法在很多领域有着举足轻重的作用,比如在石油、医学以及航天等领域,因而这门课也适用于很多工科专业的学生。但是给数学系的学生讲解这门课,常常使人面临两难的境地:一方面数学专业的学生专业课多数侧重于数学理论,如果把偏微分方程数值解完全上成理论课,倒是迎合了学生一贯的思维模式,但应用方面就会欠缺;如果大量介绍具体应用,由于缺乏相关背景知识,又会使学生迷茫。所以在讲解时,需要把理论和实际合理地结合,即要把理论介绍得深入浅出,又要把某些算法与实际问题联系起来,从而使学生在学习中产生兴趣,是值得思考的问题。另一方面从以往的教学经验来看,偏微分方程数值解的教学存在诸多问题,诸如课时少,内容多;授课一直以板书为主,教学手段比较单一;以往的教学过于偏重理论,学生学过之后不知如何应用;由于课时紧张,而没有安排上机实践等问题。针对上述问题,为了提高教学质量,激发学生学习兴趣,笔者采取了一系列的教改措施,总结如下。
1 课程改革内容
1.1 改革教学内容,选择合适的教材和教辅丛书
关于偏微分方程数值的教材有很多,但一些教材在推导公式时使用的语言晦涩,公式抽象,学生难于理解。在我校的本科生授课时,选择了清华大学陆金甫编写的《偏微分方程数值解》[1]作为教材,以及北京大学出版的李治平编写的《偏微分方程数值解讲义》[2]作为教辅。《偏微分方程数值解》在内容的处理上,体现了由浅入深、循序渐进的原则;在叙述表达上,严谨精练、清晰易读,便于教学与自学。此教材充分反映偏微分方程数值解法中的核心内容。《偏微分方程数值解讲义》的作为教辅的特点是:每章之后配置了相当数量的习题,并在书后附上了大部分习题的答案或提示,更便于学生复习、巩固、理解和拓广所学的知识。
