c语言矩形槽内电位分布,槽内电位分布.doc

作业一：MATLAB仿真及电位分布的计算

姓名:孙云霞 学号:0211433 专业:电力电子与电力传动

题：一个长直接地金属矩形槽，侧壁与底面电位均为0，顶盖电位是100V，求槽内的电位分布，画出等位线。

U=100V

U=0V U=0V

U=0V

1.程序：

a=15;b=12;

w1=ones(b,a);

w1(b,:)=ones(1,a)*100;

w1(1,:)=zeros(1,a);

for i=1:b

w1(i,1)=0;

w1(i,a)=0;

end

w2=ones(b,a);

w2=w1;s=0;maxt=1;

while(maxt>0.000001)

maxt=0;

for i=2:b-1

for j=2:a-1

w2(i,j)=( w1(i,j+1)+w1(i+1,j)+w2(i-1,j)+w2(i,j-1))/4;

s=abs(w2(i,j)-w1(i,j));

if(s>maxt) maxt=s; end

end

end

w1=w2

end

axis([0,17,0,11,0,100])

subplot(1,2,2),contour(w2,15)

hold on

x=1:1:a;y=1:1:b

[xx,yy]=meshgrid(x,y);

axis([-1.5,a+2.5,-2,13])

plot([1,1,a,a,1],[1,b,b,1,1],'k')

text(a/2,0.3,'0V','fontsize',11);

text(a/2-0.5,b+0.5,'100V','fontsize',11);

text(-0.5,b/2,'0V','fontsize',11);

text(a+0.3,b/2,'0V','fontsize',11);

2.程序运行结果(电位计算结果)

3.仿真结果