电磁场和电磁波演示验证实验1

Рjingdu=10^(-10); %%求解精度Рnum=0; %%迭代次数初始化Р%%%%¸赋初值Рv1=zeros(nodenumy+1,nodenumx+1);v1(nodenumy+1,:)=ones(1,nodenumx+1)*U;v2=v1;%%前后两次迭代值Рr=zeros(nodenumy-1,nodenumx-1); %迭代差值Рd=3; %求解精度初始化可以赋值任意大于精度的值Рwhile(d>jingdu) Р num=num+1;Р v1=v2;Р for I=2:nodenumyР for J=2:nodenumxР v2(I,J)=(v1(I,J+1)+v1(I+1,J)+v1(I-1,J)+v1(I,J-1))/4;%%简单迭代Р r(I,J)=abs(v2(I,J)-v1(I,J)); Р endР endР d=max(max(r));РendР3.超松弛迭代法程序:Р for i=2:MР for j=2:M+1Р x(i,j)=j-1;Р endРend Р y=x; %赋初值Рflag=1;N=[ 0 0 0 0 0 0 0 0 0 0 0]; %迭代次数赋初值Рfor z=1:11Рwhile flag==1Р for i=2:MР for j=2:MР b=0.25*(y(i-1,j)+x(i+1,j)+y(i,j-1)+x(i,j+1));Р y(i,j)=x(i,j)+a1(z)*(b-x(i,j));Р endР end %循环迭代Р if max(abs(x-y))<werrorР flag=0; %跳出循环Р endР x=y;Р N(z)=N(z)+1; %迭代次数加1 Р endРflag=1;Рfor i=2:MР for j=2:M+1Р x(i,j)=j-1;Р endРendРy=x; %切记再赋初值Рend