用Matlab 实现信号时域图形.pptVIP

信号系统与数字信号处理 指导老师： 班级： 学生： 学号： 一 课题名称 39.已知输入信号为x(n)=u(n)，系统的单位取样响应为h(n)=n[u(n)-u(n-6)]，y1(n)=x(n)*h(n)和y2(n)=x(n)*h(n+5),试分别绘出其时域图形并写出前13个数值解； （注：为了验证误差，后面将x(n)取16个点和26个点）。 二 设计目的 熟练掌握时域离散信号的平移，和，差及卷积运算。 能够熟练运用MATLAB软件进行信号的处理及仿真。 三 题目分析 由于该输入信号为一离散信号x(n)=u(n)，所以在大于零的区域内该信号的值恒为1。 由于h(n)=n[u(n)-u(n-6)]，所以根据“竖式算法”定理，可以将x(n)=u(n)改为x(n)=u(n)-u(n-7),就可以刚好求出题目所对的前十三个数值解，而且从图形的角度分析问题取值足够。 为了说明误差问题，后面取了16个点和26个点，以分析波形的变化。 四 设计原理（理论分析原理及思路） 1 U(n)= 1, n≥0 0 n0 2 序列x1(n),x2(n)的卷积和用x1(n)*x2(n)表示。其定义为 根据离散信号的卷积定义，可求出y1(n)=x(n)*h(n)与y2(n)=x(n)*h(n+5)的波形。 根据解析法将式子改为： = u(m)m[u(n-m)-u(n-m-6)] (根据离散信号的卷积和定义式) = mu(m)u(n-m)- mu(m)u(n-m-6) - 根据等差数列求和公式得出： = u(n)- =15 结论：当取值为无限大时，该卷积的极限等于15. 根据解析法将式子改为： u(n) 计算y2(n)=x(n)*h(n+5) 的卷积和还是利用离散信号的卷积定义 注：由于只发生横向移位，所以它的极限值不变，这里再不做过程。 五 matlab源程序及其代码 1 U(n)取0到7的源程序： x1=[0:7]; x1=[1 1 1 1 1 1 1 1]; x2=[1 1 1 1 1 1 1 1]; h1=[0:5]; h1=[0 1 2 3 4 5]; h2=[-5:0]; h2=[0 1 2 3 4 5]; y1=[0:10]; y1=[0 1 3 6 10 15 15 15 15 14 12 9 5] y2=[-5:5]; y2=[0 1 3 6 10 15 15 15 15 14 12 9 5] y1=conv(x1,h1); y2=conv(x2,h2); subplot(3,2,1),stem(0:7,x1,.),axis([0 7 0 2]),ylabel(x1(n)); subplot(3,2,2),stem(0:7,x2,.),axis([0 7 0 2]),ylabel(x2(n)); subplot(3,2,3),stem(0:5,h1,.),axis([0 5 0 6]),ylabel(h1(n)); subplot(3,2,4),stem(-5:0,h2,.),axis([-5 0 0 6]),ylabel(h2(n)); subplot(3,2,5),stem(0:12,y1,.),axis([0 12 0 15]),ylabel(y1(n)); subplot(3,2,6),stem(-5:7,y2,.),axis([-5 7 0 15]),ylabel(y2(n)); 2 u(n)取0到15的源程序： x1=[0:15]; x1=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]; x2=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]; h1=[0:5]; h1=[0 1 2 3 4 5]; h2=[-5:0]; h2=[0 1 2 3 4 5]; y1=[0:10]; y1=[0 1 3 6 10 15 15 15 15 15 15 15 15 15 15 15 15 14 12 9 5]; y2=[-5:5]; y2=[0 1 3 6 10 15 15 15 15 15 15 15 15 15 15 15 15 14 12 9 5]; y1=conv(x1,h1); y2=conv(x2,h2); subplot(3,2,1),stem(0:15,x1,.),axis([0 15 0 2]),ylabel(x1(n)); subplot(3,2,2),stem(0:15,x2,.),axis([0 15 0 2]),ylabel(x2(n)); subplot(3,2,3),stem(0:5,h1,.),axis([0 5 0 6]),ylabel(h1(n)); subplot(3,2,4),stem(-5:0,h2,.),