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

(* aw%v
u  
Demo for program"RP Fiber Power": thulium-doped fiber laser, 42wC."A  
pumped at 790 nm. Across-relaxation process allows for efficient 4eL54).1O  
population of theupper laser level. _XNR um4  
*)            ！（*  *）注释语句 yi-"hT`  
`yXx[deY  
diagram shown: 1,2,3,4,5  ！指定输出图表 o+Ti$`2<O7  
; 1: "Powersvs. Position"     ！分号是注释；光纤长度对功率的影响 G:C6`uiy`  
; 2:"Variation of the Pump Power"  ！泵浦光功率变化对信号输出功率的影响 He-Ja  
; 3:"Variation of the Fiber Length"！信号输出功率vs 光纤长度的变化，仿真最佳光纤长度 cW\Y?x   
; 4:"Transverse Profiles"             ！横向分布，横坐标为半径位置 AYu'ptDNr  
; 5:"Transition Cross-sections"    ！不同波长的跃迁横截面，横坐标波长，纵坐标为横截面 ^ R^N`V   
HT=Am  
include"Units.inc"         ！读取“Units.inc”文件中内容 JKi@Kw  
9iddanQA  
include"Tm-silicate.inc"    ！读取光谱数据 vJAAAS  
gXLZ)>+A+  
; Basic fiberparameters:    ！定义基本光纤参数 $]U5  
L_f := 4 { fiberlength }      ！光纤长度 K-k.=6mS  
No_z_steps := 50 {no steps along the fiber } ！光纤步长，大括号{ }是注释，相当于备注 })Mv9~&S  
r_co := 6 um { coreradius }                ！纤芯半径 fQTA@WAr  
N_Tm := 100e24 { Tmdoping concentration }  ！纤芯Tm离子掺杂浓度 ;E? hz  
n5*{hi  
; Parameters of thechannels:                ！定义光信道 2T(,H.O  
l_p := 790 nm {pump wavelength }                ！泵浦光波长790nm +9 p`D  
dir_p := forward {pump direction (forward or backward) }   ！前向泵浦 XvkI+c  
P_pump_in := 5 {input pump power }                    ！输入泵浦功率5W c]!D`FA*K  
w_p := 50 um {radius of pump cladding }               ！包层泵浦相应的半径 50um X)iWb(@k"7  
I_p(r) := (r <=w_p) { pump intensity profile }          ！泵浦光强度分布 .6 ?>t!&W  
loss_p := 0 {parasitic losses of pump wave }           ！泵浦光寄生损耗为0 a.&#dxgW[  
VfA5r`^  
l_s := 1940 nm {signal wavelength }                   ！信号光波长1940nm NuS|X   
w_s := 7 um                          ！信号光的半径 E,D:D3O  
I_s(r) := exp(-2 *(r / w_s)^2)            ！信号光的高斯强度分布 ~eV!!38 J  
loss_s := 0                            ！信号光寄生损耗为0 +=*m! 7Mr  
E%\7Uo-  
R_oc := 0.70 {output coupler reflectivity (right side) }      ！输出耦合反射率 b~Pxgfu"  
0.BUfuuh  
; Function for defining themodel:   ！定义模型函数，一定要有calc命令，否则函数只会被定义，但不会被执行 ^Im%D(MY  
calc ogc('HqF^'  
  begin Wt)SdF=U/  
    global allow all;                   ！声明全局变量 +65~,e  
    set_fiber(L_f, No_z_steps, '');        ！光纤参数 q'Pz3/mk  
    add_ring(r_co, N_Tm); 3+
8"  
    def_ionsystem();              ！光谱数据函数 |h=+&*(:  
    pump := addinputchannel(P_pump_in, l_p,'I_p', loss_p, dir_p);  ！泵浦光信道 eFj6p
<  
    signal_fw := addinputchannel(0, l_s, 'I_s',loss_s, forward);      ！前向信号光信道 Y9@dZw%2  
    signal_bw := addinputchannel(0, l_s, 'I_s',loss_s, backward);    ！后向信号光信道 X!#i@V  
    set_R(signal_fw, 1, R_oc);                                 ！设置反射率函数 rv%ye H  
    finish_fiber();                                   D"x$^6`c}  
  end; A{8K#@!  
!G"9xrr1  
; Display someoutputs in the Output window (on the right side): ！在Output aera区域显示输出 @X|i@{<';  
show "Outputpowers:"                                   ！输出字符串Output powers： o,Tr^e$  
show"pump:     ", P_out(pump):d3:"W"  ！输出字符串pump：和计算值（格式为3个有效数字，单位W） !,6v=n[Nz  
show"signal:   ",P_out(signal_fw):d3:"W" ！输出字符串signal：和计算值（格式为3个有效数字，单位W） DNW2;i<hsz  
R0hctT1j  
\Ul.K!b7  
; ------------- fr'huvc  
diagram 1:                   !输出图表1 ZH;y>Z  
QLB1:O>  
"Powers vs.Position"          ！图表名称 # +QWi0B  
:<t{ =0G  
x: 0, L_f                      ！命令x: 定义x坐标范围 Vh#Mp!  
"position infiber (m)", @x      ！x轴标签；@x 指示这些字符串沿坐标轴放置 pg?i
F1  
y: 0, 15                      ！命令y: 定义y坐标范围 r}y[r}vk  
y2: 0, 100                    ！命令y2: 定义第二个y坐标范围 L'\/)!cEd  
frame          ！frame改变坐标系的设置 GIRSoRVsh  
legpos 600, 500  ！图行在图表窗口中的位置（相对于左上角而言） ] [HGzHA  
hx             ！平行于x方向网格 0Y6q$h>4  
hy              ！平行于y方向网格 KCJN<  
|,1bkJt  
f: P(pump, x),    ！命令f: 定义函数图；P(pump, x)函数是计算x位置处的泵浦光功率 _Vo)<--+I  
  color = red,  ！图形颜色 W) 33;E/}  
  width = 3,   ！width线条宽度 0tW<LR-}E  
  "pump"       ！相应的文本字符串标签 }#ep}h  
f: P(signal_fw, x),  ！P(signal_fw ,x) 函数是计算x位置处的前向信号光功率 ;Y@!:p-H  
  color = blue,     W8G9rB
|T  
  width = 3, |ilv|UV  
  "fw signal" tIyuzc~U  
f: P(signal_bw, x),   ！P(signal_bw ,x) 函数是计算x位置处的后向信号光功率 TDAWI_83-  
  color = blue, yNrinYw  
  style = fdashed, Vedyy\TU  
  width = 3, X/E7o92\  
  "bw signal" (@KoqwVWc  
%_b^!FR  
f: 100 * n(x, 2),    ！n(x ,2) 函数是计算x位置处激活粒子数在能级2上的占比 u(1J=h  
  yscale = 2,            ！第二个y轴的缩放比例 H:q)^$s  
  color = magenta, N;Z`%&  
  width = 3, 2/yXY_L  
  style = fdashed, /B,:<&_-  
  "n2 (%, right scale)" ggm2%|?X  
tLD~  
f: 100 * n(x, 3),          ！n(x ,3) 函数是计算x位置处激活粒子数在能级3上的占比 KMK&[E#r  
  yscale = 2, zd*3R+>U'>  
  color = red, UT
DcX  
  width = 3, 6%g
B E  
  style = fdashed, 9^ >M>f"  
  "n3 (%, right scale)" ]g;^w?9h  
)|w*/JK\Z  
:W.jNV{e\F  
; ------------- {J,6iP{>ZN  
diagram 2:                    ！输出图表2 <#+oQ>5s  
JQ)4}t  
"Variation ofthe Pump Power"
5:Yck<  
~Na=+}.q_  
x: 0, 10 x],8yR)R  
"pump inputpower (W)", @x ?
w+ QbT  
y: 0, 10 !;o\5x<'$O  
y2: 0, 100 May&@x/oMS  
frame K<@[_W+  
hx AIZW@Nq.5  
hy H+4=|mkQ  
legpos 150, 150 Sao>P[#x  
$hG;2v  
f: (set_P_in(pump, x);P_out(signal_fw)), ！set_P_in(pump,x)改变泵浦信道功率；P_out(signal_fw)输出前向信号光 v!NB~"LQ  
  step = 5, "sF Xl  
  color = blue, Yv)Bj  
  width = 3, T*h!d(  
  "signal output power (W, leftscale)",     ！相应的文本字符串标签 J0^{,eY<  
  finish set_P_in(pump, P_pump_in) Y
]C;T  
nK+lE0  
f: (set_P_in(pump,x); 100 * n_av(2)),   ！改变泵浦信号功率对能级2上激活粒子占比的影响 n?'d|h  
  yscale = 2, `IEq@Wr#$!  
  step = 5, %ZujCZn  
  color = magenta, ,f}h}  
  width = 3, =2^Vgc  
  "population of level 2 (%, rightscale)", FE/$(7rM  
  finish set_P_in(pump, P_pump_in) BG'6;64kx6  
jja{*PZ6H  
f: (set_P_in(pump,x); 100 * n_av(3)),   ！改变泵浦信号功率对能级3上激活粒子占比的影响 ZlthYuJ  
  yscale = 2, SHMl%mw  
  step = 5, IE&_!ce  
  color = red, DdBxqkh  
  width = 3, PC*m% ?+  
  "population of level 3 (%, rightscale)", ~O \}/I28  
  finish set_P_in(pump, P_pump_in) (# JMB)  
hP?7zz$*j  
vG Y!4@[  
; ------------- Pj'62
[5z  
diagram 3:                         ！输出图表3 *"1~b
Pl  
hlX>K  
"Variation ofthe Fiber Length" 31WZJm^  
`SCy<w3$+[  
x: 0.1, 5 cp$.,V  
"fiber length(m)", @x \CcmePTN#x  
y: 0, 10 e#tIk;9Xz  
"opticalpowers (W)", @y @gNpJB]V  
frame g^1r0.Sp{8  
hx _3|6ZO  
hy A:/}`  
'<TD6jBs  
f: (set_L(x);P_out(signal_fw)),     ！改变光纤长度对信号光输出功率的影响 7'Lp8  
  step = 20,             l1&5uwuF  
  color = blue, UmJUt|  
  width = 3, &GuF\wJ{7  
  "signal output" VSh!4z1  
d7E7f  
;f: (set_L(x);P_out(pump)),                     ！改变光纤长度对泵浦信号输出功率的影响 XM$GQn]B  
   step = 20, color = red, width = 3,"residual pump" pjACFVMFX  
5? Wg%@  
! set_L(L_f) {restore the original fiber length } ]GNh)  
oew|23Ytb  
A^-iHm  
; ------------- <KtBv Ip]  
diagram 4:                                  ！输出图表4 h6g:(3t6m  
6#E7!-u(-  
"TransverseProfiles" L_^`k4ct  
3!aEClRtq  
I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) +$PFHXB  
X<Z(]`i  
x: 0, 1.4 * r_co /um Vb2\/e:k  
"radialposition (µm)", @x !nwbj21%  
y: 0, 1.2 * I_max *cm^2 Rb#/qkk/  
"intensity (W/ cm&sup2;)", @y HcA[QBh  
y2: 0, 1.3 * N_Tm ha5e(Hj?  
frame V{0%xz #  
hx RK-x?ZYH'  
hy B%fU'  
?;^5ghY$  
f: N_dop(1, x * um,0),      ！掺杂浓度的径向分布 4f!dYo4L  
  yscale = 2, UeA2c_ 5  
  color = gray, .3
@Ng  
  width = 3,  //<:k8  
  maxconnect = 1, ;}D-:J-z_  
  "N_dop (right scale)" JA<~xo[Q9  
S;8.y
j-  
f: I(pump, -1, x *um, 0) * cm^2,    ！泵浦光沿光纤径向的强度分布 ugI#ZFjJWE  
  color = red, +#~O'r]%GG  
  maxconnect = 1,           ！限制图形区域高度，修正为100%的高度 !Av1Leb9$  
  width = 3, [t{](-  
  "pump" 8"i/wMP]  
F$h'p4$T  
f: I(signal_fw, -1,x * um, 0) * cm^2,  ！信号光沿光纤径向的强度分布 088C|  
  color = blue, )pS8{c)E  
  maxconnect = 1, ")ED)&e  
  width = 3, E67XPvo1+@  
  "signal" Z(HZB
  
q!&B6]  
V9T 4+  
; ------------- Fj^AWv^/  
diagram 5:                                  ！输出图表5 >$uUuiyL4  
^X&)'H  
"TransitionCross-sections" f>niFPW"  
hO6RQ0Iv@  
I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) MtUY?O.P2  
6.'$EtH  
x: 1450, 2050 y"-{$N  
"wavelength(nm)", @x IBET'!j4"  
y: 0, 0.6 vTd-x>n  
"cross-sections(1e-24 m&sup2;)", @y dF e4K"  
frame 2:'lZQ  
hx C2G
|?=  
hy JpK[&/Ct  
nLJBq)i  
f: s12_Tm(x * nm) /1e-24,      ！Tm³⁺吸收截面与波长的关系 K2HvI7$-  
  color = red, :tLbFW[  
  width = 3, X`1p'JD  
  "absorption" Cw#V`70a  
f: s21_Tm(x * nm) /1e-24,  ！Tm³⁺发射截面与波长的关系 2r;GcjezH  
  color = blue, <tuS,.  
  width = 3, u!Bk,}CE`  
  "emission" lR@& Z6lw  
~^7r?<aKc  

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