间歇采样直接转发干扰仿真

# Python 雷达对抗仿真示例
# 实验: 间歇采样直接转发干扰仿真
import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# 信号参数设置
fs = 10e6          # 采样率 10 MHz
T = 100e-6         # 脉宽 100 us
B = 5e6            # 带宽 5 MHz

# 生成 LFM 信号
t = np.arange(0, T, 1/fs)
K = B / T
s = np.exp(1j * np.pi * K * t**2)

# 时域波形
plt.figure(figsize=(10, 4))
plt.plot(t * 1e6, np.real(s))
plt.xlabel('时间 (us)')
plt.ylabel('幅度')
plt.title('LFM 信号时域波形（实部）')
plt.grid(True)
plt.tight_layout()
plt.savefig('figure_1.png', dpi=100)
plt.close()

# 频谱分析
N = len(s)
f = np.fft.fftshift(np.fft.fftfreq(N, 1/fs))
S = np.fft.fftshift(np.fft.fft(s))
plt.figure(figsize=(10, 4))
plt.plot(f / 1e6, np.abs(S))
plt.xlabel('频率 (MHz)')
plt.ylabel('幅度')
plt.title('LFM 信号频谱')
plt.grid(True)
plt.tight_layout()
plt.savefig('figure_2.png', dpi=100)
plt.close()
print('仿真完成。')