w_f_pingaoqiang.m

来自「激光晶体热透镜效应相关计算」· M 代码 · 共 7 行

%记住粘贴过来得表达式中要加“.”，可以用查找替换功能来实现，注意要把光标放在最前面
clear;
f=9:0.5:800;
wcrytal=1./4294967296.*6065549779780615.^(1./2).*(abs(-56167123./1317006./f+790889./50460)./((56167123./1317006./f-790889./50460).*(-589./870./f-1./25-(-409979./30276./f+589./870)./f)-(80693./75690./f+137./2175+137./87.*(-409979./30276./f+589./870)./f).^2).^(1./2)).^(1./2);
wc=1000.*wcrytal;        %让单位化成um
plot(f,wc)
xlabel('f(mm)'),ylabel('w(um)')