【正文】 畢業(yè)設(shè)計（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務(wù)；學校可以采用影印、縮印、數(shù)字化或其它復(fù)制手段保存論文；在不以贏利為目的前提下，學?？梢怨颊撐牡牟糠只蛉績?nèi)容。sound(ydai,fs,bits)。plot(ydai)。title(39。Amplitude39。 xlabel(39。plot(abs(Ydai))。)。Normaliaed Magnitude(dB)39。xlabel(39。hndai=fir1(Mdai,wcdai,hanning(Mdai+1))。wsdai2=pi*。sound(yd,fs,bits)。plot(yd)。title(39。Amplitude39。 xlabel(39。plot(abs(Yd))。)。Normaliaed Magnitude(dB)39。xlabel(39。hnd=fir1(Nd,wcd,hanning(Nd+1))。Btd=wsdwpd。sound(yg,fs,bits)。plot(yg)。title(39。Amplitude39。 xlabel(39。plot(abs(Yg))。)。Normaliaed Magnitude(dB)39。xlabel(39。,hanning(Ng))。Ng=N0g+mod(N0g+1,2)。wpg=pi/2。title(39。)。subplot(2,2,3)。)。subplot(2,2,2)。原始聲音波形39。)。濾波后語音信號時域波形39。ylabel(39。axis([0,25000,0,])。濾波后語音信號頻譜圖39。ylabel(39。plot(f,abs(Yk(1:16383)))。Hk=[H H]。系統(tǒng)函數(shù)幅頻特性39。ylabel(39。f=22050*(0:16383)/32768。s39。[N,wc]=cheb2ord(wp,ws,rp,rs,39。wp=2*fs*tan(wpz*pi/2)。sync39。)wavplay(y,fs,39。Z39。xlabel(39。)Z=fft(z,32768)。z39。xlabel(39。x=[rot90(x1,3),rot90(x1,3)]。title(39。)。subplot(2,1,2)。title(39。)。figure(1)subplot(2,1,1)。)wavplay(abs(yk),fs)。yk39。xlabel(39。)yk=ifft(Yk,32768)。Yk39。xlabel(39。Yk=Z.*Hk。)figure(4)H=rot90(Hk)。Hk39。plot(f,20*log10(abs(Hk)))xlabel(39。)。[b,a]=cheby2(N,rs,wc,39。rs=100。wsz=2*500/fs。wavplay(z,fs,39。被污染語音信號頻譜圖39。ylabel(39。plot(f,abs(Z(1:16385)))。被污染語音信號時域波形39。ylabel(39。plot(abs(z))。x1=*sin(2*pi*100*t)+*sin(2*pi*220*t)。)。f39。f=22050*(0:16384)/32768。)。n39。)。濾波后語音信號時域波形39。ylabel(39。axis([0,25000,0,])。濾波后語音信號頻譜圖39。ylabel(39。plot(f,abs(Yk(1:16385)))。Hk=[H H]。系統(tǒng)函數(shù)幅頻特性39。ylabel(39。f=22050*(0:16383)/32768。s39。[N,wc]=cheb2ord(wp,ws,rp,rs,39。wp=2*fs*tan(wpz*pi/2)。sync39。)wavplay(y,fs,39。Z39。xlabel(39。)Z=fft(z,32768)。Z39。xlabel(39。)Z=fft(z,32768)。z39。xlabel(39。x=[rot90(x1,3),rot90(x1,3)]。title(39。)。subplot(2,1,2)。title(39。)。figure(1)subplot(2,1,1)。在后面的時間里，我順利完成了濾波器設(shè)計的功能。(14) 語句太長的話可以再句末加...換行 清理內(nèi)存：clear幾何分布：geornd，超幾何分布：hygernd，泊松分布：poissrnd(12) 隨機數(shù)：創(chuàng)建一個元素服從均勻分布的隨機數(shù)數(shù)組：rand轉(zhuǎn)置：transpose動畫：movie誤差圖：errorbar求矩陣維數(shù)：size或者ndims解決方法：長度不一致時，將短的矩陣增零，使長度一致；矩陣類型不同時，使用rot90(xn,3)：矩陣逆時旋轉(zhuǎn)270度，使行轉(zhuǎn)換成列或列轉(zhuǎn)換成行，或者直接使用xn’直接轉(zhuǎn)換。調(diào)試一下程序，把下標為 0 或者負數(shù)的地方修正。)。 濾波后語音信號的時域圖和頻譜圖 hanning程序的調(diào)試和運行結(jié)果。[bz,az]=bilinear(b,a,fs)。)。rp=20。其設(shè)計程序如下：wpz=[2*700/fs,2*4500/fs]。 wpzl=700Hz。z=x+y。 原語音信號的時域圖和頻譜圖。在MATLAB中通過wavplay(‘filename’)讀取語音信號數(shù)據(jù)。Hk=freqz(bz,az,16384,fs)。,39。s39。ws=2*fs*tan(wsz*pi/2)。 rs=100。 被污染語音信號的時域圖和頻譜圖。其加噪信號程序如下：t=0::。其程序如下：[y,fs,nbits]=wavread(39。 低通濾波器幅頻圖。s39。[N,wc]=cheb2ord(wp,ws,rp,rs,39。wp=2*fs*tan(wpz*pi/2)。 rp=2。x1=*sin(2*pi*5500*t)+*sin(2*pi*6000*t)。G:\Windows 39。π/2時，N為奇數(shù)時，所設(shè)計的FIR數(shù)字帶通濾波器為Ⅲ型濾波器；當N為偶數(shù)時，為Ⅳ型濾波器。數(shù)字帶通濾波器的設(shè)計也分兩種情況討論，即無相移和相移為177。π/2兩種情況。由于的頻率映射關(guān)系是根據(jù)推導(dǎo)的，所以使jΩ軸每隔2π/Ts便映射到單位圓上一周，利用沖激響應(yīng)不變法設(shè)計數(shù)字濾波器時可能會導(dǎo)致上述的頻域混疊現(xiàn)象。(2) 根據(jù)轉(zhuǎn)換后的技術(shù)指標設(shè)計模擬低通濾波器G(s)。 函)在MATLAB中，實現(xiàn)切比雪夫窗的函數(shù)為chebwin，調(diào)用格式為： w=chebwin (N,r)其中r 表示切比雪夫窗函數(shù)的傅里葉變換旁瓣幅度比主瓣低rdB(其默認值為100dB)，且旁瓣是等紋波的。sflag39。實際上，w=boxcar(N)等價于w=ones(N,1)。 窗函數(shù)法就是設(shè)計FIR數(shù)字濾波器的最簡單的方法。但是這些特性是以犧牲線性相位頻率特性為代價的，即用Butterworth、切比雪夫和橢圓法設(shè)計的數(shù)字濾波器逼近理想的濾波器的幅度頻率特性，得到的濾波器往往是非線性的。(3) 脈沖響應(yīng)不變法工作原理 沖激響應(yīng)不變法遵循的準則是使數(shù)字濾波器的單位取樣響應(yīng)與參照的模擬濾波器的脈沖響應(yīng)的取樣值完全一樣，即h(n)=ha(nT),其中T為取樣周期。即橫坐標。由于三角函數(shù)計算的重復(fù)量相當大，故FFT能極大地提高運算效率。Y可按下式計算 （21）式中，N為樣本容量，為采樣間隔。 subplot(5,1,4) plot(f,y(1:fs/2+1))。 subplot(512) stem(xn) yn=xn.*x。 x=A*exp(a*t).*sin(w*t)。 t1=0::1。 subplot(5,1,4) plot(f,y(1:fs/2+1))。 subplot(512) stem(xn) yn=xn.*x。 x=A*exp(a*t).*sin(w*t)。 t1=0::1。同時可觀察到間接用FFT求得速度要比直接用卷積求得速率要快，效率要高。Z=X1.*Y。Y=fft(y,1024)。 線性卷積結(jié)果(4) FFT 由循環(huán)卷積定理可知：對于時域序列循環(huán)卷積，可先進行FFT變換，然后頻率相乘，最后對結(jié)果進行IFFT變換，即可得到時域循環(huán)卷積結(jié)果。y(n)39。)。M=M1。xm2=[fliplr(xn2) zeros(1,686)]。)。 M=685。)當N12時，圓周卷積發(fā)生混疊；當N≥20時，圓周卷積等于線性卷積。)ylabel(39。)yn(12)=sum((xm1.*xm2)39。)ylabel(39。n=N39。n39。n=N39。m39。)title(39。figure(4)stem(yn)xlabel(39。n=039。m39。x1(m)39。)xn1=[2 0 0 7 5 7 1 7 0 1 1 9]。xlabel(39。(n)的翻轉(zhuǎn)、移位、乘積、求和的過程。n=N39。n39。n=N39。m39。)title(39。figure(6)stem(yn)xlabel(39。n=N39。m39。)title(39。figure(4)stem(yn)xlabel(39。)title(39。figure(3)stem(n,xm2)xlabel(39。)ylabel(39。)xn1=[2 0 1 2 5 7 0 5 0 1 1 2]。xlabel(39。其中翻轉(zhuǎn)采用fliplr[1]，程序如下：n=7:18。這類濾波器主要有維納濾波器、卡爾曼濾波器、線性預(yù)測濾波器以及自適應(yīng)濾波器等。數(shù)字濾波器對信號進行篩選，可通過特定頻段的信號。MATLAB是MathWorks公司于1982年推出的一款功能強大、易于使用的高效數(shù)值計算和可視化軟件，它為進行算法開發(fā)、數(shù)據(jù)計算、信號分析與可視化提供了交互式應(yīng)用開發(fā)環(huán)境，主要包括基本數(shù)學計算、編程環(huán)境(M語言)、數(shù)據(jù)可視化、GUIDE等。從發(fā)展史看，電子測量儀器經(jīng)歷了由模擬儀器、智能儀器到虛擬儀器，由于計算機性能的飛速發(fā)展，已把傳統(tǒng)儀器遠遠拋到后面，并給虛擬儀器生產(chǎn)廠家不斷帶來連鍋端的技術(shù)更新速率。通過聲音播放語句進行語音播放，可觀察聲音的變化；通過圖形處理語句和FFT得到時域圖和頻譜圖，可分析得到濾波器對頻率的濾波功能。首先通過聲音處理語句得到聲音信號的時域數(shù)據(jù)，利用FFT變換可得到頻域數(shù)據(jù)，以此進行頻率分析。 TECHNOLOGY數(shù)字信號處理課程設(shè)計 語音信號的數(shù)字濾波處理(十一) 題目： ——橢圓函數(shù)（hanning窗）濾波器 語音信號的數(shù)字濾波處理(十一) 摘要 本次設(shè)計的內(nèi)容為切比雪夫及hanning低通、高通、帶通濾波器，并利用MATLAB平臺進行設(shè)計。濾波器分別用切比雪夫II型和hanning窗設(shè)計，間接法設(shè)計IIR濾波器采用雙線性變換法，濾波器設(shè)計指標由頻譜分析得到?；谟嬎銠C