【正文】 函數(shù)嗎? 答：得到A6(z)的{ki}值如下：k(6)= k(5)= k(4)= k(3)= k(2)= k(1)= 從{ki}的值可以得到傳輸函數(shù)A6(z)是穩(wěn)定的，因?yàn)榉答佅禂?shù)的平均幅值小于整體。disp(39。)disp(39。)。disp(39。答：級聯(lián)實(shí)現(xiàn)框圖： num = input(39。sos = zp2sos(z,p,k)，生成如下有限沖激響應(yīng)傳輸函數(shù)的一個級聯(lián)實(shí)現(xiàn): H1(z)=2+10z^(1)+23z^(2)+34z^(3)+31z^(4)+16 z^(5)+4z^(6)畫出級聯(lián)實(shí)現(xiàn)的框圖。 ：由此看出，都在單位圓內(nèi)。分母系數(shù) = 39。disp(roots(num4))。disp(roots(num2))。)。zplane(num3,1)。Type 1 FIR Filter39。)。)。Time index n39。)。)。Time index n39。num2 = [b 81 81 fliplr(b)]。用設(shè)計的穩(wěn)定無限沖激響應(yīng)帶通濾波器的傳輸函數(shù)的參數(shù)α和β，生成一個二階無限沖激響應(yīng)帶阻濾波器的傳輸函數(shù)HBS（z）。用MATLAB計算并畫出它們的增益響應(yīng)，驗(yàn)證設(shè)計的濾波器是否滿足指標(biāo)。)。n = 0:M1。K = 39。grid。Hz = conv(Hz,[1 1])。M = 39。plot(w/pi,g)。ylabel(39。plot(w/pi,abs(h))。Enter the filter time shift N: 39。Time index n39。y = 2*fc*sinc(2*fc*n)。)。k = n+。ylabel(39。stem(k,y)。)。Poles:39。num = [ ]。(z)在()穩(wěn)定,這意味著H[n]是絕對可和,我們看到交替和指數(shù)衰減的脈沖響應(yīng)。])。輸入樣本數(shù) L: 39。群延遲39。num = [ 1]。如下圖1這是一個帶通濾波器。Phase in radians39。)。subplot(2,1,1)plot(w/pi,abs(h))。第二篇：數(shù)字信號處理第四章介紹第四章 線性時不變離散時間系統(tǒng)的頻域分析一、傳輸函數(shù)和頻率響應(yīng) clear。plot(w/pi,tmpY, 39。plot(tmpX,tmpY, 39。plot(tmpX,tmpY, 39。figure(4)。omega /pi39。title(39。plot(w2/pi,Phase)。)。plot(w/pi,g)。wgts = max(Dp,Ds)*[1/Ds 1/Dp 1/Ds]。DEV = [Ds Dp Ds]。Fp1 = 1800。)。UPhase = unwrap(Phase)。)。% Compute the frequency response w2 = 0:pi/511:pi。omega /pi39。)。Ds = 。As = 50 dB增益和相位響應(yīng)如下：從圖中可以看出設(shè)計的濾波器滿足要求。ylabel(39。)。xlabel(39。Gain Response39。axis([0 180 5])。% Show the Numerator Coefficients disp(39。As = 50。)。UPhase = unwrap(Phase)。)。% Compute the frequency response w2 = 0:pi/511:pi。omega /pi39。)。Dp = Ds。程序：% Program Q7_23 % Use Kaiser window to design a linear phase Lowpass % FIR Digital Filter meeting the design specification given % in .% % It is clear from the statement of the question that Mitra % wants us to use()and()for this % isn39。ylabel(39。)。xlabel(39。T3 = *ones(1,length(w2))。title(39。plot(w/pi,g)。A = [1 1 0 0]。% Design spec as given in = 2*10^3。 線性相位帶通FIR濾波器的階數(shù)估算，則使用kaiserord的結(jié)果為同樣，它也有矛盾。 ，使用remezord函數(shù)的結(jié)果為N=：firpmord([2000 2500],[1 0],[ ],10000)：在這里，firpmord給了一個比凱爾更大比凱瑟更小一點(diǎn)的結(jié)果。其中模擬過渡帶寬|FpFs|和Δω的關(guān)系：Δω=2pi*|FpFs|/FT。對于每種情況，研究希爾伯特變換器的頻率響應(yīng)的擺動行為。波紋的數(shù)量與微分器的長度之間的關(guān)系——兩者成正比。二型濾波器不可能是高通濾波器，因?yàn)楸仨氃趜=1處有零點(diǎn)，意味著w=+π。最大波紋的高度與濾波器的長度之間的關(guān)系——最大波紋的高度與長度無關(guān)。寫出所產(chǎn)生的傳輸函數(shù)的準(zhǔn)確表達(dá)式。寫出所產(chǎn)生的傳輸函數(shù)的準(zhǔn)確表達(dá)式。)。% 計算增益響應(yīng)[g, w] = gain(num,den)。% 顯示傳輸函數(shù)disp(39。% 巴特沃斯帶阻濾波器的設(shè)計Ws = [ ]。指標(biāo)如下:7 kHz的抽樣率, kHz的通帶邊界頻率， kHz的阻帶邊界頻率，, dB的通帶波紋,50 dB的最小阻帶衰減。第一篇：數(shù)字信號處理第七章介紹第七章 數(shù)字濾波器設(shè)計：無限沖激響應(yīng)濾波器的階數(shù)的估計。：無限沖激響應(yīng)濾波器設(shè)計。)。disp(den)。增益，dB39。濾波器的未畸變的相位響應(yīng)及群延遲響應(yīng)如下：。濾波器的未畸變的相位響應(yīng)及群延遲響應(yīng)如下：。波紋的數(shù)量與濾波器的長度之間的關(guān)系——波紋的數(shù)量減少與長度成正比。使其成了一個線性相位FIR濾波器。波紋的數(shù)量與微分器的長度之間有什么關(guān)系，最大波紋的高度與濾波器的長度之間有什么關(guān)系? 答：幅度響應(yīng)分別如下：從中可以觀察到由于吉布斯現(xiàn)象產(chǎn)生的幅度響應(yīng)的擺動行為。num=[fliplr(b)0 b]。(b)線性相位FIR(c)：濾波器階數(shù)和采樣頻率的關(guān)系為–對于一個給定的模擬過渡帶寬，采樣頻率的增加導(dǎo)致估算階數(shù)也相應(yīng)增加，朝下一個整數(shù)取整。這種估計經(jīng)常和FIR設(shè)計的凱瑟窗一起用。所以答案不唯一，可以選擇其中一個。：增益和相位響應(yīng)如下:=，圖如下： 漢寧窗：布萊克曼窗：切比雪夫窗： 程序如下：% Program Q7_22 % Use ParksMcClellan to design a linear phase Lowpass % FIR Digital Filter meeting the design specification given % in .%Compute and plot the gain function.%Compute and plot the unwrapped phase response.%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear。% zero to one, not zero to = [0 Wp Ws 1]。figure(1)。)。T2 = *ones(1,length(w2))。grid。Phase Response39。)。欲滿足指標(biāo)，應(yīng)調(diào)節(jié)N= 用凱澤窗設(shè)計一個有限沖激響應(yīng)低通濾波器。Ds = 10^(As/20)。Numerator Coefficients are 39。xlabel(39。)。omega /pi39。figure(3)。Unwrapped Phase(rad)39。Ws = 。h = fir1(N,Wn,Ftype,Win)。grid。title(39。grid。Phase Response39。)。Ws = 。Dp = 。Numerator Coefficients are 39。xlabel(39。)。omega /pi39。figure(3)。Unwrapped Phase(rad)39。% Design spec as given in = 1500。A = [0 1 0]。A2 = [0 0 1 1 0 0]。figure(1)。Gain in dB39。Phase = angle(Hz)。)。xlabel(39。)。% vertical line at passband edge freq tmpX = ones(1,length(tmpY))*wp2。% vertical line at stopband edge freq tmpX = ones(1,length(tmpY))*ws2。% horizontal line at Dp tmpY = ones(1,length(w))*(Ds)。用kaiserord估計濾波器的階數(shù)N=73。h = freqz(num, den, w)。Amplitude39。ylabel(39。den = [ ]。w = 0:pi/511:pi。ylabel(39。L = input(39。 脈沖響應(yīng)的樣本39。（）式（）式由圖可知：這些情節(jié)由impz給生成的因果的脈沖響應(yīng)實(shí)現(xiàn)的H(z)。clf。disp(39。Hit to continue...39。k = n+。)。y = 2*fc*sinc(2*fc*n)。Time index n39。n = [:1:]。xlabel(39。N = input(39。h = freqz(y, [1], w)。)。[g,w] = gain(num,1)。title([39。for i=1:K。plot(w/pi,g,w/pi,ThreedB,t1,t2)。title([39。)。omega /pi39。 。用MATLAB計算并畫出你所設(shè)計的濾波器的增益響應(yīng)，并驗(yàn)證它確實(shí)滿足給定的條件。num1 = [b 81 fliplr(b)]。xlabel(39。Type 1 FIR Filter39。Amplitude39。xlabel(39。Type 3 FIR Filter39。Amplitude39。title(39。subplot(2,2,3)。Type 4 FIR Filter39。)。)。den = input(39。由此我們可以總結(jié)出H1(z)穩(wěn)定，H2(z)不穩(wěn)定 ：由此看出，都在單位圓內(nèi)。[z,p,k] = tf2zp(num,den)。答：Numerator coefficient vector = [3,8,12,7,2,2] Denominator coefficient vector = [16,24,24,14,5,1] sos = 0 0 級聯(lián)實(shí)現(xiàn)框圖：：畫出級聯(lián)實(shí)現(xiàn)的框圖。[r2,p2,k2] = residue(num,den)。極點(diǎn)在39。并聯(lián)II型39。disp(p2)。畫出兩種實(shí)現(xiàn)的框圖。den = input(39。k(ii)= a1(ii+1)。disp(alpha) 使用程序 (z)的 GrayMarkel級聯(lián)格型實(shí)現(xiàn)參數(shù)如下： 晶格參數(shù)和前饋乘數(shù)分別如下：對應(yīng)GrayMarkel的結(jié)構(gòu)框圖如下：使用程序P6_3，從這些格型參數(shù)可以得到傳輸函數(shù)H1(z)是穩(wěn)定的，因?yàn)樗懈裥蛥?shù)的平方值比整體的小。den = input(39。[num2,den2] = latc2tf(k,alpha)運(yùn)行結(jié)果如下： k = alpha = 0.***849 。在高度信息化的今天，語音信號處理是信息高速公路、多媒體技術(shù)、辦公自動化、現(xiàn)代通信及智能系統(tǒng)等新興領(lǐng)域應(yīng)用的核心技術(shù)之一。同時也是眾構(gòu)成思想交流和感情溝通的最主要的途徑。一方面，對聲學(xué)語音學(xué)統(tǒng)計模型的研究逐漸深入，魯棒的語音識別、基于語音段的建模方法及隱馬爾可夫模型與人工神經(jīng)網(wǎng)絡(luò)的結(jié)合成為研究的熱點(diǎn)。當(dāng)已知一個數(shù)字音源后，可以利用計算機(jī)的處理能力，用數(shù)字的方式通過計算模擬回聲效應(yīng)。AIC23通過麥克風(fēng)輸入或者立體聲音頻輸入采集模擬信號, 并把模擬信號轉(zhuǎn)化為數(shù)字信號, 存儲到DSP的內(nèi)部RAM中,以便DSP處理。收發(fā)時鐘信號CLKX1和CLKR1由A IC23的串行數(shù)據(jù)輸入時鐘BCLK提供, 并由A IC23的幀同步信號LRCIN、LRCOUT啟動串口數(shù)據(jù)傳輸。用戶可以根據(jù)電阻阻值調(diào)節(jié)增益的大小，使語音輸入輸出達(dá)到最佳效果。 修改前程序圖 修改前程序圖 頻譜分析 左聲道及右聲道波形 實(shí)驗(yàn)分析所以，從本實(shí)驗(yàn)可知當(dāng)已知一個數(shù)字音源后，可以利用計算機(jī)的處理能力，用數(shù)字的方式通過計算模擬回聲效應(yīng)。參考文獻(xiàn)[1][M].北京:中國水利水電出版社,2004.[2]王安民,陳明欣,[M].北京:清華大學(xué)出版社,2007 [3]彭啟琮,[M].成都:電子科技大學(xué)出版社,1997 [4]李宏偉,[J].解放軍理工大學(xué)學(xué)報,2001(1):41~44 [5][M].梁曉雯,裴小平,李玉虎,:清華大學(xué)出版社,2006 [6][M].北京:機(jī)械工業(yè)出版社,2003比較圖4和圖5,可以看到1200Hz以上的頻譜明顯得到了抑制。數(shù)字信號處理技術(shù)正飛速發(fā)展，它不但自成一門學(xué)科，更是以不同形式影響和滲透到其他學(xué)科：它與國民經(jīng)濟(jì)息息相關(guān)，與國防建設(shè)緊密相連；它影響或改變著我們的生產(chǎn)、生活方式，因此受到人們普遍的關(guān)注。第七單元的課程我們學(xué)習(xí)了線性相位有限脈沖響應(yīng)（FIR）數(shù)字濾波器，窗函數(shù)法設(shè)計有限脈沖響應(yīng)（FIR）數(shù)二、專業(yè)認(rèn)識和未來規(guī)劃通信工程是一門工程學(xué)科，主要是在掌握通信基本理論的基礎(chǔ)上，運(yùn)用各種工程方法對通信中的一些實(shí)際問題進(jìn)行處理。數(shù)字信號處理課程屬于專業(yè)基礎(chǔ)課，所涵蓋的內(nèi)容主要有：離散時間信號與系統(tǒng)的基本概念及描述方法，離散傅立葉變換及快速傅立葉變換，數(shù)字濾波器結(jié)構(gòu)及設(shè)計等。近年來，國外及國內(nèi)有些學(xué)校對一般電類專業(yè)該課程的教學(xué)主要強(qiáng)調(diào)應(yīng)用性學(xué)習(xí)，主要介紹數(shù)字信號處理的用途和用法，而對其深奧的理論推導(dǎo)僅做一般介紹，并給學(xué)生提供進(jìn)行實(shí)驗(yàn)的機(jī)會，以激發(fā)學(xué)生對該課程的興趣和學(xué)習(xí)主動