信号与系统-实验6-离散时间系统的 Z 域分析

一、实验目的

1、掌握 z 变换及其性质；了解常用序列的 z 变换、逆 z 变换；

2、掌握利用 MATLAB 的符号运算实现 z 变换；

3、掌握利用 MATLAB 绘制离散系统零、极点图的方法；

4、掌握利用 MATLAB 分析离散系统零、极点的方法；

5、掌握利用 MATLAB 分析离散系统频率响应的方法。

二、实验内容

1、分别求cos(kn)和  的 z 变换。

clear;
syms k n a;
f=cos(k*n);
f1=a^k;
disp("cos(kn)的z变换为");
y=ztrans(f)
disp("a^k的z变换为" );
y1=ztrans(f1)

2、系统函数为

分析并绘制出离散系统的零、极点图。

clear;
num=[0.2,0.1,0.3,0.1,0.2];
den=[1,-1.1,1.5,-0.7,0.3];
subplot(121)
zplane(num,den)

num1=[3,-5,10,0];
den1=[1,5,2,-10];
subplot(122)
zplane(num1,den1)

3、绘制 的频率响应图形。

clear;
num=[1,1];
den=[2 2 1];

[H,w]=freqz(num,den);
Hf=abs(H);  %取幅度值实部
Hx=angle(H);  %取相位值对应相位角
clf
subplot(121)
plot(w,20*log10(Hf))  %幅值变换为分贝单位
title('离散系统幅频特性曲线')
subplot(122)
plot(w,Hx)
title('离散系统相频特性曲线')

4、已知某一离散系统的系统函数为 ，试用 MATLAB 绘出该系统的零、极点分布图及幅频特性曲线，并分析该系统 的频率特性。

clear;
B=[1 5 5 1];
A=[1 0 0 0];
N=60;
[h,t]=impz(B,A,N);

subplot(211)
plot(t,h)
title("h(t)");
xlabel("N");ylabel("h(t)");
subplot(212)
zplane(B, A) 

figure
freqz(B,A);

5、已知某离散系统的系统函数为  ，试用 MATLAB 求出 该系统的单位响应h(k)。

clear;
num=[1,0,0]
den=[1,5,6]
impz(num,den)

6、已知某序列的 z 变换为 ，试用 MATLAB 求 F(z) 的逆变换。 

clear;
syms z;
F=(z^2+2*z)/((z-0.1)*(z-0.2)*(z-0.3));
X=iztrans(F)

7、已知某离散系统的系统函数为 ，试用 MATLAB 分析 该系统的频率特性，绘制其幅频及相频特性曲线。

B=[1.25 1.25];
A=[1 -0.25];
freqz(B, A)