数字的信号处理实习报告材料68046.pdf

实用标准文案 实验一 一．实验内容 设线性时不变（ LTI ）系统的冲激响应为 h(n) ，输入序列为 x(n) 1、h(n)=(0.8) n，0≤n≤4 ； x(n)=u(n)-u(n-4) n 2、h(n)=(0.8) u(n), x(n)=u(n)-u(n-4) n 3、h(n)=(0.8) u(n), x(n)=u(n) 求以上三种情况下系统的输出 y(n) ，显示输入和输出波形。 二．实验目的 1，熟悉离散卷积的实现过程 2，掌握离散卷积的计算方法 三．算法设计 离散卷积定义为 y (n) x( k )h(n k ) k 当序列为有限长时，则 n y (n) x (k)h( n k) k 0 这里我们把冲激响应 h(n) 与输入序列 x(n) 分别输入到程序中，然后调用离散卷 积函数 y=conv(x.,h) 即可得到所要求的结果。 四．程序代码、运行结果与结果对比分析 程序代码如下： x1=[1 1 1 1 ];nx1=0:3; nh1=0:4; h1=0.8.^nh1; [y1,ny1]=conv_m(x1,nx1,h1,nh1); figure(1); subplot(221);stem(nx1,x1);title(' 序列 x1'); xlabel('n');ylabel('x(n)'); subplot(222);stem(nh1,h1);title(' 序列 h1'); xlabel('n');ylabel('h(n)'); subplot(223);stem(ny1,y1);title(' 序列 y1'); xlabel('n');ylabel('y(n)'); x2=[1 1 1 1 ];nx2=0:3; nh2=0:20; h2=0.8.^nh2; [y2,ny2]=conv_m(x2,nx2,h2,nh2); 精彩文档 实用标准文案 figure(2); subplot(221);stem(nx2,x2);title(' 序列 x2'); xlabel('n');ylabel('x(n)'); subplot(222);stem(nh2,h2);title(' 序列 h2'); xlabel('n');ylabel('h(n)'); subplot(223);stem(ny2,y2);title(' 序列 y2'); xlabel('n');ylabel('y(n)'); x3=ones(1,20);nx3=0:19; nh3=0:20; h3=0.8.^nh3; [y3,ny3]=conv_m(x3,nx3,h3,nh3); figure(3); subplot(221);stem(nx3,x3);title(' 序列 x3'); xlabel('n');ylabel('x(n)'); subplot(222);stem(nh3,h3);title(' 序列 h3'); xlabel('n');ylabel('h(n)'); subplot(223);stem(ny3,y3);title(' 序列 y3'); xlabel('n');ylabel('y(n)'); function [y,ny]=conv_m(x,nx,h,nh) nyb=nx(1)+nh(1);nye=nx(length(x))+nh(length(h)); ny=[nyb:nye]; y=conv(x,h); end 运行结果如图 1、图 2、图 3 所示：