瑞利信道的多普勒谱的原理与MATLAB仿真

一、问题

最近为了完成老师布置的作业，就通过博客这种方式来记录一下自己的解决问题的过程，问题如下：

(1) 无线信道中的多普勒谱有一种经典谱(classic spectrum)，请解释产生这种谱形状的机理；
(2) 请用Simulink或者m语言，产生一条单径瑞利信道，其多普勒谱为经典谱，其中移动速率为120km/h。并请描述：自己用计算机产生的经典谱是正确的依据。

       通过阅读与分析知网以及CSDN的相关文献后，对于多普勒经典谱的产生原理，我利用Clarke理论模型，来推导出瑞利信道的经典谱的功率谱密度，解决第一问，下面是Clarke理论模型的推导过程。

二、Clarke模型

1）计算平面波的通频带发射信号$\tilde{x}\left(t\right)$

       Clarke模型假设有N个具有任意相位的平面波，每个平面波以任意的方向到达移动台，且所有平面波的平均功率相同。设移动台的速度为$v$，平面波沿$x-y$ 平面上的水平方向到达，平面波与移动台的运动方向之间的夹角即到达角度（Angle of Arrival，AoA）为$\theta$，如图2.1所示。

图2.1: 移动台与平面波的到达角度的示意图

其中设平面波的基带发射信号为$x\left(t\right)$，载波的波长为$\lambda$，载波频率为$f_c$，则通频带的发射信号$\tilde{x}\left(t\right)$为
$$\tilde{x}\left(t\right)=\text{Re}\left[x\left(t\right)e^{j2\pi f_{c}t}\right]\text{\hspace{1em}(2.1)}$$

2）计算移动台的基带接收信号$y\left(t\right)$

       由于移动台的相对运动，到达移动台的平面波都会经历多普勒频移，且每条传播路径的多普勒频移都不同。设移动台位于有$I$条传播路径的散射信道，则通频带接收信号$\tilde{y}\left(t\right)$可表示为
$$\begin{gathered} \tilde{y}\left(t\right)=\text{Re}\left[\sum _{i=1}^I{C}_i{e}^{j2\pi \left(f_c+f_i\right)\left(t-{\tau }_i\right)}x\left(t-{\tau }_i\right)\right] \\ =\text{Re}\left[y\left(t\right)e^{j2\pi f_{c}t}\right]\text{\hspace{1em}(2.2)} \end{gathered}$$

其中$C_i,f_i,{\tau }_i$分别表示第$i$条传播路径的信道增益、多普勒频移、时延。
       而移动台的基带接收信号$y\left(t\right)$可表示为
$$y\left(t\right)=\sum _{i=1}^I{C}_i{e}^{-j{\varphi }_i\left(t\right)}x\left(t-{\tau }_i\right)\text{\hspace{1em}(2.3)}$$

其中 ϕ i ( t ) = 2 π { ( f c + f i ) τ i − f i t i } \phi _i\left( t \right) =2\pi \left\{ \left( f_c+f_i \right) \tau _i-f_it_i \right\} ϕi​(t)=2π{(fc​+fi​)τi​−fi​ti​}。
       第$i$条传播路径的多普勒频移可表示为$f_i=f_m\cos {\theta }_i=\frac{v}{\lambda }\cos {\theta }_i$，其中$f_m,{\theta }_i$分别表示最大多普勒频移、第$i$条传播路径的到达角度。第$i$条传播路径的到达角度服从$\left[-\pi ,\pi \right]$的均匀分布。

3）利用脉冲响应表示通频带接收信号$\tilde{y}\left(t\right)$

        基带信道可建模为线性时变滤波器，当路径的时延差远小于采样周期$T_s$时，路径时延${\tau }_i$
可近似为$\hat{\tau }$，则该滤波器的复基带脉冲响应为
$$h\left(t,\tau \right)=\sum _{i=1}^I{C}_i{e}^{-j{\varphi }_i\left(t\right)}\delta \left(t-{\tau }_i\right)\approx h\left(t\right)\delta \left(t-\hat{\tau }\right)\text{\hspace{1em}(2.4)}$$

其中$h\left(t\right)={\sum }_{i=1}^I{C}_i{e}^{-j{\varphi }_i\left(t\right)}$。
        假设基带发射信号$x\left(t\right)=1$，将脉冲响应代入式(2.2)，则通频带接收信号可进一步表示为
y ~ ( t ) = Re [ y ( t ) e j 2 π f c t ] = Re [ { h I ( t ) + j h Q ( t ) } e j 2 π f c t ] = h I ( t ) cos ⁡ 2 π f c t − h Q ( t ) sin ⁡ 2 π f c t (2.5) \tilde{y}\left( t \right) =\text{Re}\left[ y\left( t \right) e^{j2\pi f_ct} \right]=\text{Re}\left[ \left\{ h_I\left( t \right) +jh_Q\left( t \right) \right\} e^{j2\pi f_ct} \right]\\=h_I\left( t \right) \cos 2\pi f_ct-h_Q\left( t \right) \sin 2\pi f_ct \tag{2.5} y~​(t)=Re[y(t)ej2πfc​t]=Re[{hI​(t)+jhQ​(t)}ej2πfc​t]=hI​(t)cos2πfc​t−hQ​(t)sin2πfc​t(2.5)

其中$h_I\left(t\right)$和$h_Q\left(t\right)$分别为$h\left(t\right)$的同相分量和正交分量，可分别表示为
$$h_I\left(t\right)=\sum _{i=1}^I{C}_i\cos {\varphi }_i\left(t\right)\text{\hspace{1em}(2.6)}$$

$$h_Q\left(t\right)=-\sum _{i=1}^I{C}_i\sin {\varphi }_i\left(t\right)\text{\hspace{1em}(2.7)}$$
由中心极限定理可得，当传播路径数$I$足够大时，同相分量和正交分量可近似为高斯随机变量，因此接收信号的幅度$\left|\tilde{y}\left(t\right)\right|=\sqrt{h_I^2\left(t\right)+h_Q^2\left(t\right)}$服从瑞利分布。

4）计算瑞利信道的功率谱密度

       由于自相关函数和功率谱密度是一对傅里叶变换，对通频带接收信号$\tilde{y}\left(t\right)$的自相关函数做傅里叶变换，由此得到瑞利衰落的功率谱密度（PSD） 为
$$S_{\bar{y}\bar{y}}\left(t\right)=\{\begin{array}{l}\frac{{\Omega }_p}{4\pi f_m}\frac{1}{\sqrt{1-{\left(\frac{f-f_c}{f_m}\right)}^2}},\left|f-f_c\right|\le f_m\\ 0,\text{其他}\end{array}\text{\hspace{1em}(2.8)}$$

其中${\Omega }_p=E\left\{h_I^2\left(t\right)\right\}+E\left\{h_Q^2\left(t\right)\right\}={\sum }_{i=1}^I{C}_i^2$，式(2.8)即为经典多普勒谱的功率谱密度。

三、Jakes模型

       根据Clarke理论模型对瑞利衰落的推导，可得接收信号的幅度服从瑞利分布，且其同相分量与正交分量都是服从高斯分布的随机变量，而相位是服从$\left[-\pi ,\pi \right]$的均匀分布。
       因此Jakes模型采用复正弦波合成的方法，根据中心极限定理，当正弦波的数量足够多时，经过叠加后得到的接收信号的幅度近似服从瑞利分布，由此产生了具有经典谱的单径瑞利信道，在物理上近似实现了Clarke理论模型。下面是Jakes模型的具体步骤。

1）设置经过多普勒频移的平面波

       设有N个平面波，定义$N_0=\left(\frac{N}{2}-1\right)/2$，其中限定$N/2$为一个奇数。设经过最大多普勒频移$f_m$的平面波的频率为${\omega }_d=2\pi f_m$，初始相位为${\varphi }_N$。
       设$N_0$个经过多普勒频移的平面波，其中每一个平面波的到达角度为${\theta }_n=\frac{2\pi n}{N}$，频率为${\omega }_n={\omega }_d\cos {\theta }_n$，初始相位为${\varphi }_n,n=1,2,\cdots ,N_0$。为了使瑞利衰落的相位服从均匀分布，初始相位可设置为
$$\begin{gathered} {\varphi }_N=0 \\ {\varphi }_n=\frac{\pi n}{N_0+1},n=1,2,\cdots ,N_0\text{\hspace{1em}(2.9)} \end{gathered}$$

2）合成复正弦波
       将上述的$N_0$个经过多普勒频移的平面波的复振荡器的输出求和，然后与经过最大多普勒频移$f_m$
的平面波的复振荡器的输出相加，如图3.1所示。在复振荡器的总和中，实部$h_I\left(t\right)$和虚部$h_Q\left(t\right)$分别表示为
$$\begin{gathered} h_I\left(t\right)=2\sum _{n=1}^{N_0}\left(\cos {\varphi }_n\cos {\omega }_{n}t\right)+\sqrt{2}\cos {\varphi }_N\cos {\omega }_{d}t \\ =\left[\cos {\varphi }_1\cos {\varphi }_n\cdots \cos {\varphi }_{N_0}\cos {\varphi }_N\right]\cdot \left[\begin{array}{l}2\cos {\omega }_{1}t\\ 2\cos {\omega }_{2}t\\ ⋮\\ 2\cos {\omega }_{N_0}t\\ \sqrt{2}\cos {\omega }_{N}t\end{array}\right]\text{\hspace{1em}(2.10)} \end{gathered}$$
$$\begin{gathered} h_Q\left(t\right)=2\sum _{n=1}^{N_0}\left(\sin {\varphi }_n\cos {\omega }_{n}t\right)+\sqrt{2}\sin {\varphi }_N\cos {\omega }_{d}t \\ =\left[\sin {\varphi }_1\sin {\varphi }_n\cdots \sin {\varphi }_{N_0}\sin {\varphi }_N\right]\cdot \left[\begin{array}{l}2\cos {\omega }_{1}t\\ 2\cos {\omega }_{2}t\\ ⋮\\ 2\cos {\omega }_{N_0}t\\ \sqrt{2}\cos {\omega }_{N}t\end{array}\right]\text{\hspace{1em}(2.11)} \end{gathered}$$

图3.1: Jakes仿真模型的复正弦波合成过程
因此Jakes仿真模型的复输出可以表示为
h ( t ) = E 0 2 N 0 + 1 { h I ( t ) + j h Q ( t ) } (2.12) h\left( t \right) =\frac{E_0}{\sqrt{2N_0+1}}\left\{ h_I\left( t \right) +jh_Q\left( t \right) \right\} \tag{2.12} h(t)=2N0​+1 ​E0​​{hI​(t)+jhQ​(t)}(2.12)

其中$E_0$为衰落信号的平均幅度。由中心极限定理可得，为了使经过叠加后得到的接收信号的幅度近似服从瑞利分布，当经过多普勒频移的平面波的数量必须足够多，一般取$N_0=8$。

四、Jakes模型的MATLAB仿真

       为了完成第二问的仿真任务，我通过MATLAB仿真软件对Jakes模型进行了仿真。根据题设的移动速率$v=120km/h$的条件，由最大多普勒频移$f_m$公式可得
$$f_m=\frac{v}{\lambda }=\frac{v{f}_c}{c}\text{\hspace{1em}(2.13)}$$

其中$c=3\times 10^{8}m/s$为光速，$f_c$为载波频率。对于载波频率$f_c$我设置为2GHz，代入式(2.13)可得最大多普勒频移$f_m=222.2Hz$，将数据代入Jakes仿真模型中，可得MATLAB输出的Jakes模型的幅度在时域的特点、幅度与相位的分布、经典模型与仿真模型的自相关函数、多普勒谱如下图所示。

图4.1: Jakes模型的幅度在时域的特点        由时域图可得，Jakes模型的幅度随着时间的增加会发生衰减。

图4.2: 幅度与相位的分布        由分布图可得，在Jakes模型中接收信号的幅度基本服从瑞利分布，相位则近似服从均匀分布。

图4.3: 经典模型与仿真模型的自相关函数、多普勒谱

       这里的经典模型采用第一类0阶贝塞尔函数的方式得到经典模型的自相关函数，然后对自相关函数做傅里叶变换得到经典模型的功率谱密度，从MATLAB仿真出的经典谱与经典模型的经典谱的图形对比以及自相关函数的对比图可以看出，计算机仿真产生的经典谱与理论推导的经典谱差别不大，仿真与经典的自相关函数仍差距不大，同时我使用均方根误差RMSE进一步对两个图形的多普勒谱的幅度数据进行了对比，RMSE的公式如下：
$$RMSE=\sqrt{\frac{1}{N}\sum _{t=1}^N{\left(observe{d}_t-predicte{d}_t\right)}^2}$$

经过MATLAB仿真后输出的RMSE约为3.3902，可得仿真数据的精确度较高。

五、Jakes模型的Simulink仿真

       为了更进一步地仿真单径瑞利信道，本次我也采用了Simulink做出了如下图所示的仿真框图。

图5.1: 单径瑞利信道的仿真框图

为了更快地仿真完，设置采样频率为100kHz，仿真时间为20s。然后设置最大多普勒频移$f_m=222.2Hz$。注意在单径瑞利信道模块中开启多普勒频谱的可视化，最后连接好频谱分析仪，运行Simulink后，得到如下结果。

图5.2: 实际样本与理论多普勒谱的对比图

由对比图可得，实际测量的样本接近于理论多普勒谱（经典谱）。然后观察频谱分析仪得到的结果。

图5.3: 频谱分析仪得到的频谱图

由频谱图可得，两个峰值对应的频率约为205Hz、-205Hz，其绝对值接近于设置的最大多普勒频移$f_m=222.2Hz$。综上，可以判断出我用计算机产生的经典谱是正确的。

六、总结

       感谢大家能阅读到这里！这是我时隔一年多又一次开始动笔写技术博客了，哈哈感觉还是生疏了，对于Markdown的有些格式问题弄了挺久的，自己也查阅了很多不懂知识，参考了很多资料，感谢前辈们的帮助！自己还是写得不太好，还请读者们多多指教，不过真的感觉写博客还是很有趣的！如果这篇博客能够帮助到大家，我自己就很开心了！
       这次通过写老师布置的作业，学会了在网上查找资料，然后阅读分析，消化吸收相关的内容，基本上Clarke理论模型、Jakes仿真模型搞懂了，然后也温习了一遍MATLAB，写了写仿真代码，用轻薄本跑了跑程序，哈哈每跑一次风扇都要转起来！也熟悉了Simulink的使用，希望自己之后能继续加油。雄关漫道真如铁，而今迈步从头越！

这里是本文的参考文献：
[1]:matlab实现单径瑞利信道仿真中经典多普勒谱（Clarke模型、Jakes模型）
[2]:MIMO-OFDM无线通信技术及MATLAB实现
[3]:均方误差(MSE)和均方根误差(RMSE)和平均绝对误差(MAE)
[4]:瑞利信道建模 matlab程序原理到实现