矩形波导matlab求解.docxVIP

矩形波导matlab求解题目： 1、已知宽度为50mm, 高度为34mm的矩形波导管，长度为5m。忽略边缘效应，该波导管中的场满足方程：其中，m, n为波沿x轴和y轴的模式画出：1. m=1,n=0 2. m=2, n=3 时波导管中的 场分布1、求解定解问题2、用Matlab（C语言/Java, etc.）编程，画出求解区域中的场二维分布图，横坐标为x轴/y轴/z轴，纵坐标为电位值求解方程对于TM波，Hz=0,Ez满足如下关系：Ez=Eoz(x, y)e(-jβz)在直角坐标系中,上式可写作应用分离变量法, 令Eoz(x, y)=X(x)Y(y)代入上式, 并除以X(x)Y(y), 得：要使上式成立, 上式左边每项必须均为常数, 设分 和可得于是, Eoz(x, y)的通解为：Eoz(x,y)=(A1coskxx+A2sinkxx)(B1coskyy+B2sinkyy)其中, A1A2B1B2为待定系数, 由边界条件确定。Ez应满足的边界条件为：代入可得于是方程的基本解为截止波长λc=2*pi/(((m*pi/a)^2+(n*pi/b)^2)^0.5)截止频率fc=c/λc 无尽空间中波长λ0=3*10^8/f（f为工作频率）波导波长λg=l0/((1-(3*10^8/(lc*f))^2)^0.5)相位常数B=2*pi/lgw=B*(3*10^8)由激励源强度决定可得TM波其他横向用matlab画出场分布在matlab中新建m文件function n(ao,bo,d,Emn,f,m,n,t)%画矩形波导场结构 所有计算单位为米 a0,b0输入为毫米 %f工作频率/波长%lg 波导波长 %lc TM11模截止波长%a b 波导尺寸%d 采样精度 %t t时刻的场结构图a=ao/1000;b=bo/1000;lc=2*pi/(((m*pi/a)^2+(n*pi/b)^2)^0.5); l0=3*10^8/f;u=1/(36*pi)*10^(-9);lg=l0/((1-(3*10^8/(lc*f))^2)^0.5);c=lg;B=2*pi/lg;w=B*(3*10^8);x=0:a/d:a;y=0:b/d:b;z=0:c/d:c;[x1,y1,z1]=meshgrid(x,y,z);Ex=B./((m*pi./a)^2+(n*pi./b)^2).*(m*pi./a).*Emn.*cos(m*pi./a.*x1).*sin(n*pi./b.*y1).*sin(w*t-B.*z1);Ey=B./((m*pi./a)^2+(n*pi./b)^2).*(n*pi./b).*Emn.*sin(m*pi./a.*x1).*cos(n*pi./b.*y1).*sin(w*t-B.*z1);Ez=Emn.*sin(m*pi./a.*x1).*sin(n*pi./b.*y1).*cos(w*t-z1.*B);quiver3(z1,x1,y1,Ez,Ex,Ey,b);hold on;x2=x1-0.05;y2=y1-0.05;z2=z1-0.05;Hz=zeros(size(z2));Hy=-w.*u./((m*pi./a)^2+(n*pi./b)^2).*(m*pi./a).*Emn.*cos(m*pi./a.*x2).*sin(n*pi./b.*y2).*sin(w*t-B.*z2);Hx=-w.*u./((m*pi./a)^2+(n*pi./b)^2).*(n*pi./b).*Emn.*sin(m*pi./a.*x2).*cos(n*pi./b.*y2).*sin(w*t-B.*z2);quiver3(z2,x2,y2,Hz,Hx,Hy,r);xlabel(传输方向);ylabel(波导宽边a);zlabel(波导窄边b);hold off;End当m=1,n=1时在命令框输入n(50,34,10,1,5.5*10^9,1,1,0)场分布如图当m=2,n=3时输入如下命令n(50,34,10,1,5*10^10,2,3,0)