【正文】 % 求所給系統(tǒng)的階躍響應(yīng)；plot(t,y,“k”)。掌握并熟練使用有關(guān)離散系統(tǒng)分析的 MATLAB 調(diào)用函數(shù)及格式，以深入理解離散時間系統(tǒng)的頻率特性。= =234。 229。通過本次實驗，我對使用MATLAB快速設(shè)計濾波器的流程更為熟練，同時，也由衷地感嘆MATLAB的功能強(qiáng)大之處，它讓我們在設(shè)計時能節(jié)約大量的時間。title(39。plot(w,20*log10(abs(hk(k))))。h(n)39。subplot(1,2,1)。ws=pi*。)title(39。,m,c,39。% 凱塞窗 m=1:200。plot(m,d)subplot(4,2,5)m=200。漢寧窗39。b=bartlett(m)%三角窗 m=1:200。例如，要設(shè)計線性相位低通特性，可選擇h(n)=h(N1n)一類，而不能選h(n)=h(N1n)一類。五、思考題用雙線性變換設(shè)計數(shù)字濾波器的過程中,下面變換公式的T值的取值，對設(shè)計結(jié)果是否有影響？為什么？21z1s= T1+z1雖然采用雙線性變換法設(shè)計數(shù)字濾波器不會產(chǎn)生混疊現(xiàn)象，T得取值可以任選。)。Xk=fft(x,N)。)。mag=abs(H)。y(n)39。subplot(2,2,2)。b2=[A,2*A,A]。ylabel(39。.39。axis([0,1,100,5])。數(shù)字濾波器幅值響應(yīng)曲線39。[Nd,wdc]=buttord(wpz,wsz,rp,rs)。)。wk=2*pi*fk。%預(yù)畸變校正轉(zhuǎn)換指標(biāo) [N,wc]=buttord(wp,ws,rp,rs,39。以下是通過matlab實現(xiàn)的巴特沃斯低通IIR濾波器的設(shè)計 %1 T=1。在實驗過程中，學(xué)習(xí)如何使用matlab程序語言解決問題，是很有價值的。q=n*2*pi/N。q=n*2*pi/N。q=n*2*pi/N。q=n*2*pi/N。)。y1=[xn(1:1:M),zeros(1,NM)]。x(n)39。n1=[0:1:N1]。)。stem(wk,abs(Xk),39。結(jié)果分析：雖然周期序列不存在FT,但是一個周期序列可以利用其DFS系數(shù)X(k)表示它的頻譜分布規(guī)律，從以上各頻譜圖可以看出，隨著補(bǔ)零點數(shù)的增加，周期序列的諧波次數(shù)越來越多，其頻譜的包絡(luò)線越來越平滑連續(xù)，更能反映幅度值隨時間的變化。xlabel(39。xtide(n)39。magXk=abs([Xk(N/2+1:N)Xk(1:N/2+1)])。Xk=xn*WNnk。N=55。第三篇：數(shù)字信號處理實驗報告數(shù)字信號處理 實驗報告實驗一序列的傅立葉變換一、實驗?zāi)康?DFT算法的原理。知道輸入、輸出語句怎么形成。n=0:1:N1。M=(N1)/2。n=0:1:N1。M=(N1)/2。n=0:1:N1。M=(N1)/2。subplot(2,1,1)。(選做)(2)fc=,d=1dB,fr=,At=25dB,T=1ms；分別用脈沖響應(yīng)不變法及雙線性變換法設(shè)計一巴特沃思數(shù)字低通濾波器，觀察所設(shè)計數(shù)字濾波器的幅頻特性曲線，記錄帶寬和衰減量，檢查是否滿足要求。f = w/pi*fs/2。[h1,w]=freqz(num1,den1)。plot(f,abs(h2))。fs=4000。[B,A]=butter(3, OmegaC,39。)。)。a=。stem(abs(Xa),39。stem(n,xa,39。7 其他(1)觀察高斯序列的時域和幅頻特性，固定信號xa(n)中參數(shù)p=8，改變q的值，使q分別等于2，4，8，觀察它們的時域和幅頻特性，了解當(dāng)q取不同值時，對信號序列的時域幅頻特性的影響；固定q=8，改變p，使p分別等于8，13，14，觀察參數(shù)p變化對信號序列的時域及幅頻特性的影響，觀察p等于多少時，會發(fā)生明顯的泄漏現(xiàn)象，混疊是否也隨之出現(xiàn)？記錄實驗中觀察到的現(xiàn)象，繪出相應(yīng)的時域序列和幅頻特性曲線。4n239。xc(n)=237。n163。(3)了解應(yīng)用FFT進(jìn)行信號頻譜分析過程中可能出現(xiàn)的問題，以便在實際中正確應(yīng)用FFT。N=50。),ylabel(39。t163。15n=0:1:15。并用stem語句畫出A、B、C、D、E、F、G。FIR數(shù)字濾波器的設(shè)計出來的是兩種窗的圖形，通過兩種窗的比較，我了解了他們各自的特點，幅頻和相頻特性。矩形窗設(shè)計振幅譜39。xlabel(39。xlabel(39。%矩形窗函數(shù)設(shè)計的頻響 magH=abs(H)。Wc=*pi。用矩形窗設(shè)計一個線性相位高通濾波器。w2=。freqz(b,1,512)。（3）根據(jù)過渡帶寬度和阻帶最小衰減，確定窗函數(shù)類型和窗口長度N；（4）求FIR濾波器單位脈沖響應(yīng)；（5）分析幅頻特性，若不滿足要求，可適當(dāng)改變窗函數(shù)形式或長度N，重復(fù)上述設(shè)計過程，以得到滿意的結(jié)果。實驗七 FIR數(shù)字濾波器的設(shè)計一、實驗?zāi)康恼莆沼么昂瘮?shù)設(shè)計FIR數(shù)字濾波的原理及其設(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型模擬濾波器原型使用沖激不變轉(zhuǎn)換方法設(shè)計數(shù)字濾波器，要求參數(shù)為通帶截止頻率Wp=；通帶最大衰減Ap=1dB；阻帶截止頻率Ws=；阻帶最小衰減As=35dB。一般來說，在要求時域沖激響應(yīng)能模仿模擬濾波器的場合，一般使用沖激響應(yīng)不變法。)。xlabel(39。title(39。angx=angle(x).*180/pi。（一）程序b=[,]。離散系統(tǒng)中z變化建立了時域函數(shù)h（n）與z域函數(shù)H（z）之間的對應(yīng)關(guān)系。)。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)。在本實驗中，將學(xué)習(xí)利用MATLAB計算離散時間信號的DTFT和DFT，并加深對其相互關(guān)系的理解。stem(n,x)。title(39。grid on。%學(xué)號01 x=exp((2+j*w)*n)。title(39。%學(xué)號01 phi=0。stem(n,x)。學(xué)號后兩位248。單位階躍序列39。分析:在上圖的基礎(chǔ)上平移了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)生出各種離散時間信號。常見的離散信號如下：（1）單位沖激序列δ（n）如果δ（n）在時間軸上延遲了k個單位，得到δ（nk），即長度為N的單位沖激序列δ（n）可以通過下面的MATLAB命令獲得?？梢宰C明，只有當(dāng)2π/w為有理數(shù)時，正弦序列具有周期性?！?并繪出其圖及d(n”學(xué)號后兩位0n185。單位沖激序列39。title(39。、u(n“學(xué)號后兩位”)及0np0238。)。單位階躍序列39。z=xy。編制程序產(chǎn)生正弦序列x(n)=cos(2pn)、x(n)=cos231。phi=0。分析:該序列具有周期性,且輸出為余弦信號.(2)N=5。余弦信號39。n=0::N1。N=3。subplot(2,1,2)stem(n,imag(x)),title(39。a=1。n=0:N1。實驗三 離散時間信號的頻域分析一、實驗?zāi)康男盘柕念l域分析是信號處理中一種有效的工具。３－１clear。subplot(221)。)。xlabel(39。)。ylabel(39。plot(w/pi,imagX)。title(39。Xk=fft(xn,N)。figure N=64。stem(k,abs(Xk))。它的作用是將H（z）的有理分式表示式轉(zhuǎn)換為零極點增益形式。一種形式為[H,w]= reqz(B,A,N)其中，B與A分別表示H（z）分子和分母多項式的系數(shù)向量；N為正整數(shù)，默認(rèn)值為512；返回值w包含[0，π]范圍內(nèi)的N個頻率等分點；返回值H則是離散時間系統(tǒng)頻率響應(yīng)在0~π范圍內(nèi)N個頻率處的值。x22=x2+(x2==0)*eps。)。)。)。IIR濾波器的階數(shù)就等于所選的模擬原型濾波器的階數(shù)，所以其階數(shù)確定主要是在模擬原型濾波器中進(jìn)行的。（1）butter。%阻帶最小衰減Fs=1000。249。[Z,P,K]=cheby1(N,rp,Wn,39。plot(Q*Fs/(2*pi),abs(P))?；贐utterworth型模擬濾波器原型使用雙線性變換方法設(shè)計數(shù)字濾波器的，要求參數(shù)為截止頻率Wp=；通帶最大衰減Ap=1dB；阻帶截止頻率Ws=；阻帶最小衰減AS=40dB。ws1=ws*ts。[b,a]=lp2lp(Bap,Aap,Wn)。)ylabel(39。ylabel(39。設(shè)計時，要根據(jù)阻帶的最小衰減和過渡帶寬度來選擇恰當(dāng)?shù)拇昂瘮?shù)類型和窗口長度N。w2=。程序2:萊克曼窗設(shè)計N=48。布萊克曼窗39。w163。%矩形窗函數(shù)設(shè)計的系統(tǒng)脈沖響應(yīng) w=(0:1:500)*2*pi/500。矩形窗設(shè)計h(n)39。)。)。在實驗的過程中，我對數(shù)字濾波器的整個過程有了很好的理解和掌握。（5）觀察離散系統(tǒng)的頻率響應(yīng)。a=(+3*i)*n。15(4)繪出下列時間函數(shù)的圖形，對x軸、y軸以及圖形上方均須加上適當(dāng)?shù)臉?biāo)注： a)x(t)=sin(2pt)0163。)。2z1+z2)/(+)，求出并繪制H(z)的幅(7)計算序列{821 2 3}和序列{2 313}的離散卷積，并作圖表示卷積結(jié)果。h=filter(a,b,x1)。e239。n163。34163。n163。q=2。)。)。)。)。脈沖響應(yīng)不變法： fc=1000。f = w/pi*fs/2。[num2,den2]=bilinear(B,A,fs)。s39。)。plot(f,abs(h2),39。N=20。a)在上面所得的幾幅圖中，在截止頻率兩邊可以觀察到幅頻響應(yīng)的擺動行為。win=boxcar(N)。plot(w,angle(H))。win=boxcar(N)。plot(w,angle(H))。win=boxcar(N)。plot(w,angle(H))。我在這次實驗中的收獲很大，接觸了很多新