【正文】 11)*(2*pi/N)]/pi。)。y0=fft(yy)。)。%阻帶截至頻率 Rp=1。subplot(3,2,5)。%顯示標表示的意義 grid。)。d:39。頻率圖中，女生高 27頻成分較多。在以后的學習中還需要深入了解這方面的內容。離散序列通常用x（n）來表示，自變量必須是整數(shù)。與連續(xù)的正弦信號不同，正弦序列的自變量n必須為整數(shù)。1n=0編制程序產(chǎn)生單位沖激序列d(n)=237。title(39。stem(n,y)。0編制程序產(chǎn)生單位階躍序列u(n)=237。單位階躍序列39。title(39。y=[zeros(1,NM+1)ones(1,NM)]。246。w=2*pi。grid on。title(39。phi=0。分析:該序列具有周期性,、編制程序產(chǎn)生復正弦序列x(n)=e(2+j學號后兩位）n，并繪出其圖形。grid on。n=0:N1。(2)N=10。grid on。X(e)=+(n)=cosjwjw++++（２）計算32點序列5pn16，0≦n≦31的32點和64點DFT，分別繪出幅度譜圖形，并繪出該序列的DTFT圖形。imagX=imag(H)。magnitude39。grid。相位 part39。)。subplot(224)。)。k=0:1:N1。)。subplot(2,1,2)。在MATLAB中系統(tǒng)函數(shù)的零極點可通過函數(shù)roots得到，也可借助函數(shù)tf2zp得到，tf2zp的語句格式為[Z,P,K]=tf2zp(B,A)其中，B與A分別表示H（z）的分子與分母多項式的系數(shù)向量。MATLAB提供了求離散時間系統(tǒng)頻響特性的函數(shù)freqz，調用freqz的格式主要有兩種。x2=*exp(1*j*w)+*exp(2*j*w)+*exp(3*j*w)+*exp(4*j*w)+*exp(5*j*w)。零極點圖39。振幅39。相位部分39。二、實驗原理在IIR濾波器的設計中，常用的方法是：先根據(jù)設計要求尋找一個合適的模擬原型濾波器，然后根據(jù)一定的準則將此模擬原型濾波器轉換為數(shù)字濾波器。MATLAB提供了一組標準的數(shù)字濾波器設計函數(shù)，大大簡化了濾波器的設計工程。%通帶最大衰減 rs=35。209。)。figure(2)。)。wp1=wp*ts。[Bap,Aap]=zp2tf(Z,P,K)。頻率39。)。對此，采取兩項措施：一是將h（n）截短；二是將其右移。w1=。grid on。title(39。p0163。h=(1).^n.*sin((i).*Wc)./((i).*pi)。title(39。h39。H振幅39。對這種語言環(huán)境也有了新的了解。通常這些信號處理的過程要滿足實時且快速高效的要求，隨著DSP技術的發(fā)展，以DSP為內核的設備越來越多，為語音信號的處理提供了良好的平臺。 DSP系統(tǒng)特點和設計基本原則DSP（digital signal processor）是一種獨特的微處理器，是以數(shù)字信號來處理大量信息的器件。另一方面，為了語音識別實用化的需要，講者自適應、聽覺模型、快速搜索識別算法以及進一步的語言模型的研究等課題倍受關注。簡單的講，可以在原聲音流中疊加延遲一段時間后的聲流，實現(xiàn)回聲效果。當DSP完成對音頻數(shù)據(jù)的處理以后, AIC23再把數(shù)字信號轉化為模擬信號, 這樣就能夠在立體聲輸出端或者耳機輸出端聽到聲音。DX1和DR1分別與A IC23 的D IN 和DOUT 相連, 從而完成VC5509與AIC23間的數(shù)字信號通信。從而實現(xiàn)良好的模擬語音信號輸入與模擬信號的輸出。簡單的講，可以在原聲音流中疊加延遲一段時間后的聲流，實現(xiàn)回聲效果。[7]江濤,[J].電子技術用,2002,28(7):70~72[8]TexasInstrumentsIncorporated:TMS320VC5402Datasheet,2001致謝在本次課程設計的即將完成之際，筆者的心情無法平靜，本文的完成既是筆者孜孜不倦努力的結果，更是指導老師樊洪斌老師親切關懷和悉心指導的結果。（3）掌握簡單的繪圖命令。15n=0:1:15。stem(x2)c)x(n)=3cos(+)+2sin(+)0163。plot(t,x, 39。cos39。x1=[1 zeros(1,N1)]。(np)q239。0238。0163。0238。n=0:15。title(39。title(39。stem(n,xb,39。stem(abs(Xb),39。(3)熟悉巴特沃思濾波器、切比雪夫濾波器和橢圓濾波器的頻率特性。[num1,den1]=impinvar(B,A,fs)。OmegaC=2*pi*fc。Xb(k)=FFT[xb(n)]的波形 39。xb波形39。f=。subplot(2,3,4)。subplot(2,3,1)。n163。7 其他d)反三角波序列236。n239。20163。(2)應用FFT對典型信號進行頻譜分析。n50y(n)+(n1)(n2)=x(n)2x(n1)實驗過程與結果（含實驗程序、運行的數(shù)據(jù)結果和圖形）； clear all。t39。10sb)x(t)=cos(100pt)sin(pt)0163。n163。輸入A=[1 2 3 4]，B=[3,4,5,6]，求C=A+B，D=AB，E=A.*B,F=A./B,G=A.^B。另外，還要感謝這段時間來陪我一起努力同學，感謝我們這個小團隊，感謝每一個在學習和生活中所有給予我關心、支持和幫助的老師和同學們，幾年來我們一起學習、一起玩耍，共同度過了太多的美好時光。同時，在本次實驗中我了解了ICETEK – VC5509 – A板上語音codec芯片TLV320AIC23的設計和程序控制原理，并進一步掌握了數(shù)字回聲產(chǎn)生原理、編程及其參數(shù)選擇、控制，以及了解了VC5509DSP擴展存儲器的編程使用方法。這時可從耳機中聽到帶數(shù)字回聲道語音放送。模擬信號的輸入輸出電路如圖所示。MODE接數(shù)字地, 表示利用I2 C控制接口對AIC23傳輸控制數(shù)據(jù)。利用VC5509 的片上外設I2C(InterIntegrated Circuit, 內部集成電路)模塊配置AIC23 的內部寄存器。聲音遇到較遠物體產(chǎn)生的反射會比遇到較近的反射波晚些到達聲源位置，所以回聲和原聲的延遲隨反射物體的距離大小改變。它的強大數(shù)據(jù)處理能力和高運行速度，是最值得稱道的兩大特色。語音通信的研究涉及到人工智能、數(shù)字信號處理、微型計算機技術、語言聲學、語言學等許多領域，所以說語音的通信是一個多學科的綜合研究領域，其研究成果具有重要的學術價值。FIR數(shù)字濾波器的設計出來的是兩種窗的圖形，通過兩種窗的比較，我了解了他們各自的特點，幅頻和相頻特性。矩形窗設計振幅譜39。xlabel(39。xlabel(39。%矩形窗函數(shù)設計的頻響 magH=abs(H)。Wc=*pi。用矩形窗設計一個線性相位高通濾波器。w2=。freqz(b,1,512)。（3）根據(jù)過渡帶寬度和阻帶最小衰減，確定窗函數(shù)類型和窗口長度N；（4）求FIR濾波器單位脈沖響應；（5）分析幅頻特性，若不滿足要求，可適當改變窗函數(shù)形式或長度N，重復上述設計過程，以得到滿意的結果。實驗七 FIR數(shù)字濾波器的設計一、實驗目的掌握用窗函數(shù)設計FIR數(shù)字濾波的原理及其設計步驟；熟悉線性相位數(shù)字濾波器的特性。plot(W/pi,20*log10(abs(H)))。subplot(2,1,1)。[N,Wn]=buttord(wp2,ws2,rp,rs,39。rp=1。頻率/Hz39。[H,W]=zp2tf(Z,P,K)。179。201?；赾hebyshev I型模擬濾波器原型使用沖激不變轉換方法設計數(shù)字濾波器，要求參數(shù)為通帶截止頻率Wp=；通帶最大衰減Ap=1dB；阻帶截止頻率Ws=；阻帶最小衰減As=35dB。一般來說，在要求時域沖激響應能模仿模擬濾波器的場合，一般使用沖激響應不變法。)。xlabel(39。title(39。angx=angle(x).*180/pi。（一）程序b=[,]。離散系統(tǒng)中z變化建立了時域函數(shù)h（n）與z域函數(shù)H（z）之間的對應關系。)。k=0:1:N1。subplot(2,1,2)。axis([0 11 ])。frequency in pi unit39。title(39。plot(w/pi,realX)。ylabel(39。)。xlabel(39。w=[0:500]*pi/500 H=freqz(x,y,w)。在本實驗中，將學習利用MATLAB計算離散時間信號的DTFT和DFT，并加深對其相互關系的理解。stem(n,x)。title(39。grid on。%學號01 x=exp((2+j*w)*n)。title(39。%學號01 phi=0。stem(n,x)。學號后兩位248。單位階躍序列39。分析:在上圖的基礎上平移了1個單位.(3)N=6。M=1。n=(N1):N1。grid on。(2)N=6。程序:(1)N=4。n=0：N1 x=a.^n stem(n,x)單邊指數(shù)序列n的取值范圍為n=0。（3）矩形序列矩形序列有一個重要的參數(shù)，就是序列的寬度N。而要研究離散時間信號，首先需要產(chǎn)生出各種離散時間信號。在設計的過程中遇到了很多問題，剛剛開始時曾天真的認為只要把以前的程序改了參數(shù)就可以用了，可是問題沒有我想象中的那么簡單，單純的搬程序是不能解決問題的。，把高頻濾除，聲音變得比較低沉。%顯示標表示的意義 grid。低通濾波器濾波后信號的頻譜圖39。)。%butter函數(shù)（num為分子系數(shù)den為分母系數(shù)）[h,w]=freqz(num,den)。grid。疊加信號的頻譜圖39。ylabel(39。LineWidth39。%將兩個聲音疊加：wavplay(yy,fs)。fs=8000。%快速傅立葉變換subplot(3,2,2)。%信號的時域波形 title(39。%對語音信號進行采樣wavplay(x1,fs)。pause。double39。%錄音時間 fprintf(39。頻率w/pi39。%顯示縱坐標的表示意義 grid。%頻帶寬度 subplot(2,2,3)。%加網(wǎng)格[y,fs,bits]=wavread(39。)。幅值/A39。f=[(n11)*(2*pi/N)]/pi。)。)。%暫停命令 fprintf(39。%sound(x, Fs, bits)。%顯示標表示的意義 title(39。%在四行兩列的第八個窗口顯示圖形plot(f, abs(fftwave(n1)))。頻率w/pi39。ys1= fftfilt(h,x)。wn1=(ws1+wc1)/2。%加網(wǎng)格wavplay(ys,8000)。)。)。ylabel(39。,kaiser(n+2,beta))。beta=。幅值/A39。)。%信號的時域波形title(39。band39。%加網(wǎng)格3．1低通濾波器濾波 fs=8000。)。%加網(wǎng)格y0=fft(x)。plot(x)。)。頻率（HZ）39。[Bz,Az]=bilinear(B,A,Fs1)。wap=2*tan(wp/2)。fs=[1000,3200]。頻率（HZ）39。[Bz,Az]=bilinear(B,A,Fs1)。10[N,wc]=buttord(wap,was,rp,rs,39。rp=1。ylabel(39。[h,w]=freqz(Bz,Az,512,Fs1*22050)。[N,wc]=buttord(wap,was,rp,rs,39。Fs=22050。)。C:39。)。%原始信號的頻譜圖 title(39。%加標題 ylabel(39。%采樣間隔 n1=1:N/2。)。x=wavrecord(time*Fs,Fs,39。%抽樣頻率 time=3。a239。?Xa(jW)=229。dT(t)=頻域描述：利用傅氏變換的性質，時域相乘頻域卷積，若n=165。 DF線性相位的概念，即線性相位對h(n)、H(w)及零點的約束，了解四種FIR DF的頻響特點。近年來，國外及國內有些學校對一般電類專業(yè)該課程的教學主要強調應用性學習，主要介紹數(shù)字信號處理的用途和用法，而對其深奧的理論推導僅做一般介紹，并給學生提供進行實驗的機會，以激發(fā)學生對該課程的興趣和學習主動性。目前幾乎所有的高等院校都在電子工程類、信息工程類、通信工程類、電子技術類、自動控制類、電氣工程類、機電工程類、計算機科學類等工科電類及其他相關專業(yè)的本科生中開設了該門課程。對于我們通信專業(yè)，我覺得是個很好的專業(yè)，現(xiàn)在這個專業(yè)很熱門，這個專業(yè)以后就業(yè)的方向也很多，就業(yè)面很廣。這學期數(shù)字信號處理所含有的具體內容如下：第一單元的課程我們深刻理解到時域離散信號和時域離散系統(tǒng)性質和特點；時域離散信號和時域離散系統(tǒng)時域分析方法；模擬信號的數(shù)字處理方法。第二單元的課程我們理解了時域離散信號（序列）的傅立葉變換，時域離散信號Z變換，時域離散系統(tǒng)的頻域分析。我們畢業(yè)以后工作，可以進入設備制造商、運營商、專有服務提供商以及銀行等領域工作。隨著計算機技術、微電子技術、數(shù)字信號處理理論和方法的發(fā)展，半個世紀以來，尤其是最近的三十來年里，數(shù)字信號處理的方法和應用得到了飛躍式的發(fā)展，數(shù)字信號處理的地位和作用變得越來越重要。對該課程的改革思想主要是課程內容要適應數(shù)字信號處理技術的發(fā)展現(xiàn)狀，淡化枯燥的數(shù)學推導，輔助以現(xiàn)代化教學手段，并開設相應的實驗課。，F(xiàn)IR濾波器利用函數(shù)fftfilt對信號進行濾波。229。Xa(jWjk)=229。238。%錄音時間 fprintf(39。double39。pause。f=[(n11)*(2*pi/N)]/pi。幅值/A39。原始信號的頻譜圖39。%顯示縱坐標表示的意義 title(39。, 7[4000,24000])。ylabel(39。rp=1。s39。plot(w,abs(h))。耗損（dB）39。rs=100。s39。figure(1)。)。Fs=8000。was=2*tan(ws/2)。figure(4)。)。%播放原始信號wavplay(x,fs)。%錄制信號的時域波形title(39。%快速傅立葉變換 subplot(4,2,3)。%顯示橫坐標表示的意義 ylabel(39。beta=。,kaiser(n+1,beta))。低通濾波后信號的時域波形39。%顯示標表示的意義 grid。)。ws=2*5000/fs。[h1,w]=freqz(h)。幅值/A39。%顯示