dconv.m

1． 该光盘包含了“信号与系统分析及MATLAB实现”一书的第七章所有程序； 2． 为了方便读者阅读程序

function [f,k]=dconv(f1,f2,k1,k2)  
%The function of compute  f=f1*f2
% f:  卷积和序列f(k)对应的非零样值向量
% k：  序列f(k)的对应序号向量
% f1: 序列f1(k)非零样值向量
% f2:  序列f2(k)的非零样值向量
% k1:  序列f1(k)的对应序号向量
% k2:  序列f2(k)的对应序号向量
f=conv(f1,f2) 									%计算序列f1与f2的卷积和f
k0=k1(1)+k2(1);									%计算序列f非零样值的起点位置
k3=length(f1)+length(f2)-2;						%计算卷积和f的非零样值的宽度
k=k0:k0+k3                                      %确定卷积和f非零样值的序号向量
subplot(2,2,1)
stem(k1,f1)										%在子图1绘序列f1(k)时域波形图
title('f1(k)')
xlabel('k')
ylabel('f1(k)')
subplot(2,2,2)
stem(k2,f2)										%在图2绘序列f2(k)时波形图
title('f1(k)')
xlabel('k')
ylabel('f2(k)')
subplot(2,2,3)
stem(k,f);										%在子图3绘序列f(k)的波形图
title('f(k)f1(k)与f2(k)的卷积和f(k)')
xlabel('k')
ylabel('f(k)')
h=get(gca,'position');								
h(3)=2.5*h(3);
set(gca,'position',h)					%将第三个子图的横坐标范围扩为原来的2.5倍