山东大学第七次数字信号处理实验报告.docxVIP

数字信号处理实验报告实验七、综合练习题已知信号。现对该信号进行频谱分析，请用Matlab编程实现：1、选择合适的抽样频率和记录时间，对进行抽样，要求抽样得到的信号是周期的，且能分辨上述三个信号。2、画出上述信号的频谱。3、设计两个巴特沃斯滤波器，将分离为,其中获得, 获得。4、用升余弦窗设计两个滤波器，完成3的任务。注：所有结果画图表示。实验内容及实验结果：1、选择合适的抽样频率和记录时间，对进行抽样，要求抽样得到的信号是周期的，且能分辨上述三个信号。2、画出上述信号的频谱。以上两问的实验程序：clc;clear all;close all;t=0:0.00001:0.4;x=cos(2*pi*100*t)+0.5*cos(2*pi*105*t)+0.75*cos(2*pi*280*t); %原信号subplot(2,2,1);plot(t,x);grid on;title(原信号x(t));t3=1/1000;L=50; n1=L/t3;D=2*pi/(n1*t3);n2=0:n1-1;x2=cos(2*pi*100*t3*n2)+0.5*cos(2*pi*105*t3*n2)+0.75*cos(2*pi*280*t3*n2);k1=floor(-(n1-1)/2:(n1-1)/2);Xw1=t3*fftshift(fft(x2));subplot(2,2,2);stem(k1*D/(2*pi),abs(Xw1));grid on;title(原信号频谱X(f));t10=1/1000;t2=0:t10:0.4;x3=cos(2*pi*100*t2)+0.5*cos(2*pi*105*t2)+0.75*cos(2*pi*280*t2); %抽样后的信号subplot(2,2,3);stem(t2,x3,b);grid on;title(抽样信号);hold on;plot(t,x)L10=50; n10=L10/t10;D10=2*pi/(n10*t10);n11=0:n10-1;x10=cos(2*pi*100*n11*t10)+0.5*cos(2*pi*105*n11*t10)+0.75*cos(2*pi*280*n11*t10);k2=floor(-(n10-1)/2:(n10-1)/2);Xw11=t10*fftshift(fft(x10));subplot(2,2,4);stem(k2*D10,abs(Xw11));grid on;title(抽样信号频谱)以上两问的实验结果：3、设计两个巴特沃斯滤波器，将分离为,其中获得, 获得。实验程序：clc;clear all;close all;load Butterworth1 %巴特沃斯低通滤波器load Butterworth2 %巴特沃斯带通滤波器fs=1000;N1=128;n3=0:N1-1;t3=n3/fs;xt=cos(2*pi*100*t3)+0.5*cos(2*pi*105*t3)+0.75*cos(2*pi*280*t3);[b3,a3]=sos2tf(SOS,G);y=filter(b3,a3,xt);[Yw,w]=freqz(y,1,1024,fs);[b4,a4]=sos2tf(SOS1,G1);y1=filter(b4,a4,xt);[Yw1,w1]=freqz(y1,1,1024,fs);subplot(221)[Y3,w3]=freqz(b3,a3,1024);plot(w3,20*log10(abs(Y3)),b,linewidth,1);title(巴特沃斯低通滤波幅度响应(dB));ylabel(20log|H(e^jw)|(dB));grid on;subplot(222);stem(w,abs(Yw),b,.);grid;hold onstem(-w,abs(Yw),b,.);grid on;xlim([-600,600]);subplot(223)[Y4,w4]=freqz(b4,a4,1024);plot(w4,20*log10(abs(Y4)),b,linewidth,1);title(巴特沃斯带通滤波幅度响应(dB));ylabel(20log|H(e^jw)|(dB));grid on;subplot(224);stem(w1,abs(Yw1),b,.);grid;hold onstem(-w1,abs(Yw1),b,.);grid on;xlim([-600,600]);滤波器设计：巴特沃斯低通滤波器：巴特沃斯带通滤波器：实验结果：4、用升余弦窗设计两个滤波器，完成3的任务。实验程序：clc;clear all;close all;load Hamming1 %海明窗低通滤波器load Hamming2%海明