华科电磁场matlab仿真作业.docxVIP

电磁场作业

电气1202XXXU201200000

作业一

0V

0V

100V

100V

ΓΓΓΓΓΓ

程序框图

程序

clear;

col=61;%第一行点数

row=col;%行数

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

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

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

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

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

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

A(i,j)=100;

end

end%初始化电势完毕

temp=A;

forn=1:500%迭代次数

fori=2:row-1

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

forj=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

forj=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

forj=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));

holdon;

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

quiver(Gx,Gy);

计算机绘图

作业二

0V100V

0V

100V

1.程序框图

2.程序

clear;

col=101;

row=101;

A=zeros(row,col);

fori=11:101

A(i,11)=100;

end

forj=11:101

A(11,j)=100;

end

fori=2:10

A(i,col)=10*i-10;

end

forj=2:10

A(row,j)=10*j-10;

end

A1=A;

x=1;t=0;

whilex0.1

x=0;

fori=2:100

forj=2:100

if(i1i11)||(j1j11)

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

t=abs(A(i,j)-A1(i,j));

if(tx)

x=t;

end

end

end

end

A=A1;

end

X=1:1:row;

Y=1:1:col;

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

figure(1);

surf(rot90(A,4));

figure(2);

contour(rot90(A,4));

holdon;

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

quiver(Gx,Gy);

3.计算机绘图

作业三

1.程序框图

2.程序

clear;

N=41;%第一行点数

h=0.2/(N-1);%步长

row=(N+1)/2;%行数

E=N:-1:N-row+1;%每一行的终止点

S=1:1:N-row+1;%每一行的起始点

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

A(1,:)=ones(1,E(1))*100;%初始化电势

fori=1:ro