[工学]偏微分方程数值解第9次课.pptVIP

1 二阶线性偏微分方程 1 二阶线性偏微分方程 2 差商代替偏导建立差分格式 3 积分插值原理构造差分格式 微分方程转化为积分方程 差分格式构造 区域离散 方程离散 边界处理 差分格式相容收敛和稳定性讨论 求解结果 微分方程转化为积分方程 差分格式构造 区域离散 方程离散 第二、三边界处理 4 椭圆型差分方程极值原理 极值原理: 定理2.4.1 比较定理: 定理2.4.2 最大模估计定理: 推论2.4.3 5 五点差分格式的实用分析 相容性 收敛性 稳定性 6 计算机实现 五点差分格式的系数矩阵 算法 程序 function u=possion(f,f1,f2,f3,f4,a,b,h,max1) % SOD solve possion problem. % input: f(x,y) right-side function; % f1(x,y),f2(x,y),f3(x,y),f4(x,y)-boundary condtion % (down-upper-left-right) % a,b the interval; h- step size, max1 interative %output: u the solution on the meshgrid % initialize parameters and u m=fix(a/h)+1; % horizontal direction, m-1 part, m point n=fix(a/h)+1; %vertical direction, n-1 part u=zeros(m,n); % boundary condition u(1,1:n)=feval(f3,0,0:h:(n-1)*h); %left-side u(m,1:n)=feval(f4,a,0:h:(n-1)*h); %right-side u(1:m,1)=feval(f1,0:h:(m-1)*h,0); % down-side u(1:m,n)=feval(f2,0:h:(m-1)*h,b); % upper-side u(1,1)=(u(1,2)+u(2,1))/2; %left-down u(m,1)=(u(m-1,1)+u(m,2))/2; %right-down u(1,n)=(u(1,n-1)+u(2,n))/2; % left-upper u(m,n)=(u(m-1,n)+u(m,n-1))/2; % right-upper % lefr-side function ff=zeros(m,n); for j=1:n for i=1:m ff(i,j)=feval(f,i*h,j*h); end end %SOR parameter w=4/(2+sqrt(4-(cos(pi/(n-1))+cos(pi/(m-1)))^2)); %refine approximations and sweep operator err=0;cnt=0; while((err0.0001)(cnt=max1)) for j=2:n-1 for i=2:m-1 relx=w*(h^2*f(i,j)+u(i,j+1)+u(i,j-1)+u(i+1,j)+u(i-1,j)-4*u(i,j))/4; u(i,j)=u(i,j)+relx; if(err=abs(relx)) err=abs(relx); end end end cnt=cnt+1; end u=flipud(u); Guass-Seidel迭代法求解代数方程Au=f 将A分解 G-S迭代法 超松弛迭代法求解代数方程Au=f A大型稀疏矩阵 利用G-S迭代法 将A分解 确定出 然后通过正参数加速迭代 将其代入上式, 消去 逐次超松弛迭代法(SOD) Successive Over Relaxation Method 分量形式 例：解线性方程组 解: SOR 迭代格式 对于矩形区域possion方程第一边值问题差分格式 边界值的平均值 系数矩阵的特征值. * * * * 椭圆型差分方程数值解法 2 差商代替偏导建立差分格式 3 积分插值原理构造差分格式 4 椭圆型差分方程极值原理 5 五点差分格式的实用分析 6 计算机实现 考虑考虑两个自变量(x, y)的二阶线性偏微分方程： 椭圆型 抛物型 双曲型