【正文】 關(guān)系。)。s39。x1=*sin(2*pi*100*t)+*sin(2*pi*220*t)。)。其加噪信號(hào)程序如下：t=0::。rs=50。 原語音信號(hào)的時(shí)域圖和頻譜圖 帶通濾波器的幅頻圖和相頻圖 低通濾波器的幅頻圖和相頻圖5 總結(jié)在這次課程設(shè)計(jì)中，我遇到了許多問題，這些問題是我們初學(xué)者經(jīng)常會(huì)遇到的問題。(6) 變量顯示方式更改：format long/short/bank散點(diǎn)圖：scatter/scatter3矩陣范數(shù)：norm，條件數(shù)：cond，初等變換：rref，特征值：eig/eigs創(chuàng)建一個(gè)元素服從正態(tài)分布的隨機(jī)數(shù)數(shù)組：randnylabel(39。原語音信號(hào)頻譜圖39。f=22050*(0:16384)/32768。f39。ws=2*fs*tan(wsz*pi/2)。Hk39。Yk39。)wavplay(abs(yk),fs)。subplot(2,1,2)。xlabel(39。Z39。[N,wc]=cheb2ord(wp,ws,rp,rs,39。)。)。(3) 切比雪夫帶通濾波器完整程序：[y,fs,nbits]=wavread(39。plot(f,abs(Y(1:16385)))。n39。)。s39。)figure(4)H=rot90(Hk)。)yk=ifft(Yk,32768)。subplot(2,2,1)。plot(x2)。wcg=(wpg+wsg)/2/pi。yg=fftfilt(hng,x2)。加噪聲聲音原始波形39。[magd,radd]=freqz(hnd)。title(39。title(39。Normalized Frequency(rad)39。Frequency（hz）39。畢業(yè)設(shè)計(jì)（論文）原創(chuàng)性聲明和使用授權(quán)說明原創(chuàng)性聲明本人鄭重承諾：所呈交的畢業(yè)設(shè)計(jì)（論文），是我個(gè)人在指導(dǎo)教師的指導(dǎo)下進(jìn)行的研究工作及取得的成果。：任務(wù)書、開題報(bào)告、外文譯文、譯文原文（復(fù)印件）。title(39。title(39。[magdai,raddai]=freqz(hndai)。加噪聲聲音原始波形39。yd=fftfilt(hnd,x2)。N0d=ceil(*pi/Btd)。)。)。wsg=pi/4。t=1/fs。)wavplay(abs(yk),fs)。Yk39。Hk39。ws=2*fs*tan(wsz*pi/2)。f39。y=y(1:20191,:)z=x+y。原語音信號(hào)時(shí)域波形39。)。f39。f39。wp=2*fs*tan(wpz*pi/2)。xlabel(39。x=[rot90(x1,3),rot90(x1,3)]。title(39。yk39。xlabel(39。plot(f,20*log10(abs(Hk)))xlabel(39。)。f=22050*(0:16384)/32768。)。ylabel(39。plot(abs(y))。(2) 輸入輸出文件：save/load，在屏幕上顯示文件：type在MATLAB中通過wavplay(‘filename’)讀取語音信號(hào)數(shù)據(jù)。wsz=[2*500/fs,2*4700/fs]。其程序如下：[y,fs,nbits]=wavread(39。rp=10。G:\Windows 39。ws=2*fs*tan(wsz*pi/2)。4 設(shè)計(jì)程序的調(diào)試和運(yùn)行結(jié)果 切比雪夫低通程序的調(diào)試和運(yùn)行結(jié)果 。在MATLAB中，實(shí)現(xiàn)巴特里特窗的函數(shù)為bartlett，調(diào)用格式為： w=bartlett (N)在MATLAB中，實(shí)現(xiàn)凱塞窗的函數(shù)為kaiser，調(diào)用格式為： w=kaiser (N,beta)其中beta為窗函數(shù)的參數(shù)β。它在設(shè)計(jì)FIR數(shù)字濾波器中有很重要的作用，正確地選擇窗函數(shù)可以提高設(shè)計(jì)數(shù)字濾波器的性能，或者在滿足設(shè)計(jì)要求的情況下，減小FIR數(shù)字濾波器的階次。 IIR設(shè)計(jì)原理(1) 切比雪夫?yàn)V波器原理切比雪夫?yàn)V波器的幅頻特性具有等波紋特性。 stem(yn)。 subplot(5,1,3) stem(yn)。z=ifft(Z,1024)figure(2)stem(z)axis([0,700,1,])。figure(1)stem(yn)xlabel(39。yn=zeros(1,687)。x1(Nm)39。)ylabel(39。)ylabel(39。程序如下：n=0:11。)ylabel(39。)ylabel(39。m39。 M=17。目前已經(jīng)有許多較成熟的頻譜分析軟件，如SpectraLAB、RSAVu、dBFA等。最后將被污染信號(hào)通過設(shè)計(jì)的濾波器，實(shí)現(xiàn)濾波功能，得到濾波后的語音信號(hào)。但如果有用信號(hào)和無用信號(hào)，或有用信號(hào)和噪聲的頻譜相互重疊，那么經(jīng)典濾波器則不能實(shí)現(xiàn)理想的濾波性能。figure(2)stem(n,xm1)xlabel(39。)ylabel(39。)ylabel(39。y(n)39。m39。)ylabel(39。y(n)39。)title(39。)。G:\Windows XP 菜單命令39。 a=sqrt(2)*n*pi。 a=sqrt(2)*n*pi。但上式的計(jì)算效率很低，因?yàn)橛写罅康闹笖?shù)(等價(jià)于三角函數(shù))運(yùn)算，故實(shí)際中多采用快速Fourier變換(FFT)。有限長單位沖激響應(yīng)(FIR)數(shù)字濾波器具有以下優(yōu)良的特點(diǎn)：(1) 可在設(shè)計(jì)任意幅度頻率特性濾波器的同時(shí)，保證精確、嚴(yán)格的線性相位特性。當(dāng)sflag＝symmetric時(shí)，為對(duì)稱采樣；當(dāng)sflag＝periodic時(shí)，為周期采樣，此時(shí)hann函數(shù)計(jì)算N+1個(gè)點(diǎn)的窗，但是僅返回前N個(gè)點(diǎn)。若所設(shè)計(jì)的數(shù)字濾波器是低通的，那么上述設(shè)計(jì)工作可以結(jié)束，若所設(shè)計(jì)是高通、帶通或帶阻濾波器，那么還需進(jìn)行以下步驟。當(dāng)相移為177。設(shè)計(jì)指標(biāo)：wpz=4500Hz。[bz,az]=bilinear(b,a,fs)。 wsz=500Hz。)。設(shè)計(jì)指標(biāo)：wpz=4500Hz。s39。2. 如果 a 是函數(shù)，則這是因?yàn)閙atlab在所有已添加的路徑中都無法找到該函數(shù)對(duì)應(yīng)的m文件而導(dǎo)致的，對(duì)此，把該m文件移動(dòng)到當(dāng)前路徑下再運(yùn)行便可[7]。剛開始自己編的程序運(yùn)行就錯(cuò)誤很多，當(dāng)時(shí)都沒信心能把這個(gè)課程設(shè)計(jì)完成，但什么都要自己面對(duì)的。)Y=fft(y,32768)。figure(2)subplot(2,1,1)。)。title(39。)。Hk=[H fliplr(H)]。subplot(2,1,1)。plot(abs(y))。ylabel(39。)。sync39。,39。H=Hk39。plot(abs(yk))。xlabel(39。Y39。title(39。)。[bz,az]=bilinear(b,a,fs)。f=22050*(0:16382)/32768。n39。X1=fft(x1)。mag2=abs(X2)。subplot(2,2,1)。濾波后信號(hào)頻譜39。加噪聲聲音濾波后波形39。)。)。wpdai2=pi*。title(39。subplot(2,2,3)。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫的成果作品。本人完全意識(shí)到本聲明的法律后果由本人承擔(dān)。title(39。)。wsdai2=pi*。Amplitude39。Normaliaed Magnitude(dB)39。sound(yg,fs,bits)。 xlabel(39。xlabel(39。title(39。subplot(2,2,2)。ylabel(39。plot(f,abs(Yk(1:16383)))。f=22050*(0:16383)/32768。sync39。)Z=fft(z,32768)。title(39。)。xlabel(39。Yk=Z.*Hk。)。wavplay(z,fs,39。被污染語音信號(hào)時(shí)域波形39。)。n39。axis([0,25000,0,])。Hk=[H H]。s39。)wavplay(y,fs,39。Z39。xlabel(39。subplot(2,1,2)。在后面的時(shí)間里，我順利完成了濾波器設(shè)計(jì)的功能。求矩陣維數(shù)：size或者ndims解決方法：長度不一致時(shí)，將短的矩陣增零，使長度一致；矩陣類型不同時(shí)，使用rot90(xn,3)：矩陣逆時(shí)旋轉(zhuǎn)270度，使行轉(zhuǎn)換成列或列轉(zhuǎn)換成行，或者直接使用xn’直接轉(zhuǎn)換。[bz,az]=bilinear(b,a,fs)。 wpzl=700Hz。Hk=freqz(bz,az,16384,fs)。 rs=100。 低通濾波器幅頻圖。 rp=2。數(shù)字帶通濾波器的設(shè)計(jì)也分兩種情況討論，即無相移和相移為177。由于的頻率映射關(guān)系是根據(jù)推導(dǎo)的，所以使jΩ軸每隔2π/Ts便映射到單位圓上一周，利用沖激響應(yīng)不變法設(shè)計(jì)數(shù)字濾波器時(shí)可能會(huì)導(dǎo)致上述的頻域混疊現(xiàn)象。sflag39。由于三角函數(shù)計(jì)算的重復(fù)量相當(dāng)大，故FFT能極大地提高運(yùn)算效率。 x=A*exp(a*t).*sin(w*t)。 x=A*exp(a*t).*sin(w*t)。Y=fft(y,1024)。M=M1。)當(dāng)N12時(shí)，圓周卷積發(fā)生混疊；當(dāng)N≥20時(shí)，圓周卷積等于線性卷積。n=N39。)title(39。x1(m)39。n=N39。)title(39。)title(39。)ylabel(39。這類濾波器主要有維納濾波器、卡爾曼濾波器、線性預(yù)測濾波器以及自適應(yīng)濾波器等。通過聲音播放語句進(jìn)行語音播放，可觀察聲音的變化；通過圖形處理語句和FFT得到時(shí)域圖和頻譜圖，可分析得到濾波器對(duì)頻率的濾波功能?；谟?jì)算機(jī)軟硬件平臺(tái)的虛擬儀器可代替?zhèn)鹘y(tǒng)的測量儀器，如示波器、邏輯分析儀、信號(hào)發(fā)生器、頻譜分析儀等。動(dòng)態(tài)演示兩個(gè)序列進(jìn)行線性卷積x1(n)﹡x2(n)的翻轉(zhuǎn)、移位、乘積、求和的過程。xm2=[fliplr(xn2) zeros(1,18)]。figure(5)stem(n,xm3)xlabel(39。xlabel(39。動(dòng)態(tài)演示兩個(gè)序列進(jìn)行圓周卷積[2, 3]x1(n)⊙x2figure(3)stem(n,xm2)xlabel(39。figure()stem(n,xm1)xlabel(39。m39。因?yàn)橹苯佑肍FT進(jìn)行1024點(diǎn)卷積大于兩序列長度，所以可用線性卷積替代圓周卷積，其程序如下：n=7:686。yn(687)=sum((xm1.*xm3)39。X1=rot90(X,3)。 xn=ones(1,length(t))。 xn=ones(1,length(t))。根據(jù)采樣定理，只有頻率不超過的信號(hào)才能被正確采集，即Fourier變換的結(jié)果中頻率大于的部分是不正確的部分，故不在頻譜圖中顯示。 (3) 允許設(shè)計(jì)多通帶(多阻帶)系統(tǒng)。sflag39。當(dāng)相移為177。wsz=2*4700/fs。在MATLAB中通過wavplay(‘filename’)讀取語音信號(hào)數(shù)據(jù)。wp=2*fs*tan(wpz*pi/2)。 濾波后語音信號(hào)的時(shí)域圖和頻譜圖 切比雪夫帶通程序的調(diào)試和運(yùn)行結(jié)果。 rs=50。其濾波程序如下：Yk=Z.*Hk。MATLAB中常用語句：(1) 學(xué)會(huì)用help