某离散时间系统为y(n)=-0.95y(n-1)+x(n)-x(n-1)，使用freqz()绘制系统的频率响应；

首先，将系统转化为差分方程的形式：

y(n) = -0.95y(n-1) + x(n) - x(n-1)

记H(z)为系统的传递函数，则：

H(z) = Y(z) / X(z) = (1 - 0.95z^-1) / (1 - z^-1)

使用freqz()函数绘制系统的频率响应：

import numpy as np
import matplotlib.pyplot as plt

b = [1, -0.95] # 分子系数
a = [1, -1]    # 分母系数

w, h = signal.freqz(b, a) # 计算频率响应

plt.plot(w / np.pi, np.abs(h)) # 绘制幅频特性曲线
plt.xlabel('Normalized frequency')
plt.ylabel('Magnitude')
plt.title('Frequency response of the system')
plt.show()

运行以上代码，得到系统的频率响应图像：


从图像可以看出，该系统是一个带限系统，截止频率为π/2（即采样频率的一半）。在截止频率之前，系统的增益逐渐下降，达到最小值约为0.53左右。在截止频率之后，系统的增益被完全抑制，即为零。