1、利用MATLA昉真多普勒效應(yīng)某某某摘 要：分析多普勒效應(yīng)特性,建立數(shù)學(xué)模型,利用MATLAB軟件對其進(jìn)行仿真試驗,進(jìn)行定量分析,根據(jù)仿真試驗結(jié)果繪制出聽者接收到的信號的頻率變化曲線以及用信號處理工具 箱函數(shù)spectrogram繪制的時間頻率圖,并生成相應(yīng)的聲音.關(guān)鍵詞：Doppler effect MATLAB/ Simulink0研究背景多普勒效應(yīng)是由生在德國的奧地利物理學(xué)家多普勒(Johann Doppler 1802 1853)發(fā)現(xiàn)的.1845年,荷蘭氣象學(xué)家巴依斯巴洛Ballot)測得了聲音的多普勒效應(yīng).一輛汽車在我們身旁急馳而過,車上喇叭的音調(diào)有一個從高到低的忽然變化；站在鐵路旁邊

2、聽火車的汽笛聲也能夠發(fā)現(xiàn),火車迅速迎面而來時音調(diào)較靜止時為高,而火車迅速離去時那么音調(diào)較靜止時為低.這是日常生活中的一個多普勒效應(yīng)的例子.在天文、通信等領(lǐng)域還有眾多的例子.當(dāng)波源或觀察者相對于媒質(zhì)運動時,或者說波源和觀察者有相對運動時,觀察者接受到的震動頻率與波源震動頻率不同的現(xiàn)象,稱為多普勒效應(yīng).對于多普勒效應(yīng)的討論,一般僅限于聲源和聽者在同一直線上運動的情況.當(dāng)聲源和聽者不在同一直線上運動時,接收頻率變化比擬復(fù)雜,聽者接收到的信號波形方程也難以用解 析式表示.MATLA具有強(qiáng)大的數(shù)值計算和仿真功能以及圖形技術(shù).本文試圖從MATLAB編程的角度出發(fā),應(yīng)用MATLAB的Simulink仿真試驗

3、方法,建立仿真的試驗環(huán)境,對聲源和聽者不在同一直線上運動的情況下產(chǎn)生的多普勒效應(yīng)特性進(jìn)行分析,產(chǎn)生極好的模擬,實現(xiàn)多普勒效應(yīng)的驗證,繪制出聽者接收到的信號的頻率變化曲線以及用信號處理工具箱函數(shù) spectrogram繪制的時間頻率圖,并生成相應(yīng)的聲音.1根本原理以下公式描述了多普勒效應(yīng)現(xiàn)象的各個物理量之間的定量關(guān)系:(1-1)ff01 cos /Vs其中,f0是聲源發(fā)出的聲音的頻率；V是聽者與聲源的相對運動速度；e為速度矢量與聲源和聽者的連線夾角；vs為聲音在空氣中傳播的速度,f是聽者聽到的聲音頻率.2理論模型多普勒效應(yīng)-設(shè)聲源距離聽者的水平距離為x0 ,以V的速度沿水平方向向聽者直線駛來,其

4、軌跡與聽者的最小距離為垂直距離聲源的頻率為f0,試求聽者接受到的信號的 頻率變化曲線以及用信號處理工具箱函數(shù) ectrogram繪制的時間頻率圖,并生成相應(yīng)的 尸日.3仿真試驗下面是一段用 MATLAB的M文件依據(jù)公式1-1原理編寫的程序1-1.它會產(chǎn)生一 個名為的表現(xiàn)多普勒效應(yīng)的聲音文件,描述火車向一個距離鐵路 30米垂直距離、距火車150 水平距離的聽者開來時他聽到的聲音.將 u1送入圖3-1所示的Simulink仿真系統(tǒng), 再一次聽到該聲音,并且看到用頻譜儀表現(xiàn)的聲音頻率隨時間的變化的情況.表3-1和表3-2給出了仿真系統(tǒng)中兩個模塊的主要參數(shù)設(shè)置.圖3-1顯示聲音信號多普勒效應(yīng)的仿真系統(tǒng)

5、模型表3-1 From Multimedia File 波形文件的主要參數(shù)模塊名稱 From Multimedia File 位置 DSP System toolbox / Signal Processing Sources參數(shù)名稱參數(shù)值File name文件名Samples per output frame輸出每幀中的樣值數(shù)1024表3-2 To Audio Device 揚聲器的主要參數(shù)模塊名稱 To Audio Device 位置 DSP System toolbox / Signal Processing Sinks參數(shù)名稱參數(shù)值DeviceDefaultQueue durations

6、econds隊列延遲時間1Automatically determine Buffer size 自動確定緩沖區(qū)長度選中程序1-1%多普勒效應(yīng)程序x0=150;v=50;y0=30;c=330;f0=2000;fs=8000;t=0:1/fs:6;x_t=-x0+v.*t;y=y0;r=sqrt(x_t.A2+y.A2);costheta=x_t./r;f=f0./(1+v.*costheta./c);fmin=min(f);fmax=max(f);%m%速 m/s%m%聲音在空氣中的速度m/s%產(chǎn)首的頻率Hz%采樣頻率%規(guī)定t的范圍和步進(jìn)%火車與觀察者的水平距離%火車與觀察者的垂直距離%火車

7、與觀察者間距離%聲速矢量與火車和聽者間連線夾角的余弦%聽者接收到的頻率voltage=(f-fmin)./(fmax-fmin)*2-1;signal=*vco(voltage,fmin fmax,fs);%歸一化調(diào)頻電壓在-1 +1之間1%壓控震蕩器,由輸入電壓限制輸出信號頻率u=sin(2*pi*f0*t);%聲源發(fā)出的信號sound(u,fs);pause(5); sound(signal,fs); wavwrite(signal,fs,''); figure(1);plot(t,f);xlabel('Time/s');ylabel('Freq/H

8、z');figure(2);%播放聲源發(fā)出的信號2%暫彳亭3%»放聽者接收到的信號4%音頻數(shù)據(jù)存盤5%做出聽者接收到的信號的頻率變化曲線spectrogram(signal,kaiser(256,5),220,512,fs,'yaxis');%做出時間頻率圖 6axis(0 6,1000,spectro3000)程序難點分析1在計算機(jī)中所能存儲的數(shù)字序列,也即模擬信號必須通過取樣和量化后,變成相應(yīng) 的數(shù)字信號,才能被計算機(jī)存儲和處理.對音頻信號來說,實現(xiàn)模擬音頻信號與數(shù)字音 頻信號之間的轉(zhuǎn)化模塊就是聲卡,MATLAB可以方便地對聲卡進(jìn)行諸如采樣頻率等輸入/輸出

9、參數(shù)的配置.2播放火車靜止時的汽笛聲.使用命令"sound將數(shù)字序列以設(shè)定的采樣速率輸出到聲卡,通過聲卡轉(zhuǎn)化為模擬音頻信號.“sound用法是：sound y,F.；其中：y為取值范圍必須在 卜1, ,1區(qū)間的n行1列的數(shù)字序列單聲道輸出.Fs 為設(shè)定的采用速率,一般聲卡支持5000Hz到441000Hz.3此處暫停必不可少,且時間要足夠,以便再翻開聲音系統(tǒng),這個量與計算機(jī)硬件有 關(guān).4播放靜止的聽者聽到運動的火車的汽笛聲.5 MATLAB可以記錄的音頻信號直接保存為“ *.wav 格式.利用命令 “wavwritey,Fs,'Filename'就可以將向量 y存儲為

10、取樣頻率為 Fs的wav音頻 文件.6另一個函數(shù)發(fā)生器是壓控振蕩器VCQ,它由輸入向量來決定某個震蕩頻率的信號.此處,繪制生成信號的光譜圖.頻域測量儀器頻譜儀 Spectrum Scope計算機(jī)仿真中的頻譜儀應(yīng)用的是數(shù)字信號處理中的快速傅里葉變換FFT技術(shù),它能完成數(shù)據(jù)流從時域到頻域的變換.將時域的數(shù)據(jù)流取出一段,即FFT size快速傅里葉變換的長度確定為N,通常要求N是2的哥次方.同時需要設(shè)置相應(yīng)長度的Buffer Size修存器.N的大小即時窗的長短決定了頻譜儀的分辨率,N越大,頻率的分辨率就越高,但相應(yīng)的計算時間也愈長.希望所研究的譜線內(nèi)容出現(xiàn)在頻譜儀顯示窗的中間局部,能看到在頻率軸上

11、譜線的低端和高端的情況,以便于觀察和分析. 要做到這一點,將輸入信號的采樣頻率取為期望的頻率 顯示窗最大值的兩倍即可.注意,頻譜儀的采樣頻率與被測信號的采樣頻率要一致.4仿真結(jié)果及分析1仿真結(jié)果圖4-1所示：M文件程序運行結(jié)束時繪出的聽者接收到信號的頻率變化曲線.圖4-2所示：用信號處理工具箱函數(shù)spectrogram繪制的時間頻率圖.圖4-3所示：某個時刻仿真系統(tǒng)中的頻譜儀上顯示的功率譜.FijLOLts 2EiL» Edit i*w Ss：統(tǒng)七 1«U RF哂"g Kindr*心 I?裝/*0|國 a 圖 4-2 時間頻率圖(Swept-Frequency W

12、aveforms ), ilopp/Specrru* ScoreFi;匕Ax. 門hhLn Yhiduw 31kl 了1o o c _u o D o3 2 1 12 7C22,533.54Frequency (kHz)迎 00.5Franr.e: 3Cjp p蛆p=b*xn'll6wT=0sF,arrg: 24Ffecunc/ (kHi;dcpp/SpCijf：TUM Scope底回區(qū)Pi t 史Arte Chun«J s 此in&H Htlp宣柳齡片X33201. J -1J-2J-33-41-由b.i7DT=3sT=6s結(jié)果分析傅里葉變換是把各種形式的信號用正弦信

13、號表示,因此非正弦信號進(jìn)行傅里葉變換,會得到與原信號頻率不同的成分一一都是原信號頻率的整數(shù)倍.這些高頻信號是用來修飾頻 率與原信號相同的正弦信號,使之趨近于原信號的.所以說,頻譜上頻率最低的一個峰 往往是幅度上最高的,就是原信號頻率.由上可知,聽者接收到的信號的頻率先高于原來的汽笛聲,后低于原來的汽笛聲.5結(jié)束語在MATLAB/Simulink下,從數(shù)學(xué)模型到計算機(jī)仿真模型的轉(zhuǎn)換非常容易.本文利用 Simulink結(jié)合M文件編程的方法,使其同時具備圖形界面的直觀性和字符界面的強(qiáng)大功能. 事實上,所有Simulink的模塊以及系統(tǒng)構(gòu)建、仿真參數(shù)等均可通過編程語句實現(xiàn).針對聲源和聽者不在同一直線上運動情況,利用MATLAB