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