电子信息实验杭六.pdfVIP

实验六卷积分步演示

一、实验目的

1.掌握卷积的原理及其计算过程；

2.进一步熟悉MATLAB程序语言与程序调试。

二、实验内容

xnNhnM

1.自行编写线性卷积函数：任意给定两组离散序列()（点）和()（

点），自行编写线性卷积函数，输出并画出卷积结果：

∞

yn=xkhn−k

()()()

k=−∞

并将计算结果与MATLAB自带线性卷积函数conv进行对比。

（提示：可能用到圆周移位函数circshift和数值翻转函数fliplr）

自定义函数：

EditorWindow

functionyn=convv(xn,hn)

hn=hn;

n=length(hn)+length(xn)-1;

xn(length(xn)+1:n)=0;hn(length(hn)+1:n)=0;

yn=ones(1,n);

fori=1:n+1

a=0;

fork=1:i-1

a=a+xn(k)*hn(i-k);

end

yn(i)=a;

end

yn=yn(2:n+1);

对比：

EditorWindow

clc,clear;

xn=1:10;hn=1:21;

subplot(2,1,1)

y1=convv(xn,hn)

stem(y1)

subplot(2,1,2)

y3=conv(xn,hn)

stem(y3)

2.自行编写圆周卷积函数：任意给定两组离散序列xnxn

1()和2()，自行编写

N点圆周卷积函数，输出并画出N点圆周卷积

N−1

xnxnRxmxnm

1()⊗N2()=N1()2((−)modN)，

m=0

其中n是整数，k表示k整除N后的余数，Rzn表示对序列zn

()modNN()()

取主值区间。（提示：可能用到圆周移位函数circshift和数值翻转函数fliplr）

M函数：

EditorWindow

functiony=cirshift(x,m,N)

x=[xzeros(1,N-length(x))];

n=[0:1:N-1];

n=mod(n-m,N);

y=x(n+1);

functiony=circonvt(x1,x2,N)

x1=[x1