3

我正在尝试使用 python 绘制奈奎斯特图。但是,我无法正确处理。

这是我正在使用的代码:

import control
import matplotlib.pyplot as plt

#Creating a transfer function G
s = control.TransferFunction.s

K0 = 10
T = 10
G = K0 / (s**2 * (T*s + 1))

control.nyquist(G)
plt.grid(True)
plt.title('Nyquist Diagram of G(s)')
plt.xlabel('Re(s)')
plt.ylabel('Im(s)')
plt.show()

代码输出此图:

奈奎斯特图输出图像

但情节应该是这样的(来自 wolfram):

wolfram 奈奎斯特图

为什么控制包中的奈奎斯特图是错误的?我怎样才能正确?

非常感谢。

4

1 回答 1

0

您必须通过 matplotlib 参数进行调整。您正在使用的 nyquist 函数将 **kwars 传递给 matplotlib。

你可以看到

control.nyquist

是一个别名

# Function aliases
bode = bode_plot
nyquist = nyquist_plot

有这个输入:

def nyquist_plot(syslist, omega=None, Plot=True, color=None,
             labelFreq=0, *args, **kwargs):
"""
Nyquist plot for a system

Plots a Nyquist plot for the system over a (optional) frequency range.

Parameters
----------
syslist : list of LTI
    List of linear input/output systems (single system is OK)
omega : freq_range
    Range of frequencies (list or bounds) in rad/sec
Plot : boolean
    If True, plot magnitude
color : string
    Used to specify the color of the plot
labelFreq : int
    Label every nth frequency on the plot
*args
    Additional arguments for :func:`matplotlib.plot` (color, linestyle, etc)
**kwargs:
    Additional keywords (passed to `matplotlib`)

Returns
-------
real : array
    real part of the frequency response array
imag : array
    imaginary part of the frequency response array
freq : array
    frequencies
于 2020-10-19T22:04:53.597 回答