数字数信号处理上机实习精选报告.docVIP

仅供个人参考 专题一 离散卷积的计算 一、实验内容 设线性时不变（ LTI ）系统的冲激响应为 h(n) ，输入序列为 x(n) 1、h(n)=(0.8) n，0≤n≤4； x(n)=u(n)-u(n-4) For personal use only in study and research; not for commercial use 2、h(n)=(0.8) nu(n), x(n)=u(n)-u(n-4) 3、h(n)=(0.8) n u(n), x(n)=u(n) 求以上三种情况下系统的输出 y(n) 。 For personal use only in study and research; not for commercial use 二、实验目的 1、掌握离散卷积计算机实现。 2、进一步对离散信号卷积算法的理解 。 三、原理及算法 概要 算法：把冲激 响应 h(n) 与输 入序列 x(n) 分别输入到程序中，然后调用离散卷积函数 y=conv(x.,h) 即可得 到所要求的结果。 原理： 离散卷积定义为 y (n) x( k)h(n k) k 当序列为有限长时，则 n y(n) x(k )h(n k) k 0 四．理论计算 1、h(n)=(0.8) n，0≤n≤4； x(n)=u(n)-u(n-4) 不得用于商业用途 仅供个人参考 y(n) x(n) h( n) x(m) h(n m) m (a) 当 n0 时，y( n)=0 (b) 当 0 n 3时, n y(n) (0.8) m 0 n (c) 当 4 n 7 时， 4 y(n) (0.8) m n 3 n (d) 当 n7 时，y( n)=0 理论结果与上图实验结果图中所示吻合。 2、h(n)=(0.8) nu(n), x(n)=u(n)-u(n-4) y(n) x(n) h( n) x(m) h(n m) m (a) 当 n0 时，y( n)=0 n n (b) 当 0 n 3时, y(n) (0.8) m 0 n n (c) 当 4 n 20 时， y( n) (0.8) m n 3 20 n (d) 当 21 n 23 时， y(n) (0.8) m n 3 (e) 当 n23 时，y( n)=0 在卷积序列中，因为 n 趋于无穷，在实验中无法实现，而 0.8 的 20 次幂基本接近于 0，所以这里只取到 n 为 20 就可以了。由以上可以看出，理论结果与上图实验结果图中所 示吻合。 3、h(n)=(0.8) nu(n), x(n)=u(n) y(n) x(n) h( n) x(m) h(n m) m (a) 当 n0 时，y( n)=0 n n (b) 当 0 n 70 时 , y(n) (0.8) m 0 不得用于商业用途 仅供个人参考 (c) 当 71 n 140 时， 70 y( n) (0.8) m n 70 n (d) 当 n140 时，y( n)=0 在卷积序列中，因为 n 趋于无穷，在实验中无法实现，而 0.8 的 70 次幂的叠加基本接 近于最大值，所以这里只取到 n 为 70 就可以了。由以上可以看出，理论结果与上图实验结 果图中所示吻合。 五．程序 x1=[1 1 1 1 ];nx1=0:3; h1=[1 0.8 0.64 0.8^3 0.8^4];nh1=0:4; y1=conv(x1,h1); subplot(3,3,1);stem(nx1,x1);title( 序列 x1); xlabel(n);ylabel(x1(n)); subplot(3,3,2);stem(nh1,h1);title( 序列 h1); xlabel(n);ylabel(h1(n)); subplot(3,3,3);stem(y1);title( 序列 y1); xlabel(n);ylabel(y1(n)); x2=[1 1 1 1];nx2=0:3; nh2=0:1:20; h2=(0.8).^nh2; y2=conv(x2,h2); subplot(3,3,4);stem(x2);title( 序列 x2); xlabel(n);ylabel(x2(n)); subplot(3,3,5);stem(h2);title( 序列 h2); xlabel(n);ylabel(h2(n)); subplot(3,3,6);stem(y2);title( 序列 y2); xlabel(n);ylabel(y2(n)) nx3=0:1:20; x3=1.^nx3; nh3=0:1:20; h3=(0.8).^nh3; y3=conv(x3,h3); subplot(3,3,7);stem(nx3,x3);title( 序列 x3); xlabel(n);