数字信号处理第二章实验报告..docxVIP

实 验 报 告课 程：数字信号处理专业班级：学生姓名：学 号：年月 日2.1对M=2，运行上述程序，生成输入x[n]=s1[n]+s2[n]的输出信号。输入x[n]的哪个分量被该离散时间系统抑制？% 程序 P2_1% 一个M点滑动平均滤波器的仿真% 产生输入信号n = 0:100;s1 = cos(2*pi*0.05*n); % 一个低频正弦s2 = cos(2*pi*0.47*n); % 一个高频正弦x = s1+s2;% M点滑动平均滤波器的实现M = input(滤波器所需的长度 = );num = ones(1,M);y = filter(num,1,x)/M;clf;subplot(2,2,1);plot(n, s1);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(低频正弦);subplot(2,2,2);plot(n, s2);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(高频正弦);subplot(2,2,3);plot(n, x);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(输入信号);subplot(2,2,4);plot(n, y);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(输出信号); axis;图形显示如下：答：输入部分的高频成分成分被抑制了。2.3对滤波器长度M和正弦信号s1[n]和s2[n]的频率取其他值，运行程序P2.1,算出结果。n = 0:100;s1=cos(2*pi*0.02*n);s2=cos(2*pi*0.46*n);x = s1+s2;% M点滑动平均滤波器的实现M = input(滤波器所需的长度 = );num = ones(1,M);y = filter(num,1,x)/M;clf;figure,subplot(2,2,1);plot(n, s1);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(低频正弦);subplot(2,2,2);plot(n, s2);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(高频正弦);subplot(2,2,3);plot(n, x);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(输入信号);subplot(2,2,4);plot(n, y);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(输出信号); axis;num =[1,-ones(1,M-1)];y = filter(num,1,x)/M;figure,subplot(2,2,1);plot(n, s1);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(低频正弦);subplot(2,2,2);plot(n, s2);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(高频正弦);subplot(2,2,3);plot(n, x);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(输入信号);subplot(2,2,4);plot(n, y);axis([0, 100, -2, 2]);xlabel(时间序号n); ylabel(振幅);title(输出信号); axis;图形显示如下：答：运行结果如下图，可以看出输出信号保留了输入信号x[n]的高频分量，即保留了s2[n]分量，低频部分s1[n]被抑制了。2.5用不同频率的正弦信号作为输入信号，计算每个输入信号的输出信号。输出信号是如何受到输入信号频率的影响的？从数学上对你的结论加以证明。% 程序P2_2% 产生一个正弦输入信号clf;n = 0:200;f=input(Please input the value of f:)x = cos(2*pi*f*n);% 计算输出信号x1 = [x 0 0]; % x1[n] = x[n+1] x2 = [0 x 0]; % x2[n] = x[n]x3 = [0 0 x]; % x3[n] = x[n-1]y = x2.*x2-x1.*x3;y = y(2:20