fenbu.m

电磁学代码

function dwq2
prompt={'请输入上边界的电压值','请输入下边界的电压值','请输入左边界的电压值','请输入右边界的电压值'};
def={'','','',''};
Title='电压的边界值';
line=1;
a=inputdlg(prompt,Title,line,def);%绘制对话框
for k=1:4
b(k)=str2num(a{k});
end
prompt={'请输入横向网格数(大于零的整数)','请输入纵向网格数(大于零的整数)'};
def={'',''};
Title='电压';
line=1;
c=inputdlg(prompt,Title,line,def);%绘制对话框
m=str2num(c{1});%横向网格数
n=str2num(c{2});%纵向网格数
hx=m+1;hy=n+1;
v1=ones(hy,hx);
v1(hy,:)=b(1);%上边界的电压值
v1(1,:)=b(2);%下边界的电压值
v1(:,1)=b(3);%左边界的电压值
v1(:,hx)=b(4);%右边界的电压值
v2=v1;t=0;
k=0;%k为迭代次数
s=1;%s为误差值
while s>10^(-6)
    k=k+1; s=0;
    for i=2:hy-1
        for j=2:hx-1
            v2(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 t>s
                s=t;
            end
        end
    end
v1=v2;
end
k
mesh(v1)
colormap(jet)
title('电势分布图')