我已经从包含每个频率值的幅度和相位值的CSV文件中构建了系统的波特图。现在我需要确定静态误差常数的值，增益和相位裕度以及截止频率。由于传递函数未知，因此必须通过图形分析获得这些值。
我从csv文件中得到的图表如下所示
x1c 0d1x的数据
这是我用来绘制波特图的代码：

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt

data = pd.read_csv('example_P3.csv')
data
magnitude = data['Mag']
phase = data['Phase']
w_rad = data['W']

w_hz = w_rad/(2*np.pi)

plt.figure(figsize=(10, 6))
plt.subplot(2, 1, 1)
plt.semilogx(w_hz, magnitude, 'b')
plt.title('Diagrama de Bode - Magnitud')
plt.ylabel('Magnitud (dB)')

plt.subplot(2, 1, 2)
plt.semilogx(w_hz, phase, 'r')
plt.title('Diagrama de Bode - Fase')
plt.xlabel('Frecuencia (Hz)')
plt.ylabel('Fase (grados)')

plt.tight_layout()
plt.show()

字符串
为了确定相位和增益裕量以及截止频率，我在代码中执行了以下几行：

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt

data = pd.read_csv('example_P3.csv')
magnitude = data['Mag']
phase = data['Phase']
w_rad = data['W']

w_hz = w_rad/(2*np.pi)

plt.figure(figsize=(10, 6))
plt.subplot(2, 1, 1)
plt.semilogx(w_hz, magnitude, 'b')
plt.title('Diagrama de Bode - Magnitud')
plt.ylabel('Magnitud (dB)')

plt.subplot(2, 1, 2)
plt.semilogx(w_hz, phase, 'r')
plt.title('Diagrama de Bode - Fase')
plt.xlabel('Frecuencia (Hz)')
plt.ylabel('Fase (grados)')

# Calculate gain margin and phase margin
gain_margin_freq = w_hz[np.argmax(magnitude > 0)]  # Frequency where magnitude crosses 0 dB
phase_margin_phase = phase[np.argmin(np.abs(magnitude - 0))]  # Phase at the frequency of gain margin

# Calculate cutoff frequency and bandwidth
cutoff_freq = w_hz[np.argmax(magnitude > -3)]  # Frequency where magnitude is -3 dB
bandwidth = w_hz[np.argmax(magnitude > -3)] - w_deg[np.argmax(magnitude == np.max(magnitude))]  # Width at -3 dB

# Plot vertical lines and text annotations for gain and phase margins, cutoff frequency, and bandwidth
plt.subplot(2, 1, 1)
plt.axvline(gain_margin_freq, color='g', linestyle='--', label=f'Gain Margin: {gain_margin_freq:.2f} Hz')
plt.axvline(cutoff_freq, color='y', linestyle='--', label=f'Cutoff Frequency: {cutoff_freq:.2f} Hz')
plt.legend()

plt.subplot(2, 1, 2)
plt.axvline(gain_margin_freq, color='g', linestyle='--', label=f'Gain Margin: {gain_margin_freq:.2f} Hz')
plt.axhline(phase_margin_phase, color='m', linestyle='--', label=f'Phase Margin: {phase_margin_phase:.2f} degrees')
plt.legend()

plt.tight_layout()
plt.show()

型
对于截止频率，我所做的是尝试定位幅度的最大值，并找到3 dB下降的位置，以在图上标记它。至于裕度，对于增益，我寻找幅度跨越0 dB的频率，对于相位，我寻找相位跨越-180度的相位值。对于截止频率，我所做的是尝试定位幅度的最大值，并找到3 dB下降的位置，以在图上标记它。至于裕度，对于增益，我寻找幅度穿过0 dB的频率和相位，我寻找相位与-180度交叉处的相位值。这种方法似乎不合适，因为，如下图所示，截止频率的位置不正确。

的
这是我第一次做这种类型的事情，所以我完全找不到解决问题的方法，这导致我向做过这种类型问题的人寻求帮助，他可以指导我解决它。