雷达脉冲信号的 时域分析的原始代码，输出PDW描述字

以下是雷达脉冲信号时域分析的Python代码，可以输出PDW描述字：

import numpy as np

# 定义雷达脉冲信号
fs = 1e9  # 采样率
T = 1/fs  # 采样时间
t = np.arange(0, 2e-6, T)  # 时域范围
fc = 10e6  # 载频频率
f0 = 1e6  # 调频斜率
B = 50e6  # 带宽
s = np.exp(1j * 2 * np.pi * (fc * t + 0.5 * f0 * t**2)) * np.cos(np.pi * B * t)**2

# 计算FFT
N = len(s)
freq = np.fft.fftfreq(N, T)
S = np.fft.fft(s)

# 定义PDW描述字格式
class PDW:
    def __init__(self, time, freq, amplitude):
        self.time = time
        self.freq = freq
        self.amplitude = amplitude

# 提取PDW描述字
pdw_list = []
for i in range(1, N-1):
    if np.abs(S[i]) > np.abs(S[i-1]) and np.abs(S[i]) > np.abs(S[i+1]):
        time = t[i]
        freq = freq[i]
        amplitude = np.abs(S[i])
        pdw = PDW(time, freq, amplitude)
        pdw_list.append(pdw)

# 输出PDW描述字
for pdw in pdw_list:
    print("time:{:.6f}s, freq:{:.2f}MHz, amplitude:{:.2f}".format(pdw.time, pdw.freq/1e6, pdw.amplitude))

该代码生成一个带有调频斜率的线性调频信号，然后计算其FFT，并提取PDW描述字。输出结果包括每个PDW的时间、频率和幅度。