资源描述:
《第5章时域离散系统的网络结构》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、成绩:《数字信号处理»作业与上机实验(第二章)班级:学号:姓名:任课老师:完成时间:信息与通信工程学院2014—2015学年第1学期第5章时域离散系统的网络结构1、对信号x(n)二{1,2,3,4,5,6,7,8},n二0,1,2....7,应用如下系统进行处理:y(n)+0.ly(n-1)+0.2y(n-2)=x(n)+0.6x(/?-1)4-0.7x(n一2)分别用迭代法按差分方程直接编程、按照直接型网络结构编程、调用MATLAB的filter函数实现该信号的处理,求得y(n)。1)用三种方法求初始条件为0时的
2、y(n);(1)迭代法:①Matlab代码:n=0:7;x=[1,2,3,4,5,6,7,8];y(l)=x(l);y(2)=x(2)+0.6*x(l)-0.1*y(1);form=3:8y(m)=x(m)+0.6*x(m-l)4-0.7*x(m-2)-0.1*y(m-l)-0.2*y(m-2);endstem(nzy);title(1信号y(n)*);ylabel(M畐值冷;xlabel(*n*);①处理后的信号y(n)波形见图一:图一信号y(n)波形图③处理后信号y(n)的数值:7.9845y=1.00002.
3、50004.45006.25509.750611.528013.2971(2)按照直接型网络结构编程:©Matlab代码:n=0:7;x=[1,2,3Z4Z5,6,7,8];ml=0;m2=0;form=l:8w(m)=x(m)-0-l*ml-0•2*m2;y(m)=w(m)+ml*0.6+0.7*m2;m2=ml;ml=w(m);endstem(n,y);title(1信号y(n)1);ylabel('l
4、®值');xlabel(1n*);②处理后的信号y(n)波形见图二:fs号y(n)图二信号y(n)波形图③处
5、理后信号y(n)的数值:y二1.00002.50004.45006.25507.98459.750611.528013.2971(3)调用MATLAB的filter函数:①Matlab代码:n=0:7;x=[1,2,3,4,5,6,7,8];a=[lz0.1z0.2];b=[lz0.6z0.7];y=filter(b,azx);stem(nzy);title(*信号y(n)*);ylabel(■幅值冷;xlabel('n*);②处理后的信号y(n)波形见图三:信号y(n)图三信号y(n)波形图③处理后信号y(n)的
6、数值:y=1.00002.50004.45006.25507.98459.750611.528013.29712)当初始条件为:y(-l)=-0.l,y(-2)=0.2,x(-1)=0.3,x(-2)=0.8时,用三种方法求y(n);(1)迭代法:©Matlab代码:n=0:7;x=[1,2Z3,5Z6r7,8];y(l)=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(l)+0.7*0.3-0.1*y(l)-0.2*(-0.1);form=3:8y(
7、m)=x(m)+0.6*x(m-1)+0.7*x(m-2)-0.l*y(m-1)-0.2*y(m-2);endstem(n,y);title(*信号y(n)*);ylabel(*幅值T;xlabel();②处理后的信号y(n)波形见图四:信号y(n)图四信号y(n)波形图③处理后信号y(n)的数值:y=1.71002.65904.29216.23908.01779.750411.521413.2978(2)按照直接型网络结构编程:©Matlab代码:n=0:7;x=[1,2,3,4,5,6,7,8];v0=x;vl
8、(l)=0.3;v2(1)=0.8;wl(1)=-0.1;w2(1)=0.2;form=l:8wO(m)=-0.l*wl(m)-0.2*w2(m)+v0(m)+0.6*vl(m)+0.7*v2(m);y(m)=wO(m);v2(m+1)=vl(m);w2(m+1)=wl(m);vl(m+1)=vO(m);wl(m+1)=wO(m);endstem(nzy);title(*信号y(n)*);ylabel(*幅值T;xlabel('n');②处理后的信号y(n)波形见图二:fS号y(n)图二信号y(n)波形图③处理后信
9、号y(n)的数值:8.0177y二1.71002.65904.29216.23909.750411.521413.2978(3)调用MATLAB的filter函数:①Matlab代码:n=0:7;x=[l,2r3,4Z5,6,1,8];a=[lz0.1z0.2];b=[lz0.6,0.7];xs=[0.3Z0.8];ys=[-0.1,0.2];xi=filt