python生成LFM信号及其频谱

115 阅读 0 评论 76 点赞

我是靠谱客的博主多情夏天，这篇文章主要介绍python生成LFM信号及其频谱，现在分享给大家，希望可以做个参考。

生成时宽为10us的基带LFM信号
带宽为30MHz
第一部分为实信号第二部分为复信号
fft得到频谱，画图展示幅度谱

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
import numpy as np
import matplotlib.pyplot as plt

########################
T = 10e-6 
BW = 30e6
K = BW/T # chirp rate
Fs_LFM = 2*BW
Ts = 1/Fs_LFM
N = T/Ts
t = np.linspace(-T/2,T/2,N)
st_LFM = np.cos(np.pi*K*np.power(t,2))
plt.plot(t*1e6,st_LFM)
plt.xlabel('Time(us)')
plt.ylabel('Amplitude(W)')
plt.title('LFM signal in time domain')
plt.show()
freq = np.linspace(-Fs_LFM/2,Fs_LFM/2,N)
sf_LFM = np.fft.fftshift(np.fft.fft(st_LFM))
plt.plot(freq*1e-6,abs(sf_LFM.real))
plt.xlabel('Frequency(MHz)')
plt.ylabel('Amplitude(W)')
plt.title('Spectrum amplitude of real LFM signal')
plt.show()

########################
st_complex = np.exp(1j*np.pi*K*np.power(t,2))
sf_complex = np.fft.fftshift(np.fft.fft(st_complex))
plt.plot(freq*1e-6,abs(sf_complex))
plt.xlabel('Frequency(MHz)')
plt.ylabel('Amplitude(W)')
plt.title('Spectrum amplitude of complex LFM signal')
plt.show()