[原创]RP Fiber Power仿真设计掺铥光纤激光器代码详解 [复制链接]

(* aCG rS{  
Demo for program"RP Fiber Power": thulium-doped fiber laser, UQl?
_[G  
pumped at 790 nm. Across-relaxation process allows for efficient ?v2_7x&  
population of theupper laser level. 5|H;%T3_  
*)            ！（*  *）注释语句 ';I}6N  
46C%at M0}  
diagram shown: 1,2,3,4,5  ！指定输出图表 t<!m4Yd|#  
; 1: "Powersvs. Position"     ！分号是注释；光纤长度对功率的影响 |bnd92fvks  
; 2:"Variation of the Pump Power"  ！泵浦光功率变化对信号输出功率的影响 Z<vz%7w  
; 3:"Variation of the Fiber Length"！信号输出功率vs 光纤长度的变化，仿真最佳光纤长度 Q0
J1"*P0  
; 4:"Transverse Profiles"             ！横向分布，横坐标为半径位置 Px_8lB/;  
; 5:"Transition Cross-sections"    ！不同波长的跃迁横截面，横坐标波长，纵坐标为横截面 Wmbc `XC  
K"!rj.Da  
include"Units.inc"         ！读取“Units.inc”文件中内容 u+)!C*ho  
i2qN 0?n  
include"Tm-silicate.inc"    ！读取光谱数据 l#0
zHBc  

Thggas,  
; Basic fiberparameters:    ！定义基本光纤参数 yhH2b:nY(9  
L_f := 4 { fiberlength }      ！光纤长度 <-`bWz=+  
No_z_steps := 50 {no steps along the fiber } ！光纤步长，大括号{ }是注释，相当于备注 ioi/`iQR  
r_co := 6 um { coreradius }                ！纤芯半径 nDrRK  
N_Tm := 100e24 { Tmdoping concentration }  ！纤芯Tm离子掺杂浓度  ^*P?gG  
voZaJ2ho/O  
; Parameters of thechannels:                ！定义光信道 X_PzK'#m  
l_p := 790 nm {pump wavelength }                ！泵浦光波长790nm AH,?B*zGj  
dir_p := forward {pump direction (forward or backward) }   ！前向泵浦 Zr}>>aIJ]k  
P_pump_in := 5 {input pump power }                    ！输入泵浦功率5W =9&2udV1
 
w_p := 50 um {radius of pump cladding }               ！包层泵浦相应的半径 50um %`~4rf"7  
I_p(r) := (r <=w_p) { pump intensity profile }          ！泵浦光强度分布 ]D6<6OB  
loss_p := 0 {parasitic losses of pump wave }           ！泵浦光寄生损耗为0 7] y3<t  
I]ej ]46K  
l_s := 1940 nm {signal wavelength }                   ！信号光波长1940nm SrA6}kS  
w_s := 7 um                          ！信号光的半径 {tVA(&\<  
I_s(r) := exp(-2 *(r / w_s)^2)            ！信号光的高斯强度分布 {A}T^q!m]  
loss_s := 0                            ！信号光寄生损耗为0 q
;SD+%tI
 
S%RxYJ(  
R_oc := 0.70 {output coupler reflectivity (right side) }      ！输出耦合反射率 U'jmgHq  
l"8g9z  
; Function for defining themodel:   ！定义模型函数，一定要有calc命令，否则函数只会被定义，但不会被执行 B~o3Z  
calc |mhKD#:  
  begin pll5m7[  
    global allow all;                   ！声明全局变量 ygTfQtN  
    set_fiber(L_f, No_z_steps, '');        ！光纤参数 iu{y.}?  
    add_ring(r_co, N_Tm); DvWBvs,  
    def_ionsystem();              ！光谱数据函数 ,$]m1|t@z  
    pump := addinputchannel(P_pump_in, l_p,'I_p', loss_p, dir_p);  ！泵浦光信道 K`=U5vG^  
    signal_fw := addinputchannel(0, l_s, 'I_s',loss_s, forward);      ！前向信号光信道 z{.&sr>+v  
    signal_bw := addinputchannel(0, l_s, 'I_s',loss_s, backward);    ！后向信号光信道 (fc_V[(m"  
    set_R(signal_fw, 1, R_oc);                                 ！设置反射率函数 qA}l[:F+#  
    finish_fiber();                                   kU:Q&[/jzH  
  end; FUI/ A>  
"P5,p"k:)  
; Display someoutputs in the Output window (on the right side): ！在Output aera区域显示输出 Lh+7z>1  
show "Outputpowers:"                                   ！输出字符串Output powers： pDC`Fi  
show"pump:     ", P_out(pump):d3:"W"  ！输出字符串pump：和计算值（格式为3个有效数字，单位W） &*Z)[Bl  
show"signal:   ",P_out(signal_fw):d3:"W" ！输出字符串signal：和计算值（格式为3个有效数字，单位W） 7\dt<VV  
:4:U\k;QwA  
w^~s4Q_>>  
; ------------- .:-*89c  
diagram 1:                   !输出图表1 V\Lh(zPt  
y=GD
uU%  
"Powers vs.Position"          ！图表名称 B \V;{:  
Ow0~sFz  
x: 0, L_f                      ！命令x: 定义x坐标范围 _Wp, z`  
"position infiber (m)", @x      ！x轴标签；@x 指示这些字符串沿坐标轴放置 Mt4`~`6  
y: 0, 15                      ！命令y: 定义y坐标范围 Ornm3%p+e  
y2: 0, 100                    ！命令y2: 定义第二个y坐标范围 H2_6m5[&,  
frame          ！frame改变坐标系的设置 &S^a_L:  
legpos 600, 500  ！图行在图表窗口中的位置（相对于左上角而言） (!ud"A|ab4  
hx             ！平行于x方向网格 4YMUkwh  
hy              ！平行于y方向网格 B)DtJf  
<t*3w  
f: P(pump, x),    ！命令f: 定义函数图；P(pump, x)函数是计算x位置处的泵浦光功率 Si;eBPFH  
  color = red,  ！图形颜色 `yXHb  
  width = 3,   ！width线条宽度 En+`ZcA\z  
  "pump"       ！相应的文本字符串标签 !;PKx]/&  
f: P(signal_fw, x),  ！P(signal_fw ,x) 函数是计算x位置处的前向信号光功率 )q+;+J`>  
  color = blue,     ?cr^.LV|h^  
  width = 3, ;;17 #T2  
  "fw signal" }8'bXG+  
f: P(signal_bw, x),   ！P(signal_bw ,x) 函数是计算x位置处的后向信号光功率 oF
Jx8XU  
  color = blue, JOq&(AZe  
  style = fdashed, M-f; ,>  
  width = 3, 0o!Egq_  
  "bw signal" _nW#Cl~  
WGMEZx  
f: 100 * n(x, 2),    ！n(x ,2) 函数是计算x位置处激活粒子数在能级2上的占比 nrZZkQNI  
  yscale = 2,            ！第二个y轴的缩放比例 a/ Z\h{*  
  color = magenta, 4q^'MZm1  
  width = 3, x_KJCU  
  style = fdashed, #.W<[KZf  
  "n2 (%, right scale)" X[iQ%Y$/n  
]$WwPDZ  
f: 100 * n(x, 3),          ！n(x ,3) 函数是计算x位置处激活粒子数在能级3上的占比 .37Jrh0Iv  
  yscale = 2, )1!*N)$  
  color = red, OM7EmMa;  
  width = 3, 64-;| k4F  
  style = fdashed, Ut^ {4_EC  
  "n3 (%, right scale)" eB*0})  
8d"Ff  
+/~;y{G..z  
; ------------- >FM2T<.;  
diagram 2:                    ！输出图表2 p{4nWeH?B  
:HQQ8uQfb  
"Variation ofthe Pump Power" u IGeSd5B  
W~FM^xR?p  
x: 0, 10 g O ;oM?|  
"pump inputpower (W)", @x ^8eu+E.{  
y: 0, 10 RW04>oxVn  
y2: 0, 100 /V@9!  
frame yd4\%%]  
hx '9d] B^)F  
hy {:0TiOP5x  
legpos 150, 150 }grel5lq  
TmIw?#q^  
f: (set_P_in(pump, x);P_out(signal_fw)), ！set_P_in(pump,x)改变泵浦信道功率；P_out(signal_fw)输出前向信号光 @is!VzE  
  step = 5, ' (1`iQ;  
  color = blue, K/Qo~  
  width = 3, %RS8zN  
  "signal output power (W, leftscale)",     ！相应的文本字符串标签 V 0M&D,  
  finish set_P_in(pump, P_pump_in) DYFfq  
cg]>*lH  
f: (set_P_in(pump,x); 100 * n_av(2)),   ！改变泵浦信号功率对能级2上激活粒子占比的影响 ~ '/Yp8(  
  yscale = 2, yz7X7mAo  
  step = 5, 5:sk&0:@U  
  color = magenta, Qaeg3f3F3  
  width = 3, L
aCVI  
  "population of level 2 (%, rightscale)", $"[5]{'J  
  finish set_P_in(pump, P_pump_in) &<RK=e'*x  
|`,AAa  
f: (set_P_in(pump,x); 100 * n_av(3)),   ！改变泵浦信号功率对能级3上激活粒子占比的影响 GLE/ 1  
  yscale = 2, J70r`   
  step = 5, 9O2??N7f  
  color = red, JRYCM}C]  
  width = 3, J':x]_;  
  "population of level 3 (%, rightscale)", l1I
\kh
S  
  finish set_P_in(pump, P_pump_in) <tx`#,  
x>TIQU=\  
A$o7<Hx  
; ------------- !l5&>1?  
diagram 3:                         ！输出图表3 @b{I0+li"/  
'?)<e^  
"Variation ofthe Fiber Length" %&}gt+L(M  
&r0
U9J  
x: 0.1, 5 P?/Mrz   
"fiber length(m)", @x eB2a1<S&@  
y: 0, 10 _y5J]Yu`j  
"opticalpowers (W)", @y K$kI%eGZA  
frame X(M|T]`b:  
hx s)9d\{  
hy >\4"k4d}  
we}G%09L  
f: (set_L(x);P_out(signal_fw)),     ！改变光纤长度对信号光输出功率的影响 -a\[`JHi  
  step = 20,             %ki^XB86  
  color = blue, _?rL7oTv  
  width = 3, !&#
5*  
  "signal output" ]g
jB%R[.m  
8'|
_O  
;f: (set_L(x);P_out(pump)),                     ！改变光纤长度对泵浦信号输出功率的影响 hA:RVeS{  
   step = 20, color = red, width = 3,"residual pump" ly( LMr  
Zt/4|&w  
! set_L(L_f) {restore the original fiber length } 2`P=ekF]  
WlW7b.2.  
, G9{:  
; ------------- C7Ny-rj}IA  
diagram 4:                                  ！输出图表4 0f9*=c  
RcpKv;=iB  
"TransverseProfiles" hmp!|Q[)  
x.kIzI5  
I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) %Fp1c K  
XpIl-o&re  
x: 0, 1.4 * r_co /um "(+p1  
"radialposition (µm)", @x `BzjDI:a  
y: 0, 1.2 * I_max *cm^2 \$W\[s4I  
"intensity (W/ cm&sup2;)", @y 05s{Z.aK  
y2: 0, 1.3 * N_Tm Q/]t$  
frame ~ya@ YP]';  
hx ')zf8>,  
hy d:3OC&  
sg%Ptp  
f: N_dop(1, x * um,0),      ！掺杂浓度的径向分布 t~_bquGk  
  yscale = 2, y134m  
  color = gray, we&D
"V  
  width = 3, s_#6^_  
  maxconnect = 1, +nz0ZQ9 a  
  "N_dop (right scale)" Ex-?[Hq  

'n/L1Fn  
f: I(pump, -1, x *um, 0) * cm^2,    ！泵浦光沿光纤径向的强度分布 m5 l,Lxj  
  color = red, .1YiNmW=  
  maxconnect = 1,           ！限制图形区域高度，修正为100%的高度 %4^NX@1jV  
  width = 3, wu~?P`  
  "pump" <~}NxY\5  
(-%1z_@Y  
f: I(signal_fw, -1,x * um, 0) * cm^2,  ！信号光沿光纤径向的强度分布 Dks"(0g  
  color = blue, Uu0  
  maxconnect = 1, LO@o`JF  
  width = 3, l]~mB~  
  "signal" '#NDR:J"  
$ #bWh  
QC<O=<$Q[  
; ------------- Jxsch\  
diagram 5:                                  ！输出图表5 BPd *@l  
"
5@\"L  
"TransitionCross-sections" ,b9!\OWDF  
")i_{C,b^  
I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) (w1$m8`=  
MeDlsO  
x: 1450, 2050 &~V6g(9  
"wavelength(nm)", @x +o|I@7f  
y: 0, 0.6 #b>D^=NV>)  
"cross-sections(1e-24 m&sup2;)", @y zbAyYMtEk  
frame W.nr&yiQ  
hx mWTV)z57  
hy UO4
z~  
#k|f%!-Vo  
f: s12_Tm(x * nm) /1e-24,      ！Tm³⁺吸收截面与波长的关系  ?)2;W  
  color = red, 4uE|$  
  width = 3, O"9Or3w  
  "absorption" go?}M]c%7  
f: s21_Tm(x * nm) /1e-24,  ！Tm³⁺发射截面与波长的关系 ;4k/h/o1#  
  color = blue, hxkwT  
  width = 3, #L+ZH
s~  
  "emission" 85vyt/.,k  
Bk@&k}0  

感谢，视频上有点看不清楚