pingao_longcrystal_3w.m

激光器典型腔计算

%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')
clear;
l3=1:0.5:60;
wquarter=5./4294967296.*89397865824936929.^(1./2).*(abs(-16345./1682-20./29.*l3-(234./725-2./25.*l3).*l3)./(97734+10643.*l3-841.*l3.^2).^(1./2)).^(1./2);
wr=5./4294967296.*89397865824936929.^(1./2).*(abs(-7886./841-17./29.*l3-(309./725-2./25.*l3).*l3)./(97734+10643.*l3-841.*l3.^2).^(1./2)).^(1./2);
wl=5./4294967296.*89397865824936929.^(1./2).*(abs(-(9./725-2./25.*l3).*l3-l3-358./29)./(97734+10643.*l3-841.*l3.^2).^(1./2)).^(1./2);

wquarter=1000.*wquarter;         %让单位化成um
wr=1000.*wr;
wl=1000.*wl;        %让单位化成um

plot(l3,wquarter,'g')
hold on;
plot(l3,wr,'b')
hold on;
plot(l3,wl,'y')
xlabel('l3(mm)'),ylabel('w(um)')