import adi

import adi是一个在Python中导入ADI公司的模块的命令。ADI(Analog Devices Inc.)是一家以模拟与数字混合信号处理技术为核心的半导体公司。通过使用import adi命令，可以在Python中使用ADI提供的模块来进行相关的开发工作。该模块提供了一系列功能，包括数据采集、信号处理、通信等。使用import adi命令后，可以根据具体的需要选择相应的ADI模块并进行调用。例如，import adi.ad9361可以导入ADI的AD9361射频前端模块，使得我们可以在Python中通过该模块来进行射频信号的采集与处理。

相关问题

在下面代码中加入接收端的收到的复信号的时域图，再加一个隔直流的代码，用python import time import adi import matplotlib.pyplot as plt import numpy as np from scipy import signal # Create radio sdr = adi.Pluto() # Configure properties sdr.rx_rf_bandwidth = 4000000 sdr.rx_lo = 1900000000 sdr.tx_lo = 2000000000 sdr.tx_cyclic_buffer = True sdr.tx_hardwaregain_chan0 = -30 sdr.gain_control_mode_chan0 = "slow_attack" # Read properties print("RX LO %s" % (sdr.rx_lo)) print(sdr.sample_rate) # Create a sinewave waveform fs = int(sdr.sample_rate) N = 1024 fc = int(3000000 / (fs / N)) * (fs / N) ts = 1 / float(fs) t = np.arange(0, N * ts, ts) i = np.cos(2 * np.pi * t * fc) * 2 ** 14 q = np.sin(2 * np.pi * t * fc) * 2 ** 14 i=i+2**15 q=q+2**15 iq = i + 1j * q # Send data sdr.tx(iq) # Collect data for r in range(20): x = sdr.rx() print(x) y=abs(x) print(y) print('------------------------') plt.figure(0) plt.plot(y) f, Pxx_den = signal.periodogram(x, fs) #plt.clf() # plt.figure(1) plt.semilogy(f, Pxx_den) plt.ylim([1e-7, 1e4]) plt.xlabel("frequency [Hz]") plt.ylabel("PSD [V**2/Hz]") plt.draw() plt.pause(0.05) time.sleep(0.1) plt.show()

以下是加入接收端的收到的复信号的时域图和隔直流的代码的修改后的完整代码：

import time
import adi
import matplotlib.pyplot as plt
import numpy as np
from scipy import signal

# Create radio
sdr = adi.Pluto()

# Configure properties
sdr.rx_rf_bandwidth = 4000000
sdr.rx_lo = 1900000000
sdr.tx_lo = 2000000000
sdr.tx_cyclic_buffer = True
sdr.tx_hardwaregain_chan0 = -30
sdr.gain_control_mode_chan0 = "slow_attack"

# Read properties
print("RX LO %s" % (sdr.rx_lo))
print(sdr.sample_rate)

# Create a sinewave waveform
fs = int(sdr.sample_rate)
N = 1024
fc = int(3000000 / (fs / N)) * (fs / N)
ts = 1 / float(fs)
t = np.arange(0, N * ts, ts)
i = np.cos(2 * np.pi * t * fc) * 2 ** 14
q = np.sin(2 * np.pi * t * fc) * 2 ** 14
i = i + 2 ** 15
q = q + 2 ** 15
iq = i + 1j * q

# Send data
sdr.tx(iq)

# Collect data
for r in range(20):
    x = sdr.rx()
    print(x)
    y = abs(x)
    y_dc = y - np.mean(y)  # remove DC component
    print(y_dc)
    print('------------------------')
    plt.figure(0)
    plt.subplot(211)
    plt.plot(y_dc)
    plt.xlabel("Time [samples]")
    plt.ylabel("Amplitude")
    plt.subplot(212)
    f, Pxx_den = signal.periodogram(x, fs)
    plt.semilogy(f, Pxx_den)
    plt.ylim([1e-7, 1e4])
    plt.xlabel("Frequency [Hz]")
    plt.ylabel("PSD [V**2/Hz]")
    plt.draw()
    plt.pause(0.05)
    time.sleep(0.1)

plt.show()

在这个代码中，我们加入了一个隔直流的代码，用于移除接收到的信号的DC分量。这个代码为：

y_dc = y - np.mean(y)  # remove DC component

我们在时域图上画出了处理后的信号，并在频域图上画出了信号的功率谱密度（PSD）。

注意：这个代码只能在连接了PlutoSDR设备并安装了必要的包后才能运行。

详细解释以下Python代码：import numpy as np import adi import matplotlib.pyplot as plt sample_rate = 1e6 # Hz center_freq = 915e6 # Hz num_samps = 100000 # number of samples per call to rx() sdr = adi.Pluto("ip:192.168.2.1") sdr.sample_rate = int(sample_rate) # Config Tx sdr.tx_rf_bandwidth = int(sample_rate) # filter cutoff, just set it to the same as sample rate sdr.tx_lo = int(center_freq) sdr.tx_hardwaregain_chan0 = -50 # Increase to increase tx power, valid range is -90 to 0 dB # Config Rx sdr.rx_lo = int(center_freq) sdr.rx_rf_bandwidth = int(sample_rate) sdr.rx_buffer_size = num_samps sdr.gain_control_mode_chan0 = 'manual' sdr.rx_hardwaregain_chan0 = 0.0 # dB, increase to increase the receive gain, but be careful not to saturate the ADC # Create transmit waveform (QPSK, 16 samples per symbol) num_symbols = 1000 x_int = np.random.randint(0, 4, num_symbols) # 0 to 3 x_degrees = x_int*360/4.0 + 45 # 45, 135, 225, 315 degrees x_radians = x_degrees*np.pi/180.0 # sin() and cos() takes in radians x_symbols = np.cos(x_radians) + 1j*np.sin(x_radians) # this produces our QPSK complex symbols samples = np.repeat(x_symbols, 16) # 16 samples per symbol (rectangular pulses) samples *= 2**14 # The PlutoSDR expects samples to be between -2^14 and +2^14, not -1 and +1 like some SDRs # Start the transmitter sdr.tx_cyclic_buffer = True # Enable cyclic buffers sdr.tx(samples) # start transmitting # Clear buffer just to be safe for i in range (0, 10): raw_data = sdr.rx() # Receive samples rx_samples = sdr.rx() print(rx_samples) # Stop transmitting sdr.tx_destroy_buffer() # Calculate power spectral density (frequency domain version of signal) psd = np.abs(np.fft.fftshift(np.fft.fft(rx_samples)))**2 psd_dB = 10*np.log10(psd) f = np.linspace(sample_rate/-2, sample_rate/2, len(psd)) # Plot time domain plt.figure(0) plt.plot(np.real(rx_samples[::100])) plt.plot(np.imag(rx_samples[::100])) plt.xlabel("Time") # Plot freq domain plt.figure(1) plt.plot(f/1e6, psd_dB) plt.xlabel("Frequency [MHz]") plt.ylabel("PSD") plt.show()，并分析该代码中QPSK信号的功率谱密度图的特点

这段Python代码的作用是使用ADI Pluto SDR设备生成并传输一个QPSK信号，并将接收到的信号进行功率谱密度分析。下面是对代码的注释：

import numpy as np
import adi
import matplotlib.pyplot as plt

# 设置采样率、中心频率和采样点数
sample_rate = 1e6 # Hz
center_freq = 915e6 # Hz
num_samps = 100000 # number of samples per call to rx()

# 连接ADI Pluto SDR设备
sdr = adi.Pluto("ip:192.168.2.1")
sdr.sample_rate = int(sample_rate)

# 配置发送端的参数
sdr.tx_rf_bandwidth = int(sample_rate) # filter cutoff, just set it to the same as sample rate
sdr.tx_lo = int(center_freq)
sdr.tx_hardwaregain_chan0 = -50 # Increase to increase tx power, valid range is -90 to 0 dB

# 配置接收端的参数
sdr.rx_lo = int(center_freq)
sdr.rx_rf_bandwidth = int(sample_rate)
sdr.rx_buffer_size = num_samps
sdr.gain_control_mode_chan0 = 'manual'
sdr.rx_hardwaregain_chan0 = 0.0 # dB, increase to increase the receive gain, but be careful not to saturate the ADC

# 创建发送的QPSK信号
num_symbols = 1000
x_int = np.random.randint(0, 4, num_symbols) # 0 to 3
x_degrees = x_int*360/4.0 + 45 # 45, 135, 225, 315 degrees
x_radians = x_degrees*np.pi/180.0 # sin() and cos() takes in radians
x_symbols = np.cos(x_radians) + 1j*np.sin(x_radians) # this produces our QPSK complex symbols
samples = np.repeat(x_symbols, 16) # 16 samples per symbol (rectangular pulses)
samples *= 2**14 # The PlutoSDR expects samples to be between -2^14 and +2^14, not -1 and +1 like some SDRs

# 启动发送端并发送信号
sdr.tx_cyclic_buffer = True # Enable cyclic buffers
sdr.tx(samples) # start transmitting

# 接收接收端的信号
for i in range (0, 10):
    raw_data = sdr.rx()  # Receive samples
rx_samples = sdr.rx()
print(rx_samples)

# 停止发送端
sdr.tx_destroy_buffer()

# 计算接收到的信号的功率谱密度
psd = np.abs(np.fft.fftshift(np.fft.fft(rx_samples)))**2
psd_dB = 10*np.log10(psd)
f = np.linspace(sample_rate/-2, sample_rate/2, len(psd))

# 绘制时域图
plt.figure(0)
plt.plot(np.real(rx_samples[::100]))
plt.plot(np.imag(rx_samples[::100]))
plt.xlabel("Time")

# 绘制频域图
plt.figure(1)
plt.plot(f/1e6, psd_dB)
plt.xlabel("Frequency [MHz]")
plt.ylabel("PSD")
plt.show()

以上代码生成了一个随机QPSK信号，通过ADI Pluto SDR设备将其传输，并使用Pluto SDR设备接收该信号。接收到的信号进行了功率谱密度分析，并绘制了频域图。

QPSK信号的功率谱密度图的特点是，其频谱表现为四个簇，每个簇对应QPSK信号的一个符号。每个簇的带宽约为基带信号的带宽，且由于使用矩形脉冲，每个簇的带宽之间有一定的重叠。此外，功率谱密度图中还可以看到一些其他频率分量，这些分量可能是由于接收信号中存在其他干扰或噪声导致的。