setaspectratio函数_有哪些堪称「数学物理建模软件技巧及神器」，却鲜为人知？...

现在，已经有不少软件可以作物理的动画，进行过程模拟。本教程将结合Mathematica的动画功能，模拟大学物理中两个机械波源所发出波的干涉,以及Matlab在物理学中的应用。

通常，物理教科书上对波的干涉都只讲同频率的两个波源所发出的波的干涉，其实，那只是一个特殊情况。波的干涉是波的叠加性的一种表现，不论两个波源的频率如何，叠加都是存在的，因此，干涉也是存在的。只不过干涉现象的观察手段受到限制，人们只能观察到“稳定的干涉现象”，而同频率的条件可以造成稳定的干涉。如果波源的频率不同，干涉现象就不是稳定的，那么，这时的干涉现象会是怎样的？当然，对此问题的回答可以从干涉的数学表达式上得到，但不够直观。下面，就以在平面上传播的波(如水波)为例，将干涉的数学表达式转化成动画，模拟一般的干涉现象，并观察其特征。

1.波的叠加公式

设在水面上相距为2α的两个地方分别放置半径均为r0的两个柱体，它们各自受到垂直方向的余弦振动源的驱动，周期性地上下起伏，因而引起两个波动向周围传播，并在水面上互相叠加，形成干涉。波源的布置如图1-1所示。图1-1 水面上两个柱形波源

模拟干涉成功的要点是合理设置波的表达式。设左右两个波源的振动圆频率分别为ω1，和ω2。平面上任一点的位置矢量r={x，y}两个波源的中心分别是P1={-a，0}和P2={a，0}。在区域

波的振幅为1，各点都作同相位的简谐振动，波的表达式为

而在此区域以外，则按普通的同心圆方式发散传播。这时波的表达式为

其中，

是波速。在公式(1-2)中，已经考虑了波在传播中的相位延迟，以及传播中能量的守恒，波的衰减问题暂不考虑。

两列波叠加后的波幅表达式为

2.模拟相同频率波的干涉

在某一时刻t，p(r，t)是一个二变量的函数，用函数Plot3D[]可以画出其联合分布，模拟水波的空间波形(如果观察光波的干涉，需要波幅的平方)。该函数里有些参数需要重新设置，以使波形更加逼真，其中，观察的视角参数ViewPoint很关键，注意参数ViewPoint的使用。为了演示波形函数中各参数作用,先来看演示波形的Mathematic示例代码:

Animate[Plot3D[

Sin[Sqrt[x^2 + y^2] + t*2*Pi], {x, -8 Pi, 8 Pi}, {y, -8 Pi, 8 Pi},

PlotRange -> 10, PlotPoints -> 50, AspectRatio -> 1,

Boxed -> False], {t, 2, 0}, AnimationRunning -> False]

代码运行以后,动画在启动时,不自动启动,所以将 AnimationRunning 的值,设置为False,具体功能说明如下图所示:

波形示例如下图所示:

Mathematic操作如下: