水立方

平板除中心矩形区域外修饰为超疏水表面,往平板上滴水……
平板除中心矩形区域外修饰为超疏水表面,往平板上滴水……
用算出来的一个结果做为初值,即continuation 方法。下面用具体的例子作为说明。例子来自Solving ODEs with Matlab。
方程为:
\begin{equation*} \begin{split} f'''-R[(f')^2-ff'']+RA=&0\\ R''+Rfh'+1=&0\\ \theta ''+P_ef\theta '=&0 \end{split} \end{equation*}
边界条件为:\begin{equation*} \begin{split} f(0)=f'(0)=&0\\ f(1)=f'(1)=&1\\ h(0)=h(1)=&0\\ \theta(0)=&0\\ \theta(1)=&1 \end{split} \end{equation*}
方程描述流体描述流体流过竖直管道的问题,其中 $R$ 为常数,在模型中为雷诺数。常数$P_e=0.7R$,在模型中为佩克莱特数。$A$ 为我们要计算的未知常数。
给定 $R=100, 1000, 10000$,分别解方程。$R=100$时,很容易得到结果,$R=1000, 10000$时程序不收敛,可以把$R=100$时的结果作为初值代入,则顺利得到结果。
下面给出两个程序。
很酷,但这是科学吗?
一个质点同时做两个简谐振动,两个简谐振动的方向在同一直线上,但频率略有不同,假设振幅相同,初相位都是0。假设两个简谐振动分别为
\begin{equation*} x_1=A\cos (\omega_1 t) \end{equation*}
\begin{equation*} x_2=A\cos (\omega_2 t) \end{equation*}
则合振动为
\begin{equation*} x=x_1+x_2=A\cos (\omega_1 t)+A\cos (\omega_2 t)=2A\cos\left(\frac{\omega_1-\omega_2 }{2} t\right)\cos\left(\frac{\omega_1+\omega_2 }{2}t \right) \end{equation*}
如果$latex \omega_1 \approx \omega_2$,于是$\omega_1-\omega_2 \ll \omega_1+\omega_2$,所以合振动可看作是振幅缓变的近似简谐振动。合振动振幅周期变化的现象叫拍,合振幅每变化一个周期叫做一拍。
下面通过声音演示一下拍现象。声音来自密歇根大学音乐物理讲义
两频率分别为100Hz和120Hz的声音,如下图所示:
两声音叠加,可听到拍:
[audio wav="http://joyfulphysics.scholarnet.cn/wp-content/uploads/2016/04/beatdemo.wav"][/audio]
上述文件中依次听到100Hz和120Hz的声音之后才是拍的声音。
下面是从Youtube上下载的演示视频。
测试博客的一些设置。
如果您看到这篇文章,表示您的 blog 已经安装成功.
一 | 二 | 三 | 四 | 五 | 六 | 日 |
---|---|---|---|---|---|---|
  |   |   |   | 1 | 2 | 3 |
4 | 5 | 6 | 7 | 8 | 9 | 10 |
11 | 12 | 13 | 14 | 15 | 16 | 17 |
18 | 19 | 20 | 21 | 22 | 23 | 24 |
25 | 26 | 27 | 28 |   |   |   |