1139阅读
0回复

[求助]请教：论坛上下载的一个波导程序

UID ：33555

注册：2009-05-25
登录：2010-04-25
发帖：56
等级：仿真一级

0楼发表于: 2010-04-13 09:07:25

最近在论坛上下载了一个波导的程序，有些地方不太明白——%变量IM、JN，IS、JS，I0、J0的作用是什么？？，希望高手指点！谢谢！
c=2.99792458e8; ;
pi=3.14159267;
N=20000;
IE=60;
JE=30;
IB=IE+1;
JB=JE+1;
a=7.112e-3;
b=3.556e-3;
dl=a/IE;
IM=4;
JN=2;
IS=fix(a/dl/(IM+JN-1));
JS=fix(b/dl/(IM+JN-1));
I0=5;
J0=4;            %变量IM、JN，IS、JS，I0、J0的作用是什么？？
s=1/2;
dt=dl*s/c
T=15*dt;
t0=140*dt;
muxx=ones(IB,JE);
muyy=ones(IE,JB);
epszz=ones(IB,JB);

hx=zeros(IB,JE);
hy=zeros(IE,JB);

ez=zeros(IB,JB);

EEZ=zeros(1,N);
x=zeros(1,N);

for n=1:N;

for ii=2:IE;
     for jj=2:JE;
        ez(ii,jj)=ez(ii,jj)+( s/epszz(ii,jj) ).*( hy(ii,jj)-hy(ii-1,jj)+...
                                                         hx(ii,jj-1)-hx(ii,jj) );
     end;
   end;
   EEZ(n)=ez(I0,J0);


   hx(IS,JS)=exp(-((n*dt-t0)/T)^2);
   hy(IS,JS)=exp(-((n*dt-t0)/T)^2);

   for ii=2:IE;
     for jj=1:JE;
             hx(ii,jj)=hx(ii,jj)-(s/muxx(ii,jj)).*( ez(ii,jj+1)-ez(ii,jj) );
     end;
   end;
   for ii=1:IE;
     for jj=2:JE;
             hy(ii,jj)=hy(ii,jj)+(s/muyy(ii,jj)).*( ez(ii+1,jj)-ez(ii,jj) );
     end;
   end;
HHX(n)=hx(I0,J0);
HHY(n)=hy(I0,J0);
mesh(ez)
drawnow;
end;
hold on
x=fft(HHX);
x=abs(x);
n=1:N;
f=n/(dt*N);
fc=c/(2*a)
% f=f/fc;
plot(f(1:N/40),x(1:N/40));
xlabel('Frequency(Hz)/fc-TE10');
ylabel('Magnitude');
%calculate the cutoff frequency
position=zeros(1,20);%record the position of cut-off frequency
freq=zeros(1,20);%record the data of cut-off frequency
m=1;
for i=1:N/40-2;%find the maxium of amplitude
    x1=x(i);
    x2=x(i+1);
    x3=x(i+2);
    xx1=(x ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

发帖回复