Python 频谱分析仪‘最大保持’功能
Python 中“Max Hold”的实现
如何在Moku的频谱分析仪上使用 Python 创建“最大保持”
Moku的频谱分析仪在 iPad 和 Windows 应用程序上都包含一个数学通道。此数学通道具有实用的“最大保持”功能。
附加的 Python 脚本部署了频谱分析仪,扫描频率并将数据捕获到 Python 中。然后,我们在 Python 中计算并显示“最大保持”轨迹。
#
# Moku example: Plotting Spectrum Analyzer - Max Hold
#
# This example demonstrates how you can configure the Spectrum Analyzer
# instrument and plot ad Max Hold of its spectrum data in real-time.
#
# (c) 2023 Liquid Pty. Ltd.
#
from moku.instruments import SpectrumAnalyzer
import matplotlib.pyplot as plt
# Launch Spectrum Analyzer and connect to your device via IP
i = SpectrumAnalyzer('192.168.###.###', force_connect=True)
try:
# Set spectrum analyzer configuration
i.set_defaults()
i.set_span(70e6, 130e6)
i.set_rbw('Auto') # Auto-mode
# Configure ADC inputs
i.set_frontend(1, "50Ohm", "AC", "400mVpp")
i.set_frontend(2, "50Ohm", "AC", "400mVpp")
# Set up basic plot configurations
line1, = plt.plot([])
line2, = plt.plot([])
plt.ion()
plt.show()
plt.grid()
plt.ylim([-200, 100])
plt.autoscale(axis='x', tight=True)
# Get an initial frame of data to set any frame-specific plot parameters
frame = i.get_data()
# Format the x-axis as a frequency scale
ax = plt.gca()
counter = 0
# Get and update the plot with new data
while True:
frame = i.get_data()
plt.pause(0.0001)
fc1 = frame['ch1']
spectrum_length = len(fc1)
# Set the frame data for each channel plot
line1.set_ydata(fc1)
if counter == 0:
current_max = fc1
else:
for x in range(0,spectrum_length):
if fc1[x]>current_max[x]:
current_max[x] = fc1[x]
line2.set_ydata(current_max)
# Frequency axis shouldn't change, but to be sure
line1.set_xdata(frame['frequency'])
line2.set_xdata(frame['frequency'])
# Ensure the frequency axis is a tight fit
ax.relim()
ax.autoscale_view()
# Redraw the lines
plt.draw()
counter = counter + 1
except Exception as e:
print(f'Exception occurred: {e}')
finally:
# Close the connection to the Moku device
# This ensures network resources are released correctly
i.relinquish_ownership()