w_f_2.m

激光晶体热透镜效应相关计算

%记住粘贴过来得表达式中要加“.”，可以用查找替换功能来实现，注意要把光标放在最前面
clear;
f=8:0.5:800;
wcrytal=1./4294967296.*6065549779780615.^(1./2).*(abs(-878085209./16462575./f+198643./12615)./((878085209./16462575./f-198643./12615).*(-1516./2175./f-1./25-(-2432369./189225./f+1516./2175)./f)-(273638./189225./f+361./4350+361./174.*(-2432369./189225./f+1516./2175)./f).^2).^(1./2)).^(1./2)
wc=1000*wcrytal;        %让单位化成um
plot(f,wc)
xlabel('f(mm)'),ylabel('w(um)')