激光物理学第一次作业

global interval L k a lmlm=632.8e-9 %波长L=100*lm    %腔长a=25*lm     %腔镜n=300;k=2*pi/lm; interval=1000;x=linspace(-a,a,interval);u_=ones(1,interval);for m=1:n    for mm=1:interval        u0(mm)=fun(x(mm),u_);    end;    u_=u0/max(abs(u0));endsubplot(2,1,1)plot(x,abs(u0)/abs(u0(interval/2)))xlabel('x');ylabel('相对振幅');angle_u0=angle(u0)/pi*180;angle_u0=angle_u0-angle_u0(interval/2);subplot(2,1,2)plot(x,angle_u0)xlabel('x');ylabel('相对相位/(°C)');function y=fun(x,u)global interval L k a lmx_=linspace(-a,a,interval);step_length=2*a/(interval-1);y=sqrt(1i/(L*lm)*exp(-1i*k*L))*sum(exp(-1i*k/2/L*(-x_+x).^2).*u)*step_length;