电磁场与电磁波性学习—有限差分方法技术方案.pptx

研究型学习报告 利用有限差分法解决静电场第一类边值问题 马啸理论分析 理论分析 有限差分法——拉普拉斯方程的离散化 将区域D划分为正方形网格，在节点处代入导数 公式，可得 将其带入拉普拉斯方程，可得方程的差分格式为 理论分析 差分方程组的解——Jacobi迭代法 分量格式 向量格式 其中 理论分析 差分方程组的解——Gauss-Seidel 迭代法 分量形式 向量形式 其中 理论分析 差分方程组的解——Successive Over Relaxation 迭代法 分量形式 向量形式 其中 参考文献给出，最佳松弛因子的计算公式为 理论分析 计算机仿真 问题描述 如图所示，有一长方形的导体槽，a = 20，b = 5，设槽的长度为无限长，槽上有两块与槽绝缘的盖板，电位分别为100V和 50V，其它板电位为0V，求槽内的电位分布. 计算机仿真 获取系数矩阵 %get_coeeficient_matrix_function %输入未知数矩阵的行数m及列数n,根据有限差分法返回方程的系数矩阵及常数列向量 function [A,b]=get_coefficient_matrix_function(m,n) B=zeros(m*n,m*n); b=zeros(m*n,1); rup=50;rlw=100;rlf=0;rrt=0; i=1; j=1; k=(i-1)*n+j;%将矩阵个数的未知数转换为自然数排列 B(k,k+1)=1; B(k,k+n)=1; b(k,1)=rup+rlf; for j=2:n-1 k=(i-1)*n+j; B(k,k-1)=1;B(k,k+1)=1;B(k,k+n)=1; b(k,1)=rup; end j=n; k=(i-1)*n+j; B(k,k-1)=1;B(k,k+n)=1; b(k,1)=rup+rlf; %第一行完毕 计算机仿真 获取系数矩阵 B=0.25.*B; %得到Jacobi迭代矩阵 for i=1:m*n B(i,i)=-1; end A=-B; %得到系数矩阵   b=b/4; %A=(A-D)+D D=zeros(m*n,m*n);%A所对应的对角阵 for i=1:m*n D(i,i)=A(i,i); end f=D*b; b=f; 计算机仿真 Gauss-Seidel迭代法 h=input(所期望步长:); col_h=20/h+1;col_x=col_h-2; line_h=5/h+1;line_x=line_h-2; u=zeros(line_h,col_h); %创建符合步长要求的矩阵存放各标志点电位值   for j=2:col_h-1 for i=2:line_h-1 u(i,j)=50+50/line_h*(i-1); %初始条件按内部线性变化确定 end %初始条件对结果影响不大，但影响循环次数 end %完成差分方程初始条件的设定 [A,b]=get_coefficient_matrix_function(line_x,col_x);%使用自定义函数获得方程组系数矩阵及常数列向量 [mat,flag]=chol(A);%检查A是否正正定 if flag==0 disp(收敛) else disp(不收敛) end %若A正定，可以使用Gauss-Seidel迭代法求解方程 计算机仿真 while flag==0 %在迭代法收敛的条件下进入循环 while (t1e-6) n=n+1; t=0; %必须将误差控制项置零 for i=2:line_h-1; for j=2:col_h-1; v(i,j)=(u(i,j+1)+u(i+1,j)+v(i-1,j)+v(i,j-1))/4; %Seidel迭代法 if tabs(v(i,j)-u(i,j)) t=abs(v(i,j)-u(i,j));%用无穷大范数作为误差控制 end end end u=v; end end 计算机仿真 Successive Over Relaxation 迭代法 w_sor