loading
A=1;B=0;C=-2/R1;D=1; q1=(A*qm+B)/(C*qm+D); z0=real(q1);f=imag(q1); w0=sqrt(f*wl/pi); z=0:L1;
w=w0*sqrt(1+(z+z0).^2/f^2); plot(z,w),hold on
A=1;B=L1;C=0;D=1; q2=(A*q1+B)/(C*q1+D);
A=1;B=0;C=-1/F;D=1; q3=(A*q2+B)/(C*q2+D); z0=real(q3);f=imag(q3); w0=sqrt(f*wl/pi); z=0:L2;
w=w0*sqrt(1+(z+z0).^2/f^2); plot(z+L1,w),hold on
A=1;B=L2;C=0;D=1; q4=(A*q3+B)/(C*q3+D);
A=1;B=0;C=-2/R2;D=1; q5=(A*q4+B)/(C*q4+D); z=0:L2;
w=w0*sqrt(1+(z+z0).^2/f^2); plot(z+L1,w),hold on
A=1;B=L2;C=0;D=1; q6=(A*q5+B)/(C*q5+D);
A=1;B=0;C=-1/F;D=1; q7=(A*q6+B)/(C*q6+D); z0=real(q1);f=imag(q1); w0=sqrt(f*wl/pi); z=0:L1;
w=w0*sqrt(1+(z+z0).^2/f^2); plot(z,w),hold on
A=1;B=L1;C=0;D=1; q8=(A*q7+B)/(C*q7+D);
title('光斑半径与z的关系');xlabel('透镜光焦度D（/mm）');ylabel('光束半径')
运行结果：
题目四
取不同的l1值和R1值，计算谐振腔的稳定性，输出镜M2和透镜上的模式半径与光焦度1/F的关系。
由于此处matlab程序与前雷同，只显示结果。 （1） 取L1=200,R1=1000;
（2） 取L1=200,R1=-1000;