利用Matlab实现矩形波导电磁场分布图的绘制
利 用 M a t l a b 实 现 矩 形 波 导 电 磁 场 分 布 图 的 绘 制 附 源 程 序 通过 Matl ab 计算并绘出任意时刻金属矩形波导的主模 TE10 模的电磁场分布图。波导尺寸、工作频率及时刻均由外部给定。 A. 矩形波导中传输的主模为 TE10 模。设金属波导尺寸为 a*b TE10 模的截止波长为2*a 。其电磁场分量可推导表示如下 000s i n s i nc o ss i n s i nxzyaH H x t zaH H t zaE H x t za 1-1 上式中各参量如下 2212cgcgacfc 1-2 B. 用Matl ab 画电磁力线的步骤 1. 由外部给定的波导尺寸、工作频率参照1-2 式计算得到参量。 2. 由外部给定的绘图精度分别确定电场和磁场的坐标点。按照公式1-1 计算得到电场、磁场的分量。 3. 用qui ver3 函数绘制磁场分布。允许图像叠加。 4. 用qui ver3 函数绘制电场分布。不允许图像叠加。 C. 三维的电力磁力线分布效果图 图1 图2 C. 附程序清单 rectwavestrct1(22. 86, 10. 16, 6, 1, 9. 84*109, 0. 03); %mai n functi on rectwavestrct1(ao, bo, d, H0, f, t) % 画矩形波导场结构 所有计算单位为米 输入为毫米 %f l 0 工作频率/ 波长 %l g 波导波长 %l c TE10 模截止波长 %a b 波导尺寸 %c 传输方向 这里取为波导波长 %d 采样精度 %t t 时刻的场结构图 a=ao/1000; b=bo/1000; l c=2*a; %TE10 截止频率 l 0=3*108/f; u=4*pi *10(-7); i f(l 0l c) return; el se cl f; l g=l 0/(1-(l 0/l c)2)0. 5); c=l g; B=2*pi /l g; w=B/(3*108); x=0: a/d: a; y=0: b/d: b; z=0: c/d: c; x1, y1, z1 =meshgri d(x, y, z); %mesh(x1, y1, z1); hx=-B. *a. *H0. *si n(pi . /a. *x1). *si n(w*t-B. *z1). /pi ; hz=H0. *cos(pi . /a. *x1). *cos(w*t-z1. *B); hy=zeros(si ze(y1); qui ver3(z1, x1, y1, hz, hx, hy, b ); hol d on; x2=x1-0. 001; y2=y1-0. 001; z2=z1-0. 001; ex=zeros(si ze(x2); ey=w. *u. *a. *H0. *si n(pi . /a. *x2). *si n(w*t-B. *z2). /pi ; ez=zeros(si ze(z2); qui ver3(z2, x2, y2, ez, ex, ey, r ); xl abel ( 传输方向 ); yl abel ( 波导宽边a ); zl abel ( 波导窄边b ); hol d off; end %- - - - - - - - - - - - - - - - - - - - - - - - - - End Code- - - - - - - - - - - - - -