ch5作业.doc

成 绩： 《数字信号处理》 作业与上机实验 （第五章） 班 级： 信工11-1BF 学 号： 14113901118 姓 名： 谢清宝 任课老师： 李宏民 完成时间： 2014-5-8 信息与通信工程学院 2013—2014学年第 2 学期 第5章 时域离散系统的网络结构 1、对信号x(n)={1,2,3,4,5,6,7,8}，n=0，1，2....7,应用如下系统进行处理： 分别用迭代法按差分方程直接编程、按照直接型网络结构编程、调用MATLAB的filter函数实现该信号的处理，求得y(n)。 用三种方法求初始条件为0时的y(n); 用迭代法差分方程编程如下所示: n=0:7;x=1:8; y(1)=x(1); y(2)=x(2)+0.6*x(1)-0.1*y(1); for m=3:1:8 y(m)=x(m)+0.6*x(m-1)+0.7*x(m-2)-0.1*y(m-1)-0.2*y(m-2); end plot(n,y); 经过差分方程处理后的结果如下图1.1.1所示: 用直接型网络结构编程代码如下： M1=0; M2=0; for m=0:7 w1=x(m+1)-0.1*M1-0.2*M2; w2=w1; y(m+1)=w2+0.6*M1+0.7*M2; M2=M1; M1=w2; end plot(n,y); 直接型处理后的结果如图1.1.2所示 用filter函数写代码如下所示： x=[1,2,3,4,5,6,7,8]; S=[1,0.6,0.7,1,0.1,0.2]; G=1; [B,A]=sos2tf(S,G); y=filter(B,A,x); plot(n,y); 结果如图1.1.3所示: 2)当初始条件为：y(-1)=-0.1,y(-2)=0.2,x(-1)=0.3,x(-2)=0.8时，用三种方法求初始条件为0时的y(n); 有初始条件下,求初始条件为0时的y(n),用差分方程编程的代码如下: n=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 k=3:1:8 y(k)=x(k)+0.6*x(k-1)+0.7*x(k-2)-0.1*y(k-1)-0.2*y(k-2); end plot(n,y); 差分方程处理的结果如下图1.2.1所示: 有初始条件下,求初始条件为0时的y(n),直接型编程代码如下所示: M1=-0.1;M2=0.2;N1=0.3;N2=0.8; for m=0:7 w1=x(m+1)-0.1*M1-0.2*M2; w2=w1; y(m+1)=w2+0.6*N1+0.7*N2; M2=M1; M1=y(m+1); N2=N1; N1=x(m+1); end plot(n,y); 直接型处理后的结果如下图1.2.2所示: 有初始条件下,求初始条件为0时的y(n),用filter函数写的代码如下所示: B=[1,0.6,0.7]; A=[1,0.1,0.2]; xi=[0.3 0.8]; yi=[-0.1 0.2]; xn=filtic(B,A,yi,xi); y=filter(B,A,x,xn); plot(n,y); 用filter函数处理的结果如下图1.2.3所示: 2、x(n)为受到工业高频干扰的心电信号：x(n)=[-4，-2，0，-4，-6，-4，-2，-4，-6，-6，-4，-4，-6，-6，-2，6，12，8，0，-16，-38，-60，-84，-90，-66，-32，-4，-2，-4，8，12，12，10，6，6，6，4，0，0，0，0，0，-2，-4，0，0，0，-2，-2，0，0，-2，-2，-2，-2，0] 按直接型与级联型网络结构编程实现对心电信号的处理。 直接型结构编程代码如下所示: n=0:55; x=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0]; S=[1 2 1 1 -1.2686 0.7051 ;1 2 1 1 -1.01