%本程序用來模擬平面電磁波在空間中的傳播

%假設電磁波沿著X軸正向傳播，電場具有Z方向的分量Ez，磁場具有Y方向上的分量Hy

%由于程序中的電場分量與磁場分量均與兩個參數-時間和坐標有關，因此為了使

%兩個參數的變化對電磁場分量值的影響產生明顯改變，將時間的單位設為納秒

%反映電磁場性質的其他參數均設為常數，即自由空間中的值

clear

close?all

u0=4*pi*1e-7;??????????%自由空間中的磁導率

e0=1e-9/（36*pi）;???????%自由空間中的電介質常數

Z0=（u0/e0）^0.5;????????%自由空間中的波阻抗

f=1e8;?????????????????%電磁波的頻率

w=2*pi*f;?????????????

k=w*（u0*e0）^0.5;???????%波數

phi_E=0;???????????????%初始相位設為0

phi_H=0;

EE=20;??????????????????%電場幅度

HH=EE/Z0;

x=0:0.1:20;????????????%傳播方向上的采樣點

m0=zeros（size（x））;

gifname=‘mag_motion.gif‘;

figure

for?t=0:1:100???????????%為了消除波數與頻率之間的數量級帶來的影響，時間單位為ns

????Ez=EE*cos（k*x-w*t*1e-9+phi_E）;??%電場強度值

????Hy=HH*cos（k*x-w*t*1e-9+phi_H）;??%磁場強度值

????

????plot3（xm0Ez‘b‘‘LineWidth‘2）;??%繪制電場傳播曲線

????hold?on

????plot3（xHym0‘r‘‘LineWidth‘2）;??%繪制磁場傳播曲線

????hold?off

????xlabel（‘傳播方向‘）

????ylabel（‘磁場Hy‘）

????zlabel（‘電場Ez‘）?

????title（[‘平面電磁波傳播示意圖‘‘t=‘num2str（t）‘ns‘]‘fontsize‘14）

????set（gca‘fontsize‘12）

????drawnow

????frame=getframe（1）;%獲取當前圖像

????im=frame2im（frame）;

????[imindcm]=rgb2ind（im500）;%格式轉換

????if?t==0;

????????imwrite（imindcmgifname‘gif‘）;%創建一個gif文件

????else

????????imwrite（imindcmgifname‘gif‘‘WriteMode‘‘append‘‘DelayTime‘0.1）;

????????%向gif文件中添加一張圖片

????end

%??end

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????2303284??2017-11-16?15:13??mag_motion.gif

?????文件???????1674??2017-11-16?15:12??eleW3.m

-----------?---------??----------?-----??----

??????????????2304958????????????????????2