数字信号处理实验matlab.docVIP

已知信号，分别用FS=5000、1000对信号进行采样，请画出信号采样后的波形以及频谱 fs1=5000;fs2=1000;n=1024;tp=50/1000; t1=-2*pi:1/fs1:2*pi; t2=-2*pi:1/fs2:2*pi; xt1=exp(-1*abs(t1)); xw1=fft(xt1,n); xt2=exp(-1*abs(t2)); xw2=fft(xt2,n); figure(1); subplot(2,2,1);plot(t1,xt1); subplot(2,2,2); k=0:n-1;wk=2*k/tp;plot(wk,abs(xw1)); axis([0,fs1,0,1.2*max(abs(xw1))]) subplot(2,2,3);plot(t2,xt2); subplot(2,2,4); k=0:n-1;wk=2*k/tp;plot(wk,abs(xw2)); axis([0,fs2,0,1.2*max(abs(xw2))]) 2、已知x(n)=cos(0.48πn)+cos(0.52πn)，n：0-----9 （1）对信号做10点DFT，并画出信号频谱图 （2）为使信号频谱曲线光滑，该对原始信号怎样处理，请画出修改后的频谱图 （3）对原始信号做100点DFT，画出信号频谱图比较各图，分析结果 t=0:9; xn=cos(0.48*pi*t)+cos(0.52*pi*t); N=length(xn); n=0:N-1;k=n;nk=n*k; wn=exp(-1*j*2*pi/N); wnk=wn.^nk; xk=xn*wnk;%10点dft； figure(2) subplot(3,2,1);plot(t,xn); subplot(3,2,2);stem(n,abs(xk),.); %----------------------------------------------------- t=0:0.1*pi:10*pi; xn=cos(0.48*pi*t)+cos(0.52*pi*t); N=length(xn); n=0:N-1;k=n;nk=n*k; wn=exp(-1*j*2*pi/N); wnk=wn.^nk; xk=xn*wnk;subplot(3,2,3);plot(t,xn); subplot(3,2,4);stem(n,abs(xk),.); %------------------------------------------------------- t=0:99; xn=cos(0.48*pi*t)+cos(0.52*pi*t); N=length(xn); n=0:N-1;k=n;nk=n*k; wn=exp(-j*2*pi/N); wnk=wn.^nk; xk=xn*wnk;%100点dft； subplot(3,2,5);plot(t,xn); subplot(3,2,6);stem(n,abs(xk),.);设计一个通带截止频率为500 Hz -800 Hz，阻带截止频率为600 Hz -700 Hz，通带内波纹的最大衰减为0.1 dB，阻带内波纹的最小衰减为50 dB，采样频率为2000Hz的带通椭圆数字滤波器，画出其幅频响应图。 figure(3) fpl=500;fpu=800; fsl=600;fsu=700; fs=2000; wp=2*[fpl,fpu]/fs; ws=2*[fsl,fsu]/fs; rp=0.1;rs=50; [N,wc]=ellipord(wp,ws,rp,rs); [b,a]=ellip(N,rp,rs,wc); freqz(b,a);综合题： 人体心电图测量信号在测量的过程中会经常受到工业高频干扰，所以必须经过低通滤波器后，才能判断心脏功能的有用信息。下面是一组实际心电图信号的采样样本，其中存在高频干扰，请通过MATLAB程序，以为输入序列，滤除其中的额干扰成分。 X(n)={-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38, -60,-84,-90,-66,-32,-4,-2,8,12,12,10,6,6,4,0,0,0,0,0,-2,-2,0,0,-2, -2,-2,-2,0} 实验内容包括： 绘制原数据图形， 设计巴特沃斯低通滤波器并绘制其幅频响应曲线； 用设计的滤波器对原数据进行滤波 绘制滤波后的数据图 figure(4) xn1=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6]; xn2=[12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,8,12];