数字信号处理实验三 卷积

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

3-1

x = [1,1,2,3,5,8]; nx=length(x)

h = [5,7,8,3,3,6,4,2]; nh=length(h)

L = nx+nh-1

n = 0:(L-1);

y = conv(x,h)

stem(n,y,'*');

title('x 与h 的卷积结果')

xlabel('n');ylabel('y(n)');grid;

0246

81012020

40

60

80

100

120

140

x 与h 的卷积结果

n y (n )

3-2

function [y,ny] = convwthn(x,nx,h,nh)

ny1 = nx(1)+nh(1);

ny2 = nx(end)+nh(end);

y = conv(x,h);ny = [ny1:ny2];

3-3

x = [3,11,7,0,-1,4,2]; nx = -3:3

h = [2,3,0,-5,2,1]; nh = -1:4

[y,ny] = convwthn(x,nx,h,nh);

stem(y,ny,'*');title('y 的波形');

xlabel('n');ylabel('y(n)');grid;

-60-40-200

204060-4-2

2

4

6

8

y 的波形

n y (n )

3-4

nx = -10:10;

nh = -10:10;

x = 5+3*cos(0.2*pi*nx);

ay = [1,-0.5,0.25,0];

ax = [1,2,0,1];

m = impseq(0,-10,10); %输入一个单位冲激函数

h = filter(ax,ay,m); %用filter 函数求系统的系统函数

[y,ny] = convwthn(x,nx,h,nh) %利用卷积的方法求系统的输出 stem(ny,y,'*');title('系统输出y(n)');

xlabel('n');ylabel('y(n)');grid;

-20-15-10-50

5101520-50

5

10

15

20

25

30

35

40

45

系统输出y(n)

n y (n )

相关文档
最新文档