偏微分方程数值解法抛物型方程差分法.ppt
天马****23
亲,该文档总共17页,到这已经超出免费预览范围,如果喜欢就直接下载吧~
相关资料
偏微分方程数值解法(抛物型方程差分法).ppt
《偏微分方程数值解法》7——抛物型方程差分法2x取取无穷大范数定义简单显式差分格式设齐次方程定理若H=A-1B为正规矩阵,即HH*=H*H,则条件简单显式差分格式矩阵形式过渡矩阵的谱半径简单隐式差分格式矩阵形式过渡矩阵的谱半径C-N格式矩阵形式过渡矩阵的谱半径数值实验题用三种差分格式求初边值问题数值解数值计算实验T=input('inputT:=');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
偏微分方程数值解法抛物型方程差分法.ppt
《偏微分方程数值解法》7——抛物型方程差分法2x取取无穷大范数定义简单显式差分格式设齐次方程定理若H=A-1B为正规矩阵,即HH*=H*H,则条件简单显式差分格式矩阵形式过渡矩阵的谱半径简单隐式差分格式矩阵形式过渡矩阵的谱半径C-N格式矩阵形式过渡矩阵的谱半径数值实验题用三种差分格式求初边值问题数值解数值计算实验T=input('inputT:=');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
抛物型方程有限差分法.doc
抛物型方程有限差分法抛物方程差分法的构造在空间方向上与椭圆方程类似,在时间方向上用一阶差商代替代替一阶微商。然后在时间方向上逐层求解。特别当空间维数较高时,可以使用局部一维格式大大降低计算量。1.简单差分法考虑一维模型热传导方程(1.1),其中为常数。是给定的连续函数。(1.1)的定解问题分两类:第一,初值问题(Cauchy问题):求足够光滑的函数,满足方程(1.1)和初始条件:(1.2),第二,初边值问题(也称混合问题):求足够光滑的函数,满足方程(1.1)和初始条件:,及边值条件,假定和在相应的区域光
大连理工大学高等数值分析抛物型方程有限差分法.doc
抛物型方程有限差分法抛物方程差分法的构造在空间方向上与椭圆方程类似,在时间方向上用一阶差商代替代替一阶微商。然后在时间方向上逐层求解。特别当空间维数较高时,可以使用局部一维格式大大降低计算量。1.简单差分法考虑一维模型热传导方程(1.1),其中为常数。是给定的连续函数。(1.1)的定解问题分两类:第一,初值问题(Cauchy问题):求足够光滑的函数,满足方程(1.1)和初始条件:(1.2),第二,初边值问题(也称混合问题):求足够光滑的函数,满足方程(1.1)和初始条件:,及边值条件,假定和在相应的区域光
偏微分方程数值解法.doc
个人收集整理勿做商业用途个人收集整理勿做商业用途个人收集整理勿做商业用途一、问题用有限元方法求下面方程的数值解inonin二、问题分析第一步利用Green公式,求出方程的变分形式变分形式为:求,使得(*)以及.第二步对空间进行离散,得出半离散格式对区域进行剖分,构造节点基函数,得出有限元子空间:,则(*)的Galerkin逼近为:,求,使得(**)以及,为初始条件在中的逼近,设为在中的插值.则,有,=,代人(**)即可得到一常微分方程组。第三步进一步对时间进行离散,得到全离散的逼近格式对用差分格式.为此把