超松弛迭代法求解接地金属槽内电位分布

1、实验一用超松弛迭代法求解接地金属槽内电位分布1、 实验内容：试用超松弛迭代法求解接地金属槽内电位的分布。已知：，给定边值如图所示。给定初值：误差范围：计算迭代次数，分布。一.实验思路由边界条件用泊松方程的五点差分格式求得中央点的点位。再以所得点及边界再次利用泊松方程的五点差分格式求出另四个点，依照此方法求出其余点的电位分布。用最佳收敛因子的经验公式计算收敛因子。利用超松弛迭代法进行差分方程的求解，当遇到边界是采用边界值或者边界差分格式。直到所有节点电位满足误差条件。二实验设计原理：有限差分法有限差分法（Finite Differential Method）是基于差分原理的一种数值计算法。其基本

2、思想：将场域离散为许多小网格，应用差分原理，将求解连续函数的泊松方程的问题换为求解网格节点上的差分方程组的问题。编程时将边值编入到程序中，这样可以省略输入，从而直接输入迭代因子进行求解，可以减少编程的难度。迭代时所用公式是和书上一样，为aij=bij+w/4*(bi+1j+bij+1+aij-1+ai-1j-4*bij);其中a代表k+1，而b代表k。启动输出开始菜单（边值都已经给定）输入迭代因子w迭代次数n=0n+ 开始循环迭代函数判断相邻二次差值是否小于给定值输出n,电位aij终止三、程序运行界面及结果四源程序代码#include#include#includeusing namespac

3、e std;classoverrei /over-relaxation iterative methodprivate:intm,n;doublex,e;double *p,*q;public:int k;overrei(int m0,int n0,double e0)inti;e=e0;k=0;m=m0;n=n0;p=new double *m;for(i=0;im;i+)pi=new doublen;q=new double *m;/迭代因子求解for(i=0;im;i+)qi=new doublen;if(m=n)x=2/(1+sin(3.141592654/(m-1);elsex=2-

4、3.141592654*sqrt(2)*sqrt(1/(m-1)*(m-1)+1/(n-1)*(n-1);cout最佳收敛因子:xendl;void Initialization(); /赋边界条件void Cal(); /计算void Diedai(); /迭代函数void Show(); /输出部分;void overrei:Initialization() /赋边界条件inti,j; for(i=0;im;i+) /边界条件for(j=0;jn;j+)if(i=0)pij=100; qij=100;elsepij=0;qij=0;cout初始点位:endl;Show();void ove

5、rrei:Cal() /计算inti,j;int c=1;while(1) c=1;for(i=0;im;i+)for(j=0;je|(qij-pij)e) /相邻两次迭代误差是否小于1e-5c=0;break;if(c=0) break;if(c=1 & k!=0) break;Diedai();void overrei:Diedai() /迭代函数inti,j;double y=x/4;if(k%2)for(i=1;im-1;i+)for(j=1;jn-1;j+)qij=pij+y*(qi-1j+qij-1+pi+1j+pij+1-4*pij);elsefor(i=1;im-1;i+)fo

6、r(j=1;jn-1;j+)pij=qij+y*(pi-1j+pij-1+qi+1j+qij+1-4*qij);k+;void overrei:Show() /输出部分inti,j;for(i=0;im;i+)for(j=0;jn;j+)coutsetw(12)setprecision(6)fixedqijends;coutendl;int main()cout *endl;cout 超松弛迭代法求解接地金属槽内电位的分布endl;cout *endl;overrei A(5,5,1e-5);A.Initialization();A.Cal();cout电位分布:endl;A.Show();c

7、out迭代次数:A.kendl;return 0;实验二按对称场差分格式求解电位的分布一.实验思路只计算一半的区域，对另一半进行对称性计算，减小计算量。选取二维数组，同样是两个数组a4141、b4141,分别用于表示迭代前后的取值。由此对第20列进行特殊处理。即第21列是和第20列相等的，因为在第21列上标注了/x=0，即要求每次迭后都要求u1i20=u1i19,其余则同前。当计算出了前21列的电位值时，根据对称性可得到右半边的电位值，即要求u1i40-j=u1ij;，直接赋值。四周的电位值依旧保持不变，不能够参与迭代。并且中间对称列上的电位值为零，故也不参与迭代。由题意可知，只有第1列到第2

8、0列参与迭代。在进行精度计算时，用数组b的值与相应数组a里的值进行比较，没迭代一次都要进行比较。如果在误差允许的范围内，则继续迭代；当超过误差标准的时候则停止迭代。依旧按照比较得到最少收敛次数的方法在程序中计算出最佳收敛因子。然后对对称列右边的点进行对称赋值。最后输出输出最佳迭代因子、迭代的次数和迭代后各相应点的电位值。二.源程序代码#include#include#includevoidmain()doubleu14141,u24141;intflag=10000,i,j,e,n;doublea,d;/a最佳加速收敛因子for(j=0;j41;j+)u10j=100;u140j=0;for(

9、i=1;i40;i+)u1i0=0;u1i40=0;for(i=1;i40;i+)/内部点赋初值for(j=1;j21;j+)u1ij=2.5*(j-1);cout迭代前左半区域各电位点上的初始值:endl;for(i=0;i41;i+)for(j=0;j21;j+)coutu1ij;/输出迭代初值coutendl;coutendln;for(d=1;d2;d+=0.01)/取迭代次数最少的加速收敛因子n=0;for(j=0;j41;j+)u10j=100;u140j=0;for(i=1;i40;i+)u1i0=0;u1i40=0;for(i=1;i40;i+)for(j=1;j21;j+)u

10、1ij=2.5*(j-1);do/迭代x次for(i=0;i41;i+)for(j=0;j22;j+)u2ij=u1ij;for(i=1;i40;i+)for(j=1;j20;j+)u1ij=u2ij+(d/4)*(u2i+1j+u2ij+1+u1i-1j+u1ij-1-4*u2ij);u1i20=u1i19;for(i=1;i40;i+)/判断精度for(j=1;j1e-5)e=1;break;elsee=0;if(e=1)break;n+;/迭代次数while(e);if(nflag)flag=n;a=d;for(i=1;i40;i+)/对其对称点赋值for(j=1;j21;j+)u1i4

11、0-j=u1ij;cout迭代后各节点上最终电位近似值:endl;for(i=0;i41;i+)for(j=0;j41;j+)coutsetiosflags(ios:fixed)setprecision(5)setw(8)u1ij;/对齐输出coutendl;coutendlendl;cout收敛因子:aendl;coutendl;cout迭代次数:flagendl;coutendl; 三.程序运行界面及结果附录：程序一：#include#includevoid main()double m55,n55;int N=0,b=1;inti,j;double e=0.00001;double a=

12、2/(1+sin(3.1415926/4);for(i=0;i=4;i+)for(j=0;j=0;j-)for(i=0;i=4;i+)coutmij=mijt;coutendl;while(b=1)b=0;N=N+1;for(i=1;i=3;i+)for(j=1;j=3;j+)mij=mij+a*(mi-1j+mij-1+mi+1j+mij+1-4*mij)/4;for(i=1;i=3;i+)for(j=1;j=e)b=1;nij=mij;for(j=4;j=0;j-)for(i=0;i=4;i+)coutmij=mijt;coutendl;cuntN=Nendl;程序二：#include#i

13、nclude#includeusing namespace std;classoverrei /over-relaxation iterative methodprivate:intm,n;doublex,e;double *p,*q;public:int k;overrei(int m0,int n0,double e0)inti;e=e0;k=0;m=m0;n=n0;p=new double *m;for(i=0;im;i+)pi=new doublen;q=new double *m;/迭代因子求解for(i=0;im;i+)qi=new doublen;if(m=n)x=2/(1+si

14、n(3.141592653/(m-1);elsex=2-3.141592653*sqrt(2)*sqrt(1/(m-1)*(m-1)+1/(n-1)*(n-1);cout最佳收敛因子:xendl;void Initialization(); /赋边界条件void Cal(); /计算void Diedai(); /迭代函数void Show(); /输出部分;void overrei:Initialization()inti,j; /边界条件for(i=0;im;i+)for(j=0;jn;j+)if(i=0)pij=100; qij=100;elsepij=0;qij=0;cout初始点位:

15、endl;Show();voidoverrei:Cal()inti,j;int c=1;while(1) c=1;for(i=0;im;i+)for(j=0;je|(qij-pij)e)c=0;break;if(c=0) break;if(c=1 & k!=0) break;Diedai();voidoverrei:Diedai()inti,j;double y=x/4;if(k%2)for(i=1;im-1;i+)for(j=1;jn-1;j+)qij=pij+y*(qi-1j+qij-1+pi+1j+pij+1-4*pij);elsefor(i=1;im-1;i+)for(j=1;jn-1

16、;j+)pij=qij+y*(pi-1j+pij-1+qi+1j+qij+1-4*qij);k+;voidoverrei:Show()inti,j;for(i=0;im;i+)for(j=0;jn;j+)coutsetw(12)setprecision(6)fixedqijends;coutendl;int main()cout *endl;cout 超松弛迭代法求解接地金属槽内电位的分布endl;cout *endl;overrei A(5,5,1e-6);A.Initialization();A.Cal();cout电位分布:endl;A.Show();cout迭代次数:A.kendl;r

17、eturn 0;程序三#include#includeinti,j,N=0;/N为迭代次数doubleu155;/定义电位doubleu255;doublee=0.00001;/误差范围doublea=1.5;/加速收敛因子intfinish();/迭代结束函数intNresult();/求解迭代次数voidprint();voidmain()if(i=4)for(j=0;j5;j+)u1ij=100;/边界值elseu1ij=0;N=Nresult();print();voidprint()cout迭代次数为：Nendl;for(j=0;j5;j+)coutu24jt;coutendl;fo

18、r(j=0;j5;j+)coutu23jt;coutendl;for(j=0;j5;j+)coutu22jt;coutendl;for(j=0;j5;j+)coutu21jt;coutendl;for(j=0;j5;j+)coutu20jt;coutendl;intfinish()doubleb;b=u2ij-u1ij;if(-eb)&(be)return1;elsereturn0;intNresult()doN+;for(i=1;i4;i+)for(j=1;j4;j+)u1ij=u2ij;for(j=1;j4;j+) u21j=u11j+a/4*(u20j+u21j-1+u12j+u11j+

19、1-4*u11j);u1ij=u2ij;for(j=1;j4;j+) u22j=u12j+a/4*(u21j+u22j-1+u13j+u12j+1-4*u12j);u1ij=u2ij;for(j=1;j4;j+) u23j=u13j+a/4*(u22j+u23j-1+u14j+u13j+1-4*u13j);u1ij=u2ij;while(finish()=0);returnN; 程序四#include#include#includeusingnamespacestd;classcscddprivate:intm,n;doublex,e;double*p,*q;public:intk;cscdd

20、(intm0,intn0,doublee0)inti;e=e0;k=0;m=m0;n=n0;p=newdouble*m;for(i=0;im;i+)pi=newdoublen;q=newdouble*m;for(i=0;im;i+)qi=newdoublen;if(m=n) x=2/(1+sin(3.141592654/(m-1);elsex=2-3.141592654*sqrt(2)*sqrt(1/(m-1)*(m-1)+1/(n-1)*(n-1);cout最佳收敛因子:xendl;voidbjtj();/赋边界条件voidjishuan();/计算voidDiedai();/迭代函数voidsuchu();/输出部分;voidcscdd:bjtj()/赋边界条件inti,j;for(i=0;im;i+