test_w_l3.m

来自「激光器典型腔计算」· M 代码 · 共 11 行

%syms l3;
%s=solve('5/8589934592*527702830840913505^(1/2)*(abs(44436/21025+5846/1576875*l3)/(-118750+2575*l3-4*l3^2)^(1/2))^(1/2)=25e-3','l3')
l3=51:0.5:800;
wsesam=25/4294967296*78852147137147995.^(1/2)./(-35766250+723775.*l3-169.*l3.^2).^(1/4);
wcrytal=25/8589934592*1055405661681827010^(1/2).*(abs(-238/2523+496/63075.*l3)./(-35766250+723775.*l3-169.*l3.^2).^(1/2)).^(1/2);
ws=1000*wsesam;         %让单位化成um
wc=1000*wcrytal;        %让单位化成um
plot(l3,ws,':')
hold on;
plot(l3,wc)
xlabel('l3(mm)'),ylabel('w(um)')