间序号存入向量ny中，利用stem(ny,y)画出y[n]。
代码：
图形：
21.510.50-0.5-1-1.5-2-101234567
2． 考虑两个有限长序列h[n]和x[n]用MATLAB向量h和x表示，其相应的时间
序号由nh=[a:b]和nx=[c:d]给出。调用y=conv(h,x)将会在向量y中得到
y[n]?x[n]?h[n]的正确序列值，但是必须要确定对应的一组时间标号向量
ny。为了帮助构造向量ny，现考虑序列
h[n]??[n?a]??[n?b]和
x[n]??[n?c]??[n?d]，用解析法求卷积y[n]?x[n]?h[n]。根据所得结果，确
定利用a，b，c和d表示的ny应该是什么。为了验证结果，证实当
a?0,b?N?1,c?0和d?M?1时，y[n]的长度是
M+N-1。
答：
ny=a+c：b+d； 因为当a?0,b?N?1,c?0和d?M?1时，ny=0：M+N-2；因此y[n]的长度是M+N-1
3．考虑由下式给出的输入x[n]和单位冲激响应h[n]
?1?x[n]???u[n?2] ?2?h[n]?u[n?2]n
如果想用conv计算
0?n?24y[n]?x[n]?h[n]，就必须处理x[n]和h[n]的无限长问题。将
的x[n]的值存入向量x，将0?n?14的的h[n]值存入向量h中，再将调用
函数conv(h,x)的结果存入向量y中。因为已经将h[n]和x[n]截断了，要论证conv的输出只有一部分是真实的。试标明在输出中哪些值是真实的，哪些值不是真实的。求参数a,b,c和d的值，以使得nx=[a:b]和nh=[c:d]，并由2的答案构成y的正确时间序号。利用stem画出代码如下：
function x=heaviside3(n); k=length(n); for j=1:k; if n(j)>-1; x(j)=1; else x(j)=0; end; end;
y[n]并指出y[n]中哪些值是真实的，哪些值不真实。
n1=0:24; n2=0:14;
x=(1/2).^n.*heaviside3(n-2); h=heaviside3(n2+2); y=conv(x,h) stem(0:38,y);
结果分析：输出之中，前面15个是真实的，后面的都不真实，a=2，b=26，c=-2，d=12，图中y[n]中前面15个是真实的，后面的都不真实 中等题
对于这些练习将研究一种称之为块卷积的方法，这一方法经常用于音乐或语音处理系统的数字滤波器的事实实现中，因为这是希望有较短的处理延时。这一方法特别在用一个相对较短的滤波器处理一个很长的输入序列时最为有用。将输入序列分成一些很短的段，其中每一段都能用相当少的延时单独进行处理。卷积的线性特性能保证所有各段的输出叠加就等于整个序列与滤波器单位冲激响应的卷积。
例如：假设有限长单位冲激响应h[n]的滤波器仅在0?n?P?1内为非零，输入
?序列x[n]的长度比P大很多。现将x[n]分成长度为L的一些段x[n]??r?0xr[n?rL]，
式中L?P，且
?x[n?rL] 0?n?L-1xr[n]???0 其余n
3． 对于h[n]?(0.9)n(u[n]?u[n?10])和x[n]?cos(n2)sin(2?n5)，直接利用conv计算
0?n?99内的y[n]?h[n]?x[n]，并用stem画出y[n]。
代码：
图形：
5．设L?50，现将x[n]分成两个序列。计算y0[n]?h[n]?x0[n]和y1[n]?h[n]?x1[n]，这里x0[n]和x1[n]分别是x[n]的前50个和后50个样本。输出y[n]的形式给出
y[n]?y0[n]?y1[n?k]。
求出合适的k值并注意y0[n]和y1[n]都是长度为L?P?1。当y0[n]和y1[n]相加在一
起时，一般一定有一个两者都不为零的区域。正是这个原因，这种块卷积的方法称为重叠相加法。用这种方法计算y[n]，并画出0?n?99内的y[n]，所得结果与4求得的一样吗？ 代码： n=0:99;
h=(0.9).^n.*(heaviside3(n)-heaviside3(n-10)); n1=0:49;
x0=cos(n1.^2).*sin(2*pi.*n1./5); n2=50:99;
x1=cos(n2.^2).*sin(2*pi.*n2./5); y0=conv(x0,h); y1=conv(x1,h); y=zeros(1,199); for j=1:109; if j<50; y(j)=y0(j); else if j<59;
y(j)=y0(j)+y1(j-49); else
y(j)=y1(j-49); end;
end; a=length(y); stem(0:a-1,y); axis([0 100 -2 3]); 图形：