西南交通大学MATLAB实习任务17.doc

任 务：_ 17 （任务号_组号）03___ 课程名称： MATLAB 　 课程代码： 9990267 实习周数：　　　　　　1　　　 　　　 学　　分： 1 实习单位：　　　物理学院　　　　　　　 　　　　 实习地点：　　图书馆机房　　　　　　　　 实习时间：　　　15:00 至 17：00 实习的目的和任务： 实习目的 编程绘制波导内的TE波空间分布； 掌握quiver（）、pcolor（）绘图语句的使用。 掌握动画制作语句的使用。 实习任务 建立矩形波导内的TE波电场磁场方程（包括Ex，Ey，Hx，Hy，Hz）。 利用MATLAB的函数，绘制TE10模和TE11模电场和磁场在xoy、yoz、xoz平面上的场分布（共3*3*2-4=14个图，不用画XOZ和YOZ面的电场矢量图）。 生成TE10和TE11模式电场Ey分布动态图（TE10_Ey.avi,TE11_Ey.avi，并在动画中写入作者学号的和姓名）。 提示： 波导内电场（x分量）：Ex=1j*w*mu*Hm/kc^2*(n*pi/b)*cos(m*pi/a.*x).*sin(n*pi/b*y).*exp(-1j*kz.*z); 作二维色图：subplot(2,1,1); pcolor(x,y,Ex); title(Ex) shading interp 作矢量图：quiver(x,y,Ex,Ey) 参考图： t=T/4 TE21模式-XOY面（3个图） 场大小 场矢量 二、实习主要内容：（仿真结果及分析，含源程序＋注释＋图） %TE10在XOY平面上的分布图形 H0=1;m=1;n=0;mu0=pi*4.e-7;eps0=8.85e-12;mu=mu0;eps=eps0; %参数设定 f=1e8;Tp=1/f;w=2*pi*f;v=1./sqrt(eps*mu);lambda=v/f; a=2*lambda; b=1*lambda; c=4*lambda; k=w/v;kc=sqrt((m*pi/a)^2+(n*pi/b)^2); kz=sqrt(k^2-kc^2); figure(1) y=b/4; for t=0:Tp/100:2*Tp %进入循环 z=0; %设定x,y,z的值 [x,y]=meshgrid(0:a/100:a,0:b/5:b); ejwt=exp(1j*w*t); Ey=-1j*w*mu*H0/kc^2*(m*pi/a)*sin(m*pi/a.*x).*cos(n*pi/b.*y).*exp(-1j*kz.*z);%公式计算 Ey=real(Ey*exp(1j*w*t)); subplot(2,1,1); %绘制电场分量Ey在XOY面的分布 pcolor(x,y,Ey) %绘制图形 shading interp title(Ey of XOY) % axis([0 a 0 c -1*(w*mu*Hm/kc^2)*(m*pi/a) (w*mu*Hm/kc^2)*(m*pi/a)]) Ex=1j*w*mu*H0/kc^2*(n*pi/b)*cos(m*pi/a.*x).*sin(n*pi/b.*y).*exp(-1j*kz.*z); Ex=real(Ex*exp(1j*w*t)); subplot(2,1,2) %绘制电场分量Ey在平面的分布 pcolor(x,y,Ex) shading interp title(Ex of XOY) hold on figure %绘制第二个图形 Hx=1j*kz*H0/kc^2*m*pi/a*sin(m*pi/a*x).*cos(n*pi/b.*y).*exp(-1.j*kz*z); %磁场公式1 Hx=real(Hx*exp(1j*w*t)); subplot(3,1,1) %绘制磁场分量Hx在XOY面的分布 pcolor(x,y,Hx) shading interp title(Hx of XOY) Hy=1j*kz*H0/kc^2*n*pi/b