周期矩形脉冲信号的频谱python

周期矩形脉冲信号的频谱可以用傅里叶变换求得。以下是Python代码实现：

import numpy as np
import matplotlib.pyplot as plt

# 周期矩形脉冲信号的参数
T = 1    # 周期
D = 0.5  # 占空比

# 构造周期矩形脉冲信号
t = np.linspace(0, 5*T, 1000)
x = np.zeros_like(t)
x[(t % T) < D*T] = 1

# 计算频谱
X = np.fft.fft(x)
freq = np.fft.fftfreq(len(x), t[1]-t[0])

# 绘制时域图像
plt.subplot(2, 1, 1)
plt.plot(t, x)
plt.xlabel('Time')
plt.ylabel('Amplitude')
plt.title('Periodic Rectangular Pulse Signal')

# 绘制频域图像
plt.subplot(2, 1, 2)
plt.plot(freq, abs(X))
plt.xlabel('Frequency')
plt.ylabel('Magnitude')
plt.title('Frequency Spectrum')
plt.tight_layout()
plt.show()

运行代码后，会输出周期矩形脉冲信号的时域图像和频域图像。其中，时域图像展示了周期矩形脉冲信号的波形，频域图像展示了该信号的频谱。