删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

啁啾脉冲放大激光系统中展宽器色散的解析算法

本站小编 Free考研考试/2021-12-29

摘要:时域色散精密控制是超短激光产生及其应用中的关键技术之一, 它通过控制各波长的光程产生相对延迟从而改变脉冲宽度. 展宽器是啁啾脉冲放大激光系统中对激光脉冲展宽的装置, 基于光线追迹法研究光线在展宽器中的传输路径, 可计算飞秒脉冲中各波长的光程, 进而计算脉冲展宽量并应用于系统设计. 由于展宽器的光程表达式复杂, 直接对其求导获得色散表达式较困难, 目前只能采用数值导数获得近似解, 这在计算过程中会引入误差, 不利于激光系统精确设计和优化. 本文介绍了一种易实现的求解展宽器色散的解析算法, 通过归纳展宽器光程表达式特点, 引入四个基元函数, 将光程表达式分解和反复代换, 可得到高阶色散的精确解析值. 本文首先对Martinez型展宽器重新光线追迹, 获得与Offner型展宽器一致的相位表达式, 其次通过解析算法获得了两种展宽器的精确高阶色散值, 最后将解析算法与数值算法的结果进行了比较. 该解析算法对于啁啾脉冲放大系统的参数设计具有实用价值.
关键词: 色散补偿/
啁啾脉冲放大/
脉冲展宽/
高阶色散

English Abstract


--> --> -->
飞秒激光啁啾脉冲放大(CPA)是结合激光脉冲宽度操控与激光能量放大以实现超快超强激光输出的技术. 基于这种技术, 目前已经可以获得峰值功率10 PW以上[1]的超强激光脉冲(1 PW = 1015 W), 为强激光与物质相互作用的研究提供了强大的技术工具. 峰值功率为太瓦量级的啁啾脉冲放大激光系统相比前者可以实现较高的脉冲重复频率, 在先进工业制造和生物医学领域具有优异的应用价值. 如在激光加工领域, 飞秒激光加工相比于传统的机械加工具有非接触的优点, 相比于长脉冲激光加工可以直接将介质等离子体化实现高精度的冷加工; 在生物医学领域, 飞秒激光可进行多光子成像获得比传统光学成像更深的穿透深度、超衍射极限的空间分辨率及更低的背景荧光噪声[2], 也可以用于外科手术和靶向治疗. 当前的超快激光啁啾脉冲放大技术正向着产生更高的功率密度、更洁净的激光脉冲和以新激光材料为基础打造更低造价的新型激光系统等方向发展, 这些都需要可靠的光学理论模型作为支撑. 本文在理论上提出了一种超快激光啁啾脉冲放大系统中展宽器色散的解析计算方法, 对于精确计算展宽器色散、优化系统结构和设计新型激光系统具有实用价值.
超快激光啁啾脉冲放大的原理是: 在飞秒激光的种子光进入光放大器前, 先对激光脉冲进行时域展宽, 再将展宽后的啁啾脉冲导入光放大器进行能量放大, 最终对脉冲宽度进行压缩. 这种方法可以避免直接将飞秒脉冲进行能量放大造成的光学器件损伤, 因而能获得更高的脉冲能量和相对优良的空间和时域光束质量. 啁啾脉冲放大技术的关键步骤之一就是激光脉冲的展宽和压缩, 是通过控制飞秒脉冲的时域色散来实现的. 利用色散装置使脉冲中不同波长的光通过的光程不同, 可以实现光脉冲宽度的变化. 不过这一过程引入的高阶色散与激光放大器里的高阶材料色散不能完全抵消会导致脉冲出现边峰, 使得脉冲畸变或不够洁净, 影响其后续应用, 因而精确分析啁啾脉冲放大系统展宽和压缩器提供的高阶色散值对于优化系统就非常重要, 这时需要对脉冲色散装置进行光线追迹, 分析光脉冲的相位及其各阶色散.
常用的脉冲展宽压缩装置有棱镜对、光栅对、棱栅对、块材料、光纤、Offner型、Martinez型或Barty型展宽器等介质或光学系统结构. 对于块材料和光纤等固体介质, 可以直接使用相位表达式$\varPhi = KL$代入材料折射率方程并对光频率ω求导进行色散分析[3], 由于可以较方便地获得色散的解析表达式, 本文不再叙述过程. 对于复杂的展宽压缩光学系统, 则需要进行光线追迹进行光程计算. 早在1969年, Treacy[4]就使用光线追迹对平行光栅对的相位和各阶色散进行了详细分析; 之后Fork等[5]、Martinez等[6]和张留洋等[7]对棱镜对的色散进行了计算; Kane和Squier[8]、Durfee等[9]、郭淑艳等[10]、徐婷婷等[11]、苏娟等[12]和Bucht等[13]对棱栅对的色散进行了研究; Cheriaux等[14]、Zhang等[15]、蒋捷等[16]、Wang等[17]、田金荣等[18]对Offner型展宽器进行了光线追迹与计算; Martinez[19]、Zhang等[20]、王勇等[21]、宋晏蓉等[22]和Su等[23]对Martinez型展宽器进行了计算分析; 孙大睿等[24]、Backus等[25]、Zhang等[26]对啁啾脉冲放大系统的色散优化给出了建议; 除此之外, 研究人员也分析了有限的光束宽度、光栅失谐、光栅拼接、非匹配光栅、曲面光栅、光谱宽度、面型平整度和球差等因素对展宽器或压缩器中脉冲的影响, 在这些方面有较多的文献报道, 本文不再展开介绍. 这些针对色散装置的研究为啁啾脉冲放大激光系统设计提供了很好的理论支持. 然而, 对于以上的展宽压缩光学系统装置, 由于相位表达式复杂, 难以直接推导出高阶解析导数表达式, 目前主要是通过设定自变量步长计算函数值再进行差商来进行数值求导的方法, 而这样会给系统引入误差, 不利于系统参数的设计和优化, 尤其是高阶数值导数值的误差会逐级累积增大; 另外步长越大, 使得偏离误差越大, 而步长越小, 计算量增大, 还会引入振荡误差, 也难以快速、精确地获得色散值, 不利于展宽器参数设计和色散的精密控制. 本文提出了一种解析算法, 并对Offner和Martinez型展宽器的各阶色散进行了解析计算, 也对数值算法和解析算法的结果进行了比较. 这种算法对于棱镜对、棱栅对等具有复杂相位表达式的解析导数求解也具有普遍适用性, 对于优化设计展宽器和精确补偿色散具有实用意义.
Offner展宽器和Martinez展宽器是啁啾脉冲放大激光系统中最常用的两种展宽器. 由于以前的研究人员针对展宽器的光线追迹过程各有描述, 有的为了简化计算也会使用近似, 从而使得色散计算过程和获得的展宽器相位表达式在表述和精度上都存在差异. 为了便于计算和比较, 本文首先对两种展宽器的相位表达式进行统一.
Offner型展宽器的结构和光路如图1所示, 本文直接使用文献中的光线追迹结果作为进行计算的相位Φ表达式[18]:
图 1 Offner型展宽器的结构与光路图
Figure1. Scheme of Offner stretcher.

$\begin{split} & {\varPhi _{{\rm{offner}}}} = {\varPhi _{{\rm{Lo}}}} - {\varPhi _{{\rm{Go}}}} = k{L_{\rm{o}}}= \\ \, & \frac{{2{\text{π}}}}{\lambda }(A{B_{\rm{o}}} + B{C_{\rm{o}}} + C{D_{\rm{o}}} + D{E_{\rm{o}}} + E{F_{\rm{o}}}) - \frac{{2{\text{π}}}}{d}E{A_{\rm{o}}}, \end{split}$
其中${\varPhi _{{\rm{Lo}}}}$是总光程引入的相位, ${\varPhi _{{\rm{Go}}}}$是光栅衍射相位修正因子, ABo, BCo, CDo, DEo, EFo, EAo是展宽器中各段光路的表达式, 与文献[18]表述相同, 引用如下:
$\begin{split}A{B_{\rm{o}}} =\,& \frac{{R\sin ({\theta _1} - \theta _2)}}{{\sin {\theta _1}}}, \\B{C_{\rm{o}}} =\,& C{D_{\rm{o}}} = \frac{{R\sin ({\theta _3} - {\theta _2})}}{{\sin {\theta _3}}}, \\D{E_{\rm{o}}} =\, & \frac{{R\sin \left( {{\theta _1} + 2{\theta _3} - 3{\theta _2}} \right)}}{{\sin \left( {{\theta _1} + 2{\theta _3} - 4{\theta _2}} \right)}} \\ \,&- \frac{{\left[ {\left( {R \!- \!L} \right)\sin \left( {{\theta _1} \!+\! 2{\theta _3}\! -\! 4{\theta _2}} \right) \!+\! R\sin {\theta _2}} \right]\sin {\theta _0}}}{{\sin \left( {{\theta _0} \!+\! {\theta _1} \!+\! 2{\theta _3} \!- \!4{\theta _2}} \right)\sin \left( {{\theta _1}\! + \!2{\theta _3}\! -\!4{\theta _2}} \right)}}, \end{split}$
$\begin{split}E{A_{\rm{o}}} =\,& \frac{{\left( {R - L} \right)\sin \left( {{\theta _1} + 2{\theta _3} - 4{\theta _2}} \right) + R\sin {\theta _2}}}{{\sin \left( {{\theta _0} + {\theta _1} + 2{\theta _3} - 4{\theta _2}} \right)}},\\E{F_{\rm{o}}} =\,& E{A_{\rm{o}}}\sin {\theta _5},\end{split}$
式中的各角度如下所示, 其中θ0为光栅与OO'的夹角, 是个初始设定值, 可以决定衍射光束在凹面镜上的位置, 此处设定为激光中心波长λ0在衍射后位于凹面镜的中心, 以最大程度地利用凹面镜实际使用面积. θ1θ5分别为光线与水平线、凹面镜法线、凸面镜法线、凹面镜法线、光栅法线和光栅的夹角.
$\begin{split}{\theta _0} =& \frac{{\text{π}}}{2} - \arcsin \left( {\frac{{{\lambda _0}}}{d} - \sin {\gamma _0}} \right),\\{\theta _1} = &\frac{{\text{π}}}{{\rm{2}}} - {\theta _0} - \arcsin \left( {\frac{\lambda }{d} - \sin {\gamma _0}} \right), \\{\theta _2} =& {\theta _4} = \arcsin \left( {\frac{{R - L}}{R}\sin {\theta _1}} \right), \\{\theta _3} =& \arcsin \left( {2\sin {\theta _2}} \right), \\{\theta _5} =& \arcsin \left[ {\frac{\lambda }{d} - \cos \left( {{\theta _0} + {\theta _1} + 2{\theta _3} - 4{\theta _2}} \right)} \right].\end{split}$
对Martinez型展宽器进行光线追迹可以得到相位Φ的表达式[19-22], 为了统一描述和追迹方法, 本文对Martinez型展宽器重新进行光线追迹并获得与Offner型展宽器一致的相位表达式为:
$\begin{split} & {\varPhi _{{\rm{Martinez}}}} = {\varPhi _{{\rm{LM}}}} - {\varPhi _{{\rm{GM}}}} = k{L_{\rm{M}}} \\ = \, & \frac{{2{\text{π}}}}{\lambda }(A{B_{\rm{M}}} + B{C_{\rm{M}}} + C{D_{\rm{M}}} + D{E_{\rm{M}}} + E{F_{\rm{M}}}) \\ & - \frac{{2{\text{π}}}}{d}E{A_{\rm{M}}}, \end{split}$
式中的各参数代表的意义与上述Offner型展宽器中的参数意义相同, 具体的值由下面过程给出. 用于追迹的Martinez型展宽器结构及光路如图2所示, 其中的红色实线为传输光, 深黄色的虚线为光栅, 黑色的实线曲线和直线分别为凹面反射镜和平面反射镜, 黑色的虚线曲线为凹面镜相对于平面镜的镜像, 其余黑色虚线直线为辅助线. 光线从从A点入射到光栅, 经凹面镜B点、平面镜C、凹面镜D和光栅E点后出射; θ0为光栅与OO'的夹角, 选取为此时激光中心波长λ0的衍射光线为AO'; θ1θ6分别为光线与水平线、凹面镜法线、平面镜法线、凹面镜法线、光栅法线和光栅的夹角; 图中凹面镜圆心位于O点, 圆半径为R; 平面镜位于R/2处; 光栅凹面镜间距O'AL. BD'延长线交OO'A'点; DE'延长线交OO'G点; AFEF. 由几何关系和正弦定理, 可得各段光路的光程表达式为:
图 2 Martinez型展宽器的结构与光路图
Figure2. Scheme of Martinez Stretcher.

$\begin{split}A{B_{\rm{M}}} =\,& \frac{{R\sin ({\theta _1} - \theta {}_2)}}{{\sin {\theta _1}}}, \\B{C_{\rm{M}}} +\,& C{D_{\rm{M}}} = A'D{'_{\rm{M}}} - A'{B_{\rm{M}}} \\ =\, & \frac{{R\sin ({\theta _3} + {\theta _4})}}{{\sin {\theta _3}}} - \frac{{A{B_{\rm{M}}}\sin {\theta _1}}}{{\sin {\theta _3}}}, \\D{E_{\rm{M}}} = \,& D{G_{\rm{M}}} - E{G_{\rm{M}}} \\ =\, & D{G_{\rm{M}}} - \frac{{A{G_{\rm{M}}}\sin {\theta _0}}}{{\sin \left( {\dfrac{{\text{π}}}{2} + {\theta _5}} \right)}} = \frac{{R\sin ({\theta _3} + {\theta _4})}}{{\sin ({\theta _3} + 2{\theta _4})}} \\ & - \left[ {R - L - \frac{{R\sin {\theta _4}}}{{\sin ({\theta _3} + 2{\theta _4})}}} \right]\frac{{\sin {\theta _0}}}{{\sin \left( {\dfrac{{\text{π}}}{2} + {\theta _5}} \right)}}, \\ E{A_{\rm{M}}} =\,& \frac{{A{G_{\rm{M}}}\sin ({\theta _3} + 2{\theta _4})}}{{\sin \left( {\dfrac{{\text{π}}}{2} + {\theta _5}} \right)}} \\=\, & \left[ {R - L - \dfrac{{R\sin {\theta _4}}}{{\sin ({\theta _3} + 2{\theta _4})}}} \right]\frac{{\sin ({\theta _3} + 2{\theta _4})}}{{\sin \left( {\dfrac{{\text{π}}}{2} + {\theta _5}} \right)}}, \\ E{F_{\rm{M}}} =\,& E{A_{\rm{M}}}\cos {\theta _6}, \end{split}$
式中的角度值为:
$\begin{split}{\theta _0} =\,& \frac{{\text{π}}}{2} - \arcsin \left( {\frac{{{\lambda _0}}}{d} - \sin {\gamma _0}} \right), \\{\theta _1} =\,& \frac{{\text{π}}}{2} - {\theta _0} - \arcsin \left( {\frac{\lambda }{d} - \sin {\gamma _0}} \right), \\{\theta _2} =\,& \arcsin \left( {\frac{{R - L}}{R}\sin {\theta _1}} \right), \\{\theta _3} =\,& 2{\theta _2} - {\theta _1}, \\{\theta _4} =\,& \arcsin \left[ {\left( {\frac{{\sin {\theta _2}}}{{\sin {\theta _3}}} - 1} \right)\sin {\theta _3}} \right], \end{split}$
$\begin{split} & {\theta _5} = {{\text{π}}}/{{\rm{2}}} - ({\theta _0} + {\theta _3} + 2{\theta _4}),\\ &{\theta _6} = {{\text{π}}}/{{\rm{2}}} - \arcsin \left( { {\lambda }/{d} - \sin {\theta _5}} \right). \end{split}$
至此, 完成了对两种展宽器的相位表达式的统一, 并给出了对应项的参数, 由于展宽器结构上的差异, 总的相位表达式虽然相同, 但具体的参数会有不同. 通过表达式也可以看出, 直接计算获得光通过两种展宽器后的相位Φ关于光频率ω的各阶解析导数是比较困难的.
通过分解各段光程表达式并循环代换的方式来获得解析计算结果. 观察相位式中的各段光路表达式可以发现, 所有式子都由四个初等函数构成, 即y = sin(x), y = arcsin(x), y = A/B, y = A + B. 为了求解解析解, 只需求解这四个初等函数的导数, 然后把所有表达式依次代入这4个函数即可. 以y = sin(x)为例, 参数x, y都是光波长λ或频率ω的导数, 该初等函数的一阶到四阶导数使用简单的求导规则就可以获得:
$\left\{ {\begin{aligned}& {y = \sin x}, \\& {\dot y = \cos x \cdot \dot x}, \\& {\ddot y = - \sin x \cdot {{\dot x}^2} + \cos x \cdot \ddot x}, \\& {\dddot y = - \cos x \cdot {{\dot x}^3} - 3\sin x \cdot \dot x \cdot \ddot x + \cos x \cdot \dddot x}, \\& {y^{(4)}} = \sin x \cdot {{\dot x}^4} - 6\cos x \cdot {{\dot x}^2} \cdot \ddot x - 3\sin x \cdot {{\ddot x}^2} \\ & \quad\quad\;\; - 4\sin x \cdot \dot x \cdot \dddot x + \cos x \cdot {x^{(4)}} .\end{aligned}} \right.$
那么只要输入参数x的各阶导数值, 代入(3)式就可以获得参数y = sin(x)的各阶解析导数值. 同样的道理, 可得y = arcsin(x), y = A/B, y = A + B这三个初等函数的解析导数表达式(见附录). 只要已知输入参数xAB的导数值, 就可以获得这三个初等函数的对应的解析值. 使用软件编程, 将以上表达式创建4个自定义函数, 保存为函数文件后, 对于各段光路反复调用这4个函数就能完成相位的求解.
以求解Martinez型展宽器的${\theta _1}$为例, 为了最终求解相位关于角频率的导数, 将追迹波长λ变换为角频率ω, 得${\theta _1} = \dfrac{{\text{π}}}{2} - {\theta _0} - \arcsin \left( {\dfrac{{2{\text{π}}{\rm{c}}}}{{\omega d}} - \sin {\gamma _0}} \right)$, 式中的${\theta _0}$是光栅与对称轴的夹角, ${\gamma _0}$是光线在光栅上的入射角, d是光栅常数, c是光速, 2和π是常数, 各阶导数为0; 角频率ω是自变量, 其一阶导数为1, 高阶导数为0. 所以参照上述初等函数发现, $\dfrac{{2{\text{π}}{\rm{c}}}}{{\omega d}} = \dfrac{{2{\text{π}}{\rm{c}}}}{d\omega} = A/B$, 其中此处$A = \dfrac{{2{\text{π}}c}}{d}$, 各阶导数为0, $B = \omega $, 一阶导数为1, 各阶导数为0, 代入到A/B的解析导数式中, 可得$\dfrac{{2{\text{π}}c}}{{\omega d}}$的各阶解析导数值, 再代入A + B的解析导数式中求得$\dfrac{{2{\text{π}}c}}{{\omega d}} - \sin {\gamma _0}$的各阶导数值, 继续代入arcsin(x)的解析导数表达式中求得$\arcsin \left( {\dfrac{{2{\text{π}}c}}{{\omega d}} - \sin {\gamma _0}} \right)$的各阶导数值, 再代入A + B的解析表达式中求解得到$\dfrac{{\text{π}} }{2} - {\theta _0} - \arcsin \left( {\dfrac{{2{\text{π}}c}}{{\omega d}} - \sin {\gamma _0}} \right)$的各阶解析导数值, 这样就获得了${\theta _1}$的各阶解析导数值. 使用同样的方法, 依次得到${\theta _2}\text{---}{\theta _6}$、光程ABMEFM以及最终的相位Φ的各阶解析导数值. 虽然该过程需要把表达式分解成很多项依次代换, 但是使用程序设计很容易实现, 所以代换过程的程序的编写和运算都非常快, 并且计算结果是解析值, 比通常的数值运算更为精确.
基于上述算法, 对Offner型展宽器和Martinez型展宽器进行了编程计算. 使用的两种展宽器的计算参数如表1所列.
计算参数Offner型
展宽器
Martinez型
展宽器
凹面镜曲率半径R/mm12101210
光栅凹面镜间距L/mm1005400
激光光栅入射角γ0/(°)2222
光栅刻线密度/mm–112001200
展宽程数22


表1用于计算的Offner型与Martinez型展宽器参数
Table1.Structural parameters of Offner and Martinez stretchers.

采用这个参数时, 计算可得对于800 nm激光波长, 两种展宽器产生的二阶色散都为1.6378 × 106 fs2. 使用高斯成像公式可得, Offner型展宽器的等效光栅对间距为$ 2(R-L)\cos\theta_0 $, Martinez型展宽器的为$ 2(R/2-L)\cos\theta_0 $, 此处两种展宽器参数对应的等效光栅对间距值都为332.407 mm.
使用该解析方法对两种展宽器进行色散计算, 结果如图3所示. 红线给出的是Offner型展宽器的一阶到四阶色散, 其中图3(a)为群延迟, 图3(b)为群延迟色散或二阶色散, 图3(c)图3(d)分别为三阶色散和四阶色散. 同理, 图3(a)图3(d)的蓝线分别给出的是Martinez型展宽器的一阶到四阶色散图. 可见采用该参数时, 两种展宽器提供的色散大概是相同的, 而Offner型展宽器在光谱边带处提供更大的色散量.
图 3 Offner型与Martinez型展宽器700?900 nm的解析计算色散曲线 (a) GD; (b) GDD; (c) TOD; (d) FOD
Figure3. Analytical dispersion results of Offner and Martinez stretcher: (a) GD; (b) GDD; (c) TOD; (d) FOD.

接下来对两种展宽器的解析算法与数值算法的计算结果进行比较. 使用差商数值算法并选择数值计算步长分别为0.1和5.0 nm, 此时获得的两种展宽器的二阶色散值与解析色散值由图4给出, 其中图4(a), (c)为Offner型展宽器的结果, 图4(b), (d)为Martinez型展宽器的结果. 从图4(a), (b)中可以看出, 随着步长数值的增大, 二阶色散值的数值计算误差也相对于解析计算值增大, 所以数值计算的步长不宜过大. 图4(c), (d)为两种展宽器的四阶色散值, 其中红色点线的0.1 nm步长的值表现出显著振荡, 这是由于步长过小使得数值导数的分母过小, 引入了数值差分误差, 因而数值计算的步长也不宜过小. 本文通过多次模拟发现, 采用数值差商计算展宽器色散过程中, 为避免数值导数差分误差和累积误差过大, 仅进行二阶色散的计算时, 最佳步长可以选为0.0006 nm左右, 若进行四阶色散的计算, 最佳步长在0.3 nm左右. 通过对比数值求解与解析计算结果, 可以发现本文的解析算法计算可以避免数值计算的两种误差问题, 尤其在进行高阶色散的计算时表现出更好的实用性.
图 4 (a) Offner型展宽器的二阶色散值; (b) Martinez型展宽器的二阶色散值; (c) Offner型展宽器的四阶色散值; (d) Martinez型展宽器的四阶色散值. 图中蓝色实线为解析计算值, 红色点线是步长为0.1 nm的数值计算值, 黑色虚线是步长为5 nm的数值计算值
Figure4. (a) GDD of Offner stretcher; (b) GDD of Martinez stretcher; (c) FOD of Offner stretcher; (d) FOD of Martinez stretcher. Blue solid lines in the figure are analytical results; red dotted lines are numerical results with step of 0.1 nm; black dash lines are numerical results with step of 5 nm.

介绍了一种啁啾脉冲放大激光中展宽器色散的解析计算方法. 基于光线追迹, 本文首先给出了与Offner型展宽器光线追迹结果相统一的Martinez型展宽器追迹公式, 进而使用初等函数循环代换的方式获得了两种展宽器的解析计算结果, 接下来本文对两种展宽器的数值计算结果与解析计算结果进行了比较, 最终获得了两种展宽器一阶到四阶色散的解析计算结果. 可以看出, 该方法得到的解析计算结果更为精确, 对于啁啾脉冲放大激光系统中的展宽器设计及色散优化提供了实用的工具; 另外本算法可以不依赖于特定程序设计语言的内置函数, 可以较好地进行推广, 能做成手机应用程序或可视化应用程序进行使用. 当然, 接下来研究的是综合整个啁啾脉冲放大系统的色散进行优化, 并结合色散测量装置FROG或SPIDER进行更为精确的色散控制, 以期待获得更好的色散补偿结果.
感谢中国科学院物理研究所张伟博士在啁啾脉冲放大激光系统实验参数方面的有益讨论和技术支持.
附录 用于展宽器光线追迹解析算法的四个初等函数的四阶导数式
y = sin(x)函数的一至四阶全导数:
$\small \left\{ {\begin{aligned} & {y = \sin x,} \\& {\dot y = \cos x \cdot \dot x,} \\& {\ddot y = - \sin x \cdot {{\dot x}^2} + \cos x \cdot \ddot x,} \\& {\dddot y = - \cos x \cdot {{\dot x}^3} - 3\sin x \cdot \dot x \cdot \ddot x + \cos x \cdot \dddot x,} \\& {y^{(4)}} = \sin x \cdot {{\dot x}^4} - 6\cos x \cdot {{\dot x}^2} \cdot \ddot x - 3\sin x \cdot {{\ddot x}^2} \\& \quad\quad\;\; - 4\sin x \cdot \dot x \cdot \dddot x + \cos x \cdot {x^{(4)}}.\end{aligned}} \right.\tag{A1}$
y = arcsin(x)函数的一至四阶全导数:
$\small \left\{ {\begin{aligned}& {y = \arcsin x,} \\ & {\dot y = {{(1 - {x^2})}^{ - \frac{1}{2}}} \cdot \dot x,} \\& {\ddot y = x{{(1 - {x^2})}^{ - \frac{3}{2}}} \cdot {{\dot x}^2} + {{(1 - {x^2})}^{ - \frac{1}{2}}} \cdot \ddot x,} \\& \dddot y = {{(1 - {x^2})}^{ - \frac{3}{2}}} \cdot {{\dot x}^3} + 3{x^2} \cdot {{(1 - {x^2})}^{ - \frac{5}{2}}} \cdot {{\dot x}^3} \\& \quad\;\; + 3x{{(1 - {x^2})}^{ - \frac{3}{2}}} \cdot \dot x \cdot \ddot x + {{(1 - {x^2})}^{ - \frac{1}{2}}} \cdot \dddot x, \\& {y^{(4)}} = 9x{(1 - {x^2})^{ - \frac{5}{2}}} \cdot {{\dot x}^4} + 15{x^3}{(1 - {x^2})^{ - \frac{7}{2}}} \cdot {{\dot x}^4} \\& \quad\quad\;\; + 18{x^2}{(1 - {x^2})^{ - \frac{5}{2}}} \cdot {{\dot x}^2} \cdot \ddot x \\ & \quad\quad\;\; + 4x{(1 - {x^2})^{ - \frac{3}{2}}} \cdot \dddot x \cdot \dot x + 6{(1 - {x^2})^{ - \frac{3}{2}}} \cdot {{\dot x}^2} \cdot \ddot x \\& \quad\quad\;\; + {(1 - {x^2})^{ - \frac{1}{2}}} \cdot {x^{(4)}} + 3x{(1 - {x^2})^{ - \frac{3}{2}}} \cdot {{\dot x}^2}.\end{aligned}} \right.\tag{A2}$
y = A/B函数的一至四阶全导数:
$\small \left\{ {\begin{aligned}& y = A{B^{ - 1}}, \\& {\dot y = \dot A{B^{ - 1}} - {B^{ - 2}}A\dot B,} \\& {\ddot y = \ddot A{B^{ - 1}} - 2{B^{ - 2}}\dot A\dot B + 2{B^{ - 3}}A{{\dot B}^2} - {B^{ - 2}}A\ddot B,} \\& \dddot y = \dddot AB - 3{B^{ - 2}}\ddot A\dot B - 3{B^{ - 2}}\dot A\ddot B \\ & \quad\;\; + 6{B^{ - 3}}\dot A{{\dot B}^2} + 6{B^{ - 3}}A\ddot B\dot B - 6{B^{ - 4}}A{{\dot B}^3} - {B^{ - 2}}A\dddot B, \\& {y^{(4)}} = {A^{(4)}}{B^{ - 1}} - 4{B^{ - 2}}\dddot A\dot B - 4{B^{ - 2}}\dot A\dddot B \\ & \quad\quad\;\; - 6{B^{ - 2}}\ddot A\ddot B + 24{B^{ - 3}}\dot A\dot B\ddot B - 36{B^{ - 4}}A{{\dot B}^2}\ddot B \\ & \quad\quad\;\;+ 12{B^{ - 3}}\ddot A{{\dot B}^2} - 24{B^{ - 4}}\dot A{{\dot B}^3} + 8{B^{ - 3}}A\dot B\dddot B \\ & \quad\quad\;\;+ 6{B^{ - 3}}A{{\ddot B}^2} + 24{B^{ - 5}}A{{\dot B}^4} - {B^{ - 2}}A{B^{(4)}}. \end{aligned}} \right. \tag{A3}$
y = A + B函数的一至四阶全导数:
$\small\left\{ {\begin{aligned}& {y = A + B,} \\& {\dot y = \dot A + \dot B,} \\& {\ddot y = \ddot A + \ddot B,} \\& {\dddot y = \dddot A + \dddot B,} \\& {{y^{(4)}} = {A^{(4)}} + {B^{(4)}}.} \end{aligned}} \right. \tag{A4}$

相关话题/计算 激光 系统 过程 设计

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • Parity-time对称性对电注入半导体激光器的模式控制
    摘要:非厄米的描述对于开放系统有重要意义,满足parity-time对称性的哈密顿量,其参数在一定范围内可以使能量具有实的本征值.本文通过模拟,研究了损耗大小以及结构对称性对条形波导中的parity-time对称性的影响,并通过实验发现了电注入条件下由parity-time对称破缺导致的脊条波导模式 ...
    本站小编 Free考研考试 2021-12-29
  • 激光等离子体光丝中太赫兹频谱的调控
    摘要:研究了双色激光场激发空气成丝产生太赫兹辐射频谱的变化规律.实验观察到随驱动光功率和光丝长度增加,太赫兹光谱主要发生红移的现象.分析表明,由于等离子体密度的增加,太赫兹辐射的趋肤深度减小,等离子体吸收主导了红移的发生.在光丝足够短的条件下,趋肤深度远大于光丝长度,从而产生等离子体振荡主导的太赫兹 ...
    本站小编 Free考研考试 2021-12-29
  • 强激光间接驱动材料动态破碎过程的实验技术研究
    摘要:强激光驱动加载已成为冲击波作用下材料动态破碎过程研究的一种有效手段.采用间接驱动方式,设计合适的腔型进行物理实验研究,可实现更大且更均匀的冲击加载一维区.采用数值模拟和物理实验方法,研究强激光间接驱动材料动态破碎过程的实验技术.首先,利用IRAD程序设计适用于开展动态破碎过程研究的半柱腔,其直 ...
    本站小编 Free考研考试 2021-12-29
  • 天线方向系数的一类计算逼近方法
    摘要:天线的方向系数是天线的核心性能指标之一,准确计算方向系数是高性能天线应用的核心要求.本文基于平面近场测试理论、实测数据和快速傅里叶变换算法,系统阐述基于近场测试来数值计算天线方向系数的原理,并进行深入的误差分析.本文选择一种方向图函数和方向系数已知的被测天线,来检验所讨论的误差评估方案.评估分 ...
    本站小编 Free考研考试 2021-12-29
  • 双层耦合非对称反应扩散系统中的超点阵斑图
    摘要:通过线性耦合Brusselator模型和Lengyel-Epstein模型,数值研究了双层耦合非对称反应扩散系统中图灵模之间的相互作用以及斑图的形成机理.模拟结果表明,合适的波数比以及相同的对称性是两个图灵模之间达到空间共振的必要条件,而耦合强度则直接影响了图灵斑图的振幅大小.为了保证对称性相 ...
    本站小编 Free考研考试 2021-12-29
  • 纳秒脉冲激光诱导空气等离子体的近红外辐射特性
    摘要:开展纳秒激光诱导空气等离子体近红外辐射特性的实验研究,对波长为532nm的脉冲ns激光诱导产生的空气等离子体的近红外光谱进行测量.结果表明:空气等离子体的近红外辐射在光谱范围为1100—2400nm内由连续谱和线状谱组成,光谱指认表明线谱主要来源于N,O原子的中性原子谱和氮分子的振动光谱.通过 ...
    本站小编 Free考研考试 2021-12-29
  • 离散可积系统: 多维相容性
    摘要:对比已有完善而系统理论的微分方程领域,差分方程理论尚处于发展之中.近年来离散可积理论的进展,带来了差分方程理论的革命.多维相容性是伴随离散可积系统研究出现的新的概念,作为对离散可积性的一种理解,提供了构造离散可积系统的B?cklund变换、Lax对和精确解的工具.本文旨在综述多维相容性的概念及 ...
    本站小编 Free考研考试 2021-12-29
  • 可积谐振系统中的极端波事件研究进展
    摘要:从微观角度上讲,单个极端异常波事件可视为可积模型方程的时空局域有理函数解.本文主要讨论了三类典型的可积谐振相互作用模型(即长波短波谐振方程,三波谐振相互作用方程,非线性薛定谔和麦克斯韦-布洛赫方程)的基阶Peregrine异常波解及其相关研究进展;明确指出了这些基阶异常波解形式具有普适性,可推 ...
    本站小编 Free考研考试 2021-12-29
  • 可积系统多孤子解的全反演对称表达式
    摘要:多孤子解是非线性数学物理系统的基本激发模式.文献中存在各种类型的表达式,如广田(Hirota)形式,朗斯基(Wronskian)或双朗斯基形式和法夫(Phaffian)形式.最近在多地系统的研究中,我们发现使用一种全新但等价的形式具有极为简洁和方便的优点.本文主要综述多种类型可积非线性系统的多 ...
    本站小编 Free考研考试 2021-12-29
  • 光纤激光器中包层功率剥离器散热性能的优化
    摘要:光纤激光器系统中的包层功率剥离器在去除残余抽运光和高阶激光时,由于光热转换会产生大量的热能,所以将热能高效的耗散成为当前的研究热点.本文对国内外现有的五种剥离器进行了仿真研究与对比,发现用高折胶法制作剥离器时,改变填胶孔的形状,可以有效地增大热源与传热介质间的表面积-体积比,从而降低剥离器工作 ...
    本站小编 Free考研考试 2021-12-29