《二维场的数值计算实验报告》由会员分享，可在线阅读，更多相关《二维场的数值计算实验报告（5页珍藏版）》请在装配图网上搜索。
1、二维场的数值计算实验报告一实验目的通过采用两种不同的计算方法，上机计算二维场的电位分布，了解数值法与 解析法的优缺点和区别，掌握超松弛迭代方法。二实验内容1根据分离变量法求得的二维场的用级数表示的解，编写程序，上机计算二 维场中各点的电位；2根据有限差分法原理，编写超松弛迭代法程序，上机计算二维场中各点的 电位。三实验题目如图所示横截面为矩形的无限长直导体槽，内填空气。已I IIT IT(0, y ) = V (a, y ) = V(x,0 ) = V (x, b ) = V34o1. 根据分离变量法求得的二维场的用级数表示的解，编写程序， 上机计算二维场中各点的电位；2. 根据有限差分法原理
2、，编写超松弛迭代法程序，上机计算二维场中各点的电 位。四.实验设计过程1 分离变量法：在利用解析法求解边界问题时，一般都是利用麦克斯韦方程，在已知边界 值的情况下，结合条件球场电位函数。结果是一个复杂的级数函数表达式。虽然 不能具体体现场的分布，但可以清楚地算出场中的每一个数值，非常的精确。根 据上述题目，电位分布与z轴无关，并且槽内点电荷密度P二0，由此可以知道 其满足的拉普拉斯方程及边界条件如下：d 2 d 2+ dx 2dy 2(x, y)二 0 ,y=0(x, y )= Vy=b3通过对第一个方程进行分离变量法得到方程的通解为：山,y )=(C sin k x + C cos k x)
3、D sinh k y + D cosh k y)1x2x1y2y再将后四个方程边界条件代入求得电位函数:n为奇数1n兀.n兀smhy sm -n sinh 巴 b a aa2.有限差分法在题目中,有限差分法与分离变量法的最主要区别是对于（1）中的拉普拉斯方程进行 近似计算，将场域划分为网格节点，将网格节点电位作为未知数的差分近似替代 该点的偏导数，从而把待求得偏微分方程的边值问题转化为一组相应的差分方程 问题。最后利用超松弛迭代法计算出数值。由于点电荷密度为0则拉普拉斯方程的有限差分方程形式为具体步骤如下：1. 网格设置2. 区域设置3. 设置边界条件4. 应用模式设置选择（静电学）5方程参数
4、设定五实验程序先假设V、V、V、V的数值分别为0, 0, 100V, 0, a和b的数值为10和5cm。12341分离变量法： cleara=0;X,Y=meshgrid(0:0.05:2, 2:0.05:4);for n=1b=4*100/pi; a=b*sinh(n*pi*Y/2).*sin(n*pi*X/2)./n*sinh(n*pi*4/2); a=a+a;enda=a;figurecontour(a,20);2有限差分法：clearhx=40;hy=40;vl=ones(hy, hx);vl(hy,:)=ones(l, hx)*100;vl(l,:)=zeros(l, hx);for
5、 i=l:hyvl(i, 1)=0;v1(i, hx)=0;endv2=v1; maxt=1; t=0;k=0while(maxt1e-9)k=k+1;maxt=0;for i=2: hy-1for j=2:hx-1v2(i, j)=(v1(i, j+1)+v1(i+1, j)+v2(i-1, j)+v2(i, j-1)/4; t=abs(v2(i, j)-v1(i, j);if(tmaxt)maxt=t;endendendv1=v2;endcontour(v2, 20)hold onx=1:1:hx;y=1:1:hy;Gx, Gy=gradient(- v2, 10, 10); quiver
6、(xx, yy, Gx, Gy, 8, r) hold off六实验结果1分离变量法实验图形七实验结果分析比较和总结对比上述两种方法求解思路以及所得到的实验图形可知，虽然它们的解题思路不同，但是在通过程序进行画图后，他们的电位分布图是一样的，都无一例外 地真实准确地展示了静态场的电位分布。但是，他们也有不同的地方，分离变量 法求得的是函数，然后通过离散来展示场的分布；而有限差分法得到的解本身是 离散的，通过编程不断迭代来得到数值以接近分离变量法得到的函数。求解静态场分布的方法是多样的，但最终借助于MATLAB这个程序使我们 的结果趋于一致，可见我们求解问题并不一定只局限于一个方法，多样的方法最 终都会走向一个结果，我们需要综合这些方法，在其中摸索发现，以期有更好的 、p r r.前进。
二维场的数值计算实验报告
