信号与系统分析2022春季作业-参考答案:第十三次作业

Posted 卓晴

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了信号与系统分析2022春季作业-参考答案:第十三次作业相关的知识,希望对你有一定的参考价值。

作业题目 目 录
Contents
参考答案 系统频率特性 离散傅里叶变换 离散傅里叶变换性质

 

§00 业题目

  作业要求链接: 信号与系统2022春季作业-第十三次次作业 : https://zhuoqing.blog.csdn.net/article/details/124978410

 

§13 考答案


1.1 系统频率特性

1.1.1 确定系统的幅频特性

  (1) 系统函数具有两个极点,分别位于 p 1 = − 2 , p 2 = − 3 p_1 = - 2,p_2 = - 3 p1=2,p2=3 。改系统可以看成两个低通滤波器串联。
H 1 ( s ) = 1 ( s + 2 ) ( s + 3 ) ,     R e [ s ] > − 2 H_1 \\left( s \\right) = 1 \\over \\left( s + 2 \\right)\\left( s + 3 \\right),\\,\\,\\,\\mathop\\rm Re\\nolimits \\left[ s \\right] > - 2 H1(s)=(s+2)(s+3)1,Re[s]>2

▲ 图1.1.1 系统的零极点分布

  下面使用Python中的系统仿真软件包给出系统的频率特性。

▲ 图1.1.2 系统的幅频特性和相频特性

from headm import *
from scipy.signal           import bode
from scipy.signal           import TransferFunction

sys = TransferFunction([1],[1,5,6])
w,mag,phase = bode(sys)

plt.subplot(2,1,1)
plt.semilogx(w, mag)

plt.xlabel("Frequency")
plt.ylabel("Magnitude")
plt.grid(True)
plt.tight_layout()

plt.subplot(2,1,2)
plt.semilogx(w, phase)

plt.xlabel("Frequency")
plt.ylabel("Phase")
plt.grid(True)
plt.tight_layout()
plt.show()

  (2) 系统具有两个零点和两个极点。具体分布如下图所示:

H 2 ( s ) = s 2 s 2 + 2 s + 1 ,     R e [ s ] > − 1 H_2 \\left( s \\right) = s^2 \\over s^2 + 2s + 1,\\,\\,\\,\\mathop\\rm Re\\nolimits \\left[ s \\right] > - 1 H2(s)=s2+2s+1s2,Re[s]>1

▲ 图1.1.3 系统的零极点分布

  这个系统可以看成两个高通滤波串联。 下面给出系统的幅频特性和相频特性。

▲ 图1.1.4 系统的频率特性

sys = TransferFunction([1,0,0],[1,2,1])
w,mag,phase = bode(sys)

  (3) 系统具有左右对称的零点和极点。 符合全通特性 的定义。

H 3 ( s ) = s 2 − s + 1 s 2 + s + 1 ,    R e [ s ] > − 1 2 H_3 \\left( s \\right) = s^2 - s + 1 \\over s^2 + s + 1,\\,\\,\\mathop\\rm Re\\nolimits \\left[ s \\right] > - 1 \\over 2 H3(s)=s2+s+1s2s+1,Re[s]>21

▲ 图1.1.5 系统的零极点分布

▲ 图1.1.6 系统的频率特性

from scipy.signal           import bode
from scipy.signal           import TransferFunction

sys = TransferFunction([1,-1,1],[1,1,1])
w,mag,phase = bode(sys)

plt.subplot(2,1,1)
plt.semilogx(w, mag)

plt.xlabel("Frequency")
plt.ylabel("Magnitude")
plt.axis([min(w), max(w), -1, 1])
plt.grid(True)
plt.tight_layout()

plt.subplot(2,1,2)
plt.semilogx(w, phase)

plt.xlabel("Frequency")
plt.ylabel("Phase")
plt.grid(True)
plt.tight_layout()
plt.show()

1.1.2 确定系统的频率特性

  (1) H ( z ) = 2 z z − 0.6 H\\left( z \\right) = 2z \\over z - 0.6 H(z)=z0.62z

  系统函数具有一个位于原点的零点,以及位于0.6处的极点。当 z = e j ω z = e^j\\omega z=ejω 在单位圆上从 1 到 -1 移动时,对应的极点复矢量逐步增加,造成幅频特性下降。 因此对应的是低通滤波器特性。该系统的零点与极点都位于单位圆内,所以相片呈现为最小相位系统, 根据复矢量几何特性可以分析处相频特性先减小、后增加。

▲ 图1.1.7 离散时间系统的零极点分布

▲ 图1.1.8 系统的幅频特性和相频特性

from headm import *

def hs(w):
    z = exp(1j*w)

    return 2*z/(z-0.6)

w = linspace(0, 2*pi, 500)

hw = hs(w)
hwabs = abs(hw)
hwangle = angle(hw)

plt.subplot(2,1,1)
plt.plot(w, hwabs, label='Amplitude')

plt.xlabel("Frequency")
plt.ylabel("Amplitude")
plt.grid(True)
plt.tight_layout()

plt.subplot(2,1,2)
plt.plot(w, hwangle, label='Angle')

plt.xlabel("Frequency")
plt.ylabel("Phase")
plt.grid(True)
plt.tight_layout()
plt.show()

  (2) 系统具有一个二阶零点和一对共轭极点,分布如下图所示。根据极点所处在的位置, 在 ω = π / 2 \\omega = \\pi /2 ω=π/2 会存在峰值。 根据零点所处的位置, 在 ω = − π \\omega = - \\pi ω=π 幅值处在低谷。 所以系统的幅频特性整体上符合“带通特性”。

  该系统的零点位于单位圆外面, 所以系统不属于最小相位系统, 而是最大相位系统, 相角处于一直下降趋势。

H ( z ) = ( 0.96 + z − 1 ) 2 0.36 z − 2 + 1 H\\left( z \\right) = \\left( 0.96 + z^ - 1 \\right)^2 \\over 0.36z^ - 2 + 1 H(z)=0.36z2+1(0.96+z1)2

▲ 图1.1.9 离散时间系统的零极点分布

  下面是使用Python绘制的系统的幅频特性和相频特性。从图中可以看出,由于极点距离单位圆不是特别近,所以带通特性不明显。单位圆外的零点距离单位圆比较近,整体上,该系统的幅频特性更接近于低通滤波器特性。

  绘制的相频特性由于需要把相角限制在 ± π \\pm \\pi ±π 之内,所以中间具有两次跳变,但从变化趋势上来看, 相位始终是下降的。

▲ 图1.1.10 系统的幅频特性和相频特性

1.1.3 绘制幅频特性和相频特性

◎ 求解: 为了给出 H ( s ) H\\left( s \\right) H(s) 频率特性比较准确的绘制, 将题目中的零极点位置按照比例定义它们的取值,然后通过Python程序绘制出对应的幅频特性与相频特性。 坐标采用 x-轴对数,y-轴对数方式进行绘制,也就是绘制出系统的波特图。

  (a) 这是一个单个极点的 H ( s ) H\\left( s \\right) H(s) ,假设 H ( s ) = 1 s + 1 H\\left( s \\right) = 1 \\over s + 1 H以上是关于信号与系统分析2022春季作业-参考答案:第十三次作业的主要内容,如果未能解决你的问题,请参考以下文章

2021年春季学期-信号与系统-第十三次作业参考答案-第十小题

2021年春季学期-信号与系统-第十三次作业参考答案-第十一小题

2021年春季学期-信号与系统-第十三次作业参考答案-第七小题

2021年春季学期-信号与系统-第十三次作业参考答案-第四小题

2021年春季学期-信号与系统-第十三次作业参考答案-第五小题

2021年春季学期-信号与系统-第十三次作业参考答案-第八小题