平行平面腔模的迭代解法-MATLAB.pdf

《平行平面腔模的迭代解法-MATLAB.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、平行平面腔模的迭代解法——MATLAB实现clear;clc;globallamdaLakNlamda=632.810^(-9);L=100lamda;a=25lamda;k=2pi/lamda;N=300;%将积分区间分成N小段turn=input('渡越次数=');x=linspace(-a,a,N+1);%x轴上等间距插入N+1个点u=ones(1,N+1);%初始条件fors=1:turnu=calc(u);endA=abs(u);%相对振幅phase=angle(u)/pi180;%转换成角度phase=phase-phase((

2、N+2)/2);%相对于中心位置处的相位差%---------------画图部分---------------%subplot(211),plot(x,A,'b');xlabel('x/m'),ylabel('相对振幅'),title('相对振幅分布');boxon,gridon;subplot(212),plot(x,phase,'r');xlabel('x/m'),ylabel('相对相位/度'),title('相对相位分布');boxon,gridon;%------------函数部分------------%unction[u2

3、]=calc(u1)globallamdaLakNu2=zeros(1,N+1);fors=1:N+1fort=0:Nu2(s)=u2(s)+sqrt(1i/lamda/Lexp(-1ikL))(2a/N)exp(-1ik/2/L(2a/N(s-1t))^2)u1(t+1);endendu2=u2/max(abs(u2));%归一化end%----------------渡越次数=300------------------%12