【正文】 也認(rèn)識到自己的不足，以后要進(jìn)一步學(xué)習(xí)。在以后的學(xué)習(xí)中還需要深入了解這方面的內(nèi)容。再利用所學(xué)的操作，發(fā)現(xiàn)所寫的程序還是沒有能夠運(yùn)行，通過不斷地調(diào)試，運(yùn)行，最終得出了需要的結(jié)果。在設(shè)計的過程中遇到了很多問題，剛剛開始時曾天真的認(rèn)為只要把以前的程序改了參數(shù)就可以用了，可是問題沒有我想象中的那么簡單，單純的搬程序是不能解決問題的。原因可能為噪聲與信號頻率相近無法濾除。頻率圖中，女生高 27頻成分較多。通過帶通濾波器后，聲音比較適中。，把高頻濾除，聲音變得比較低沉。d:39。d:39。%播放濾波后信號 grid。%顯示標(biāo)表示的意義 grid。幅值/A39。)。%標(biāo)題 xlabel(39。低通濾波器濾波后信號的頻譜圖39。plot(f, abs(fftwave(n1)))。%顯示標(biāo)表示的意義 grid。幅值/A39。)。%信號的時域波形title(39。subplot(3,2,5)。%信號送入濾波器濾波，ys為輸出 fftwave=fft(ys)。%butter函數(shù)（num為分子系數(shù)den為分母系數(shù)）[h,w]=freqz(num,den)。%最小衰減[N,Wn]=buttord(Wp,Ws,Rp,Rs)。%阻帶截至頻率 Rp=1。%采樣頻率 Wp=2*1000/fs。grid。幅值/db39。)。xlabel(39。疊加信號的頻譜圖39。plot(f,abs(y0(n1)))。y0=fft(yy)。)。ylabel(39。時間/t39。)。%信號的時域波形title(39。LineWidth39。subplot(3,2,3)。f=[(n11)*(2*pi/N)]/pi。%返回采樣點數(shù) df=fs/N。%將兩個聲音疊加：wavplay(yy,fs)。)。%對錄入信號進(jìn)行采樣[x1,fs,bits]=wavread(39。 C: 39。fs=8000。噪聲信號的頻譜圖39。)。%噪聲信號的頻譜圖 ylabel(39。%快速傅立葉變換subplot(3,2,2)。)。ylabel(39。)。%信號的時域波形 title(39。subplot(3,2,1)。f=[(n11)*(2*pi/N)]/pi。%返回采樣點數(shù) df=fs/N。%對語音信號進(jìn)行采樣wavplay(x1,fs)。 C: 39。)。%按任意鍵播放語音信號 wavwrite(x,Fs,39。pause。按Enter鍵回放錄音39。)。%錄制語音信號fprintf(39。double39。)。%暫停命令 fprintf(39。,time)。%錄音時間 fprintf(39。%加網(wǎng)格分析：就時域圖看，男生的時域圖中振幅比女生的高，對于頻譜圖女生的高頻成分比較多Fs=8000。)。%顯示橫坐標(biāo)表示的意義 ylabel(39。頻率w/pi39。)。%原始信號的頻譜圖 title(39。%快速傅立葉變換subplot(2,2,4)。%顯示縱坐標(biāo)的表示意義 grid。幅值/A39。)。%錄制信號的時域波形title(39。%頻帶寬度 subplot(2,2,3)。%采樣間隔 n1=1:N/2。%播放原始信號 N=length(y)。)。%加網(wǎng)格[y,fs,bits]=wavread(39。)。%顯示橫坐標(biāo)表示的意義 ylabel(39。頻率w/pi39。)。%原始信號的頻譜圖 title(39。%快速傅立葉變換 subplot(2,2,2)。%顯示縱坐標(biāo)的表示意義 grid。幅值/A39。)。%錄制信號的時域波形title(39。subplot(2,2,1)。f=[(n11)*(2*pi/N)]/pi。%返回采樣點數(shù) df=fs/N。%播放原始信號,解決落后半拍wavplay(x,fs)。 C: 39。)。%按任意鍵播放語音信號wavwrite(x,Fs,39。pause。按Enter鍵回放錄音39。)。%錄制語音信號 fprintf(39。double39。)。%暫停命令 fprintf(39。,time)。%錄音時間 fprintf(39。 錄制一段異性的聲音進(jìn)行頻譜分析Fs=8000。%sound(x, Fs, bits)。其中低通的濾波后與原來聲音沒有很大的變化，其它兩個都又明顯的變化sound(xlow,Fs,bits)。%網(wǎng)格wavplay(ys1,8000)。)。%顯示標(biāo)表示的意義 title(39。幅值/A39。)。xlabel(39。%在四行兩列的第八個窗口顯示圖形plot(f, abs(fftwave(n1)))。%網(wǎng)格fftwave=fft(ys1)。)。ylabel(39。頻率w/pi39。)。%繪制后信號的時域的圖形 title(39。subplot(4,2,6)。ys1= fftfilt(h,x)。,kaiser(n+1,beta))。h=fir1(n,wn/pi,39。wn=[wn1 wn2]。wn1=(ws1+wc1)/2。ws2=2*pi*3000/fs。wc2=2*pi*3200/fs。beta=。%加網(wǎng)格wavplay(ys,8000)。)。%顯示標(biāo)表示的意義title(39。幅值/A39。)。xlabel(39。%在四行兩列的第七個窗口顯示圖形 plot(f,abs(fftwave(n1)))。%網(wǎng)格fftwave=fft(ys)。)。%標(biāo)題 ylabel(39。高通濾波后信號的時域波形39。)。ylabel(39。頻率w/pi39。%在四行兩列的第五個窗口顯示圖形 plot(ys)。ys=fftfilt(h,x)。,kaiser(n+2,beta))。h=fir1(n,wn/pi, 39。wn=(ws+wc)/2。wc=2*4800/fs。beta=。%加網(wǎng)格wavplay(ys,8000)。)。%顯示標(biāo)表示的意義title(39。幅值/A39。)。%繪制模值 xlabel(39。%網(wǎng)格subplot(4,2,4)。)。ylabel(39。頻率w/pi39。)。%信號的時域波形title(39。subplot(4,2,2)。%信號送入濾波器濾波，ys為輸出 fftwave=fft(ys)。[h1,w]=freqz(h)。band39。n=ceil(*pi /width)。width=wswc。wc=2*pi*1000/fs。%加網(wǎng)格3．1低通濾波器濾波 fs=8000。)。%顯示縱坐標(biāo)表示的意義 title(39。幅值 39。)。%加標(biāo)題 xlabel(39。原始信號的頻譜圖39。plot(f,abs(y0(n1)))。%加網(wǎng)格y0=fft(x)。)。%加標(biāo)題 ylabel(39。原始信號的時域波形39。plot(x)。%頻帶寬度 figure(4)。%采樣間隔 n1=1:N/2。%播放原始信號 N=length(x)。)。123．將聲音信號送入濾波器濾波x=wavread(39。)。ylabel(39。頻率（HZ）39。)。title(39。[h,w]=freqz(Bz,Az,512,Fs1*1000)。[Bz,Az]=bilinear(B,A,Fs1)。s39。)。[N,wc]=buttord(wap,was,rp,rs,39。wap=2*tan(wp/2)。T=1。wp=2*pi*fp/Fs。rp=1。fs=[1000,3200]。MATLAB程序如下。)。ylabel(39。頻率（HZ）39。)。title(39。[h,w]=freqz(Bz,Az,512,Fs1*22050)。[Bz,Az]=bilinear(B,A,Fs1)。s39。high39。)。10[N,wc]=buttord(wap,was,rp,rs,39。wap=2*tan(wp/2)。T=1。wp=2*pi*fp/Fs。rp=1。fs=5000。9 巴特沃思低通濾波器MATLAB程序如下。)。ylabel(39。頻率（HZ）39。)。title(39。[h,w]=freqz(Bz,Az,512,Fs1*22050)。[Bz,Az]=bilinear(B,A,Fs1)。s39。)。[N,wc]=buttord(wap,was,rp,rs,39。wap=2*tan(wp/2)。ws=2*pi*fs/Fs。rs=100。Fs=22050。%低通濾波fp=1000。grid。振幅39。)。%標(biāo)題 xlabel(39。截取后的聲音時域圖形39。%從4000點截取到24000結(jié)束 plot(x)。C:39。這樣就可以解決這個問題。%加標(biāo)題grid。原始信號的頻譜圖39。)。%顯示橫坐標(biāo)表示的意義 ylabel(39。頻率w/pi39。)。%原始信號的頻譜圖 title(39。subplot(2,1,2)。%加網(wǎng)格y0=fft(x)。)。%加標(biāo)題 ylabel(39。原始信號的時域波形39。plot(x)。%頻帶寬度 figure(2)。%采樣間隔 n1=1:N/2。%存儲語音信號N=length(x)。C:39。%暫停命令wavplay(x,Fs)。)。%文字提示 fprintf(39。錄音結(jié)束39。)。x=wavrecord(time*Fs,Fs,39。錄音中......39。%文字提示 pause。按Enter鍵錄音%ds39。%抽樣頻率 time=3。2fh四．具體實現(xiàn)在MATLAB中用wavrecord、wavread、wavplay、wavwrite對聲音進(jìn)行錄制、讀取、回放、存儲。奈奎斯特采樣定理：要使實信號采樣后能夠不失真還原，采樣頻率必須大于信號最高頻率的兩倍：Ws179。0WWs/2W179。a239。X(jW)239。?(jW)與X(jW)的關(guān)系：理想抽樣信號的頻譜是連續(xù)信號頻譜的Xaa周期延拓,重復(fù)周期為Ws(采樣角頻率)。Xa(jWjkWs)Tk=165。?Xa(jW)=229。dT(t)則有?(jW)=1X(jW)*D(jW)XaaT2p1165。xXa(jW)171。d(tnT)165。dT(t)=頻域描述：利用傅氏變換的性質(zhì)，時域相乘頻域卷積，若n=165。165。xa(nT)d(tnT)xn=165。連續(xù)信號經(jīng)理想抽樣后時域、頻域發(fā)生的變化（理想抽樣信號與連續(xù)信號頻譜之間的關(guān)系）理想抽樣信號能否代表原始信號、如何不失真地還原信號即由離散信號恢復(fù)連續(xù)信號的條件（抽樣定理）理想采樣過程描述： 時域描述：?a(t)=xa(t)dT(t)=229。 DF線性相位的概念，即線性相位對h(n)、H(w)及零點的約束，了解四種FIR DF的頻響特點。二．設(shè)計內(nèi)容錄制一段個人自己的語音信號，并對錄制的信號進(jìn)行采樣；畫出采樣后語音信號的時域波形和頻譜圖；給定濾波器的性能指標(biāo)，采用窗函數(shù)法和雙線性變換法設(shè)計濾波器，并畫出濾波器的頻率響應(yīng)；然后用自己設(shè)計的濾波器對采集的信號進(jìn)行濾波，畫出濾波后信號的時域波形和頻譜，并對濾波前后的信號進(jìn)行對比，分析信號的變化；回放語音信號；換一個與你性別相異的人錄制同樣一段語音內(nèi)容，分析兩段內(nèi)容相同的語音信號頻譜之間有什么特點；再錄制一段同樣長時間的背景噪聲疊加到你的語音信號中，分析疊加前后信號頻譜的變化，設(shè)計一個合適的濾波器，能夠把該噪聲濾除；三．設(shè)計原理，利用函數(shù)wavrecord()，wavwrite(),wavread(),wavplay()對語音信號進(jìn)行錄制，存儲，讀取，回放。使用zplane獲得G(z)的零極點分布圖：G(z)全通和的分解：G(z)的功率補(bǔ)充傳輸函數(shù)H(z)的表達(dá)式如下： 生成一個五階因果有界實低通橢圓傳輸函數(shù)G(z)的全通和的分解。所得到的傳輸函數(shù)和式()給出的傳輸函數(shù)相同嗎? 答：運(yùn)行結(jié)果： k = *** 0.***293 alpha = *** 0.***485 。實現(xiàn)式()給出的傳輸函數(shù)。)。% here is the lattice/ladder realization from the transfer f: [k,alpha] = tf2latc(num,den)% now check inversion disp(39。num = num/den(1)。Denominator coefficient vector = 39。)。所得到的傳輸函數(shù)和式()給出的傳輸函數(shù)相同嗎？ 答：程序如下： format long num = input(39。用該程序?qū)崿F(xiàn)式()所示的傳輸函數(shù)。 使用程序 (z)的 GrayMarkel級聯(lián)格型實現(xiàn)參數(shù)如下：對應(yīng)GrayMarkel的結(jié)構(gòu)框圖如下：使用程序P6_3，從這些格型參數(shù)可以得到傳輸函數(shù)H2(z)是穩(wěn)定的，因為所有格型參數(shù)的平方值比整體的小。)。disp(k)disp(39。格型參數(shù)是39。a1(1:ii+1)=(a1(1:ii+1)k(ii)*a1(ii+1:1:1))/(1k(ii)*k(ii))。for ii = N:1:1, alpha(N+2ii:N+1)= alpha(N+2ii:N+1)alpha(Nii+1)*a1(2:ii+1)。a1 = den/den(1)。N = length(den)1。分母系數(shù)向量 = 39。)。：無限沖激響應(yīng)傳輸函數(shù)的GaryMarkel實現(xiàn)num = input(39。顯示實現(xiàn)的框圖。答：并聯(lián)I型框圖：并聯(lián)II型框圖：：全通傳輸函數(shù)的實現(xiàn):As(z)是一個穩(wěn)定的傳輸函數(shù)嗎? 答：運(yùn)行結(jié)果：k(5)= k(4)= k(3)= k(2)= k(1)= 2從{ki}的值我們可以得到傳輸函數(shù)A5(z)是穩(wěn)定的，因為對所有的1:A6(z)足一個穩(wěn)定的傳輸函數(shù)嗎? 答：得到A6(z)的{ki}值如下：k(6)= k(5)= k(4)= k(3)= k(2)= k(1)= 從{ki}的值可以得到傳輸函數(shù)A6(z)是穩(wěn)定的，因為反饋系數(shù)的平均幅值小于整體。答：并聯(lián)I型框圖：并聯(lián)II型框圖：()所示因果無限沖激響應(yīng)傳輸函數(shù)的兩種不同并聯(lián)形式實現(xiàn)。()所示因果無限沖激響應(yīng)傳輸函數(shù)的兩種不同并聯(lián)形式實現(xiàn)。)。disp(39。)。disp(39。)。)disp(39。disp(39。)。disp(39。)。disp(39。)。)disp(39。disp(39。[r1,p1,k1] = residuez(num,den)。分母系數(shù)分量 = 39。)。答：級聯(lián)實現(xiàn)框圖： num = input(39。只用四個乘法器生成級聯(lián)框圖：：級聯(lián)和并聯(lián)實現(xiàn)： 畫出級聯(lián)實現(xiàn)的框圖。顯示新的級聯(lián)結(jié)構(gòu)的框圖。，生成如下有限沖激響應(yīng)傳輸函數(shù)