海洋波動方程典型例題傅氏解的可視化實驗

張春玲　胡松

摘? 要：波動被稱為海洋的節拍，是一種常見且非常重要的海洋現象。其數理方程的傅氏解以復變函數的級數為基礎，物理意義不易被學生理解。本文借助于海洋數據處理中常用的Matlab編程軟件，選取波動方程有界混合問題的典型例題，以直觀地方式再現其傅氏解的物理意義。

關鍵詞：波動方程;可視化;分離變量;海洋現象

引言

海洋波動是海水運動的重要形式之一，從海表面至海洋深處，均存在尺度不同的各種波動現象。按照成因可以分為風浪、涌浪、地震波、長重力波、潮波、行星波等;按照發生的位置可以分為表面波、內波和邊緣波;按照波長相對于水深的大小，可以分為淺水波和深水波;按照波形是否傳播又可以分為前進波和駐波[1] 。但不管什么形式的波動，都是能量傳播的表現，某一點的波動必然帶動其鄰近質點，通過周圍物質的周期運動來傳送或傳播能量，故波動的特征是隨時間和空間具有周期性變化[2] 。在數學上，通常利用波動方程來描述這一運動規律，并借助于傅里葉解法（也稱分離變量法）來求解有界波動問題。但實際海洋中的波動嚴格來說都不是真正的周期運動，而是一種十分復雜的現象。其傅里葉解的形式為無窮級數，既無限多個正弦波或余弦波的疊加，這也從理論上刻畫了這一特征[3] 。對于這個傅氏解的物理意義，單從公式上，學生很難理解，更不易與實際海洋波動現象聯系在一起。隨著計算機技術的發展，很多軟件及仿真技術用于數理方程的可視化研究，其中海洋數據處理軟件Matlab就不失為一種較好地工具[4-6]。本文以一個有界波的混合問題為例，利用Matlab對其傅氏解進行可視化處理，以便更直觀地解釋其物理意義。

1 典型例題的傅氏解

一維有界波動滿足的混合問題為： ，利用分離變量法求得的解為? ，表示實際波形是由一系列頻率成倍增長、相位不同、振幅也不同的駐波疊加而成的。假設研究區域為 ，考慮兩端為固定邊界（邊界上的波動振幅為零），具有一定初始位移，初始速度為零，在不受外力作用下的自由振動。其中，當 ，初試位移為 而當 ，初位移為 。由題目的條件，該波動滿足方程組 ，分離變量法求得的方程的解是? 。

2 傅氏解的可視化實驗

首先給定參數波速a，初始位移折線的斜率h及研究區域長度l，然后以一定時間步長和空間步長將研究區域和波動時間離散化，在此基礎上，把所有波動進行疊加。實驗代碼如下：

實驗得到的不同波數所對應的波形如圖1所示。由圖可以看出，波數確定的單個波波形均為正弦波，理論解具有明顯的周期性。而這些波疊加在一起就使得海面隨時間推移，在不同的位置點上表現出不同的波動現象（圖2）。波數為10時，由初試擾動引起的波動能量在擾動點達到最大值，而后向兩邊傳播（n=30），在到達邊界后（n=100），前進波與反射波能量抵消（n=102），對于某個固定的位置點，其波動形式，在反射波與前進波疊加前后，表現為相反的振幅大?。╪=180;n=200）。

3 結論

傅里葉解法是求解海洋中運動方程的基本解法，其以無窮級數的形式給出了理論解。通過matlab的可視化功能，將其理論解進行直觀地展示，可以更形象地揭示波動方程的運動情況。而且，可視化實驗程序中，可以自由調節各個參數的大小，如波速、波數、時間步長和空間步長等，對應于不同參數，研究區域各個點的波動情況也一目了然。

參考文獻

[1] 馮士筰，李鳳歧，李少菁.海洋科學導論[M].1999.

[2] D.A.羅斯，b.A，羅斯，et al.海洋學導論[M].科學出版社，1984.

[3] 姜禮尚.數學物理方程講義.第3版[M].2007.

[4] 彭芳麟.數學物理方程的MATLAB解法與可視化[M].2004.

[5] 郝玉華.一維弦振動方程的可視化處理[J].鹽城工學院學報（自然科學版），2006（04）：16-19.

[6] 徐彬，XuBin.Matlab在復變函數與積分變換課堂教學中的應用[J].湖北理工學院學報，2016，32（3）：68-72.5（2）：40-40.

作者簡介：張春玲，（1981-），女，山東，講師，博士，海洋科學，主要從事海洋數據分析研究。

胡松，（1978-），男，廣西，教授，博士，海洋科學學院副院長，分管教學。

基金項目：上海海洋大學教學改革與研究項目（A1-2005- 20-300348）;國家自然科學基金項目“合成孔徑雷達圖像中白冠覆蓋率的特征及參數化”（41606196）。