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

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

《华科电磁场matlab仿真作业.docx》由会员分享，可在线阅读，更多相关《华科电磁场matlab仿真作业.docx（12页珍藏版）》请在冰点文库上搜索。

华科电磁场matlab仿真作业

电磁场作业

电气1202XXXU201200000

一．作业一

ΓΓΓΓΓΓ

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矩正存储每点电势

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）;

3.计算机绘图

二．作业二

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;

whilex>0.1

x=0;

fori=2:

100

forj=2:

100

if（i>1&&i<11）||（j>1&&j<11）

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（t>x）

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:

row

A（i,S（i））=0;

A（i,E（i））=0;

end%初始化电势完毕

temp=A;

forn=1:

600%迭代次数

fori=2:

row

forj=S（i）:

E（i）

ifj==S（i）||j==E（i）

temp（i,j）=0;

else

temp（i,j）=（3*A（i-1,j）+3*A（i+1,j）+A（i,j-1）+A（i,j+1））/8;

end

end

end

A=temp;

end

X=row:

-1:

1;

Y=N:

-1:

1;

[X,Y]=meshgrid（X,Y）;

figure

（1）;

surf（rot90（A,2））;

figure

（2）;

contour（rot90（A,2））;

3.计算机绘图