利用MATLAB实现循环卷积.pdfVIP

利用MATLAB实现循环卷积

还是那句话，这么详细地去讲解循环卷积的计算过程目

的并不是为了的让大家去手工

计算循环卷积，只是为了让大家对“循环卷积”有更清楚

的认识罢了。要计算的话让计算

机去做就好了。

tab中没有现成的计尊循环卷积的函数，需要自己编一

个，这里给大家准备了一个，

参见附件中的：circonv.m文件。

下面看一个循环卷积的实例。

两个长度为8的序列：

Xn=[22220000]

Yn=[01234567]

计算一下这两个序列的循环卷积。

»=8;

»n=0:l:7;

»Xn=[22220000];

»Yn=[O1234567];

»Sn=circonv(Xn,Yn,);

»subplot(31l);stem(n,Xn);

»subplot(312);stem(n,Yn);

»subplot(313);stem(n,Sn);

Matlab卷积函数conv2

功能：是二维卷积运算函数(与convmtx2相似)。

如果a和b是两个离散变量nl和n2的函数，则关于a

和b的二维卷积运算数学公式如下：

S8

以〃「〃2)=II左2)b(〃i一九1，〃2一42)

AI=coZ22=~o

用法：

C=conv2(A,B)

C=conv2(hcol,hrow,A)

,,

C二conv2(...,shape)

C=conv2(A,B)计算数组A和B的卷积。如果一个数组

描述了一个二维FIR滤波器，则另一个数组被二维滤波。当

A的大小为［ma.na］,B的大小为［mb,nb］时，C的大小为

5

［ma+mb-1,mb+nb-1］。shape见下表

参数值含义

full默值，返叵全部二维卷积值。

same返回与A大小相同卷积值的中间部分

valid当all(size(A)=size(B)),C的大小为m［a+mb-l,mb+nb-l］；否则，C返回口。在n维卷

积运算中,C的大小为max(size(A)-size(B)+l,O)

例子:

s=U21;000;-1-2-1];

A=zeros(lO);

A(3:7,3:7)=ones(5);

H=conv2(A,s);

mesh(H)

Matlab的帮助：

conv2

2-Dconvolution

Syntax

C=conv2(A,B)

C=conv2(hcoI,hrow,A)

C=conv2(.../shape*)

Description

C=conv2(A,B)computesthetwo-dimensional

convolutionofmatricesAandB.Ifoneo