數(shù)字信號處理實驗三

1Z2、掌握并熟練使用有關(guān)離散系統(tǒng)分析的調(diào)用函數(shù)及格式，以深入理解離散時間LSIH(z)H(z)=B(z)/A(z),aH(z)a

b(1)b(2)z1b(3)z2...

bb

1a(2)z1a(3)z2...

1)z中nb和na分別是H(Z)分子與分母多項式的階次，在有關(guān)的系統(tǒng)分析的文件a[a(1),a(2),...,a(na例1 0.0018360.007344z10.011016z20.007374z3H(z)

（y（nx=ones(100)x(n)=1,n=1~100；t=1:100;%t用于后面的繪圖；b=[0.001836,0.007344,0.011016,0.007374,0.001836];%形成向量b；a=[1,-3.0544,3.8291,-2.2925,.55075];%a；%求所給系統(tǒng)的輸出,plot(t,x,'r.',t,y,'k-');gridon;%x(n)（綠色）y(n)（黑色）% h=impz(b，a，N)Nh(n)n10開始。2h(n)。 %stem(t,h,'.');grid3H(ejwNN2的整次冪；w是返回頻率軸坐標向量，供繪圖用；0~Fs/2。31H(ejw。clear %ylabel('AmplitudeFreq.Res.')plot(w,Hphase);gridon;ylabel('PhaseFreq.Res.')調(diào)用格式：z 或zzp41FIRH(z)11.7z11.53z20.648z3的極零圖。zne(b,a);%求并畫出所給系統(tǒng)的極零圖，該系統(tǒng)為IIR系統(tǒng)b=[1-1.71.53-0.68];zne(b,a);%求并畫出第二個系統(tǒng)的極零圖，該系統(tǒng)為FIR系H(z)0.2z20.8，編程并繪出系統(tǒng)的單位階躍響應(yīng)y(n),頻率響應(yīng)H(ejwx=ones(100)x(n)=1,n=1~100；t=1:100;%t用于后面的繪圖；b=[0,0,-0.2];%形成向量b；a=[1,0,0.8];%%plot(t,x,'r.',t,y,'k-');gridon;%x(n)（綠色）y(n)（黑色）% %ylabel('AmplitudeFreq.Res.')plot(w,Hphase);gridon;ylabel('PhaseFreq.Res.') 實驗四FFT DFT的基本性質(zhì)。FFTFFTFFT對連續(xù)信號和時域離散信號進行譜分析的方法，了解可能出現(xiàn)的分析誤差及 DFTDFT復(fù)習(xí)FFT算法原理與編程思想。熟悉FFT算法的實現(xiàn) DFTfft格式（1：y=fft(x) 計算x的FFT變換y。當x為矩陣，計算x中每一列信號的離散傅氏變換。當x的長度為2的冪時，采用基2算法，否則采用基算法。格式（2：y=fft(x,n) 計算x的n點FFT，當x長度大于n時，截斷x，否則補零。Plot線性繪圖函數(shù)。stem：繪制離散序列圖。subplot：多坐標設(shè)置與定位當前坐標系。n 0n 4n0 04 0n 4n （4）xn 4（5）

4nsin8三.實驗內(nèi)容6

N8,16x6

N8,16（參考）%fftn1=0:3;x1=[1,1,1,1];subplot(221);stem(n1,x1);title('x1序列');k1=0:7;y11=fft(x1,8);magy11=[abs(y11)];subplot(222);stem(k1,magy11);title('x1的8點k2=0:15;y12=fft(x1,16);magy12=[abs(y12)];subplot(224);stem(k2,magy12);title('x116n2=0:7;x2=[1,2,3,4,4,3,2,1];subplot(221);stem(n2,x2);title('x2序列');k1=0:7;y21=fft(x2,8);magy21=abs(y21);subplot(222);stem(k1,magy21);title('x2的8點FFT');k2=0:15;y22=fft(x2,16);magy22=abs(y22);subplot(224);stem(k2,magy22);title('x2的16點n3=0:7;x3=[4,3,2,1,1,2,3,4];subplot(221);stem(n3,x3);title('x3序列');k1=0:7;y31=fft(x3,8);magy31=abs(y31);subplot(222);stem(k1,magy31);title('x3的8點FFT');k2=0:15;y32=fft(x3,16);magy32=abs(y32);subplot(224);stem(k2,magy32);title('x3的16點n41=0:7;x41=cos((pi/4)*n41);subplot(221);stem(n41,x41);title('x48點序列');k1=0:7;y41=fft(x41,8);magy41=abs(y41);subplot(222);stem(k1,magy41);title('x4的8點n42=0:15;x42=cos((pi/4)*n42);subplot(223);stem(n42,x42);title('x416點序列');k2=0:15;y42=fft(x42,16);magy42=abs(y42);subplot(224);stem(k2,magy42);title('x416點n51=0:7;x51=sin((pi/8)*n51);subplot(221);stem(n51,x51);title('x58點序列');k1=0:7;y51=fft(x51,8);magy51=abs(y51);subplot(222);stem(k1,magy51);title('x5的8點n52=0:15;x52=sin((pi/8)*n52);subplot(223);stem(n52,x52);title('x516點序列');k2=0:15;y52=fft(x52,16);magy52=abs(y52);subplot(224);stem(k2,magy52);title('x516點

x6=cos(f1*2*pi*t)+cos(f2*2*pi*t)+cos(f3*2*pi*t);subplot(221);plot(t,x6);title('x6%x618y61=fft(x6,8);magy61=abs(y61);subplot(222);stem(k61,magy61);title('x68%x6216y62=fft(x6,16);magy62=abs(y62);subplot(224);stem(k62,magy62);title('x616x7nx4nx5nN16FFTX7k

，根據(jù)DFT

figure(7);;x42=cos((pi/4)*n42);x52=sin((pi/8)*n52);n7=0:15;x7=x42+x52;subplot(221);stem(n7,x7);title('x716點序列');k7=n7;y7=fft(x7,16);magy7=abs(y7);subplot(222);stem(k7,magy7);title('x716FFT');X4=real(y7);X5=imag(y7)*i;subplot(223);stem(k7,abs(X4));title('x416點序列');subplot(224);stem(k7,abs(X5));title('x516點序列');x8nx4njx5nN16FFTX8k。當已求出X8k，根據(jù)DFT

n8=0:15;x8=x42+x52*j;subplot(221);stem(n8,abs(x8));title('x816點序列');k8=n8;y8=fft(x8,16);magy8=abs(y8);subplot(222);stem(k8,magy8);title('x816FFT'