第五章孤立波

第五章孤立波 一个轮廓清晰又光滑的水堆，犹如一个大鼓包，沿着运河一直向前推进。

第五章孤立波 第一节历史回顾第二节 KdV方程第三节正弦—高登方程第四节非线性薛定谔方程与光学孤立子

第一节历史回顾 1.一个奇特的水波 2.孤立波与孤立子

1.一个奇特的水波 罗素的发现一个奇特的水波约170年前，苏格兰海军工程师罗素 (J.Scott Russell)在一次偶然中观察到一种奇特的水波。 1844年，他的报告：“我看到两匹骏马拉着一条船沿运河迅速前进。当船突然停止时，随船一起运动的船头处的水堆并没有停止下来。它激烈地在船头翻动起来，随即突然离开船头，并以巨大的速度向前推进。一个轮廓清晰又光滑的水堆，犹如一个大鼓包，沿着运河一直向前推进在行进过程中其形状与速度没有明显变化。我骑马跟踪注视，发现它保持着起始时约 30 英尺长，1-1.5 英尺高的浪头，约以每小时8-9英里的速度前进后来，在运河的拐弯处消失了”。罗素称之为孤立波 - Solitary wave。

1.一个奇特的水波 罗素的发现水槽中的实验罗素在一长水槽的一端，用一重锤垂落入水中，反复的观察重锤激起的水浪的运动。实验结论水波移动速度 v、水的深度 d 及水波幅度 A的关系为： B为比例常数实验结果说明水波的运动速度与波幅高度有关，波幅高的速度较快，且波幅的宽度对高度之比也相对较窄。

1.一个奇特的水波 漫长的发展史 KdV方程半个世纪后，1895年，两位荷兰科学家科特维格(Kortweg)与德弗雷斯(de Vries) 认为：罗素观察到的孤立波是波动过程中非线性效应与色散现象互相平衡的结果。他们建立了KdV方程：解 sech(x)为双曲正割函数,具有钟形形状。 FPU问题又过半个多世纪，1955年，美国阿尔莫斯国家实验室，著名物理学家费米(E.Fermi)、帕斯塔(J.Pasta)和乌莱姆(Ulam) 设计了一个数值计算实验：“非线性弹簧联结的64个质点组成弦的振动”，发现初始对少数质点激发，长时间后能量几乎全部回到了初始集中在少数质点上的状态。

1.一个奇特的水波 漫长的发展史 1965年两位美国数学家，采布斯基(Zabusky)与克鲁思卡尔(Kruskal)，用计算机计算发现， FPU 问题与 KdV 方程的解直接有关。此后，人们发现，在许多物理体系中都存在KdV方程，说明孤立波是一种普遍存在的物理现象。KdV方程成为数学物理的一个基本方程孤立波方程在理论上和实验上已对孤立波巳作了大量的研究。此后发现，除KdV方程外还有其它微分方程具有孤立波解。在数学上通常认为下列非线性方程的解的性质具有孤立波特性。 (1)KdV方程 (2)正弦—高登(Sine-Gordon)方程 (3)户田(M.Toda)非线性晶格方程 (4)非线性薛定谔方程(NLSE)

2.孤立波与孤立子 孤立波在形态上孤立波是存在于自然界里的相干结构(coherent structure，或称拟序结构)。从美丽的木星上的巨型红斑到固体中的电荷密度波都属于这样的有序结构。从运动形式上相干结构与混沌运动既是相互对立的。混沌运动是非线性中奇妙的无序状态，相干结构反映了非线性系统中的惊人有序性。在尺度上：大到天文范围(木星上巨型红斑达4×108米，约地球与月亮之间的距离；泰国安达曼海面出现的孤立波约150公里宽；水面上孤立水波的尺寸在1米量级)，小到纳米(二硫化钽晶体中的电菏密度孤立波)。

2.孤立波与孤立子 孤立子计算发现，两个在空间传播的孤立波具有碰撞特性，说明： (1)孤立波非常的稳定；(2)象一个物质粒子。人们将具有碰撞特性孤立波称为“孤立子- soliton”，简称“孤子”。孤立子是由非线性场所激发的、能量不弥散的、形态上稳定的准粒子。

2.孤立波与孤立子 孤立子定义 ①向单方向传播的行波； ②分布在空间的一个小区域中； ③波动形状不随时间演变而发生变化； ④孤立波之间的相互作用具有类似粒子一样的弹性碰撞。孤立子具有一切粒子所具有的特性，如能量、动量、质量、电荷、自旋等等，也遵循一般的自然规律，如能量、动量、质量守恒定律。它又有波动性，存在于一切可以出现波动的介质里。孤立波子哪里？孤立波除存在于浅水层外，还可在水层深处。固体介质、电磁场、等离子体、生物体、以及微观粒子的波动性中都可能有孤立波存在。它是一种行波，既可以速度 v 在空间传播，又可以处于静止状态。

2.孤立波与孤立子 孤立波类型 ( i )波包型 (ii)凹陷型 (iii)扭结型 (iv)反扭结型 ( i )、(ii) 两种是在时， (iii)、(iv) 两种是在时, u(x)趋近于不同的数值

第二节 KdV 方程 1. 波动中的非线性会聚效应 2. 波动中的色散 3. KdV方程 4. KdV方程的孤立波解

1. 波动的会聚效应 浪花的形成微风吹拂，水面只掀起层层碎浪；劲风吹来，浪尖则卷起浪花。同样的情况可以出现在海滩边。远处传来的海浪越近海岸，浪头越高，终于在离海岸不远处卷起了浪花。这是因为海滩对水浪运动产生某种阻滞力，浪的较低部分受到阻滞力大，较高部分阻滞力小。因此当水浪高处前进速度大，低处前进速度小，水浪会在前进中越来越前倾，在某一时刻波前出现坍塌，卷起了浪花。当水浪的不同部分有不同行进速度时，会出现会聚效应，形成浪花。

1.波动的会聚效应 浪花的形成数学表述设介质中x 处的粒子密度 n(x,t)，由粒子守恒 dx/dt = v 如果速度v = v0是常数，方程具有行波解：n = (x - v0 t) 介质的移动速度 v0 即波速。在初始时刻介质中出现的扰动n(x,0)= F(x)，波动将在传播中保持不变。波动将以速度v0无畸变地沿 x 方向前进。如果波动的速度 v 与介质的密度 n 有关，设： n(x,t)=F [x - v(n) t] 当出现波包前沿变陡，形成波包会聚。

2.波动中的色散 平面波的相速一个频率 w为沿 x 方向传播的平面波为：等相位面运动速度代表一列平面波的传播速度—相速等相位面 f =k x- w t= const

2.波动中的色散 色散波一个波动可以看成许多平面波（谐波）w1、w2、w3… 的合成：如果所有谐波都以同一的速度行进，w1/k1=w2/k2=...=常数，是非色散波；如果每个谐波都有不同的行进速度， w/k≠常数，是色散波。色散波将在传播中因弥散而消失。

2. 波动中的色散 色散关系设一波动方程：得关系将解代入：得色散关系由色散关系可求得波的群速线性波动方程与色散关系间存在着对应关系。于是便可以在波动方程与色散关系之间建立直接的对应。根据对应关系，可以由色散关系直接构造出波动方程来。

3. KdV 方程 KdV方程的推导对于不可压缩介质，粒子数密度 n 应用粒子速度 v 来替代，即有在重力作用下水波的色散关系：(g-重力加速度，h 水深) w (k) = g k tank (k h) 级数展开近似式利用考虑导致波形坍塌的非线性效应得KdV方程

3. KdV 方程 KdV方程的孤立波解设方程的解为方程变为整理后其解是罗素观察到的水面上奇特水波

3. KdV 方程 KdV方程与同宿线将KdV方程二次积分：雅可比矩阵本征方程两个一阶方程对于 V(u) 的极小点， V”(u) >0 ，l2的实部为零的共轭复根是中心点，围绕中心点的相轨线是椭圆；对于V(u)的极大点， V”(u) <0 ，l2为符号不同的实根，因而是鞍点。方程有一个鞍点，沿着一条流形出发，绕了一圈之后又回到了鞍点，流形为同宿线

KdV方程的孤立波解-1 设方程的解为方程变为整理后积分一次二次积分整理后有其中

KdV方程的孤立波解-2 讨论：函数是三次曲线，三个零点b1, b2, b3。根据三次代数方程的解与系数的关系，设改写为写成积分 V(u)在区间[b1, b2 ]内为负；[b2, b3] 为正；其余区间为无界对于[b2, b3]，引入代换上面积分变为椭圆函数模数

KdV方程的孤立波解-3 设积分这是KdV方程的椭圆余弦波解，u(x)是一列周期为 L的行波。cn( x ) 波周期为 u(x)的周期 L为它不具局域性质，不是孤立波。

KdV方程的孤立波解-4 孤立波解 1. 当模数，即，椭圆余弦函数这是振幅十分小的余弦波解 2. 当模数，即，椭圆余弦函数双曲正割函数当V(u)中，得：这是罗素观察到的水面上奇特水波

第三节正弦—高登方程 1. 一维原子链与正弦—高登 (sine-Gordon)方程 2. 正弦－高登方程孤立波解

1. 一维原子链与sine-Gordon方程 S -G 方程一维原子链模型：一串周期地束缚在非线性弹簧上的原子。原子链的哈密顿为：设势函数，第k个原子的运动方程为： (V(y)为外场势能)

1. 一维原子链与sine-Gordon方程 S -G 方程通过将分离变量运动方程过渡到连续变量一般情况 S -G 方程 V=cosy时

2. sine-Gordon方程孤立波解 S -G 方程解方程的解取形式则得常微分方程利用特殊函数积分公式两种特殊情况 (1) ，snx~sinx (2) ，snx~thx 为椭圆函数的模数

2. sine-Gordon方程孤立波解 S -G方程解由双曲函数公式经推导得最后得解称正、反扭曲解

2. sine-Gordon方程孤立波解 S -G方程解解的性质： (1) 扭折解当时， ; 当时， (2)反扭折解当时， ; 当时， Sine—Gordon方程解可以解释的物理现象：晶格位错的传播、磁体中畴壁的运动电荷密度波、基本粒于模型、超导约瑟夫逊结的列阵构成的传输线等。

S -G方程与异宿线 S-G方程可化成方程组：与单摆方程形式相同因此有与单摆同样的相图，坐标原点 (0,0) 为一中心点；坐标（-p, 0）与（p, 0）处为鞍点。当时相应椭圆轨道；当时相应异宿线与异宿线上支对应的是扭折解波与异宿线下支对应的是反扭折解

第四节非线性薛定谔方程与光学孤立子 1. 光纤中的光脉冲压缩效应 2. 非线性薛定谔方程（NLSE）及孤立波解 3. 光学孤立子的传播特性 4. 孤立子激光器 5. 全光型孤立子通信

1. 光纤中的光脉冲压缩效应 光纤导光原理与KdV方程孤立子类似，由非线性薛定谔方程(NLSE)描述的光学孤立子是光纤的色散效应与非线性压缩效应平衡的结果。光纤是二氧化硅为主要成分的细丝，直经在数微米至数十微米，外面为包层和涂敷层。设线芯、包层和涂敷层的折射率分别为 n1、n2、n0。n1＞ n2，所以包层界面有临界全反射角。如果满足全反射条件，光线将在光纤和包层的界面上，产生全反射而向前传播。

1. 光纤中的光脉冲压缩效应 光纤的非线性效应入射进光纤光束为强激光，常看成准单色光。在强光作用下，光纤出现非线性极化。极化矢量 P 与场强 E 的关系： c(1)，c(2)，c(3)，分别为线性的与二次、三次极化率，通常由电感应矢量 D 与极化矢量 P 。忽略高次非线性效应：，介质折射率为由于 n1 光在光纤中产生相移相移量与光强有关，因而导致不同部位的光脉冲有不同相移，称自相位调制(SPM)，SPM相移引起频率移动 n1为非线性折射率 L 为光纤长度

1. 光纤中的光脉冲压缩效应 光脉冲的啁啾效应设光脉冲是准单色光： w0为中心频率对脉冲前沿，，；对脉冲后沿，，。因此脉冲前、后沿频率变化是不同的，称为“啁啾”。假设介质有负色散特性，，即群速随频率升高而增加。脉冲前沿，前沿速度比脉冲中部 ( )慢脉冲后沿，后沿速度比脉冲中部快。前沿变慢后沿加快的结果造成脉宽变窄称自相位调制压缩

2.非线性薛定谔方程(NLSE) NLSE方程导出将准单色光传播常数 k按中心频率w0处展开略去高次项考虑与场强相关的非线性项：波动方程与色散关系有对应关系：得方程： NLSE 方程

2.非线性薛定谔方程(NLSE) NLSE 的孤立波解变换到以群速 vg的运动坐标系，并用y表示光场： NLSE 方程变为：设解为： NLSE 的解是受孤立波脉冲u(ξ)调制的，即包络为孤立波的光脉冲波。

3.光学孤立子的传播特性 光学孤立子设在x=0 (x=0) 处对光纤输入一个双曲正割型脉冲波系数 A 可为任意值，当为整数 N 时为稳定孤立子，N＝1 为基本孤立子。当输入光脉冲幅度超过基本孤立子幅度时，在传输中非线性压缩超过色散，光脉冲会进一步压缩，形成 N≥2时为的高阶孤立子。高阶孤立子在传播中波形发生周期变化。对于 N＝2 的二阶孤立子解二阶孤立子以ξ＝π/2 周期发生周期性波形变化。在半周期ξ＝π/4处，在孤立子主峰的两侧，各出现有一个小峰。 N＝3的三阶孤立子在传播中变化更复杂。它在 1/4 与 3/8 周期处，两侧各有一个小峰的高大尖峰，在半周期处那个高大尖峰又分裂为两个峰。

3. 光学孤立子的传播特性 光学孤立子

4. 孤立子激光器 第一个孤立子激光器目前光学孤立子产生有用两种方法： (1) 将锁模激光器产生的超短激光脉冲输入光纤，在光纤中产生受激拉曼散射形成光学孤立子。 (2) 光纤直接接入激光器反馈回路，是激光器一部分，称孤立子激光器。实验上贝尔电话实验室的 Mollenaure等人首先观察到光纤中的孤立子。所用的试验装置包括光脉冲放大与压缩两部分，使用两个光学谐振腔：色心激光器谐振腔与光纤脉冲压缩谐振腔。两个腔耦合成复合的激光谐振腔。

5. 全光型孤立子通信 光学通信线性通信利用光纤进行通信，通信容量大大增加，费用大幅度下降，是通信事业大发展的技术基础之一。目前光纤通信技术采用低强度光脉冲的线性通信。在光纤中传布不可避免地产生色散，造成光脉冲的加宽与变形，大大影响到光信息传送质量与距离。长距离高质量传送信息，需要设置许多中继站，每隔 100km 设置一个。光学孤立波通信 1973年，Hasegawa 和 Tappert 两人，利用非线性薛锷方程首次导出在光纤的反常色散区能够形成光学孤立子，建立了光学孤立波通信理论基础。由于传播中保持稳定不变的能量与波形，孤立子通信具有诱人的美妙前景。孤立子在传播中不可能没有一点能量损失。如在传输中得到能量充称全光型孤立子光纤通信。1982年，Y.Kodama等人提出拉曼泵浦技术来实现孤立子能量补充。

5. 全光型孤立子通信 掺铒光纤孤立子通信拉曼泵浦技术：两列不同频率的光波在光纤中共同传输，当能量足够大时，频率高的光波会将其部分的能量转移给低额的光波。波长980mm的半导体激光为泵浦光，通过波分复合器与波长为1550mm的信号光一起进入掺铒光纤。两路光脉冲在掺铒光纤中同方向传播。泵浦光为强光脉冲，在传播中对铒离子产生一路激发，信号光在受到激发的光纤中传播，强度逐步增强，完成对信号光能量的补充。利用拉曼泵浦技术，现已成功进行了数万 km长距离试验。

NLSE 的孤立波解 变换到以群速 vg的运动坐标系，并用y表示光场： NLSE 方程变为：设解为：或选择参数k，使 u’ 系数为零乘以u’，并积分 H-积分常数，表体系哈密顿量

NLSE 的孤立波解 两边平方根对x 微分利用椭园积分即有 NLSE 的解是受孤立波脉冲u(ξ)调制的，即包络为孤立波的光脉冲波。 H＝0

NLSE 与同宿线 对 NLSE 右边第一项体系动能，第二项体系势能 V(u) 是 u 的四次曲线，有三个奇点：两个极小点，两个极小点间的极大点在 [u, u’] 相图上，极小点对应中心点，邻域是椭圆轨线，极大点相对应鞍点。沿任一条离开鞍点的轨线出发，在绕了一圈之后回到了鞍点，鞍点是同宿点，相应的轨线为同宿线。积分

第五章 孤立波