离散卷积的算法分析及MATLAB实现(实例程序)

摘自：张登奇,陈佳.离散卷积的算法分析及MATLAB实现[J].湖南理工学院学报(自然科学版). 2013(02)

摘要：离散卷积是信号处理的基本运算,快速卷积和分段卷积是计算离散卷积的重要算法.文章以离散线性卷积的概念为基础,介绍了计算卷积的常用方法和运算流程,列举了MATLAB实现的程序.这些程序能动态演示卷积运算的全部过程,既可帮助理解卷积运算的原理流程,也可作为教学讲解的演示工具. 1、离散卷积常规算法的动态演示程序

clear;clc;close all;format compact;

xn=[5,4,3,2,1],M=length(xn),%输入任意序列并计算长度M

hn=[1,1,1],N=length(hn),%输入任意脉冲响应并计算长度N

m=[-(M-1):M+N-2],%设置代换变量的范围以便x(m)翻转和移位

xm=[zeros(1,M-1),xn,zeros(1,N-1)],%补零以便与m对应绘图

subplot(4,1,1);stem(m,xm,'r.'),%%绘输入序列x(m)

ylabel('x(m)'), grid,title('(a) 输入序列x(m)'),%%加标签网格和标题

hm=[zeros(1,M-1),hn,zeros(1,M-1)],%补零以便与m对应绘图

subplot(4,1,2);stem(m,hm,'r.'),%%绘脉冲响应

ylabel('h(m)'),grid,title('(b) 脉冲响应h(m)'),%%加标签网格和标题

yn=zeros(1,2*M+N-2);%卷积输出初始化

for n=0:M+N-2,%逐个计算卷积输出

if n==0;xmfy=[fliplr(xn),zeros(1,M+N-2)],%实现翻转

else for k=M:-1:1; xmfy(k+n)=xmfy(k+n-1);end ; xmfy(n)=0; xmfy,%实现翻转后移位并显示

end

subplot(4,1,3);stem(m,xmfy,'b.'),%%绘制翻转移位序列

ylabel('x(n-m)'), grid,title('(c) 翻转移位x(n-m)'),%%加标签网格和标题

yn(M+n)=sum(xmfy.*hm);%计算第n位输出并与m位置对应

subplot(4,1,4);stem(m,yn,'r.'),%%绘制卷积输出序列

axis([min(m),max(m),min([0,conv(xn,hn)]),max([0,conv(xn,hn)])]),%%控制绘图坐标

ylabel('y(n)'), grid,title('(d) 卷积输出y(n)'),%%加标签网格和标题

pause,%每算完一位卷积输出后暂停,按任意键继续

end

yn, %显示线性卷积结果

2、循环卷积运算过程的动态演示程序

clc;close all;clear;format compact;

xn=[5,4,3,2,1],M=length(xn),%输入任意序列并计算长度M

hn=[1,1,1],N=length(hn),%输入任意脉冲响应并计算长度N

L=8,if Lm=0:L-1,%设置代换变量并与两圆周上的编号对应

xm=[xn,zeros(1,L-M)],%按圆周编号方向排列不足补零

subplot(4,1,1);stem(m,xm,'r.'),%%绘输入序列x(m)

ylabel('x(m)'),grid,title('(a) 输入序列x(m)'),%%加标签网格和标题