数字信号处理实验三 卷积
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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 )