朱月恒 电磁场与电磁波实验.docVIP

电磁场与电磁波 实验报告 同实验者：冯泽涵 041010113 魏仁奎 041010119 实验 二维静电场边值问题的数值仿真实验 对于静电场边值问题，某些情况下用解析法求解很困难，甚至不可能。在实际求解过程中，可以将待求解电磁场问题化为离散的数值问题，即求解待求量在一些离散点上的数值。本实验根据二维静电场所满足的微分方程和边界条件，利用有限差分法来求解无源区域内一系列离散点上的电位，从而观测在给定边界条件下无源区域静电场的电位分布。 实验目的 1、学习用有限差分法求解静电场的方法。 2、观察边界条件对场的影响。 实验原理 1、方程的离散 有限差分法对于解决任何偏微分方程都是一种可行的方法，因为所有的电磁场问题都可以表示成标量或矢量偏微分方程，所以可以利用有限差分法解决各种介 质中电磁场的空间分布。有限差分法把要解决的区域划分为有限个离散点，并将偏微分方程用一组差分方程代替。因此，这种方法是近似的。但是，如果我们把离散点取得足够密，就能够把误差减小到可以接受的程度。 2、边界条件的处理 边界条件在电磁场问题的求解中起着至关重要的作用，因为不同电磁场问题可以满足相同的偏微分方程，它们之所以有不同的解就是因为满足的边界条件不同。为求得区域内各节点上近似解，还需将原问题中的边界条件化成离散的边界条件，即确定出在边界节点处的值。 3、差分方程的求解 实验内容 1、自行给定区域尺寸 a、b 和，取 Δ x=a/ 20 ， Δ y=b /10 ，在 y=b 的边界上令f ( x)=10 ，计算各内部节点处的电位值，并绘制出区域内的电位分布（可用Matlab 中的Mesh或Surf 函数）。 程序结果： 区域内电位分布： 2、其它条件不变，在y=b 的边界上令f(x )=10sin(x)，计算并绘制出区域内的电位分布。观察两种情况下电位分布的差别，理解边界条件对场分布的决定作用。 程序结果： 区域内电位分布： 3、对上述两种情况，分别以列表的形式给出 x= 10Δ x 和 y= 5Δ y 处各节点的电位与分离变量法（无穷级数截断到第100 项）所得结果的比较。 分离变量法在y= 5Δ y的结果： 分离变量法在x= 10Δ x 的结果： 程序代码：（其中a=20,b=10） a=20;B=10;pi=3.14159; dx=a/20;dy=B/10; M=a/dx-1;N=B/dx-1; m=M; x=1:a/20:a; %f=10*sin(pi*x/a); %f=x; A=zeros(N*M,N*M); % 事先把系数矩阵各元素置0，在重新赋值前保持为0 b=zeros(N*M,1); % 事先把已知列向量各元素置0，在重新赋值前保持为0 for j=1:N for i=1:M k=(j-1)*M+i; if j~=1 j~=Ni~=1 i~=M A(k,k-1)=1.0/dx^2; A(k,k)=-2.0*(1.0/dx^2+1.0/dy^2); A(k,k+1)=1.0/dx^2; A(k,k-m)=1.0/dy^2; A(k,k+m)=1.0/dy^2; elseif i==1j~=1j~=N A(k,k)=-2.0*(1.0/dx^2+1.0/dy^2); A(k,k+1)=1.0/dx^2; A(k,k-m)=1.0/dy^2; A(k,k+m)=1.0/dy^2; elseif i==Mj~=1j~=N A(k,k-1)=1.0/dx^2; A(k,k)=-2.0*(1.0/dx^2+1.0/dy^2); A(k,k-m)=1.0/dy^2; A(k,k+m)=1.0/dy^2; elseif j==1i~=1i~=M A(k,k-1)=1.0/dx^2; A(k,k)=-2.0*(1.0/dx^2+1.0/dy^2); A(k,k+1)=1.0/dx^2; A(k,k+m)=1.0/dy^2; elseif j==Ni~=Mi~=1 A(k,k-1)=1.0/dx^2; A(k,k)=-2.0*(1.0/dx^2+1.0/dy^2); A(k,k+1)=1.0/dx^2; A(k,k-m)=1.0/dy^2; b(k)=(-1)*10/dy^2; % f(i*dx)是上边界给定的已知电位 elseif i==1j==1 A(k,k)=-2.0*(1.0/dx^2+1.0/dy^2); A(k,k+1)=1.0/dx^2; A(k,k+m)=1.0/dy^