omg.rar_Coupled Line_Coupled Microstrip_shielding matlab

FI=100; hx=29;hy=29; CX1=6;CX2=12;CY1=6;CY2=12; CX3=18;CX4=24;CY3=6;CY4=12; DX=1+CX2-CX1;DY=1+CY2-CY1; DX0=1+CX4-CX3;DY0=1+CY4-CY3; v=ones(hy,hx)*50; mesh1=ones(hy,hx)*2; v(1,:)=zeros(1,hx); mesh1(1,:)=zeros(1,hx); v(hy,:)=zeros(1,hx); mesh1(hy,:)=zeros(1,hx); v(:,1)=zeros(hy,1); mesh1(:,1)=zeros(hy,1); v(:,hx)=zeros(hy,1); mesh1(:,hx)=zeros(hy,1); v(CY1:CY2,CX1:CX2)=ones(DX,DY)*FI; mesh1(CY1:CY2,CX1:CX2)=ones(DX,DY); v(CY3:CY4,CX3:CX4)=ones(DX0,DY0)*FI; mesh1(CY3:CY4,CX3:CX4)=ones(DX0,DY0); k=0; difmax=1.0; while(difmax>1e-6) k=k+1; difmax=0.0; for i=2:hy-1 for j=2:hx-1 m=mesh1(i,j); if (m>=2) vold=v(i,j); v(i,j)=(1/4)*(v(i-1,j)+v(i,j-1)+v(i+1,j)+v(i,j+1)); dif=v(i,j)-vold; dif=abs(dif); if(dif>difmax)difmax=dif;end end end end end subplot(1,2,1),mesh(v) axis([-2,hx+3,-2,hy+3,0,100]) subplot(1,2,2),contour(v,13) hold on plot([CX1,CX1,CX2,CX2,CX1],[CY1,CY2,CY2,CY1,CY1],'k') plot([CX3,CX3,CX4,CX4,CX3],[CY3,CY4,CY4,CY3,CY3],'k') text(CX1+0.6,CY1+(CY2-CY1)/2,'U=100','fontsize',10); text(CX3+0.6,CY3+(CY4-CY3)/2,'U=100','fontsize',10); text(hx/2,hy+1,'U=0','fontsize',10); text(hx/2,0,'U=0','fontsize',10); text(-1.7,hy/2,'U=0','fontsize',10); text(hx+0.4,hy/2,'U=0','fontsize',10); hold off