17-18-离散时间信号和离散时间系统-Matlab实现-A方案.ppt

% 线性卷积 x = [1 2 3 4]; h = [1 2 3]; y = conv(h,x); n = 0:5; %画图 stem(n,y); xlabel(Time index n); ylabel(Amplitude); title(Output Obtained by Convolution); grid on; function [y,ny] = conv_m (x,nx,h,nh) % 序列y为序列x和序列h的卷积 % ny,nx,nh 分别为y,x和h的位置向量 ny0 = nx(1)+nh(1); % 卷积后位置初值的计算 nyf = nx(end) + nh(end); % 卷积后位置终值的计算 y = conv(x,h); % 卷积序列数值的计算 ny = [ny0 : nyf]; % 卷积序列位置向量的计算 % 给定输入序列 nx = [-4:2]; x = [3,-3,7,0,-1,5,2]; % 给定脉冲响应序列 nh = [-1:4]; h = [2,3,0,-5,2,1]; % 带位置序列的卷积结果 [y,ny] = conv_m (x,nx,h,nh) stem(ny,y); xlabel(Time index n