DSP實驗考試題目及參考程序.doc

1.計算x1=1,0,7,3與x2=5,7,8,0,1的序列和x。畫出x1,x2,x的圖形。（提示：x1,x2具有相同長度，位置對應(yīng)才能相加）n=0:4;x1=1 0 7 3 0;stem(n,x1)n=0:4;x2=5 7 8 0 1;stem(n,x2) n=0:4;x=x1+x2;stem(n,x); 2. 零、極點和幅度頻率響應(yīng)。clear allw=pi*freqspace(100);%確定二維頻率響應(yīng)的頻率空間，w=0:1:50*pi/50;a=0.0528 0 0.1295 0.1295 0.0797 ;b=1 -1.8107 2.4947 -1.8801 0.9537 0.2336; z,p,k=tf2zp(a,b);subplot(2,1,1);%子圖的位置：2行1列兩個圖形中的第1個位置zplane(z,p);%為零點，為極點h=freqz(a,b,w);subplot(2,1,2); stem(w,abs(h);%幅度響應(yīng)subplot(3,1,3); stem(w,angle(h);% 相頻響應(yīng)3.零、極點分布圖和幅度頻率響應(yīng)。w=pi*freqspace(100);a=2 4 6;b=1 -0.8 0.4; z,p,k=tf2zp(a,b);subplot(2,1,1); zplane(z,p) ;幅度響應(yīng)h=freqz(a,b,w);subplot(2,1,2);stem(w,abs(h); 相頻響應(yīng)subplot(3,1,3);stem(w,angle(h); 4. 分別計算16點序列 的16點和32點DFT，繪出幅度譜圖形，并繪出該序列的DTFT圖形。clear all;n1=0:15;x1=cos(5*pi*n1/16);subplot(4,1,1);plot(n1,x1);u1=fft(x1);subplot(4,1,2);stem(n1,u1);n2=0:31;x2=cos(5*pi*n2/16);subplot(4,1,3);plot(n2,x2);u2=fft(x2);subplot(4,1,4);stem(n2,u2);5. fp=0.2KHz,Ap=1dB, fr=0.3KHz,Ar=25dB,T=1ms；用脈沖響應(yīng)不變法設(shè)計一Butterworth數(shù)字低通濾波器。N,wn=buttord(wp,wr,Ap,Ar,s) 計算巴特沃斯模擬濾波器的階數(shù)N和3dB截止頻率wn。B,A = BUTTER(N,Wn,s)設(shè)計一個階數(shù)為N，頻率為Wn的低通模擬濾波器。num,den=impinvar(B,A,fs)參數(shù)b、a為模擬濾波器傳遞函數(shù)的分子系數(shù)和分母系數(shù)，fs是采樣頻率，默認(rèn)1Hz。num和den返回設(shè)計出的數(shù)字濾波器的傳遞函數(shù)的分子和分母系數(shù)。%脈沖響應(yīng)不變法fp=200;fr=300;Ap=1;Ar=25;T=1e-3;fs=1/T;wp=2*pi*fp;wr=2*pi*fr;N1,wn1=buttord(wp,wr,Ap,Ar,s);B1,A1=butter(N1,wn1,s);num1,den1=impinvar(B1,A1,fs);H1,w=freqz(num1,den1);f=w/pi*fs/2;plot(f,20*log10(abs(H1),-.)Grid6. （fft在信號分析中的運用）使用頻譜分析方法從噪聲污染的信號x(t)中鑒別出有用的信號。 x(t)=sin(2pi*100t)+sin(2pi*200t)+rand(size(t);采樣信號=1000Hz ，采樣周期為0.001s. t=0:0.001:1；提示：對x(t)進(jìn)行512點傅里葉變換，設(shè)置頻坐標(biāo)，分別觀察時/頻圖形。t=0:0.001:1;x=sin(2*pi*100*t)+sin(2*pi*200*t)+rand(size(t); subplot(2,1,1);plot(t,x);xlabel(t);ylabel(Mangitude);title( xt的圖形);F_512=fft(x,1024);kk=1000/1024*(0:511); subplot