数字信号处理实验..

（2）系统单位脉冲响应序列产生子程序。本实验要用到两种系统。 a、ha(t)=u(t)-u(t-1)产生子程序： close all; clc; clear; n=-20:20;

y1=[zeros(1,20),ones(1,21)]; y2=[zeros(1,21),ones(1,20)]; subplot(311) stem(n,y1,'fill'); xlabel('n'); ylabel('u(t)'); subplot(312) stem(n,y2,'fill'); xlabel('n');

ylabel('u(t-1)'); subplot(313); y=y1-y2; stem(n,y,'fill'); xlabel('n'); ylabel('u(t)-u(t-1)');

b、ha(t)=e-2tu(t)子程序设计： close all; clear; clc; t=0:0.01:2; a=exp(-2*t); stem(t,a); axis([0,1,0,1.5]); xlabel('t');

ylabel('exp(-2t)u(t)'); hold on;

（3）有限长序列线性卷积子程序，可以直接调用MATLAB语言中卷积函数conv。 close all; clc; clear; a=[1 0 1 2 3]; b=[0 1 3 2 1]; c=conv(a,b); m=length(c)-1; n=0:1:m; stem(n,c,'fill'); title('有限序列卷积');

grid on;

3.通过编程实现方波信号输入系统，观测输出信号，与实际卷积计算结果对比，验证实验结果的正确性。 验证方波卷积的程序如下： close all; clc; clear; Ts=0.2; t=0:Ts:2; a=rectpuls(t,2); b=conv(a,a)*Ts; m=(length(b)-1)*Ts; n=0:Ts:m; t1=0:0.01:2-0.01;