【正文】 Function Library)基本函數(shù)：求和、正弦、余弦和復(fù)數(shù)運算等；特殊函數(shù)：矩陣求逆、矩陣特征值、貝塞爾函數(shù)和快速付里葉變換等； 3）MATLAB語言：(MATLAB Language)一種高級編程語言，包括控制流的描述、函數(shù)、數(shù)據(jù)結(jié)構(gòu)、輸入輸出及面對對象編程； 4）句柄圖形：(Handle Graphics) 可以對各種圖形對象進行更為細膩的修飾和控制，建立完整的圖形界面的應(yīng)用程序。 工作空間窗口：(Workspace Window) 用于儲存各種變量和結(jié)果的空間，顯示變量的名稱、大小、字節(jié)數(shù)及數(shù)據(jù)類型，對變量進行觀察、編輯、保存和刪除。*雙擊變量名；*選擇該窗口工具欄上的打開圖標；*鼠標指向變量名，點擊鼠標右鍵，彈出選擇菜單，然后選項操作。只有在當前目錄或搜索路徑下的文件及函數(shù)可以被運行或調(diào)用。還可用cd命令在命令窗口設(shè)置當前目錄，如： cd c:\mydir 可將c盤上的mydir目錄設(shè)為當前工作目錄。啟動平臺：(Launch Pad) 幫助用戶方便地打開和調(diào)用MATLAB的各種程序、函數(shù)和幫助文件。事實上Help工具完全具備它的功能。它將數(shù)值分析、矩陣計算、科學(xué)數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多強大功能集成在一個易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計以及必須進行有效數(shù)值計算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計語言（如C、Fortran）的編輯模式，代表了當今國際科學(xué)計算軟件的先進水平。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計算方面首屈一指。MATLAB的基本數(shù)據(jù)單位是矩陣，它的指令表達式與數(shù)學(xué)、工程中常用的形式十分相似，故用MATLAB來解算問題要比用C、FORTRAN等語言完成相同的事情簡捷得多，并且MATLAB也吸收了像Maple等軟件的優(yōu)點，使MATLAB成為一個強大的數(shù)學(xué)軟件。可以直接調(diào)用，用戶也可以將自己編寫的實用程序?qū)氲組ATLAB函數(shù)庫中方便自己以后調(diào)用，此外許多的MATLAB愛好者都編寫了一些經(jīng)典的程序，用戶可以直接進行下載就可以用。39。%讀出信號，采樣率和采樣位數(shù) y=y(:,1)。%算出y的長度Y = fft(y,sigLength)。 %Y乘Y的共軛,得到功率譜halflength=floor(sigLength/2)。 %采樣頻率的橫坐標figure。xlabel(39。)。 %計算時間軸figure。xlabel(39。)。為了減少在傳輸時的耗損，人們一般是先對傳輸信號進行特殊處理，然后再傳遞。我們知道音頻信號的頻率在我們的聽覺范圍內(nèi)，音頻信號在無線傳輸?shù)倪^程中會受到各種聲音的干擾而產(chǎn)生能量消耗，因此限制了傳輸?shù)木嚯x，所以在現(xiàn)實的傳送過程中要現(xiàn)對語音信號進行調(diào)制再發(fā)送出去。A是任意振幅。1之間。圖42調(diào)制后的語音信號由圖可知，經(jīng)過調(diào)制后的語音信號頻譜都搬移到fc=20000的頻段。fc=1200。Ap=1。 wc=2*fc/fs。 [n,wn]=ellipord(wp,wc,Ap,As)。 freqz(b,a,512,fs)。因為在這一處理過程中，讓低頻信號全部通過，所以稱為低通濾波。低通濾波器(lowpassfilter)容許低頻信號通過,但減弱(或減少)頻率高于截止頻率的信號的通過。當使用在音頻應(yīng)用時,它有時被稱為高頻剪切濾波器或高音消除濾波器。 程序如下：[y,fs,bits]=wavread(39。)。fp=1000。rp=。fc=40000。wp=2*fp/fc。 %將模擬域轉(zhuǎn)化成數(shù)字域。value39。[N,Wc]=buttord(wp,ws,rp,rs)。[num,den]=butter(N,Wc)。else %選擇radiobutton2，則制做切比雪夫Ⅰ型濾波器。 %估算切比雪夫Ⅰ型濾波器的階數(shù)N和截止頻率Wc。 %求傳輸函數(shù)的分子和分母多項式的系數(shù)。 %根據(jù)傳輸函數(shù)的分子和分母多項式的系數(shù)得到模擬濾波器，并將y通過該濾波器，得到x。 %將濾波后的信號x進行回放。 %將x信號進行快速傅里葉變換。 %在axes3坐標軸上畫出x信號的波形圖。 %在axes4坐標軸上畫出頻譜圖。所以聲音稍微有些發(fā)悶、低沉，但是很接近原來的聲音。通過這次練習是我進一步了解了信號的產(chǎn)生、采樣及頻譜分析的方法。讓我感受到只有在了解課本知識的前提下，才能更好的應(yīng)用這個工具；并且熟練的應(yīng)用MATLAB也可以很好的加深我對課程的理解，方便我的思維。同時我相信，進一步加強對MATLAB的學(xué)習與研究對我今后的學(xué)習將會起到很大的幫助。39。%讀出信號，采樣率和采樣位數(shù)。%我這里假設(shè)你的聲音是雙聲道，我只取單聲道作分析，如果你想分析另外一個聲道，請改成y=y(:,2) sigLength=length(y)。 Pyy = Y.* conj(Y) / sigLength。 f=Fs*(0:halflength)/sigLength。plot(f,Pyy(1:halflength+1))。Frequency(Hz)39。 t=(0:sigLength1)/Fs。plot(t,y)。Time(s)39。附錄二：%抽取的源程序[y,fs,bs]=wavread(39。,[25000,75000])。lg=length(y)。yx=yx(1:lg/2)。figure(1)。plot(y)。原始信號39。subplot(212)。yf=yf(1:lg/2)。ycq=zeros(1,lg)。 %通過循環(huán)，每隔80個點將抽取后的值賦值為原函數(shù)的 ycq(i)=y(i)。figure(2)。plot(ycq)。抽取后的信號39。subplot(212)。ycqf=ycqf(1:lg/2)。附錄三：“河南成建學(xué)院”語譜圖源程序：frmsize = 256。39。 % 讀取語音文件x = filter([1 ],[1],x)。 frmt = floor(n/frmsize)。 % 分幀xFFT = abs(fft(x))。 % 奈奎斯特頻率以內(nèi)是無效的xFFTdB = 20*log10(xFFT+eps)。 % 能量太高的截斷xFFTdB(xFFTdB40) = 40。 % 能量整理到 [0,60] 之間，當然你可以整理到[0,255]范圍subplot(2,1,1)。 title(39。)。image(flipud(xFFTdB))。e:\39。% sound(y,fs,bits)。%采樣點數(shù)40000figure(1)。plot(y)。原始信號波形39。subplot(212)。title(39。)。%載波頻率y1=modulate(y,fc,fs,39。)。%采樣點數(shù)40000figure(2)。plot(y1)。調(diào)制后信號波形39。subplot(212)。title(39。)。y2=y+y1。%采樣點數(shù)40000figure(3)。plot(y2)。調(diào)制信號+原始信號波形39。subplot(212)。title(39。)。figure(4)。fc=1200。Ap=1。wc=2*fc/fs1。[n,wn]=ellipord(wp,wc,Ap,As)。freqz(b,a,512,fs1)。X=fft(x,40000)。subplot(211)。title(39。)。plot(abs(X))。低通濾波后信號頻譜39。 sound(x,40000)。e:\39。x=double(x)。%將幅度歸一化到[1,1]framelen=240。%設(shè)定幀長為240幀移為80amp1=10。zcr1=10。maxsilence=3。status=0。silence=0。temp2=enframe(x(2:length(x)),framelen,frameinc)。diffs=(temp1temp2)。%計算過零率amp=sum(abs(enframe(x(1:length(x)1),framelen,frameinc)),2)。amp2=min(amp2,max(amp)/8)。x2=0。 goto=0。 status=2。 count=count+1。 count=count+1。 count=0。 count=count+1。 if silencemaxsilence %靜音不夠長,沒有結(jié)束 count=count+1。 %語音長度太短,認為是噪聲 status=0。 count=0。 end end case 3, break。for i=1:length(amp)1。(zcr(i+1)) j=j+1 jiesu1(j)=i。jiesu1(i)~=1 jiesu(j)=jiesu1(i)。 if ((zcr(i))amp。 endendj=1for i=1:length(kaishi1)1 if((kaishi1(i+1)kaishi1(i))80) kaishi(j)=kaishi1(i)。x2=x1+count1。speech39。line([x1*frameinc x1*frameinc],[1,1],39。,39。)。color39。red39。subplot(312)plot(amp)。energy39。line([x1,x1],[min(amp),max(amp)],39。,39。)。color39。red39。subplot(313)plot(zcr)。zcr39。line([x1 x1],[min(zcr),max(zcr)],39。,39。)。color39。red39。