偏微分方程数值解法(抛物型方程差分法).ppt

1/17偏微分方程数值解法偏微分方程数值解法 7抛物型方程差分法抛物型方程差分法2差分格式稳定性概念差分格式稳定性概念显、隐格式稳定性分析显、隐格式稳定性分析稳定性分析的矩阵方法稳定性分析的矩阵方法2/17xt抛物型方程抛物型方程简单显式差分格式简单显式差分格式在实际应用时在实际应用时,取逐层计算形式取逐层计算形式.当初始层数据有误差当初始层数据有误差时时,误差会逐层传播误差会逐层传播,影响以后各层的解影响以后各层的解.记记 的误差为的误差为 ,设设 无误差无误差,则有则有3/17取取设初始层上设初始层上,仅有仅有 ,其它点处无误差其它点处无误差在各计算层上在各计算层上,误差传播得到控制误差传播得到控制4/17取取设初始层上设初始层上,仅有仅有 ,其它点处无误差其它点处无误差在各计算层上在各计算层上,误差传播没有得到控制误差传播没有得到控制5/17无穷大范数定义无穷大范数定义双层差分格式双层差分格式记矩阵记矩阵双层格式的矩阵形式双层格式的矩阵形式双层差分格式初值稳定概念双层差分格式初值稳定概念:任意解都满足任意解都满足其中其中 M 与与 无关无关.k k06/17简单显式差分格式简单显式差分格式稳定性分析稳定性分析,设设此时差分格式稳定此时差分格式稳定7/17设齐次方程设齐次方程系数矩阵可逆系数矩阵可逆 记记 称之为过渡矩阵称之为过渡矩阵 常系数差分格式常系数差分格式H 的谱半径的谱半径:定理定理:双层差分格式稳定的必要条件是双层差分格式稳定的必要条件是,存在与存在与 无无关的常数关的常数 M1,使得使得8/17定理定理 若若 H=A-1B 为正规矩阵为正规矩阵,即即 HH*=H*H,则条件则条件是双层差分格式按欧氏范数稳定的充分条件是双层差分格式按欧氏范数稳定的充分条件注：欧氏范数注：欧氏范数(或离散或离散L2范数范数)9/17简单显式差分格式矩阵形式简单显式差分格式矩阵形式过渡矩阵过渡矩阵特征值特征值 10/17过渡矩阵的谱半径过渡矩阵的谱半径 极值点满足极值点满足 显式差分格式稳定充分条件显式差分格式稳定充分条件.11/17简单隐式差分格式矩阵形式简单隐式差分格式矩阵形式特征值特征值 过渡矩阵过渡矩阵12/17过渡矩阵的谱半径过渡矩阵的谱半径 隐式差分格式无条件稳定隐式差分格式无条件稳定.13/17C-N 格式矩阵形式格式矩阵形式 特征值特征值 14/17过渡矩阵的谱半径过渡矩阵的谱半径 C-N 格式是无条件稳定的格式是无条件稳定的.15/17数值实验题数值实验题 用三种差分格式求用三种差分格式求初边值问题数值解初边值问题数值解 并与准确解比较并与准确解比较 显格式显格式隐格式隐格式C-N格式格式16/17数值计算实验数值计算实验显格式显格式:input T:=1error=7.9443e-006k=200 CN格式格式input T:=1error=2.6227e-006k=5017/17T=input(input T:=);h=1/10;ta=1/200;r=ta/(h*h);s=1-2*r;x=0:h:1;N=length(x);t=0;uk=sin(pi*x);II=2:N-1;k=0;%显格式计算程序显格式计算程序while tT t=t+ta;un=s*uk(II)+r*(uk(II-1)+uk(II+1);uk=0,un,0;k=k+1;endux=exp(-pi*pi*t).*sin(pi*x);error=max(abs(ux-uk)plot(x,ux,x,uk,:or)