线性调频信号(Chirp信号)公式

Chirp信号的公式

对于固定频率$f_1$的信号，表达式为：
$$r(t)=cos(2\pi \cdot f_1\cdot t)\text{\hspace{1em}(1)}$$
其中余弦信号的相位是频率对时间的积分：
$$\theta (t)=2\pi {\int }_0^t{f}_1\cdot \tau d\tau =2\pi f_{1}t\text{\hspace{1em}(2)}$$
Chirp信号的频率是随时间变化的，最简单的一种变化就是线性变化，从$t=t_0$时刻的频率$f_{start}$开始，变化到$t=t_1$时刻的频率$f_{end}$，那么任意时刻的信号频率可以表示为：
$$f(t)=\frac{f_{end}-f_{start}}{t_1-t_0}\times t+f_{start}\text{\hspace{1em}(3)}$$
对应得余弦信号的相位变为：
$$\begin{gathered} \theta (t)=2\pi {\int }_{t_0}^{t}f(\tau )d\tau =2\pi {\int }_{t_0}^t\frac{f_{end}-f_{start}}{t_1-t_0}\times \tau +f_{start}d\tau \\ =2\pi [\frac{f_{end}-f_{start}}{t_1-t_0}\cdot \frac{1}{2}(t-t_0{)}^2+f_{start}(t-t_0)]\text{\hspace{1em}(4)} \end{gathered}$$
令$t_0=0$，则上式可以化简为：
$$\theta (t)=2\pi [\frac{f_{end}-f_{start}}{t_1-t_0}\cdot \frac{1}{2}{t}^2+f_{start}\cdot t]\text{\hspace{1em}(5)}$$
对应的Chirp信号可以表示为：
$$r(t)=cos[2\pi (\frac{f_{end}-f_{start}}{t_1-t_0}\cdot \frac{1}{2}{t}^2+f_{start}\cdot t)]\text{\hspace{1em}(6)}$$
令起始时间$t_0=0$，结束时间$t_1=0.2048s$，起始频率$f_{start}=250Hz$，结束频率$f_{end}=2000Hz$，调频斜率$b=\frac{f_{end}-f_{start}}{t_1-t_0}=8545Hz/s$。
假设采样时间$t_s=1/f_s=100\mu s$,得到的离散化波形如下：

import matplotlib.pyplot as plt
import numpy as np
def chirpSignal(start_time,start_fre,end_time,end_fre,t):
    return np.cos(2*np.pi*((end_fre-start_fre)/(end_time-start_time)*0.5*(t*t-start_time*start_time)+start_fre*(t-start_time)))
t_start = 0
t_end = 0.2048
fre_start = 250
fre_end = 2000
t_s = 0.0001
y=[]
t_list = np.arange(t_start, t_end, t_s)
for t in t_list:
    y.append(chirpSignal(t_start,fre_start,t_end, fre_end, t))
plt.plot(t_list,y)

import scipy.signal
data = scipy.signal.chirp(t_list, fre_start,t_end, fre_end)
plt.plot(t_list, data)