kdv matlab,求助：Matlab关于孤子波仿真的问题

孤子这个名词首先是在物理的流体力学中提出来的。1834年在爱丁堡到格拉斯哥的运河上，苏格兰海军工程师罗素(J.Scott Russell)观察到一种奇特的水波，并在1844年发表的一份报告中描述了这种奇特水波，并称这种波为孤立波(Solitary wave)。他是这样描述的：“我看到两匹骏马拉着一条船沿运河迅速前进。当船突然停止时，随船一起运动的船头处的水堆并没有停止下来。它激烈地在船头翻动起来，随即突然离开船头，并以巨大的速度向前推进。一个轮廓清晰又光滑的水堆，犹如一个大鼓包，沿着运河一直向前推进，在行进过程中其形状与速度没有明显变化。我骑马跟踪注视，发现它保持着起始时约30英尺(1英尺＝25.4厘米)长，1-1.5英尺高的浪头，约以每小时8-9英里的速度前进。后来它的高度逐渐减低，经过约一英里(1英里＝1.609千米)的追逐后，在运河的拐弯处消失了”。

为了探究上述的水波鼓包，罗素在水槽的一端用一重锤垂落入水中，对重锤激起的水浪的运动情况进行了反复的观察，如图1.1所示。他发现这种水浪与运河中出现的奇特水波是本相同。通过实验，他还总结出水波的移动速度v、水的深度d及水波幅度A之间的关系：

公式一.jpg (2 KB, 下载次数: 0)

2009-5-21 20:47 上传

(1.1)

其中B为一某比例常数。这实验结果说明，水波的运动速度与波幅的高度有关，波幅高的速度