华科电磁场matlab仿真作业

华科电磁场m a t l a b仿真作

业

-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

电磁场作业

电气1202 XXX U201200000一．作业一

1.程序框图

2.程序

clear;

col = 61; %第一行点数

row = col; %行数

span = 0.3/(col-1); %步长

End = ones(1,col)*col; %每一行的终止点

Start = ones(1,col); %每一行的起始点

A = zeros(row,col); %A矩正存储每点电势

for i = (col-1)/3+1:(col-1)*2/3+1

for j = (col-1)/3+1:(col-1)*2/3+1

A(i,j) =100;

end

end %初始化电势完毕

temp = A;

for n= 1:500 %迭代次数

for i = 2:row-1

if ( i<((col-1)/3+1)||i>( (col-1)*2/3+1 ) )

for j = Start(i)+1:End(i)-1

temp(i,j)=(A(i-1,j) +A(i+1,j) +A(i,j-1) +A(i,j+1))/4;

end

else

for j = 2:(col-1)/3

temp(i,j)=(A(i-1,j) +A(i+1,j) +A(i,j-1) +A(i,j+1))/4;

end

for j = 2*(col-1)/3+2:col-1

temp(i,j)=(A(i-1,j) +A(i+1,j) +A(i,j-1) +A(i,j+1))/4;

end

end

A = temp;

end

end

X = row:-1:1;

Y = col:-1:1;

[X,Y] = meshgrid(X,Y);

figure(1);

surf(rot90(A,2));

figure(2);

contour(rot90(A,2));

hold on;

[Gx,Gy] = gradient(A,1,1);

quiver(Gx,Gy);

3.计算机绘图

二．作业二

1.程序框图

2.程序100V

0V