三阶非线性光束传输特性参数的差分计算.pdf

44 激《光杂志》l999年第 20卷第 5期 LASERJOURNAL(Vo1．20，No5．1999) 三阶非线性光束传输特性参数的差分计算 巍z 保 二{) 。 (华南师范大学抽理系 ，传输光学实验室 ，广州 51O631) 矩鬻 。 0fpa瑚m缫螂m嚷 薷竺耋 莹 三 苎翟墨慧墓譬最高 要方星薯 嚣 苎in： 三阶非线性光束的维传输为模型发展 ：h ： ；三ooller~- 解 一 种快速算 。 Njd。[]利‘用 墓 ： … 三阶非线性哲，射率…． 一，、、 ，一耋姜…‘篙…五‘怎’警…装…鲁一I嚣0黼0……I…～0 宝 +v +II=o (2) 确定微分矩阵。不妨假设u为区问[0，L]上二阶可导 是归 苎 。 一 J _化 引 F 篁兰 F‘、 ， ’ ! ‘’是． 的复莹 )和动量 一 由此到微分矩阵具方法菇o。‘… …’ V上，表征光束传输特性 的参数可以表示为 考虑 p≥3，：(p—1) 值多 式如下… 发散角 。： 日2～2f 一 ) (4) ．．()= ’ 0 (．t--~TI) (86) ． ～ 2．(— ) (4) ～ 。。… 三一、 一 曲率半径 H： 由r在起点{ 0)和终点(i=N一1)处插值需要额外 · ( 1·+ ·(一 ‘”， (51 的节 ·固此需要对函数 u进行延拓。根据实际情况 R 质量因子 M M=1(w 一 ) (6】 盘9 鼍套毒耋鑫差喜 8；No． 柬708譬； 辇釜 厦家高技惯性 势函数 Ⅱ 1999哥 5月 j7日收稿 《激光杂志)1999年第 z0卷第 5期 LASERJOURNAL(VM．20，No5．1999) 45 的需要 ，本文采用 对称的延拓方式，即假定 的函数 u 需要注意的是，以上关于线性算子 L的矩阵表示在形 在 x=0附近关于轴 x=0对称，函数在 x=L附近关于 式上是正确的，但在实际计算中，由于 r0=0，在计算 中 轴 x=L对称。由此得出在 (8)式中+ 出现了{暮l，：。：罟的情况。根据{箐比塔法则有 0时 ， 4 =一 —，u(x1)= “( ) ÷=辜，由此可以将计算矩阵L的过程作相应的修 iN一1时， =2L一 2N 卜2，u(z】)= “(2Ⅳ一卜2) 改，具体来说就是将 (1