偏微分方程数值解课程设计
课程设计报告
课程：偏微分方程数值解学号：
姓名：
班级：
教师：
《偏微分方程数值解》
课程设计指导书
一.课程设计的目的
1.帮助掌握偏微分方程数值解相关知识。
2.理解偏微分方程数值解差分隐格式解决自由振动方程问题的方法。
3.锻炼编写程序代码的能力。 二.设计名称
差分法求自由振动问题的周期解。 三.设计要求
1.要求写出差分隐格式的理论方法。
2.要求编写matlab 程序，画出函数图形。
3.要求写出实验总结及心得体会。 四.设计题目
用差分法求自由振动问题的周期解：
2222000,,0|0,|sin (0,)(2,)t t u u
x t t x u u x t u t u t π==???-=-∞<<∞>????
??==???
=???
要求用差分隐格式求解，其中14
θ=
。
五.设计细则 1.区域剖分：
构造上式的差分逼近，取空间步长h 和时间步长τ，用两族平行直线
??
?===±±===
,2,1,0,,
,2,1,0,n n t t j jh x x n j τ 作矩形网格。
2.离散格式：
显格式：
于网点),(n j t x 用Taylor 展式，并整理方程得：
???
??
??--++=+-++==-+-++-121121102
10102100
)1(2)(),()()1()]()([2),(n j n j n j n j n j j j j j j j j u u r u u r u x x r x x r u x u τ?????
隐格式：
上述显格式并不是绝对稳定的差分格式，为了得到绝对稳定的差分格式，用第1-n 层、
n 层、1+n 层的中心差商的权平均去逼近xx u ，得到下列差分格式：
?
???
??
???+-++--++-=+-+-++==----+-++-+++-++-]22)21(2[2),
()()1()]()([2),(2111112112111112
211102
10102100h u u u h u u u h u u u a u u u x x r x x r u x u n j n j n j n j n j n j n j n j n j n j n j n j j j j j j j j θθθττ?????其中10≤≤θ是参数。当0=θ时就是显格式，而当4
1
=θ时可以证明该格式绝对稳定。
隐格式的矩阵形式是：
???
????????
???????????=???????????????????????????????????????????????
?-+-+-+-+--+-+-+++122111121121
12222
222
2222221212121J J j n J n J n j n n z z z z z u u u u u r r r r r r r r r r r r
θθθθ
θθθθθ
θ
θθ
其中：
1
111111122]2()2)(21[(-----+-+-++-++--=n j n j n j n j n j n j n j n j j u u u u u u u u r z θθ
3.格式稳定性： 1）显格式： 显格式稳定的充分必要条件是：网格比1 当4 1 = θ时隐格式绝对稳定。 4.数值例子： 可以证明 ) (),(t x e t x u +=是波动方程222 2u u t x ??=?? 的一个解析解。 那么，为了更精确的得到误差估计，在这里选取) (),(t x e t x u +=作数值实验。 取0201x t π ≤≤≤≤，并且将时间步长10等分，空间步长100等分（即 1/10,2/100h τπ==）。这样网格比/15 r ah π τ== <，从稳定性分析可知，此时格式稳 定。 六.程序代码： 1.主函数： T=1;%取时间长度为1 b=0.5; h=2*pi/100; tao=1/10; f=inline('0','x','t');%f=0 fx1=inline('0'); fx2=inline('sin(x)'); ft1=inline('sin(t)');%此题中取ft1=ft2=sin(t)=sin(2*pi+t) ft2=inline('sin(2*pi+t)'); [X,Y,U,Z]=chfenmethed(f,fx1,fx2,ft1,ft2,T,h,tao); mesh(X,Y,U);%绘制函数图像 shading flat; xlabel('X','FontSize',14); ylabel('t','FontSize',14); zlabel('U','FontSize',14); title('函数'); figure; mesh(X,Y,Z);%绘制误差图像 xlabel('X','FontSize',14); ylabel('t','FontSize',14); zlabel('Z','FontSize',14); title('误差'); 2.编写差分隐格式函数： function [X,T,U,Z]=chfenmethed(f,fx1,fx2,ft1,ft2,T,h,tao) %u_tt-u_xx=f(x,t) 0 %u(0,t)=ft1,u(2*pi,t)=ft2,此题中ft1=ft2 %u(x,0)=fx1(x),此题中fx1=0 %u_t(x,0)=fx2(x),此题中fx2=sin(x) x=0:h:2*pi; t=0:tao:T; m=length(x); n=length(t); s=tao/h; [X,T]=meshgrid(x,t); Z=zeros(n,m); U=zeros(n,m); for i=2:m-1 U(1,i)=feval(fx1,x(i)); U(2,i)=U(1,i)+tao*feval(fx2,x(i))+tao^2/2*(1/h^2* ... (feval(fx1,x(i+1))-2*feval(fx1,x(i))+feval(fx1,x(i-1))+feval(f,x(i),0))); Z(2,i)=abs(U(2,i)-f0(x(i),t(2))); end for j=1:n U(j,1)=feval(ft1,t(j)); U(j,m)=feval(ft2,t(j)); end A=-0.5*s^2*ones(1,m-2); C=A; B=(1+s^2)*ones(1,m-2); UU=zeros(1,m-2); f1=UU; for i=3:n for j=2:m-1 UU(j-1)=f0(x(j),t(i)); f1(j-1)=0.5*s^2*U(i-2,j-1)-(1+s^2)*U(i-2,j)... +0.5*s^2*U(i-2,j+1)+2*U(i-1,j)... +tao^2*feval(f,x(j),t(i-1)); end f1(1)=f1(1)+0.5*s^2*U(i,1); f1(end)=f1(end)+0.5*s^2*U(i,m); U(i,2:m-1)=zgf(A,B,C,f1); Z(i,2:m-1)=abs(U(i,2:m-1)-UU); end 3.编写迭代函数，用以实现差分隐格式： function x=zgf(A,B,C,f) n=length(B); B1=zeros(1,n-1); Y=zeros(1,n); x1=zeros(1,n); B1(1)=C(1)/B(1); for i=2:n-1 B1(i)=C(i)/(B(i)-A(i)*B1(i-1)); end Y(1)=f(1)/B(1); for i=2:n Y(i)=(f(i)-A(i)*Y(i-1))/(B(i)-A(i)*B1(i-1)); end x1(n)=Y(n); for i=n-1:-1:1 x1(i)=Y(i)-B1(i)*x1(i+1); end x=x1; 4.编写精确函数用以求误差： function z=f0(x,t)%精确解函数 z=exp(x+t); 七.设计结果与分析： 最终得到的两张图像如下： 八.设计体会与建议： 设计成绩：教师签名： 年月日 偏微分方程数值解试题(０6B ） 参考答案与评分标准 信息与计算科学专业 一(10分）、设矩阵A 对称，定义)(),(),(2 1)(n R x x b x Ax x J ∈-=，)()(0x x J λλ?+=.若0)0('=?，则称称0x 是)(x J 的驻点(或稳定点).矩阵A 对称(不必正定)，求证0x 是)(x J 的驻点的充要条件是:0x 是方程组 b Ax =的解 解: 设n R x ∈0是)(x J 的驻点，对于任意的n R x ∈,令 ),(2),()()()(2 000x Ax x b Ax x J x x J λλλλ?+-+=+=, （3分) 0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0，则有0||||),(2000=-=--b Ax b Ax b Ax ，得到b Ax =0． (3分) 反之,若n R x ∈0满足b Ax =0,则对于任意的x ,)(),(2 1)0()1()(00x J x Ax x x J >+==+??,因此0x 是)(x J 的最小值点. (４分） 评分标准：)(λ?的展开式3分, 每问3分，推理逻辑性１分 二（10分）、 对于两点边值问题:?????==∈=+-=0 )(,0)(),()('b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ],[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的Ritz 形式和Galerkin 形式的变分方程。 解: 设}0)(),,(|{11=∈=a u b a H u u H E 为求解函数空间,检验函数空间．取),(1b a H v E ∈,乘方程两端,积分应用分部积分得到 (３分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(1 b a H v E ∈? 即变分问题的Galerkin 形式. （3分) 一．填空（1553=?分） 1．若步长趋于零时，差分方程的截断误差0→lm R ，则差分方程的解lm U 趋近于微分方 程的解lm u . 此结论_______（错或对）； 2．一阶Sobolev 空间{} )(,,),()(21 Ω∈''=ΩL f f f y x f H y x 关于内积=1),( g f _____________________是Hilbert 空间； 3．对非线性（变系数）差分格式，常用 _______系数法讨论差分格式的_______稳定性； 4．写出3 x y =在区间]2,1[上的两个一阶广义导数：_________________________________， ________________________________________； 5．隐式差分格式关于初值是无条件稳定的. 此结论_______（错或对）。 二．（13分）设有椭圆型方程边值问题 用1.0=h 作正方形网格剖分 。 （1）用五点菱形差分格式将微分方程在内点离散化； （2）用截断误差为)(2 h O 的差分法将第三边界条件离散化； （3）整理后的差分方程组为 三．（12）给定初值问题 x u t u ??=?? ， ()10,+=x x u 取时间步长1.0=τ，空间步长2.0=h 。试合理选用一阶偏心差分格式（最简显格式）, 并以此格式求出解函数),(t x u 在2.0,2.0=-=t x 处的近似值。 1．所选用的差分格式是： 2．计算所求近似值： 四．（12分）试讨论差分方程 ()h a h a r u u r u u k l k l k l k l ττ + - = -+=++++11,111 1 逼近微分方程 0=??+??x u a t u 的截断误差阶R 。 思路一：将r 带入到原式，展开后可得格式是在点（l+1/2,k+1/2）展开的。 思路二：差分格式的用到的四个点刚好是矩形区域的四个顶点，可由此构造中心点的差分格 式。 偏微分方程数值解 偏微分方程地构建科学、工程学和其他领域的数学模型的主要手段。一般情况下,这些模型都需要用数值方法去求解。本书提供了标准数值技术的简明介绍。借助抛物线型、双曲线型和椭圆型方程的一些简单例子介绍了常用的有限差分方法、有限元方法、有限体方法、修正方程分析、辛积分格式、对流扩散问题、多重网络、共轭梯度法。利用极大值原理、能量法和离散傅里叶分析清晰严格地处理了稳定性问题。本书全面讨论了这些方法的性质,并附有典型的图像结果,提供了不同难度的例子和练习。 本书可作为数学、工程学及计算机科学专业本科教材,也可供工程技术人员和应用工作者参考。 偏微分方程数值解---学习总结(2) 关于SobolveSobolve空间的几个重要定理 迹定理 : ΩΩ是 RdRd 的一个有界开子集，具有李普希茨连续边界?Ω?Ω, s>12s>12, 则 a.存在唯一的连续线性映射γ0:Hs(Ω)→Hs?12(?Ω),满足γ0v=v ∣∣?Ω,?v∈Hs(Ω)∩C0(Ωˉˉˉˉ), b.存在唯一的连续映射R0：Hs?12(?Ω)→Hs(Ω),满足γ0°R0°φ=φ,?φ∈Hs?12(?Ω).(1)(2)(1)a.存在唯一的连续线性映射γ0:Hs(Ω)→Hs?12(?Ω),满足γ0v=v|?Ω,?v∈ Hs(Ω)∩C0(Ωˉ),(2)b.存在唯一的连续映射R0：Hs?12(?Ω)→Hs(Ω),满足γ0°R0°φ=φ,?φ∈Hs?12(?Ω). 迹定理把区域内部与边界联系起来. 上面定理中边界?Ω?Ω当被它的一个子集ΣΣ代替时，结论依然成立. S=1时， γ0:H1(Ω)→H12(?Ω)?L2(?Ω)||γ0v||0,?Ω≤||γ0v||2,?Ω≤C||v||1=C(||v||0+||?v||0).γ0:H1(Ω)→H12(?Ω)? L2(?Ω)||γ0v||0,?Ω≤||γ0v||2,?Ω≤C||v||1=C(||v||0+||? v||0). 注意几个范数 ||?||k||?||0||?||1||??||0=||?||k,2=||?||L2=||?||1,2=(||?||20+||??||20)12=|?|1.(3)(4)(5)(6)(3)||?||k=||?||k,2(4)||? ||0=||?||L2(5)||?||1=||?||1,2=(||?||02+||??||02)12(6)||?? ||0=|?|1. 庞加莱不等式(Poincare inequality): 假设ΩΩ是 RdRd 的一个有界联通开子集，ΣΣ是边界?Ω?Ω的一个非空的李普希茨连续子集. 则存在一个常数 CΩ>0CΩ>0满足 ∫Ωv2(x)dx≤CΩ∫Ω|?v(x)|2dx,?v∈H1Σ(Ω),其中H1Σ(Ω)={v ∈H1(Ω),γΣv=v∣∣Σ=0}.∫Ωv2(x)dx≤CΩ∫Ω|?v(x)|2dx,?v∈HΣ1(Ω),其中HΣ1(Ω)={v∈H1(Ω),γΣv=v|Σ=0}. 偏微分方程数值解试题(06B) 参考答案与评分标准 信息与计算科学专业 一（10分）、设矩阵A 对称，定义)(),(),(2 1 )(n R x x b x Ax x J ∈-= ，)()(0x x J λλ?+=.若0)0('=?,则称称0x 是)(x J 的驻点（或稳定点）.矩阵A 对称（不必正定），求证0x 是)(x J 的驻点的充要条件是：0x 是方程组 b Ax =的解 解: 设n R x ∈0是)(x J 的驻点,对于任意的n R x ∈,令 ),(2 ),()()()(2 000x Ax x b Ax x J x x J λλλλ?+ -+=+=, (3分) 0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有 0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若 n R x ∈0满足 b Ax =0,则对于任意的 x ,)(),(2 1 )0()1()(00x J x Ax x x J >+ ==+??,因此0x 是)(x J 的最小值点. (4分) 评分标准:)(λ?的展开式3分, 每问3分,推理逻辑性1分 二（10分）、 对于两点边值问题：????? ==∈=+-=0 )(,0)() ,()(' b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ] ,[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的Ritz 形式和 Galerkin 形式的变分方程。 解: 设}0)(),,(|{11 =∈=a u b a H u u H E 为求解函数空间,检验函数空间.取),(1 b a H v E ∈,乘方程两端,积分应用分部积分得到 (3分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(1 b a H v E ∈? 即变分问题的Galerkin 形式. (3分) 令?-+=-=b a dx fu qu dx du p u f u u a u J ])([21),(),(21)(22,则变分问题的Ritz 形式 偏微分方程数值解实验报告 1、用有限元方法求下列边值问题的数值解：''()112x -y +y =2s i n ,0 1、 %Ritz Galerkin方法求解方程 function u1=Ritz(x) %定义步长 h=1/100; x=0:h:1; n=1/h; a=zeros(n-1,1); b=zeros(n,1); c=zeros(n-1,1); d=zeros(n,1); %求解Ritz方法中内点系数矩阵 for i=1:1:n-1 b(i)=(1/h+h*pi*pi/12)*2; d(i)=h*pi*pi/2*sin(pi/2*(x(i)+h))/2+h*pi*pi/2*sin(pi/2*x(i+1))/2; end %右侧导数条件边界点的计算 b(n)=(1/h+h*pi*pi/12); d(n)=h*pi*pi/2*sin(pi/2*(x(i)+h))/2; for i=1:1:n-1 a(i)=-1/h+h*pi*pi/24; c(i)=-1/h+h*pi*pi/24; end %调用追赶法 u=yy(a,b,c,d) %得到数值解向量 u1=[0,u] %对分段区间做图 plot(x,u1) %得到解析解 y1=sin(pi/2*x); hold on plot(x,y1,'o') legend('数值解','解析解') function x=yy(a,b,c,d) n=length(b); q=zeros(n,1); p=zeros(n,1); q(1)=b(1); p(1)=d(1); for i=2:1:n 偏微分方程数值解期末复习（2011硕士） 一、考题类型 本次试卷共六道题目，题型及其所占比例分别为： 填空题20%；计算题80% 二、按章节复习内容 第一章 知识点：Euler法、向前差商、向后差商、中心差商、局部截断误差、整体截断误差、相容性、收敛性、阶、稳定性、显格式、隐格式、线性多步法、第一特征多项式、第二特征多项式、稳定多项式、绝对稳定等； 要求: 会辨认差分格式, 判断线性多步法的误差和阶; 第二章 知识点：矩形网格、(正则,非正则)内点、边界点、偏向前(向后,中心)差商、五点差分格式、增设虚点法、积分插值法、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和、稳定性等； 要求: 建立椭圆型方程边值问题的差分格式, 极值原理; 第四章 知识点：最简显格式、最简隐格式、CN格式、双层加权格式、Richardson 格式、网格比、传播因子法(分离变量法) 、传播因子、传播矩阵、谱半径、von Neumann条件、跳点格式、ADI格式、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和稳定性等； 要求: 建立抛物型方程边值问题的差分格式, 计算局部截断误差; 第五章 知识点：左偏心格式、右偏心格式、中心格式、LF格式、LW格式、Wendroff 格式、跳蛙格式、特征线、CFL条件等； 要求: 建立双曲型方程边值问题的差分格式, 计算局部截断误差; 第七章 要求: 会用线性元(线性基)建立常微分方程边值问题的有限元格式 三 练习题 1、 已知显格式21131()22 n n n n u u h f f +++-=-，试证明格式是相容的，并求它的阶。 P39+P41 2、用Taylor 展开原理构造一元函数一阶导数和二阶导数的数值微分公式。 提示：向前、向后和中心差商与一阶导数间关系，二阶中心差商与二阶导数 之间的关系 课件 3、用数值微分方法或数值积分方法建立椭圆型方程 2222(,),(,),u u f x y x y x y ??--=?∈Ω?? :01,01x y Ω≤≤≤≤ 内点差分格式。 P75+课件 4、构造椭圆型方程边值问题的差分格式. P101 (4)题 5、构建一维热传导方程220,(0)u u Lu a a t x ??=-=>??的数值差分格式(显隐格式等)。 参考P132-135相关知识点 6、设有逼近热传导方程22(0)u u Lu a f a const t x ??≡-==>??的带权双层格式 ()()1111111122(1)2k k j j k k k k k k j j j j j j u u a u u u u u u h θθτ++++-+-+-??=-++--+?? 其中[0,1]θ∈，试求其截断误差。并证明当2 1212h a θτ=-时，截断误差的阶最 高阶为24()O h τ+。 P135+P165+课件 7、传播因子法证明抛物型方程22(0)u u Lu a f a const t x ??≡-==>??的最简显隐和六点CN 格式稳定性。 P156+课件 8、对一阶常系数双曲型方程的初边值问题 0,0,0,0,(,0)(),0,(0,)(),0, u u a t T x a t x u x x x u t t t T φψ???+=<≤<<∞>?????=≤<∞??=≤≤? 专业班级 姓名 学号 开课系室数学与计算科学学院 考试日期 偏微分方程数值解试卷 一（15分）、（1）简述用差分方法求解抛物型方程初边值问题的数值解的一般步骤.(2)写出近似一阶偏导数 n m x u |??的三种有限差分逼近及其误差阶,写出近似 n m x u |22 ??的差分逼近及其误差阶. 评分标准： （1） 7分，三个离散4分，其他步骤3分 （2） 8分，每个格式及误差2分。 二（15分）、（1）以抛物型方程的差分格式为例，解释差分格式的相容性,稳定性和收敛性概念,分析相容性,稳定性和收敛性与误差的关系，简述 Lax 等价性定理。(2) 简述差分格式稳定性分析的Fourier 级数法（或称为Neumann Von 方法，分离变量法）的一般步骤。 （1）8分，解释概念6分，等价关系2分 （2）7分，典型波2分，放大因子与条件3分，其他2分 三（20分）、对于边值问题 ?? ???=?=∈=??+???0 |) 1,0()1,0(),(,92 222G u G y x y u x u （1）建立该边值问题的五点差分格式（五点棱形格式又称正五点格式），推导截 断误差的阶。 （2）取3/1=h ，求边值问题的数值解（写出对应的方程组的矩阵形式并求解） （3）就取5/1=h 的情况写出对应方程组的系数矩阵（用分块矩阵表示）。 解：（1）7分，离过程与格式 第二页（共五页） 四（20分）、对于初边值问题??? ????≤≤==<<=≤<< 偏微分方程数值解试题 1、考虑一维的抛物型方程： 2200, [0,], 0t T (,), (,)(,0)() x x u u x t x u x t u u x t u u x x ππνπ?==??=∈≤≤??=== （1）导出时间离散是一阶向前Euler 格式，空间离散是二阶精度的差分格式； （2）讨论（1）中导出的格式的稳定性； （3）若时间离散为二阶精度的蛙跳格式， 11 2n n n t t u u u t t +-=?-= ?? 空间离散是二阶精度的中心差分，问所导出的格式稳定吗？为什么？ 2、考虑Poission 方程 2(,)1, (,)0, in AB and AD (,)0, in BC and CD u x y x y u n u x y -?=∈Ω ?=?= 其中Ω是图1中的梯形。 使用差分方法来离散该方程。由于梯形的对称性，可以考虑梯形的一半，如图2， 图2 从物理空间到计算区域的几何变换 图1 梯形 为了求解本问题，采用如下方法：将Ω的一半投影到正方形区域?Ω ，然后在?Ω上使用差分方法来离散该方程。在计算区域?Ω 上用N N ?个网格点，空间步长为1/(1)N ξη?=?=-。 （1）引入一个映射T 将原区域Ω（带有坐标,x y ）变换到单位正方形?Ω（带有坐标,ξη）。 同时导出在新区域上的方程和边界条件。 （2）在变换区域，使用泰勒展开导出各导数项在区域内部和边界点上的差分格式。 3、对线性对流方程0 constant >0u u a a t x ??+=??，其一阶迎风有限体积法离散格式为 1?n j u +=?n j u a t x ?-?（?n j u 1?n j u --） （1）写出0a <时的一阶迎风有限体积法的离散格式； （2）写出a 为任意符号的常数的一阶迎风有限体积法的守恒形式。 （3）使用0 u u u t x ??+=??说明一阶迎风有限体积法不是熵保持的格式。 4、对一维Poission 方程 , (0,1) (0)(1)0 x xx u xe x u u ?-=∈? ==? 将[]01，分成(1)n +等分，写出用中心差分离散上述方程的差分格式，并问： （1）该差分格式与原微分方程相容吗？为什么？ （2）该差分格式稳定吗？为什么？ （3）该差分格式是否收敛到原微分方程的解？为什么？ （4）取(1)6n +=，写出该差分格式的矩阵表示。 5、叙述二重网格方法的执行过程，并对一维常微分方程边值问题 2 25, (0,1) (0)(1)0 xx u x x x u u πππ?-=∈? ==?（sin(5)+9sin(15)） 给出限制算子和延拓算子矩阵（以细网格h ：7n =，粗网格2h ：3n =为例）。 6、对一阶波动方程 01(,0)sin(), (0,1)2(0,)(1,)u u t x u x x x u t u t π???+=???? ? =∈?? =??? （1）写出用中心差分进行空间离散，用一阶向后Euler 进行时间离散的差分格式； 1 / 7 偏微分方程数值解试题 1、考虑一维的抛物型方程： 2200, [0, ], 0t T (,), (,)(,0)() x x u u x t x u x t u u x t u u x x ππνπ?==??=∈≤≤??=== （1）导出时间离散是一阶向前Euler 格式，空间离散是二阶精度的差分格式； （2）讨论（1）中导出的格式的稳定性； （3）若时间离散为二阶精度的蛙跳格式， 11 2n n n t t u u u t t +-=?-=?? 空间离散是二阶精度的中心差分，问所导出的格式稳定吗？为什么？ 2、考虑Poission 方程 2(,)1, (,)0, in AB and AD (,)0, in BC and CD u x y x y u n u x y -?=∈Ω ?=?= 其中Ω是图1中的梯形。 使用差分方法来离散该方程。由于梯形的对称性，可以考虑梯形的一半，如图2， 图2 从物理空间到计算区域的几何变换 图1 梯形 2 / 7 为了求解本问题，采用如下方法：将Ω的一半投影到正方形区域?Ω ，然后在?Ω上使用差分方法来离散该方程。在计算区域?Ω 上用N N ?个网格点，空间步长为1/(1)N ξη?=?=-。 （1）引入一个映射T 将原区域Ω（带有坐标,x y ）变换到单位正方形?Ω（带有坐标,ξη）。 同时导出在新区域上的方程和边界条件。 （2）在变换区域，使用泰勒展开导出各导数项在区域内部和边界点上的差分格式。 3、对线性对流方程 0 constant >0u u a a t x ??+=??，其一阶迎风有限体积法离散格式为 1?n j u +=?n j u a t x ?-?（?n j u 1?n j u --） （1）写出0a <时的一阶迎风有限体积法的离散格式； （2）写出a 为任意符号的常数的一阶迎风有限体积法的守恒形式。 （3）使用0 u u u t x ??+=??说明一阶迎风有限体积法不是熵保持的格式。 4、对一维Poission 方程 , (0,1)(0)(1)0 x xx u xe x u u ?-=∈?==? 将[]01，分成(1)n +等分，写出用中心差分离散上述方程的差分格式，并问： （1）该差分格式与原微分方程相容吗？为什么？ （2）该差分格式稳定吗？为什么？ （3）该差分格式是否收敛到原微分方程的解？为什么？ （4）取(1)6n +=，写出该差分格式的矩阵表示。 5、叙述二重网格方法的执行过程，并对一维常微分方程边值问题 225, (0,1)(0)(1)0 xx u x x x u u πππ?-=∈?==?（sin(5)+9sin(15)） 给出限制算子和延拓算子矩阵（以细网格h ：7n =，粗网格2h ：3n =为例）。 6、对一阶波动方程 1(,0)sin(), (0,1)2(0,)(1,)u u t x u x x x u t u t π???+=?????=∈??=??? （1）写出用中心差分进行空间离散，用一阶向后Euler 进行时间离散的差分格式； 一、 问题 用有限元方法求下面方程的数值解 2 u u u f t ?-?+=? in (]0,T Ω? 0u = on []0,T ?Ω? ()00,u x u = in Ω 二、 问题分析 第一步 利用Green 公式，求出方程的变分形式 变分形式为：求()()21 00,;u L T H ∈Ω，使得 ()())(2 ,,,,u v u v u v f v t ???+??+= ???? ()10v H ?∈Ω （*） 以及 ()00,u x u =. 第二步 对空间进行离散，得出半离散格式 对区域Ω进行剖分，构造节点基函数，得出有限元子空间：()12,,,h NG V span ???=???，则（*）的Galerkin 逼近为： []0,t T ?∈，求()()1 0,h h u t x V H ∈?Ω，使得 ()()()()() () )(2 ,,,,h h h h h h h d u t v u t v u t v f v dt +??+= h h v V ?∈ （**） 以及()0,0h h u u =，0,h u 为初始条件0u 在h V 中的逼近，设0,h u 为0u 在h V 中的插值. 则0t ?≥，有()()1 N G h i i i u t t ξ? == ∑，0,h u =01 N G i i i ξ?=∑，代人（**）即可得到一常微分方程组. 第三步 进一步对时间进行离散，得到全离散的逼近格式 对 du dt 用差分格式.为此把[]0,T 等分为n 个小区间[]1,i i t t -，其长度1i i T t t t n -?=-= ,n t T =. 这样把求i t 时刻的近似记为i h u ，0 h u 是0u 的近似.这里对（**）采用向后的欧拉格式，即 ()()() () )(2 11 11 1 ,,,,i i i i h h h h h h h i h u u v u v u v f v t ++++-+??+ = ? h h v V ?∈ (***) i=0,1,2…,n-1. 0 h u =0,h u 由于向后欧拉格式为隐式格式且含有非线性项，故相邻两时间步之间采用牛顿迭代，即： 偏微分方程数值解 一（10分）、设矩阵A 对称正定，定义)(),(),(2 1 )(n R x x b x Ax x J ∈-= ，证明下列两个问题等价：(1)求n R x ∈0使)(min )(0x J x J n R x ∈=;(2)求下列方程组的解：b Ax = 解: 设n R x ∈0是)(x J 的最小值点,对于任意的n R x ∈,令 ),(2 ),()()()(2 000x Ax x b Ax x J x x J λλλλ?+ -+=+=, (3分) 因此0=λ是)(λ?的极小值点,0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若 n R x ∈0满足 b Ax =0,则对于任意的x , )(),(2 1 )0()1()(00x J x Ax x x J >+ ==+??,因此0x 是)(x J 的最小值点. (4分) 评分标准:)(λ?的表示式3分, 每问3分,推理逻辑性1分 二（10分）、对于两点边值问题：????? ==∈=+-=0 )(,0)() ,()(b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ] ,[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的Ritz 形式和 Galerkin 形式的变分方程。 解: 设}0)()(),,(|{11 0==∈=b u a u b a H u u H 为求解函数空间,检验函数空间.取),(10b a H v ∈,乘方程两端,积分应用分部积分得到 (3分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(10 b a H v ∈? 即变分问题的Galerkin 形式. (3分) 二、改进的Euler 方法 梯形方法的迭代公式(1.10)比Euler 方法精度高,但其计算较复杂,在应用公式(1.10)进行计算时,每迭代一次,都要重新计算函数),(y x f 的值,且还要判断何时可以终止或转下一步计算.为了控制计算量和简化计算法,通常只迭代一次就转入下一步计算.具体地说,我们先用Euler 公式求得一个初步的近似值1+n y ,称之为预测值,然后用公式(1.10)作一次迭代得1+n y ,即将1+n y 校正一次.这样建立的预测－校正方法称为改进的Euler 方法: 预测: ),,(1n n n n y x hf y y +=+ 校正 : )].,(),([2 111+++++=n n n n n n y x f y x f h y y (1.15) 这个计算公式也可以表示为 11(,), (,), 1(). 2p n n n c n n p n p c y y hf x y y y hf x y y y y ++?=+??=+?? ?=+??? 例1 取步长0.1h =，分别用Euler 方法及改进的Euler 方法求解初值问题 d (1),01, d (0) 1. y y xy x x y ?=-+≤≤???=? 解 这个初值问题的准确解为()1(21)x y x e x =--. 根据题设知 ).1(),(xy y y x f +-= (1) Euler 方法的计算式为 )],1([1.01n n n n n y x y y y +?-=+ 由1)0(0==y y , 得 ,9.0)]101(1[1.011=?+??-=y ,8019.0)]9.01.01(9.0[1.09.02=?+??-=y 这样继续计算下去，其结果列于表9.1. (2) 改进的Euler 方法的计算式为 110.1[(1)],0.1[(1)], 1(), 2p n n n n c n p n p n p c y y y x y y y y x y y y y ++?=-?+?=-?+??? ?=+??? 由1)0(0==y y ，得 x ∈R n 2 ( Ax, x) , J ( x + x) = ? (1) = ? (0) + ( Ax, x) > J ( x ) ,因此 x 是 J ( x ) 的最小值点. (4 分) 2 二（10 分）、对于两点边值问题： ? dx dx a(u , v) = ?b ( p . + q u v)dx = ?b fvdx = f (v) , ? v ∈ H 1 (a , b ) dx dx a a 偏微分方程数值解 一（10 分）、设矩阵 A 对称正定，定义 J ( x ) = 1 ( Ax , x ) - (b , x ) ( x ∈ R n ) ，证明下 2 列两个问题等价：(1)求 x ∈ R n 使 J ( x ) = min J ( x ) ;(2)求下列方程组的解：Ax = b 解: 设 x ∈ R n 是 J ( x ) 的最小值点,对于任意的 x ∈ R n ,令 ?(λ) = J ( x + λx) = J ( x ) + λ( Ax - b , x) + λ2 (3 分) 因此 λ = 0 是 ?(λ) 的极小值点 , ? ' (0) = 0 ,即对于任意的 x ∈ R n , ( Ax - b , x) = 0 ,特 0 别取 x = Ax - b ,则有 ( Ax - b , Ax - b ) =|| Ax - b || 2 = 0 ,得到 Ax = b . (3 分) 0 0 反 之 , 若 x ∈ R n 满 足 Ax = b , 则 对 于 任 意 的 x , 1 0 0 0 评分标准: ?(λ) 的表示式 3 分, 每问 3 分,推理逻辑性 1 分 ? d du ?Lu = - ( p ) + qu = f x ∈ (a, b ) ?? u (a) = 0, u (b ) = 0 其中 p ∈ C 1 ([a , b ]), p ( x ) ≥ min p ( x ) = p x ∈[a,b ] min > 0, q ∈ C ([a , b ]), q ≥ 0, f ∈ H 0 ([a , b ]) 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的 Ritz 形式和 Galerkin 形式的变分方程。 解 : 设 H 1 = {u | u ∈ H 1 (a , b ), u (a ) = u (b ) = 0} 为求解函数空间 , 检验函数空间 . 取 v ∈ H 1 (a, b ) ,乘方程两端,积分应用分部积分得到 (3 分) du dv 即变分问题的 Galerkin 形式. (3 分) 第十章 偏微分方程数值解法 偏微分方程问题，其求解十分困难。除少数特殊情况外，绝 大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。 §1 差分方法的基本概念 1.1 几类偏微分方程的定解问题 椭圆型方程：其最典型、最简单的形式是泊松（Poisson ）方程 ),(22 2 2y x f y u x u u =??+??=? 特别地，当 0),(≡y x f 时，即为拉普拉斯（Laplace ）方程，又称 为调和方程 22 22 =??+??=?y u x u u Poisson 方程的第一边值问题为 ?? ?? ?Ω ?=Γ=Ω∈=??+??Γ∈),(),(),(),(),(22 22y x y x u y x y x f y u x u y x ? 其中 Ω为以Γ为边界的有界区域，Γ为分段光滑曲线， ΓΩY 称为定解区域，),(y x f ，),(y x ?分别为Ω，Γ上的已知连 续函数。 第二类和第三类边界条件可统一表示为 ),(),(y x u u y x ?α=??? ? ??+??Γ∈n 其中n 为边界Γ的外法线方向。当0=α时为第二类边界条件， 0≠α时为第三类边界条件。 抛物型方程：其最简单的形式为一维热传导方程 2 20(0)u u a a t x ??-=>?? 方程可以有两种不同类型的定解问题： 初值问题 ?? ???+∞ <<∞-=+∞<<-∞>=??-??x x x u x t x u a t u )()0,(,00 22 ? 初边值问题 2 212 00,0(,0)()0(0,)(),(,)()0u u a t T x l t x u x x x l u t g t u l t g t t T ????-=<<< 李微分方程数值解习题解答 1-1 如果0)0(' =?，则称0x 是)(x J 的 驻点（或稳定点）.矩阵A 对称（不必正定），求证0x 是)(x J 的驻点的充要条件是：0x 是方程组 b Ax =的解 证明:由)(λ?的定义与内积的性线性性质,得 ),()),((2 1 )()(0000x x b x x x x A x x J λλλλλ?+-++=+= ),(2 ),()(2 00x Ax x b Ax x J λλ+ -+= ),(),()(0'x Ax x b Ax λλ?+-= 必要性:由0)0(' =?,得,对于任何n R x ∈,有 0),(0=-x b Ax , 由线性代数结论知, b Ax b Ax ==-00,0 充分性: 由b Ax =0,对于任何n R x ∈, 0|),(),()0(00'=+-==λλ?x Ax x b Ax 即0x 是)(x J 的驻点. §1-2 补充: 证明)(x f 的不同的广义导数几乎处处相等. 证明:设)(2I L f ∈,)(,221I L g g ∈为)(x f 的广义导数,由广义导数的定义可知,对于任意 )()(0I C x ∞∈?,有 ??-=b a b a dx x x f dx x x g )()()()(' 1?? ??-=b a b a dx x x f dx x x g )()()()('2?? 两式相减,得到 )(0)()(021I C x g g b a ∞ ∈?=-??? 由变分基本引理,21g g -几乎处处为零,即21,g g 几乎处处相等. 补充:证明),(v u a 的连续性条件(1.2.21) 证明: 设'|)(|,|)(|M x q M x p ≤≤,由Schwarz 不等式 “十二五”国家重点图书出版规划项目 信息与计算科学丛书 67 偏微分方程数值解法 陈艳萍鲁祖亮刘利斌编著 内 容 简 介 本书试图用较少的篇幅描述偏微分方程的几种数值方法. 主要内容包括：Sobolev空间初步, 椭圆边值问题的变分问题, 椭圆问题的有限差分方法, 抛物型方程的有限差分方法, 双曲型方程的有限差分方法, 椭圆型方程的有限元方法, 抛物及双曲方程的有限元方法, 椭圆型方程的混合有限元方法, 谱方法等. 本书内容丰富, 深入浅出, 尽可能地用简单的方法来描述一些理论结果, 并根据作者对有限差分、有限元、混合有限元、谱方法的理解和研究生教学要求, 全面、客观地评价各种数值计算方法，并列举一些数值计算的例子, 阐述许多新的学术观点. 本书可作为高等学校数学系高年级本科生和研究生的教材或参考书, 也可作为计算数学工作者和从事科学与工程计算的科研人员的参考书. 图书在版编目(CIP)数据 偏微分方程数值解法/陈艳萍, 鲁祖亮, 刘利斌编著. —北京：科学出版社, 2015.1 (信息与计算科学丛书67) ISBN 978-7-03-000000-0 Ⅰ. ①偏… Ⅱ. ①陈… ②鲁… ③刘… Ⅲ. ① Ⅳ.① 中国版本图书馆CIP数据核字(2014) 第000000号 责任编辑: 王丽平／责任校对: 彭涛 责任印制: 肖钦／封面设计: 陈敬 出版 北京东黄城根北街16号 邮政编码: 100717 https://www.wendangku.net/doc/0f17146662.html, 印刷 科学出版社发行 各地新华书店经销 * 2015年1月第一版开本: 720×1000 1/16 2015年1月第一次印刷印张: 14 字数: 280 000 定价: 88.00元 (如有印装质量问题, 我社负责调换) 1. 课本2p 有证明 2. 课本812,p p 有说明 3. 课本1520,p p 有说明 4. Rit2法，设n u 是u 的n 维子空间，12,...n ???是n u 的一组基底，n u 中的任一元素n u 可 表为1n n i i i u c ?==∑ ，则,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???=== -=-∑∑是12,...n c c c 的二次函数，(,)(,)i j j i a a ????=，令 () 0n j J u c ?=?，从而得到12,...n c c c 满足1 (,)(,),1,2...n i j i j i a c f j n ???===∑，通过解线性方程组，求的i c ，代入1 n n i i i u c ?==∑， 从而得到近似解n u 的过程称为Rit2法 简而言之，Rit2法：为得到偏微分方程的有穷维解，构造了一个近似解，1 n n i i i u c ?== ∑， 利用,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???===-=-∑∑确定i c ，求得近似解n u 的过程 Galerkin 法：为求得1 n n i i i u c ? == ∑形式的近似解，在系数i c 使n u 关于n V u ∈，满足(,)(,) n a u V f V =，对任 意 n V u ∈或（取 ,1j V j n ?=≤≤） 1 (,)(,),1,2...n i j i j i a c f j n ???===∑的情况下确定i c ，从而得到近似解1 n n i i i u c ?==∑的过程称 Galerkin 法为 Rit2-Galerkin 法方程： 1 (,)(,)n i j i j i a c f ???==∑ 5. 有限元法：将偏微分方程转化为变分形式，选定单元的形状，对求解域作剖分，进而构 造基函数或单元形状函数，形成有限元空间，将偏微分方程转化成了有限元方程，利用 有效的有限元方程的解法，给出偏微分方程近似解的过程称为有限元法。 6. 解：对求解区间进行网格剖分，节点01......i n a x x x x b =<<<<=得到相邻节点1,i i x x - 一、（10分）简叙偏微分方程数值解研究的内容。建立一个偏微分方程数值格式，需要研究和讨论哪些问题？一个好的数值格式应该到达哪些要求？ 二、（10分）简叙用差分方法求解偏微分方程数值解的方法和步骤。 三、（10分）请解释收敛性、稳定性、相容逼近性三个概念。它们之间有何区别和联系？ 四．对如下两点边值问题： ?? ???==<<=+-,0)1( )0(,10 ,222u u x f u dx u d 将区间[0，1]作剖分：N N x x x x <<<<=-1100 ， （1） 在上述网格剖分下建立中心差分格式； （10分） （2） 写出局部截断误差；并问当网格剖分满足什么条件 时，上述差分格式的误差阶为2； （10分） （3） 利用极值定理证明解的存在唯一性。 （10分） 四、 （1） 中心差分格式： )()()(2))()()()((2 1111u R x f x u h x u x u h x u x u h h i i i i i i i i i i i +=+---+- -+++ （2） 局部截断误差： )(3 )(2 331h O dx u d h h u R i i i i +??????--=+ 显然当剖分为等距时，误差阶为2 （3） 用极值定理证存在唯一性： 只需证明齐次（边值与右端恒零）问题只有平凡解。实际上，设i u 是齐次问题的解，则由极值定理，i u 既不能在区间内部取正的极大，也不能取负的极小，因此0=i u 。 七．已知两点边值问题：??? ??==<<=-0 )1( )0(10 ,22u u x f dx u d （A ） （1） 写出问题（A ）的虚功原理（含证明）； （25分） （2） 写出建立线性有限元方程组h h h b U A =的主要步骤； （10 分） （3） 对于线性有限元方程组h h h b U A =，设h A 产生扰动h A ?，h b 产偏微分方程数值解期末试题及标准答案
偏微分方程数值解法试题与答案
偏微分方程数值解
偏微分方程数值解试题及答案
偏微分方程数值解实验报告
偏微分方程数值解复习题(2011硕士)
偏微分方程数值解试题06B答案
偏微分方程数值解(试题)
偏微分方程数值解(试题)
偏微分方程数值解法
偏微分方程数值解试题参考答案
偏微分方程数值解例题答案
偏微分方程数值解试题参考答案
第十章-偏微分方程数值解法
偏微分方程数值习题解答
偏微分方程数值解法
偏微分方程数值解法答案
偏微分方程数值解试卷