第5章时域离散系统的网络结构程序.doc

成 绩： 《数字信号处理》 作业与上机实验 （第五章） 班 级： 13-电信-2BF 学 号： 姓 名： 任课老师： 李宏民 完成时间： 15.11.12 信息与通信工程学院 2015—2016学年 第 1学期 第5章 时域离散系统的网络结构 1、对信号x(n)={1,2,3,4,5,6,7,8}，n=0，1，2....7,应用如下系统进行处理： 分别用迭代法按差分方程直接编程、按照直接型网络结构编程、调用MATLAB的filter函数实现该信号的处理，求得y(n)。 用三种方法求初始条件为0时的y(n); 迭代法: ①代码： m=0:7; x=[1,2,3,4,5,6,7,8]; y(1)=x(1); y(2)=x(2)+0.6*x(1)-0.1*y(1); for n=3:8 y(n)=x(n)+0.6*x(n-1)+0.7*x(n-2)-0.1*y(n-1)-0.2*y(n-2); end stem(m,y); title(信号y); ylabel(幅值); xlabel(m); ②处理后的信号y波形见图一： 图一 信号y波形图 ③处理后信号y(n)的数值： y =1.0000 2.5000 4.4500 6.2550 7.9845 9.7506 11.5280 13.2971 直接型网络结构编程: ①代码： m=0:7; x=[1,2,3,4,5,6,7,8]; M1=0; M2=0; for n=1:8 w(n)=x(n)-0.1*M1-0.2*M2; y(n)=w(n)+0.6*M1+0.7*M2; M2=M1; M1=w(n); end stem(m,y); title(信号y); ylabel(幅值); xlabel(m); ②处理后的信号y波形见图二： 图二 信号y波形图 ③处理后信号y(n)的数值： y =1.0000 2.5000 4.4500 6.2550 7.9845 9.7506 11.5280 13.2971 调用MATLAB的filter函数: ①代码： m=0:7; x=[1,2,3,4,5,6,7,8]; a=[1,0.1,0.2]; b=[1,0.6,0.7]; y=filter(b,a,x); stem(n,y); title(信号y); ylabel(幅值); xlabel(m); ②处理后的信号y波形见图三： 图三 信号y波形图 ③处理后信号y(n)的数值： y =1.0000 2.5000 4.4500 6.2550 7.9845 9.7506 11.5280 13.2971 当初始条件为：y(-1)=-0.1,y(-2)=0.2,x(-1)=0.3,x(-2)=0.8时，用三种方法求y(n); (1)迭代法: ①代码： m=0:7; x=[1,2,3,4,5,6,7,8]; y(1)=x(1)+0.6*0.3+0.7*0.8-0.1*(-0.1)-0.2*0.2; y(2)=x(2)+0.6*x(1)+0.7*0.3-0.1*y(1)-0.2*(-0.1); for n=3:8 y(n)=x(n)+0.6*x(n-1)+0.7*x(n-2)-0.1*y(n-1)-0.2*y(n-2); end stem(m,y); title(信号y); ylabel(幅值); xlabel(m); ②处理后的信号y波形见图四： 图四 信号y波形图 ③处理后信号y(n)的数值： y =1.7100 2.6590 4.2921 6.2390 8.0177 9.7504 11.5214 13.2978 (2)按照直接型网络结构编程: ①代码： m=0:7; x=[1,2,3,4,5,6,7,8]; v0=x; v1(1)=0.3