【波导仿真】基于matlab圆柱形金属波导TM43模场分布【含Matlab源码 3562期】.zip

warning off; close all clear clc dbstop if error %% 参数设置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a = 1; m = 4; n = 3; E0 = 1; landa = 532*10^-9; c=3*10^8; Umn = 14.372; Kc = Umn/a; K=2*pi/landa; Kz = sqrt(K^2-Kc^2); %% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Phi = linspace(0,2*pi,100); R = linspace(0,1,100); Z = linspace(0,3,100); [r1,phi] = meshgrid(R,Phi); X1 = r1.*cos(phi); Y1 = r1.*sin(phi); z1=0.5*pi/Kz; Hr1 = -1i*m./r1.*besselj(m,Kc*r1).*sin(m*phi)*exp(-1i*Kz*z1); Hr1=abs(real(Hr1)); surf(X1,Y1,X1*0+3,Hr1) shading interp hold on [z,phi] = meshgrid(Z,Phi); r = 1; X2 = r*cos(phi); Y2 = r*sin(phi); z2 = 3*2*pi/Kz*z/3; Hr2 = -1i*m./r.*besselj(m,Kc*r).*sin(m*phi).*exp(-1i*(Kz*z2-0.5*pi)); Hr2 = abs(real(Hr2)); surf(X2,Y2,z,Hr2) shading interp colormap gray theAxes=axis; for j=1:20; Hr1 = -1i*m./r1.*besselj(m,Kc*r1).*sin(m*phi)*exp(-1i*(Kz*z1-j*pi/20)); Hr1 = abs(real(Hr1)); Hr2 = -1i*m./r.*besselj(m,Kc*r).*sin(m*phi).*exp(-1i*(Kz*z2-0.5*pi-j*pi/20)); Hr2 = abs(real(Hr2)); surf(X1,Y1,X1*0+3,Hr1); hold on surf(X2,Y2,z,Hr2); shading interp colormap gray axis(theAxes) fmat(:,j)=getframe; end movie(fmat,10)