驻定相位原理（POSP）以及线性调频信号的频谱

最近看论文遇到了驻定相位原理，问老师直接给了我一本书让我看，看半天只有一段…不是这个方向的，半路出家做毕业设计需要用到这个定理，有错误的话请不吝赐教。

一、驻定相位原理

在数字信号处理中，经常需要将一个时域信号转换到频域进行处理，比如说需要对其进行傅里叶变换
那么对于这样的一个复数信号：
$$u\left(t\right)=a\left(t\right)e^{j\phi \left(t\right)}$$
对其进行傅里叶变换时，就会得到如下式子：
$$S(\omega )={\int }_{-\infty }^{\infty }a(t)e^{j\phi (t)}{e}^{-j\omega t}dt$$
这个积分事实上是很难计算的，并且通常没有解析解。我们现在通常是对时域信号以高于奈奎斯特频率的频率对其进行采样，再进行FFT，利用计算机求得数值解，但在没有计算机的年代，人们只能通过其他方法在特定前提下对其进行估算，比如stationary phase approximation.

This method originates from the 19th century, and is due to George Gabriel Stokes and Lord Kelvin. It is closely related to Laplace’s method and the method of steepest descent, but Laplace’s contribution precedes the others.

对于以下积分式：
$$P=\int U(t)cosV(t)dt=\text{Re}\left[\int U(t)e^{jV(t)}dt\right]$$

给定如下前提，相位$V(t)$对于时间来说是捷变的，幅度$U(t)$是缓变的。对于线性调频信号来说，即
$$k\to \infty$$

也就是说，当相位随时间变化很快的时候，那么在相位变化一个周期的时间内，幅度可视为一个常数，我们对其积分就近似为零（正负抵消），于是上述积分的值就基本上由V(t)变化缓慢的点决定，即critical points(相位梯度为零的点)。这就是驻定相位原理的基本思想。

二、利用驻定相位原理计算一般调频信号的频谱

设窄带信号的复振幅为
$$\mu (t)=a(t)e^{j\phi (t)}$$

其傅立叶变换
$$S(\omega )={\int }_{-\infty }^{\infty }a(t)e^{j\phi (t)}{e}^{-j\omega t}dt$$

根据驻定相位原理，该积分在
$$\frac{d}{dt}[\omega t-\phi (t)]=0$$

处才显著不为0，用 $t_k$ 表示驻定相位点（一阶导数为0），则
$$\omega ={\varphi }^{\prime }(t_k)(\ast )$$

对相位项 $\omega t-\varphi (t)$在驻定相位点附近进行泰勒展开
$$\omega t-\phi (t)=\omega t_k-\phi (t_k)+\left[\omega -{\phi }^{\prime }(t_k)\right](t-t_k)-\frac{{\phi }^{\prime \prime }(t_k)}{2}(t-t_k{)}^2+\dots$$

忽略二阶以上的高次项（当t远离$t_k$时,根据驻定相位原理，积分接近于0，故t必定在$t_k$附近取值，即$t-t_k$为小量），并结合（*）式消去一阶项得到：
$$\omega t-\phi (t)\approx \omega t_k-\phi (t_k)-\frac{{\phi }^{\prime \prime }(t_k)}{2}(t-t_k{)}^2$$

代入信号的频域表达式
S ( ω ) = a ( t k ) exp { − j [ ω t k − φ ( t k ) ] } ∫ t k − δ t k + δ exp [ j φ ′ ′ ( t k ) 2 ( t − t k ) 2 ] d t S(\omega)=a(t_k)\text{exp}\{-j[\omega t_k-\varphi (t_k)]\}\int^{t_k+\delta}_{t_k-\delta}\text{exp}[j\frac{\varphi''(t_k)}{2}(t-t_k)^2]dt S(ω)=a(tk​)exp{−j[ωtk​−φ(tk​)]}∫tk​−δtk​+δ​exp[j2φ′′(tk​)​(t−tk​)2]dt

作变量代换
$$t-t_k=u及\frac{{\phi }^{\prime \prime }(t_k)}{2}{u}^2=\frac{\pi y^2}{2}$$

则
$$du=\sqrt{\pi }[{\phi }^{\prime \prime }(t_k){]}^{-1/2}dy$$

代入信号的频域表达式
S ( ω ) = 2 π a ( t k ) φ ′ ′ ( t k ) exp { − j [ ω t k − φ ( t k ) ] } ∫ 0 φ ′ ′ ( t k ) π δ exp ( j π y 2 2 ) d y S(\omega)=2\sqrt{\pi}\frac{a(t_k)}{\sqrt{\varphi''(t_k)}}\text{exp}\{-j[\omega t_k-\varphi (t_k)]\}\int^{\sqrt{\frac{\varphi''(t_k)}{\pi}}\delta}_{0}\text{exp}(j\frac{\pi y^2}{2})dy S(ω)=2π ​φ′′(tk​) ​a(tk​)​exp{−j[ωtk​−φ(tk​)]}∫0πφ′′(tk​)​ ​δ​exp(j2πy2​)dy

式中积分为菲涅尔积分(Fresnel Integrals可以去查菲涅尔积分表)。若积分上限较大，则菲涅尔积分趋于
$$\frac{1}{\sqrt{2}}\text{exp}(j\frac{\pi }{4})$$

则信号的频域表达式为
$$S(\omega )=\sqrt{2\pi }\frac{a(t_k)}{\sqrt{\left|{\phi }^{\prime \prime }(t_k)\right|}}\text{exp}\left[-j\left(\omega t_k-\phi (t_k)-\frac{\pi }{4}\right)\right]$$

$♣$ 驻定相位点不一定只有一个，所以上式是不是应该加个$\Sigma$?

三、特殊线性调频信号的频谱

所谓线性调频就是

设LFM信号为
$$u(t)=rect(\frac{t}{T_p})\text{exp}(j\pi k_r{t}^2)$$

其傅立叶变换为
$$\begin{gathered} U(\omega )={\int }_{-\infty }^{\infty }u(t)\text{exp}(-j\omega t)dt={\int }_{-T_p/2}^{T_p/2}\text{exp}\left[j\left(\pi k_r{t}^2-\omega t\right)\right]dt \\ =\text{exp}\left(-j\frac{{\omega }^2}{4\pi k_r}\right){\int }_{-T_p/2}^{T_p/2}\text{exp}\left[j{\left(\sqrt{\pi k_r}t-\frac{\omega }{2\sqrt{\pi k_r}}\right)}^2\right]dt \end{gathered}$$

作变量代换
$$\begin{gathered} {\left(\sqrt{\pi k_r}t-\frac{\omega }{2\sqrt{\pi k_r}}\right)}^2=\frac{\pi }{2}{x}^2\Rightarrow x=\sqrt{2k_r}t-\frac{\omega }{\pi \sqrt{2k_r}} \\ dx=\sqrt{2k_r}dt \end{gathered}$$

则
$$U(\omega )=\frac{1}{\sqrt{2k_r}}\text{exp}\left(-j\frac{{\omega }^2}{4\pi k_r}\right){\int }_{-X_1}^{X_2}\text{exp}\left(j\frac{\pi x^2}{2}\right)dx$$

其中
$$X_1=\sqrt{2k_r}\frac{T_p}{2}+\frac{\omega }{\pi \sqrt{2k_r}}=\frac{\pi k_r{T}_p+\omega }{\pi \sqrt{2k_r}}$$
$$X_2=\sqrt{2k_r}\frac{T_p}{2}-\frac{\omega }{\pi \sqrt{2k_r}}=\frac{\pi k_r{T}_p-\omega }{\pi \sqrt{2k_r}}$$

菲涅尔积分
$$C(X)={\int }_0^X\cos \left(\frac{\pi x^2}{2}\right)dx$$
$$S(X)={\int }_0^X\sin \left(\frac{\pi x^2}{2}\right)dx$$
$$C(-X)=-C(X),S(-X)=-S(X)$$
$$\underset{X\to \infty }{\lim }C(X)=\underset{X\to \infty }{\lim }S(X)=0.5$$

则矩形包络LFM信号的频谱精确解为
$$U(\omega )=\frac{1}{\sqrt{2k_r}}\text{exp}\left(-j\frac{{\omega }^2}{4\pi k_r}\right)\left[C(X_1)+jS(X_1)+C(X_2)+jS(X_2)\right]$$

其幅度谱和相位谱分别为
$$\left|U(\omega )\right|=\frac{1}{\sqrt{2k_r}}\sqrt{{\left[C(X_1)+C(X_2)\right]}^2+{\left[S(X_1)+S(X_2)\right]}^2}$$
$$\phi (\omega )=-\frac{{\omega }^2}{4\pi k_r}+\arctan \frac{S(X_1)+S(X_2)}{C(X_1)+C(X_2)}$$

在感兴趣的频率范围内，分式
$$\frac{S(X_1)+S(X_2)}{C(X_1)+C(X_2)}\approx 1$$

当时宽带宽积足够大时，有
$$\left|U(\omega )\right|=\frac{1}{\sqrt{k_r}}$$
$$\phi (\omega )=-\frac{{\omega }^2}{4\pi k_r}+\frac{\pi }{4}$$

窄带信号 $a(t)e^{\left[j\varphi (t)\right]}$ 的频谱为
$$S(\omega )=\sqrt{2\pi }\frac{a(t_k)}{\sqrt{\left|{\phi }^{\prime \prime }(t_k)\right|}}\text{exp}\left[-j\left(\omega t_k-\phi (t_k)-\frac{\pi }{4}\right)\right]$$

对矩形包络的LFM信号（把下面这些特殊信号的取值带到上面的一般结论中）
$$\begin{gathered} a(t)=rect(\frac{t}{T_p}) \\ \phi (t)=\pi k_r{t}^2\Rightarrow {\phi }^{\prime }(t)=2\pi k_{r}t\Rightarrow {\phi }^{\prime \prime }(t)=2\pi k_r \\ \frac{d}{dt}\left[\phi (t)-\omega t\right]=0\Rightarrow 2\pi k_{r}t-\omega =0\Rightarrow t_k=\frac{\omega }{2\pi k_r} \end{gathered}$$

利用驻定相位原理得到LFM信号的幅度和相位谱分别为
$$\left|S(\omega )\right|=\sqrt{2\pi }\frac{rect\left(\frac{\omega }{2\pi k_r{T}_p}\right)}{\sqrt{2\pi k_r}}=\frac{1}{\sqrt{k_r}}rect(\frac{\omega }{\Delta \omega })$$
$$\Phi (\omega )=-\omega t_k+\phi (t_k)+\frac{\pi }{4}=-\omega \frac{\omega }{2\pi k_r}+\pi k_r{\left(\frac{\omega }{2\pi k_r}\right)}^2+\frac{\pi }{4}=-\frac{{\omega }^2}{4\pi k_r}+\frac{\pi }{4}$$

四、快速POSP

  设$g(t)$是一个调频信号
$$g(t)=w(t)\mathrm{e}\mathrm{x}\mathrm{p}[j\varphi (t)]$$
其中$w(t)$是实包络，$\varphi (t)$为信号调制相位。假设与相位相比，包络为时间缓变函数。
  该信号的频谱为
$$\begin{array}{rl}G(f)& ={\int }_{-\infty }^{+\infty }g(t)\mathrm{e}\mathrm{x}\mathrm{p}[-j2\pi ft]dt\\ & ={\int }_{-\infty }^{+\infty }w(t)\mathrm{e}\mathrm{x}\mathrm{p}[j\varphi (t)-2\pi ft]dt\\ & ={\int }_{-\infty }^{+\infty }w(t)\mathrm{e}\mathrm{x}\mathrm{p}[j\theta (t)]dt\end{array}$$
最终，$G(f)$的频谱形式如下：
$$G(f)\approx C_{1}W(f)\mathrm{e}\mathrm{x}\mathrm{p}[j(\Theta (f)\pm \frac{\pi }{4})]$$
其中变量定义如下：
$♣$ $C_1$为一个通常可忽略的常数
$$C_1=\sqrt{\frac{2\pi }{|{\varphi }^{{}^{\prime \prime }}(t)|}}$$
$♣$ $W(f)$为频域包络
$$W(f)=w[t(f)]$$
$♣$ $\Theta (f)$为频域相位
$$\Theta (f)=\theta [t(f)]$$
$♣$ $t(f)$由信号时频关系给出
$$\frac{d\theta (t)}{dt}=0$$
$♣$ $\frac{\pi }{4}$的符号由${\varphi }^{{}^{\prime \prime }}(t)$的符号给出。与$C_1$类似，在绝大多数分析中都可忽略该常数相位的影响。

五、利用快速POSP求解LFM的频谱

LFM信号频谱为
$$G(f)={\int }_{-\infty }^{\infty }rect(\frac{t}{T})\text{exp}(j\pi K{t}^2)\text{exp}(-j2\pi ft)dt$$
其中包络$w(t)$为矩形，被积相位为
$$\theta (t)=\pi K{t}^2-2\pi ft$$
对其求导，并令为0，可得时频关系：
$$t=\frac{f}{K}$$
频域相位为
$$\Theta (f)=\theta (t=f/K)=\pi K(\frac{f}{K}{)}^2-2\pi f(\frac{f}{K})=-\pi \frac{f^2}{K}$$
频域包络为
$$W(f)=w(t=f/K)=\mathrm{r}\mathrm{e}\mathrm{c}\mathrm{t}(\frac{f}{|K|T})$$
忽略前面的系数和线性相位，可得频谱为：
$$G(f)=\mathrm{r}\mathrm{e}\mathrm{c}\mathrm{t}(\frac{f}{KT})\mathrm{e}\mathrm{x}\mathrm{p}(-j\pi \frac{f^2}{K})$$
可以看到与前面的结论是一致的。

附录（我的草稿）

参考文献

（1）https://wenku.baidu.com/view/0f5ef8f74693daef5ef73dfd.html
（2）https://en.wikipedia.org/wiki/Chirp_spectrum
（3）https://en.wikipedia.org/wiki/Stationary_phase_approximation
（4）合成孔径雷达成像算法与实现/IAN G.CUMMING