首页 -> 登录 -> 注册 -> 回复主题 -> 发表主题
光行天下 -> OptiSystem,PhotonDesign,Rsoft -> RP Fiber Power仿真设计掺铥光纤激光器代码详解 [点此返回论坛查看本帖完整版本] [打印本页]

小火龙果 2020-05-28 16:28

RP Fiber Power仿真设计掺铥光纤激光器代码详解

(* )%Ru#}1X6  
Demo for program"RP Fiber Power": thulium-doped fiber laser, 1;&;5  
pumped at 790 nm. Across-relaxation process allows for efficient |}[nH>  
population of theupper laser level. 4v dNMV~  
*)            !(*  *)注释语句 otr>3a*'  
&mE?y%  
diagram shown: 1,2,3,4,5  !指定输出图表 .Q>!B?)  
; 1: "Powersvs. Position"     !分号是注释;光纤长度对功率的影响 :0 ^s0l  
; 2:"Variation of the Pump Power"  !泵浦光功率变化对信号输出功率的影响 8o-bd_  
; 3:"Variation of the Fiber Length"!信号输出功率vs 光纤长度的变化,仿真最佳光纤长度 E \{<;S  
; 4:"Transverse Profiles"             !横向分布,横坐标为半径位置 q\$6F)ha3  
; 5:"Transition Cross-sections"    !不同波长的跃迁横截面,横坐标波长,纵坐标为横截面 xh0xSqDM  
t FgX\4  
include"Units.inc"         !读取“Units.inc”文件中内容 TJy4<rb  
K:r\{#9  
include"Tm-silicate.inc"    !读取光谱数据 ~ 3T,&?r  
\cuS>G  
; Basic fiberparameters:    !定义基本光纤参数 cA;js;x@  
L_f := 4 { fiberlength }      !光纤长度 D>sYPrf  
No_z_steps := 50 {no steps along the fiber } !光纤步长,大括号{ }是注释,相当于备注 RuAlB*  
r_co := 6 um { coreradius }                !纤芯半径 .ve *Vp  
N_Tm := 100e24 { Tmdoping concentration }  !纤芯Tm离子掺杂浓度 RT*5d;l0  
!.Zt[g}  
; Parameters of thechannels:                !定义光信道 w'ybbv{c  
l_p := 790 nm {pump wavelength }                !泵浦光波长790nm Z=+03  
dir_p := forward {pump direction (forward or backward) }   !前向泵浦 {G*QY%j^  
P_pump_in := 5 {input pump power }                    !输入泵浦功率5W H:S,\D?%2x  
w_p := 50 um {radius of pump cladding }               !包层泵浦相应的半径 50um 6Nh0  
I_p(r) := (r <=w_p) { pump intensity profile }          !泵浦光强度分布 BC/_:n8O  
loss_p := 0 {parasitic losses of pump wave }           !泵浦光寄生损耗为0 y79qwM.  
.~klG&>aV  
l_s := 1940 nm {signal wavelength }                   !信号光波长1940nm sR_xe}-  
w_s := 7 um                          !信号光的半径 8_we: 9A  
I_s(r) := exp(-2 *(r / w_s)^2)            !信号光的高斯强度分布 R+=a`0_S  
loss_s := 0                            !信号光寄生损耗为0 BW Uq%o,@g  
'Ud| Ex@A9  
R_oc := 0.70 {output coupler reflectivity (right side) }      !输出耦合反射率 ..KwTf  
c^.l 2Q!  
; Function for defining themodel:   !定义模型函数,一定要有calc命令,否则函数只会被定义,但不会被执行 LSd*| 3E}n  
calc g+[kde;(^  
  begin fA^Em)cs2  
    global allow all;                   !声明全局变量 5lmO:G1  
    set_fiber(L_f, No_z_steps, '');        !光纤参数 raB+,Oi$G  
    add_ring(r_co, N_Tm); 3$p#;a:=n  
    def_ionsystem();              !光谱数据函数 Yx)o:#2  
    pump := addinputchannel(P_pump_in, l_p,'I_p', loss_p, dir_p);  !泵浦光信道 /c-nE3+rn  
    signal_fw := addinputchannel(0, l_s, 'I_s',loss_s, forward);      !前向信号光信道 KCR N}`^  
    signal_bw := addinputchannel(0, l_s, 'I_s',loss_s, backward);    !后向信号光信道 rVryt<2:@r  
    set_R(signal_fw, 1, R_oc);                                 !设置反射率函数 TKI$hc3|L  
    finish_fiber();                                   RR|\- 8;  
  end; V1bh|+o9  
.v`b[4M4  
; Display someoutputs in the Output window (on the right side): !在Output aera区域显示输出 xJ(:m<z  
show "Outputpowers:"                                   !输出字符串Output powers: yTAvF\s$(  
show"pump:     ", P_out(pump):d3:"W"  !输出字符串pump:和计算值(格式为3个有效数字,单位W) ffcLuXa  
show"signal:   ",P_out(signal_fw):d3:"W" !输出字符串signal:和计算值(格式为3个有效数字,单位W) 'K@0Wp  
[*Ju3  
sH!O0WL  
; ------------- +m$5a YX  
diagram 1:                   !输出图表1 --k:a$Nt  
WjMRH+  
"Powers vs.Position"          !图表名称 79 Bg]~}Z  
J !#Zi#8sF  
x: 0, L_f                      !命令x: 定义x坐标范围 >|udWd^$3  
"position infiber (m)", @x      !x轴标签;@x 指示这些字符串沿坐标轴放置 \cySWP[  
y: 0, 15                      !命令y: 定义y坐标范围 y{=NP  
y2: 0, 100                    !命令y2: 定义第二个y坐标范围 ;#anZC;  
frame          !frame改变坐标系的设置 Plo,XU  
legpos 600, 500  !图行在图表窗口中的位置(相对于左上角而言) zQ8!rCkg4  
hx             !平行于x方向网格 3fkk [U  
hy              !平行于y方向网格 PEXq:TA  
_a'A~JY  
f: P(pump, x),    !命令f: 定义函数图;P(pump, x)函数是计算x位置处的泵浦光功率 #}yTDBt  
  color = red,  !图形颜色 KS<Jv;  
  width = 3,   !width线条宽度 ^gR+S  
  "pump"       !相应的文本字符串标签 tJD] (F  
f: P(signal_fw, x),  !P(signal_fw ,x) 函数是计算x位置处的前向信号光功率 h'5Cp(G  
  color = blue,     : d' 5O8  
  width = 3, 5vOCCW  
  "fw signal" 3'eG ;<F  
f: P(signal_bw, x),   !P(signal_bw ,x) 函数是计算x位置处的后向信号光功率 k*$3i  
  color = blue, X[hM8G  
  style = fdashed, !~ rt:Z  
  width = 3, _" N\b%CkO  
  "bw signal" }DIF%}UK\  
$:vkX   
f: 100 * n(x, 2),    !n(x ,2) 函数是计算x位置处激活粒子数在能级2上的占比 9]yW_]P  
  yscale = 2,            !第二个y轴的缩放比例 Fr:5$,At7-  
  color = magenta, I$jvXl=$  
  width = 3, .ER98  
  style = fdashed, ygViPz<J  
  "n2 (%, right scale)" VXKT\9g3A  
8A2 z 5Aa  
f: 100 * n(x, 3),          !n(x ,3) 函数是计算x位置处激活粒子数在能级3上的占比 ]7/gJ>g,  
  yscale = 2, NGTe4Crx  
  color = red, L\  j:  
  width = 3, 3LK%1+)4  
  style = fdashed, Fb_S&!  
  "n3 (%, right scale)" e3(<8]`b[  
k W-81  
8l)  
; ------------- .iv3q?8.b  
diagram 2:                    !输出图表2 ?pwE0N^  
`E>vG-9  
"Variation ofthe Pump Power" At9X]t  
?[= U%sPu=  
x: 0, 10 kX;$}7n  
"pump inputpower (W)", @x w~@[ r4W  
y: 0, 10 HeIS;gfUY  
y2: 0, 100 8LrK94  
frame P_Uutn~  
hx  ]PX}b  
hy .W.;~`EW  
legpos 150, 150 JP=ZUu  
5h^[^*A?  
f: (set_P_in(pump, x);P_out(signal_fw)), !set_P_in(pump,x)改变泵浦信道功率;P_out(signal_fw)输出前向信号光 0HibY[_PbD  
  step = 5, ]+%=@mWYs  
  color = blue, k0ItG?Cv  
  width = 3, bZ5n,KQA5  
  "signal output power (W, leftscale)",     !相应的文本字符串标签 VifmZ;S@Y  
  finish set_P_in(pump, P_pump_in) w|Qd`  
U^$E'Q-VK  
f: (set_P_in(pump,x); 100 * n_av(2)),   !改变泵浦信号功率对能级2上激活粒子占比的影响 P. V\ov7m2  
  yscale = 2, %z)EO9vtr  
  step = 5, uMiyq<  
  color = magenta, a$}6:E  
  width = 3, f=_g8+}h  
  "population of level 2 (%, rightscale)", =vEkMJ Os  
  finish set_P_in(pump, P_pump_in) ukihx?5  
>t #\&|9I  
f: (set_P_in(pump,x); 100 * n_av(3)),   !改变泵浦信号功率对能级3上激活粒子占比的影响 "$)yB  
  yscale = 2, ?qT(3C9p  
  step = 5, -c={+z "  
  color = red, A*0*sZ0  
  width = 3, GX38~pq  
  "population of level 3 (%, rightscale)", A ,<@m2  
  finish set_P_in(pump, P_pump_in) HdCk!Fv  
&?T${*~  
UrK"u{G  
; ------------- JlIS0hnv  
diagram 3:                         !输出图表3 9d\N[[Vu]R  
+by|  
"Variation ofthe Fiber Length" !5.v'K'  
!-(J-45  
x: 0.1, 5 ^5x4q  
"fiber length(m)", @x JQT4N[rEE  
y: 0, 10 >hb- 5xC  
"opticalpowers (W)", @y @ ;J|xkJ  
frame fsI`DjKi)  
hx A-0m8<  
hy _85E=  
vKkf2 7  
f: (set_L(x);P_out(signal_fw)),     !改变光纤长度对信号光输出功率的影响 =t.F2'<[Z  
  step = 20,             uI9lK  
  color = blue, (`mOB6j  
  width = 3, Sf/W9Jw  
  "signal output" 5KaSWw/  
W-XN4:,qI  
;f: (set_L(x);P_out(pump)),                     !改变光纤长度对泵浦信号输出功率的影响 C->[$HcRa  
   step = 20, color = red, width = 3,"residual pump" 8Mb$+^zU  
Y[,U_GX/R  
! set_L(L_f) {restore the original fiber length } 8;.` {'r  
!F,s"  
P^;WB*V  
; ------------- k>-'AWH^v  
diagram 4:                                  !输出图表4 UP*yeT,P,  
B>aEH b  
"TransverseProfiles" i ~P91  
R\ e#$"a5  
I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) v1K4$&{F  
u2om5e:  
x: 0, 1.4 * r_co /um w6v1 q:20  
"radialposition (µm)", @x `#<eA*^g5  
y: 0, 1.2 * I_max *cm^2 /0!$p[cjm  
"intensity (W/ cm&sup2;)", @y |vT=Nnu  
y2: 0, 1.3 * N_Tm jL2MW(d^Q  
frame !$XO U'n  
hx bY&YSlO  
hy 6(sfpK'  
(ai72#nFtb  
f: N_dop(1, x * um,0),      !掺杂浓度的径向分布 ]Y$&78u8t  
  yscale = 2, K1 6s)S'  
  color = gray, jaNH](V  
  width = 3, 5?)}F/x  
  maxconnect = 1, qG*_w RF  
  "N_dop (right scale)" WT1q15U(=  
'9Qd.q7s|b  
f: I(pump, -1, x *um, 0) * cm^2,    !泵浦光沿光纤径向的强度分布 B~e7w 4  
  color = red, ur`}v|ZY  
  maxconnect = 1,           !限制图形区域高度,修正为100%的高度 %pM :{Z  
  width = 3, eKS:7:X  
  "pump" >sB=\  
&a~L_`\'  
f: I(signal_fw, -1,x * um, 0) * cm^2,  !信号光沿光纤径向的强度分布 n *Q4G}p  
  color = blue, xQZ MCd  
  maxconnect = 1, J$<:/^t  
  width = 3, ^M"HSewo  
  "signal" 8L@UB6b\  
}]Qmt5'NI  
WMRYT"J?N]  
; ------------- [v~Uy$d\  
diagram 5:                                  !输出图表5 JO]`LF]  
Sk!v,gx  
"TransitionCross-sections" JJ*0M(GG  
Zu)i+GeG  
I_max :=maxr(I(pump, -1, 0, 0), I(signal_fw, -1, 0, 0)) ,2E`:#$  
+q 4W0  
x: 1450, 2050 {lTR/  
"wavelength(nm)", @x d#RF0,Y9  
y: 0, 0.6 5I wX\  
"cross-sections(1e-24 m&sup2;)", @y F9ZOSL 8Q  
frame EgkZ$ah  
hx zS,%msT^A  
hy 'JJ1#kKa  
/QM0.{Ypl  
f: s12_Tm(x * nm) /1e-24,      !Tm3+吸收截面与波长的关系 vK[v eFH  
  color = red, ^p'D<!6sK  
  width = 3, K[`4vsE  
  "absorption" |F8;+nAVF#  
f: s21_Tm(x * nm) /1e-24,  !Tm3+发射截面与波长的关系 pqRO[XEp2  
  color = blue, uQXs>JuD  
  width = 3, c4FOfH|  
  "emission" >Lo6='G  
W ??;4  
lileisgsz 2021-09-28 09:47
感谢,视频上有点看不清楚
查看本帖完整版本: [-- RP Fiber Power仿真设计掺铥光纤激光器代码详解 --] [-- top --]

Copyright © 2005-2025 光行天下 蜀ICP备06003254号-1 网站统计