信号与系统实验三 信号的卷积计算

p=0.01;
t1=0:p:1;x1=ones(size(t1)); %size()命令求出矩阵的行数和列数，ones()生成全1阵；
t2=0:p:2;x2=ones(size(t2));
y=conv(x1,x2); %计算序列x1与x2的卷积和y
y=y*p;%缩小幅度
t0=t1(1)+t2(1); %计算序列y非零样值的起点位置
t3=length(x1)+length(x2)-2; %计算卷积和y的非零样值的宽度
t=t0:p:(t3*p+t0); %确定卷积和y非零样值的时间向量
subplot(2,2,1)
plot(t1,x1) %在子图1绘x1⑴时域波形图
title('x1')
xlabel('t')
ylabel('x1')
subplot(2,2,2)
plot(t2,x2) %在子图2绘x2 (t)时波形图
title('x2')
xlabel('t')
ylabel('x2')
subplot(2,2,3)
plot(t,y); %画卷积y(t)的时域波形
g=get(gca,'position'); %获取坐标轴的未知属性
g(3)=2.5*g(3);
set(gca,'position',g) %将第三个子图的横坐标范围扩为原来的2.5倍
title('y(t)=x1(t)*x2(t)')
xlabel('t')
ylabel('y(t)')

n1 = 0:2;%设置自变量n1范围
n2 = 0:3;%设置自变量n2范围
n3 = 0:5;%设置自变量n3范围
x = [2,1,2];%给出序列X的对应取值
h = [0,1,2,3]; %给出序列h的对应取值
y = conv( x, h ); %求序列x和h的卷积和
subplot( 3,1,1 ); stem( n1,x ); title('signal x'); %作出离散序列x的图
subplot( 3,1,2 ); stem( n2,h ); title('signal h'); %作出离散序列h的图
subplot( 3,1,3 ); stem( n3,y ); title('signal y'); %作出离散序列y的图