实验四 离散系统分析

实验四 离散系统分析

实验目的

 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频率特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频率响应和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法

1.已知某离散LTI系统的差分方程为$y[k]-1.143y[k-1]+0.4128y[k-2]=0.0675x[k]+0.1349x[k-1]+0.0675x[k-2]$

1） 初始状态y[-1]=1,y[-2]=2,输入x[k]=u[k]，计算并画出系统的完全响应。

a=[1 -1.143 0.4128];
b=[0.0675 0.1349 0.0675];
N=50;
x=ones(1,N);
zi=filtic(b,a,[1,2])
y=filter(b,a,x,zi);
stem(y);

2） 当以下3个信号分别通过离散系统时，分别计算并画出离散系统的零状态响应：

a=[1 -1.143 0.4128];
b=[0.0675 0.1349 0.0675];
N=25;
k=1:N;
x1=cos(pi/10*k);
y1=filter(b,a,x1) ;
figure(1);
subplot(121);stem(x1);title('输入信号');
subplot(122);stem(y1);title('响应');
x2=cos(pi/5*k);
y2=filter(b,a,x2); 
figure(2);
subplot(121);stem(x2);title('输入信号');
subplot(122);stem(y2);title('响应');
x3=cos(7*pi/10*k);
y3=filter(b,a,x3); 
figure(3);
subplot(121);stem(x3);title('输入信号');
subplot(122);stem(y3);title('响应');

2.已知某因果LTI离散系统的系统函数为

1） 计算并画出系统的单位脉冲响应，计算并画出零极点分布图

den=[1 -1.035 0.8264 -0.2605 0.04033];
num=[0.03571 0.1428 0.2143 0.1428 0.03571];
figure(1);
impz(num,den,-10:30);
figure(2);
p=roots(den)
z=roots(num)
zplane(z,p);

2） 当如下信号通过系统时，计算并画出系统的零状态响应

den=[1 -1.035 0.8264 -0.2605 0.04033];
num=[0.03571 0.1428 0.2143 0.1428 0.03571];
N=30;
f=ones(1,N)+cos(pi/4*k)+cos(pi/2*k);
k=0:1:N-1;
y=filter(num,den,f);
stem(k,y);