submain0.m

用matlab数值分析激光自再现模的FOX-li的迭代分析

function submain0(handles)
global lamda L T D U3 N dover;
%lamda=0.001;
%L=1000;               %腔长
%D=4;                 %腔镜直径
k=2*pi/lamda;
%N=150;               %精度相关
%T=100;                %跌代次数
CH=ones(N,N);
for m=1:N
    CH(m,:)=(m-1)*D/N:(-1)*D/N:(m-N)*D/N;
end
CH=CH.^2;
CH=exp(-i*k*CH/2/L);
CH=sqrt(i*exp(-i*k*L)/lamda/L)*CH;
Ci=fix((N+1)/2);      %中心坐标
U10=ones(1,N);      %初始面函数
for n=1:N
    U10(n)=1+0*i;
end
U1=U10;
for n=1:T
    if dover==0
U2=U1*CH;
U2=U2./U2(Ci);
U1=U2;
U21=abs(U2);            %光强
plot(handles.axes1,U21);
set(handles.axes1,'ylim',[0 1]);
U22=angle(U2);
plot(handles.axes2,U22);               %输出复振幅
jdstr=num2str(fix(100*n/T));
set(handles.laser_jd,'string',jdstr);
drawnow;
    end
end
U3=U2;