孫麗男，張馥菊

(1.黑河學(xué)院;2.大連理工大學(xué);3.哈爾濱師范大學(xué))

0 引言

弦振動問題是數(shù)學(xué)物理方程中的經(jīng)典問題，弦振動理論在生產(chǎn)、生活中有著較為廣泛的應(yīng)用.如對于斜拉橋、高空觀光纜車運行安全和游客人身安全而言，檢測并控制鋼索工作過程中的張力是至關(guān)重要的，而此類問題最終歸結(jié)為弦振動問題的研究.另外，工廠車間進(jìn)行絡(luò)紗時，可以觀察到紗會以相當(dāng)快的速度沿著中心軸進(jìn)行退繞，達(dá)到一定條件時便會形成多節(jié)氣圈，這種現(xiàn)象也可以通過弦振動理論加以說明.弦振動是波動的一種特殊形式，因此弦振動理論在聲學(xué)中也占有舉足輕重的地位，弦振動理論與現(xiàn)代先進(jìn)的軟件和硬件相結(jié)合，在許多領(lǐng)域諸如音樂物理學(xué)、材料學(xué)和系統(tǒng)分析中也都得到了廣泛的應(yīng)用［1］.弦振動問題是由弦振動方程附加初邊值條件構(gòu)成的，弦振動方程是一類重要的偏微分方程，它關(guān)于時間和空間變量的導(dǎo)數(shù)均是二階的，在解決實際問題的過程中，遇到的弦振動問題通常既有初始條件又有邊界條件，是一類混合問題.對于一些簡單的混合弦振動問題可以通過分離變量法、疊加原理、齊次化原理及傅里葉變換等方法求其解析解或者是廣義解，但是對于一些復(fù)雜的弦振動問題，傳統(tǒng)的方法求解起來就比較困難，或者求出的解的形式相對復(fù)雜，與實際相結(jié)合較為困難.這時，就可以通過數(shù)值方法進(jìn)行求解，借助計算機(jī)進(jìn)行輔助分析，從而使求解過程變的輕松，并且能夠?qū)⑺媒Y(jié)果以圖形等形式直觀地顯示出來［3］.

文中考慮了以下兩端固定的有界弦振動問題

其中，上述問題的解析解可以利用分離變量等方法進(jìn)行求解［2］，但為了更直觀地展示上述問題解的形式，以下將采用數(shù)值方法進(jìn)行求解并利用功能強(qiáng)大的Matlab軟件進(jìn)行輔助分析，最后將解析解和數(shù)值解進(jìn)行比較.

1 數(shù)值分析

該節(jié)首先給出求解問題(1)的有效的數(shù)值方法，然后利用這些方法及Matlab軟件求出數(shù)值解，并且比較幾種數(shù)值方法的優(yōu)劣.

1.1 數(shù)值方法

求解偏微分方程的數(shù)值方法有很多，比較常見的方法就是有限元法和差分法，文獻(xiàn)［4］中利用有限元法對(1)進(jìn)行了求解.文中將利用差分法對(1)進(jìn)行求解，主要采用顯式、隱式兩種差分格式.

1.1.1 區(qū)域剖分

首先，對定解問題(1)的求解區(qū)域W={(x，t)|0＜x＜l，t＞0}進(jìn)行剖分.在區(qū)間［0，1］內(nèi)取n-1個等分點xi=ih，i=1，…，n-1將區(qū)間n等分，其中x0=0，xn=l，h為空間步長.同樣，在時間區(qū)間內(nèi)插入若干等分點tj=jτ，j=1，2，…，其中，t0=0，τ為時間步長.用兩族平行直線x=xi，1＜i＜n;t=tj，j=1，2，… 將區(qū)域W剖分成矩形網(wǎng)格，其中(xi，tj)為結(jié)點，以下將討論問題(1)在結(jié)點處的數(shù)值解.

1.1.2 方程離散

用差分方法求解偏微分方程的關(guān)鍵是對方程的離散，下面采用顯示、隱式差分格式對問題(1)中方程進(jìn)行離散，得到如下形式的差分方程.

(i)顯示差分方程

(i)隱式差分格式

其中i=1，2，…，n-1;j=1，2，…，令則上式可改寫為

1.1.3 初值條件的處理

問題(1)對應(yīng)的初值條件可進(jìn)行如下處理

1.1.4 邊值條件的處理

問題(1)對應(yīng)的邊值可進(jìn)行如下處理

1.2 數(shù)值實驗

下面將利用計算機(jī)實驗對問題(1)的數(shù)值解進(jìn)行輔助分析，并將其可視化，實驗主要采用的軟件是MATLAB軟件.

首先，根據(jù)顯示差分格式編制MATLAB程序模擬問題(1)的數(shù)值解，獲得動態(tài)結(jié)果如圖1所示.

圖1

在圖1實驗中取初始位移為

初始速度ψ(x)取為零，時間步長τ=0.0005，空間步長h=0.0024.通過對上述結(jié)果進(jìn)行觀察，發(fā)現(xiàn)初始狀態(tài)的波平均分成兩個波向相反方向傳播，在它們到達(dá)弦的端點之前，和波在無限長弦上的傳播是一致的，到達(dá)弦的端點以后，兩個波同時發(fā)生反射，出現(xiàn)了半波損失，這是有別于無限長弦的振動的，接下來，被反射回來的兩個波進(jìn)行相向運動并在弦的中間相遇，相遇后兩個波疊加繼續(xù)傳播，整個波的傳播過程清晰明了，實現(xiàn)了問題(1)的數(shù)值解的可視化，有利于進(jìn)一步研究更為復(fù)雜的弦振動問題.

對于問題(1)還可以根據(jù)隱式差分格式編制MATLAB程序模擬它的數(shù)值解，相應(yīng)程序的編碼要比顯示差分格式復(fù)雜的多，涉及到三對角方程組的求解等技術(shù)性問題，具體程序如下

clear;clc;

format short e

a=1;l=1;n=420;=0.0005;

N=input(＇請輸入運行次數(shù) N的值:＇);

h = l/n;x = zeros(n+1，1);c =a^2*τ^2/h^2;

for i=1:n

x(i+1)=i*h;

end

u=zeros(n+1，N);

foc i=181:240

u(i，1)=0.05*sin(pi*x(i)*7/l);

u(i，2)=0.05*sin(pi*x(i)*7/l)+1/2*c*0.05*(sin(pi*x(i+1)*7/l)-2*sin(pi*x(i)*7/l) + sin(pi*x(i -1)*7/l));

end

h=plot(x，u(:，1)，＇linewidth＇，3);axis(［0，1， - 0.05，0.05］);set(h，＇EraseMode＇，＇xor，＇MackecSize＇，18)

for i=1:N

set(h，＇XData＇，x，＇YData＇，u(:，1));drawnow;

B=zeros(n-1，1);A=zeros(n-2，1);R=zeros(n-2，1);S=zeros(n-1，1);E=zeros(n-1，1);

F=zeros(n-2，1);G=zeros(n-2，1);I=zeros(n-1，1);H=zeros(［n-1，n-1］);

for i=1:n-2

E(i)=-(1+c);F(i)=1/2*c;G(i)=1/2*c;I(i)=u(i，1);

end

E(n-1)=-(1+c);H=diag(E)+diag(F，1)+diag(G，-1);I(n-1)=u(n-1，1);J=H*I;

foc i=1:n-2

B(i)=1+c;A(i)=-1/2*c;C(i)=-1/2*c;S(i)=2*u(i，2)+J(i);

end

B(n-1)=1+c;S(n-1)=2*u(n-1，2)+J(n-1);u(1，2)=0;u(n+1，2)=0;S(1)=u(1，2)+J(1);

%追趕法

Y=zecos(n-1，1);X=Y;v=zecos(1，n-1);g=v;v(1)=B(1);Y(1)=S(1);A(2:n-1)=A;

foc i=2:n-1

g(i)=A(i)/v(i-1);v(i)=B(i)-g(i)*C(i-1);Y(i)=S(i)-g(i)*Y(i-1);

end

U=zecos(n-1);L=eye(n-1);

foc i=1:n-2

L(i+1，i)=g(i+1);

end

foc i=1:n-2

U(i，i)=v(i);

U(i，i+1)=C(i);

end

U(n-1，n-1)=v(n-1);X(n-1)=Y(n-1)/v(n-1);foc i=n-2:-1:1

X(i)=(Y(i)-C(i)*X(i+1))/v(i);

end

u(1:n-1，3)=X;u(:，1)=u(:，2);u(:，2)=u(:，3);

end

運行結(jié)果的動態(tài)展示如圖2所示.

圖2

利用隱式差分格式得到的數(shù)值結(jié)果與顯式差分格式總體上是一致的，二者均較好地逼近了解析解［2］.

其中的系數(shù)是:當(dāng)n≠7時，

當(dāng)n=7時另外，B=0.利用MATLABn軟件也可以清晰的比較解析解中級數(shù)項數(shù)的多少對解的精確度的影響，當(dāng)級數(shù)達(dá)到50項以上才可以有較高的精確度，另外，根據(jù)模擬結(jié)果可以觀察到解析解中級數(shù)的每一項是一個駐波，不同頻率的駐波疊加成行波，生動形象.

2 結(jié)束語

以上利用MATLAB軟件及有限差分方法求解兩端固定的有界弦振動問題，首先利用區(qū)域轉(zhuǎn)化的思想對方程進(jìn)行離散，然后編制MATLAB程序進(jìn)行求解，并將數(shù)值解可視化，使對問題(1)的解的物理意義有了更深刻的認(rèn)識，另外，通過與解析解比較，也充分說明了所選擇數(shù)值方法的可靠性和計算機(jī)輔助分析的重要性.

［1］ 李韻，郭怡文，呂郁文.基于 Matlab環(huán)境的弦振動方程的圖像與音效模擬.科協(xié)論壇［J］，2009(7):74-75.

［2］ 谷超豪，李大潛，等.數(shù)學(xué)物理方程:第三版［M］.北京:高等教育出版社，2012.

［3］ 孫麗男，張馥菊.基于MATLAB的一類輸運問題的數(shù)值分析.哈爾濱師范大學(xué)自然科學(xué)學(xué)報，2012，28(1):14-17.

［4］ 高峰，陳君若，等.有界弦振動方程的有限元方法求解［J］.機(jī)械設(shè)計，2008，25(11):15-18.