杭电数字信号处理实验1

1、信号、系统与信号处理实验n实验报告专业word可编辑姓 名： 王健学号:14072119班级：14083413上课时间：周五-六七八实验名称：离散时间系统的时域特性分析实验目的线性时不变离散系统在时域中可以通过常系数线性差分方程来描述，冲击响应序列可以刻画其时域特性。本实验通过使用matlab函数研究离散时间系统的时域特性，以加深对离散时间系统的差分方程，冲击响应和系统的线性和时不变特性的理解。实验内容与要求考虑如下差分方程描述的两个离散时间系统：系统 1: y(n)=0.5*x(n)+0.27*x(n-1)+0.77*x(n-2);系统 2: y(n)=0.45*x(n)+0.5*x(n-1

2、)+0.45*x(n-2)+0.53*y(n-1)-0.46*y(n-2);输入;x(n)= x=cos(20*pi*n/256)+cos(200*pi*n/256); 0<=n<=299;(1)编程求上述两个系统的输出，并分别画出系统的输入与输出波形(2)编程求上述两个系统的冲击响应序列，并画出其波形(3)若系统的初始状态为0;判断系统2是否为时不变？是否为线性？三、实验程序与结果(1)编程求上述两个系统的输出 ，并分别画出系统的输入与输出波形clear all;clc;close all;n=0:299;x=cos(20*pi*n/256)+cos(200*pi*n/256);

3、b1=0.5 0.27 0.77;a1=1；y1=filter(b1,a1,x);b2=0.45 0.5 0.45;a2=1 -0.53 0.46;y2=filter(b2,a2,x);subplot(3,1,1);stem(n,x);title('输入');subplot(3,1,2);stem(n,y1);title('系统 1 输出');subplot(3,1,3);stem(n,y2);title('系统 2 输出');结果：(2)编程求上述两个系统的冲击响应序列，并画出其波形n=0:299;x=cos(20*pi*n/256)+cos(

4、200*pi*n/256);b1=0.5 0.27 0.77;a1=1；y1=impz(b1,a1,30);b2=0.45 0.5 0.45;a2=1 -0.53 0.46;y2=impz(b2,a2,30);subplot(2,1,1);stem(y1);title('系统 1 冲激响应');subplot(2,1,2);stem(y2);title('系统 2 冲激响应');结果:I O力 0O色口口 O口由口00口内00 4O 口 051D 152025.系统2冲激值应(3)若系统的初始状态为0;判断系统2是否为时不变？是否为线性？%验证系统时不变n=0:

5、199;D=50;x=cos(20*pi*n/256)+cos(200*pi*n/256);xd=zeros(1,D) x;b2=0.45 0.5 0.45;a2=1 -0.53 0.46;ic=0 0;y2=filter(b2,a2,x,ic);y2d=filter(b2,a2,xd,ic);subplot(3,1,1);stem(n,y2);title('原输出');subplot(3,1,2);stem(n,y2d(1:length(y2);title('延迟后产生的输出');subplot(3,1,3);stem(n,y2d(1+D:length(y2)

6、+D)-y2);title('差值信号);结果:%验证系统线性n=0:199;q1=2;q2=-3;x1=cos(20*pi*n/256)+cos(200*pi*n/256);x2=sin(20*pi*n/256)+sin(200*pi*n/256);x=q1*x1+q2*x2;b=0.45 0.5 0.45;a=1 -0.53 0.46;y1=filter(b,a,x1);y2=filter(b,a,x2);y=filter(b,a,x);subplot(3,1,1);stem(n,q2*y2+q1*y1);title(' q2*y2+q1*y1');subplot(

7、3,1,2);stem(n,y);title('Ta*x1+b*x2');subplot(3,1,3);stem(n,y-(q2*y2+q1*y1);title('差值信号');结果四、仿真结果分析1 .两个系统的冲击响应都为有限序列 ，系统1在n大于4的时候序列都为 0,系统2 在n大于13的时候序列都为0。其都为稳定系统。2.验证系统时不变，其差值信号都为0,表明原输出和延迟D个时间之后的输出相同，说明其是一个时不变系统。验证线性得出的差值信号是一组无规律的其值为10的-15次方数量级的序列，因为10的-15次方足够小，可以说明输出 Ta*x1+b*x2和q2*y2+q1*y1 是一串相同的序列，为此 可以表明其是一个线性系统 。四、实验问题解答与体会第一次数字信号处理实验 ，虽然题目看起来简单，但是编程的时候却会有陷阱 ，加上 自己的粗心用了好久才完成 。以后做实验一定不能大意 ，一定要预习，特别