DSP实验考试题目及参考程序.doc

1.计算x1=1,0,7,3与x2=5,7,8,0,1的序列和x。画出x1,x2,x的图形。（提示：x1,x2具有相同长度，位置对应才能相加）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. 零、极点和幅度频率响应。clear allw=pi*freqspace(100);%确定二维频率响应的频率空间，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);%幅度响应subplot(3,1,3); stem(w,angle(h);% 相频响应3.零、极点分布图和幅度频率响应。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) ;幅度响应h=freqz(a,b,w);subplot(2,1,2);stem(w,abs(h); 相频响应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；用脉冲响应不变法设计一Butterworth数字低通滤波器。N,wn=buttord(wp,wr,Ap,Ar,s) 计算巴特沃斯模拟滤波器的阶数N和3dB截止频率wn。B,A = BUTTER(N,Wn,s)设计一个阶数为N，频率为Wn的低通模拟滤波器。num,den=impinvar(B,A,fs)参数b、a为模拟滤波器传递函数的分子系数和分母系数，fs是采样频率，默认1Hz。num和den返回设计出的数字滤波器的传递函数的分子和分母系数。%脉冲响应不变法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)进行512点傅里叶变换，设置频坐标，分别观察时/频图形。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