时变多径信道的仿真

更精致的排版请查看blog

1. 仿真要求

在Matlab下进行编程，设计多径衰落信道仿真模型，具体要求如下：

• 可配置的参数包括多径数目、多径时延、各多径分量平均功率、多普勒频率扩展

• 可采用正弦波叠加或者频域滤波两种方式进行实现

• 多普勒功率谱为高斯功率谱

• 需要对信道仿真模型输出的信号的特性进行分析，绘制出其实际的信道响应特性

2. 原理分析

2.1. 多径衰落的概念及特性

多径传播是移动无线信道中信号产生衰落的主要原因。复杂的无线传播环境中，电波受到外 物的阻挡产生了反射、绕射和散射，当信号到达移动台时，已经不是原来的一路信号了，而是多个路径的多个信号的叠加。因为多路径传输造成了路径距离存在远近差别，所以发送的 一路信号在发散为多个路径后，各信号到达接收机的时间和相位都不一样。各信号按各自不 同相位叠加，相位相同的信号相加，信号增强，反之衰减，信号便产生失真。这样造成叠加 的接收信号的幅度和发送的单一信号幅度的快速巨大的差异，这称之为多径衰落。

多径衰落存在两种特性：一种是接收信号在幅度上的衰减，另一种是信号在信道上传输产生的时延扩展。具体地说，空间上，当接收端在运动时，随着距离的变化接收信号的幅度会发生衰落效应，其中接收信号的短时间的剧烈变化曲线是由于本地反射造成的，运动距离的变化对信号的平均值产生影响，运动距离的影响包括地形的变化和信号在空间传输产生的能量损耗。时间上，信号传输速度是一定的，但多径传播会导致距离的不同，由此带来信号传输时间的不一样，对同一接收端，接收同一信号在不同时间到达的合成信号，信号彼此之间存在到达时间差，即时延。这就使得接收信号相对发送信号在时间轴上产生信号宽度的增加称为时延扩展。一般来说，多径效应在信号幅度上引起的变化是模拟系统关注的重点；而它产生的时延扩展现象是数字系统考虑的重点。

2.2. 多普勒频移

当移动台以固定速率 $v$ 在长度为 $d$ 的路径上运动时，收到来自远端信号 $S$ 发出的信号，如图所示。无线电波从 $S$ 出发，在移动台运动的两个端点分别被移动台接收时所走的路径差

$$\Delta d=d\cos \theta =v\Delta t\cos \theta$$

这里$\Delta t$是移动台从左端点运动到右端点所用的时间，$\theta$是两个端点与入射波的夹角。由于源端距离很远，可以假设两端点处$\theta$是相同的。所以，由路程差造成的接收信号的相位变化值为

$$\Delta \varphi =\frac{2\pi \Delta d}{\lambda }=\frac{2\pi v\Delta t}{\lambda }\cos \theta$$

由此可得出频率变化值，即多普勒频移$f_d$为

$$f_d=\frac{1}{2\pi }\frac{\Delta \varphi }{\Delta t}=\frac{v}{\lambda }\cos \theta$$

2.3. 多径衰落信道的分类

信号通过无线信道时所产生的多径时延、多普勒频移及信号包络所服从的各种分布等，由此导致了信号通过无线信道时经历了不同类型的衰落。接收的多路信号时延差异带来信号间的干扰，由此产生了时间色散和频率点上信号的衰落即频率色散，这两种色散可能产生4种衰落效应。总体来说，产生衰落需要信号、信道以及接收端的移动，缺一不可。概括起来这4种衰落效应是：由时间色散导致的平坦衰落和不同频段上衰落变化的频率选择性衰落；由频率色散造成的信号在频点上幅值减小，导致的快衰落效应和慢衰落效应。

2.3.1. 平坦衰落和频率选择性衰落的判断条件

理论分析结果表明当信道的频带宽度大于原始发送信号的频带带宽，且在整个带宽范围内信号增益恒定不变和相位线性变化，则接收端的信号就是平坦衰落下的信号。此种衰落，信道的多径特性不会改变接收信号的频谱分布。平坦衰落信道的条件可概括为：

$$B_s\ll B_c,T_s\gg {\sigma }_{\tau }$$

其中$T_s$为信号周期(信号带宽$B_s$的倒数)；${\sigma }_{\tau }$为信道的时延扩展；$B_c$为相关带宽。

信道相当于一张门，它的带宽没有需要传送的原始信号的宽，就相当于一个高个子需要通过一张很矮的门，通过后的信号不再能保持信号原有的幅度值和相位值，他们的线性变化也无法保证，信号会表现出选择性衰落，即高个子想要通过需要弯腰才行。由信道的多径传播效应特性，信号通过信道必然存在延迟，延迟与带宽成反比，因此，我们发现在这种选择性衰落发生时信道时延会比信号周期长。从窄的信道即矮门中得到的波形是经历了衰减和时延的失真的信号。上述条件转化为公式如下：

$$B_s>B_c,T_s<{\sigma }_{\tau }$$
通常若$T_s\le {\sigma }_{\tau }$，可认为这样的信道会使信号产生频率选择性衰落，在实际应用时需要根据调制方式而定。

2.3.2. 快衰落信道和慢衰落信道的判断条件

当发送端发出的原始信号的周期长度与信道的相干时间相比要大，信道多普勒扩展$B_D$同样比原始发送信号的带宽$B_s$大，信道冲激响应在单个符号时间内变化迅速，接收端收到的信号会剧烈衰变，此衰落为快衰落。所以信号经历快衰落的条件是：

$$T_s>T_c,B_s<B_D$$

当发送端原始信号的符号周期长度比信道的相干时间短很多，且多普勒频谱扩展$B_D$比基带信号的带宽$B_s$短得多时，信道冲激信号的响应变化时间远没有传送的原始信号码元的周期时间长，这种条件下的信道属于慢衰落信道。由分析得出慢衰落的条件是：

$$T_s\ll T_c,B_s\gg B_D$$

3. 模型建立

3.1. 随机过程模型

为刻画多径衰落信道提出了各种各样的模型，它们几乎都使用了随机过程来描述衰落。描述多径的模型有两类，离散多径模型，(有限数量的多径分量)和散射多径模型(多径分量的连续体)。在移动无线通信中，第一类模型通常用于移动无线信道的波形级的仿真，而第二类模型通常用在窄带调制的对流层信道。在两种情况下，信道都被建模为具有复低通等效响应$\tilde{c}(\tau ,t)$的线性时变系统。如果有N个离散的多径分量存在，则信道的输出是N个被延迟和衰减的输入信号之和。因此

$$\tilde{y}(t)=\sum _{k=1}^{N(t)}{\tilde{a}}_k(t)\tilde{x}(t-{\tau }_k(t))$$

冲激响应$\tilde{c}(\tau ,t)$为

$$\tilde{c}(\tau ,t)=\sum _{k=1}^{N(t)}{\tilde{a}}_k(t)\tilde{x}(t-{\tau }_k(t))$$
其中N(t)为多径分量的数量，而$\widetilde{a_k}(t)$和${\tau }_k(t)$是在t时刻第k条路径的复衰减和延迟。

如前所述，多径信道可能是时不变的。然而，对所研究的所有实际信道而言，信道可以刻画为时变的(衰落)。引起时变的原因有两个：

1.即使发射机和接收机是固定的，
环境也在变化。例如电离层的变化、植物的摇摆以及反射体和散射体的移动。

2.即使环境可能是静态的，发射机和接收机也是移动的。因此，在实际的多径信道中，N、$a_k$和${\tau }_k$都可能随时间随机地变化，如图所示为一个实例。
由衰落引起的接收信号随机波动可以被建模为时间t的随机过程$\tilde{c}(\tau ,t)$。如果接收信号是各路径的大量散射分量之和，则由中心极限定理可知，$\tilde{c}(\tau ,t)$表示为时间的复高斯过程。在任意时刻t，概率密度函数的实部和虚部都是高斯的。这个模型意味着对于每个$\tau$或${\tau }_k$，信号路径是由大量不能分离的分量组成，因此，$\tilde{c}(\tau ,t)$和${\tilde{a}}_k(t)$都是时间的复高斯过程。

如果$\tilde{c}(\tau ,t)$具有零均值，则包络$R=|\tilde{c}(\tau ,t)|$具有形式如下的瑞利概率密度函数

$$f_R(r)=\frac{r}{{\sigma }^2}\exp \left(-\frac{r^2}{2{\sigma }^2}\right),r>0$$
其中${\sigma }^2$是$\tilde{c}(\tau ,t)$的实部和虚部的方差。

如果$\tilde{c}(\tau ,\tau )$具有非零均值，则意味着有视距非衰落的路径(称为镜像分量)存在，包络$R=|\tilde{c}(\tau ,T)|$具有形式如下的莱斯概率密度函数

$$f_R(r)=\frac{r}{{\sigma }^2}{I}_0\left(\frac{Ar}{{\sigma }^2}\right)\exp \left(-\frac{r^2+A^2}{2{\sigma }^2}\right),r>0$$
其中A是$\tilde{c}(\tau ,t)$的非零均值，$I_0(z)$是修正的Bessel函数，定义如下

$$I_0(z)=\frac{1}{2\pi }{\int }_0^{2\pi }\exp \left(z\cos (u)\right)\mathrm{d}u$$
比值$K=(A^2)/({\sigma }^2)$叫莱斯分布因子，是非衰落分量功率对衰落分量功率的比值。$K\gg 1$表示衰落不严重，而$K\ll 1$表示存在严重的衰落。

$\tilde{c}(\tau ,t)$的概率密度函数描述了复冲激响应的瞬时值，而时变特性可以使用变量为t的随机过程适当的自相关函数或功率谱密度来建模。

3.2. 信道响应中的时变模型

要对信道的时变特性进行数学建模，可将$\tilde{c}(\tau ,t)$视为一个参数的广义平稳(wide
sense stationary, WSS)随机过程，其自相关函数为

$$R_{\tilde{c}\tilde{c}}({\tau }_1,{\tau }_2,\alpha )=E\left[{\tilde{c}}^{\ast }({\tau }_1,t)\tilde{c}({\tau }_2,t+\alpha )\right]$$
在大多数的多径信道中，假设对应于不同延迟(即不同路径)的衰或和相位偏移不相关，这种不相关散射(uncorrelated
scatterig, US)假设，有

$$R_{\tilde{c}\tilde{c}}({\tau }_1,{\tau }_2,\alpha )=R_{\tilde{c}\tilde{c}}({\tau }_1,\alpha )\delta ({\tau }_1-{\tau }_2)$$
上式体现了广义平稳和不相关散射两方面的假设，它经常被称为衰落的WSSUS模型，最早是由Bello提出的。用$R_{\tilde{c}\tilde{c}}({\tau }_1,\alpha )$来表示自相关函数，则有

$$R_{\tilde{c}\tilde{c}}({\tau }_1,\alpha )=E\left[{\tilde{c}}^{\ast }(\tau ,t)\tilde{c}(\tau ,t+\alpha )\right]$$
对自相关函数作傅里叶变换，我们可以得到衰落信道的频域模型，用功率谱密度表示为

$$S(\tau ,\lambda )=F{R}_{\tilde{c}\tilde{c}}(\tau ,\alpha )={\int }_{-\infty }^{\infty }{R}_{\tilde{c}\tilde{c}}(\tau ,\alpha )\exp (-j2\pi \lambda \alpha )\mathrm{d}\alpha$$
$S(\tau ,\lambda )$被称为信道的散射函数，是一个有两个自变量的函数，其中一个是时域变量(延迟),另一个是频域变量，即多普勒频率变量。散射函数以延迟和多普勒频率的函数形式，给出了一个信道平均功率输出的度量。

通过散射函数，我们可以得到信道的最重要的参數，它们将影响在信道上运行的通信系统的性能。我们从定义如下的"多径强度"曲线开始人手

$$p(\tau )=R_{\tilde{c}\tilde{c}}(\tau ,0)=E\left[|\tilde{c}(\tau ,t){|}^2\right]$$

上式表示作为延迟函数的平均接收功事，上式通常被称为功率-延迟曲线。可以证明，$p(\tau )$与散射函数通过下式相关联

$$p(\tau )={\int }_{-\infty }^{\infty }S(\tau ,\lambda )\mathrm{d}\lambda$$
刻画衰落的另一个有用的函数是多普勒功率谱，
它也可以通过下式由散射函数推导出来

$$S_d(\lambda )={\int }_{-\infty }^{\infty }S(\tau ,\lambda )\mathrm{d}\tau$$
这些函数的相互关系如图所示

{width=“0.8\linewidth”}

3.3. 进一步假设

• 每个多径信号延迟、衰减互相彼此独立，不相关

• 多普勒功率谱密度模型存在多个多径分量，其中每个分量的延
  迟不同，但都具有相同形状的多普勒功率谱，因此每个多径分量
  的多普勒功率谱仍然是高斯的

• 每个多径信号的复衰减因子包络瑞利分布，相位均匀：即实部高斯分布，虚部也是高斯的，
  实部虚部不相关。换言之，每路多径信号经过瑞利信道。

• 假设多径分量的数量和延迟结构的变化相比复衰减因子的的变化要缓慢，因此，在仿真期间可以认为是常数。

综合以上假设，可以得到第$k$
条路径的时域相关函数$R_k(\alpha )=E\left[{\tilde{a}}_k^{\ast }(t){\tilde{a}}_k(t+\alpha )\right]$
对$R_k(\alpha )$做傅里叶变换可以得到 $S_k(\lambda )$,即多普勒功率谱。

题目限定$S_k(\lambda )$是高斯的，则
$R_k(\alpha )$也是高斯的，因此仿真需要产生一个
时域高斯相关的随机过程的一个样本函数。可以通过图1{reference-type=“ref”
reference=“fig:gain”}来实现。

4. 仿真模型

多径衰落信道的仿真模型可以用如图所示的延时抽头线(TDL)来实现。

{width=“0.8\linewidth”}

4.1. 抽头增益过程的产生

抽头增益过程是具有高斯概率密度函数和任意功率谱密度函数的平稳随机过程。最简单的抽头增益过程模型假设它们都是不相关的复高斯过程，并具有零均值、相同功率谱密度和不同方差。在这种情况下，
如图所示，抽头增益过程可以通过对白高斯过程进行滤波来产生。

{#fig:gain width=“0.8\linewidth”}

选择滤波器传递函数，使其产生所需的多普勒功率谱密度。换句话说，选择$H(f)$使

$$S_{\tilde{g}\tilde{g}}(f)=S_d(f)=S_{\tilde{w}\tilde{w}}(f)|\tilde{H}(f){|}^2=|\tilde{H}(f){|}^2$$
其中$S_{\tilde{w}\tilde{w}}(f)$是输入白噪声过程的功率谱密度，可以将其设置为1；$S_{\tilde{g}\tilde{g}}(f)$是给定的抽头增益过程的多普勒功率谱密度，选择滤波器增益使得$\tilde{g}(t)$具有归一化功率1。图中的静态增益${\sigma }_n$说明不同抽头可以具有不同功率级或方差。如果抽头增益的功率谱密度不同，则需要使用不同的滤波器来实现不同的抽头。

5. 仿真结果

仿真设计了一个复基带信号

$$\tilde{s}(t)=\sum _{i=1}^3\exp \left(2\pi f_{i}t\right)$$
其中
$f_1=1Hz,f_2=3Hz,f_3=5Hz$,多径数量$n=10$，等作为
自定义函数$\mathrm{c}\mathrm{h}\mathrm{a}\mathrm{n}\mathrm{n}\mathrm{e}\mathrm{l}.\mathrm{m}$文件的输入。

5.1. 快衰落

n=10;%多径数量
tau_max=0.001;%最大可能多径时延
tau=[0,tau_max*rand(1,n-1)];%随机生成n个时延
Ptau=1./(tau+1);%根据时延分配功率
f1=1;f2=3;f3=5;%设计基带频率
%设计基带复信号
singal=@(tt) exp(j*2*pi*f1*tt)+exp(j*f2*2*pi*tt)+exp(j*f3*2*pi*tt);
channel(n,tau,Ptau,1,singal);%快衰落

在这种输入参数下，将会发生快衰落，仿真结果如图2{reference-type=“ref”
reference=“fig:kaui1”},3{reference-type=“ref”
reference=“fig:kaui2”},4{reference-type=“ref”
reference=“fig:kaui3”}

{#fig:kaui1 width=“0.8\linewidth”}

{#fig:kaui2 width=“0.8\linewidth”}

{#fig:kaui3 width=“0.8\linewidth”}

1. 观察图2{reference-type=“ref”
   reference=“fig:kaui1”}，可以发现，复衰减因子的功率谱是高斯的，满足题目要求。

2. 观察图3{reference-type=“ref”
   reference=“fig:kaui2”},可以发现，信道的冲激响应在信道的符号周期内变换很快

3. 观察图4{reference-type=“ref”
   reference=“fig:kaui3”},可以发现，输出信号的频谱在$f_1,f_2,f_3$处产生
   了频谱扩展，即产生了新的频谱分量

5.2. 慢衰落

channel(n,tau,Ptau,0.001,singal);%观察慢衰落

{#fig:man1 width=“0.8\linewidth”}

{#fig:man2 width=“0.8\linewidth”}

{#fig:man3 width=“0.8\linewidth”}

1. 观察图5{reference-type=“ref”
   reference=“fig:man1”}，可以发现，复衰减因子的功率谱几乎是一个冲激函数。

2. 观察图3{reference-type=“ref”
   reference=“fig:kaui2”},可以发现，信道的冲激响应在信道的符号周期内变换的很慢

3. 观察图4{reference-type=“ref”
   reference=“fig:kaui3”},可以发现，输出信号的频谱在$f_1,f_2,f_3$处几乎不产生
   频谱扩展，即不产生新的频谱分量

4. 这些是慢衰落的一些特性

5.3. 平坦衰落

channel(n,tau,Ptau,0.01,singal);%观察平坦衰落

{#fig:tan2 width=“0.8\linewidth”}

{#fig:tan3 width=“0.8\linewidth”}

1. 观察图8{reference-type=“ref”
   reference=“fig:tan2”},可以发现，信道的冲激响应在信道的符号周期内有较大变化

2. 观察图9{reference-type=“ref”
   reference=“fig:tan3”},可以发现，输出信号的频谱在$f_1,f_2,f_3$具有相同的功率，信道
   对3个频率有相等的增益。虽然频域没有发生变化，但因为信道是时变的，所以
   输出信号的幅度会有剧烈的变动。

3. 这些是平坦衰落的一些特性

5.4. 选择性衰落

tau_max=1;%最大可能多径时延
tau=[0,tau_max*rand(1,n-1)];%随机生成n个时延
Ptau=1./(tau+1);
channel(n,tau,Ptau,0.001,singal);%选择性衰落

{#fig:xuan2 width=“0.8\linewidth”}

{#fig:xuan3 width=“0.8\linewidth”}

1. 观察图11{reference-type=“ref”
   reference=“fig:xuan3”},可以发现，输入信号的频谱在$f_1,f_2,f_3$具有相同的功率，但是输出信号的
   对3个频率有不同的功率，即信道对不同频段有不同的增益。频谱发生失真，
   但是输出信号时域幅度变动不大。

2. 这些是选择性衰落的一些特性

6. 代码

clear
clf
n=10;%多径数量
tau_max=0.001;%最大可能多径时延
tau=[0,tau_max*rand(1,n-1)];%随机生成n个时延
Ptau=1./(tau+1);%根据时延分配功率
f1=1;f2=3;f3=5;%设计基带频率
%设计基带复信号
singal=@(tt) exp(j*2*pi*f1*tt)+exp(j*f2*2*pi*tt)+exp(j*f3*2*pi*tt);
% channel(n,tau,Ptau,1,singal);%快衰落
% channel(n,tau,Ptau,0.001,singal);%观察慢衰落
% channel(n,tau,Ptau,0.01,singal);%观察平坦衰落
tau_max=1;%最大可能多径时延
tau=[0,tau_max*rand(1,n-1)];%随机生成n个时延
Ptau=1./(tau+1);
channel(n,tau,Ptau,0.001,singal);%选择性衰落

function out=channel(n,tau,Ptau,Bd,signal_in)
% n   = 1;        %多径数目
% tau =  ;        %多径时延
% Ptau=  ;        %各多径分量平均功率
% Bd  =  ;        %多普勒频谱的3dB带宽
% signal_in       %输入的信号表达式，为一匿名函数


Fs=100;Ts=1/Fs;
N=2e+4;
t=(0:N-1)*1/Fs;


%%设计高斯滤波器，目的是为每一路产生相关的复衰减因子
a=sqrt(1/4*log(2)/Bd^2);%高斯滤波器参数，若Bd为3dB带宽a^2Bd^2=1/4ln2


H_gauss_fre =@(f) exp(-a^2*f.^2)/sqrt(2*pi/(2*a^2));%频域
H_gauss_time=@(t) sqrt(t)./a.*exp(-pi^2/a^2*t.^2);%时域
Gauss_t=(0:N-1)*Ts;
% sFilt = filter([0 1],[1,1], s); % Example filtered signal

%%产生n路单位复高斯白噪声
No_complex=wgn(n,N,1,'linear','complex');

%%各路的增益
K=sqrt(Ptau);

signal_out=zeros(n,N);
at=zeros(n,N);

for i=1:n
    %生成第i个路径复增益的一个样本函数
    at(i,:)=filter(K(i)*H_gauss_time(Gauss_t),1,No_complex(i,:));
    %第i路输出
    signal_out(i,:)=at(i,:).*signal_in(t-tau(i));
end

%%合并第n个路径信号
out=sum(signal_out);


%%绘图:复增益
f1=figure; 
subplot 311
plot(real(No_complex(1,:)));
title('复高斯白噪音')
% at_conv=conv(No_complex,H_gauss_time(t),'same');
[r,lags]=xcorr(at(1,:),'unbiased');
% var(at)
maxlags=max(lags);
% maxlags
index1=maxlags-100:maxlags+100;
subplot 312
plot(lags(index1)*Ts,real(r(index1)));
title('高斯滤波后的波形:复衰减因子的时域自相关')
subplot 313
Fa=abs(fftshift(fft(real(r))));
Fa_normaliz=Fa/max(Fa);
index2_shift=round(20*Bd*N/Fs);
if(round(20*Bd*N/Fs)>=N)
    index2_shift=N-1;
end
index2=N-index2_shift:N+index2_shift;
plot(lags(index2)*Fs/(2*N),Fa_normaliz(index2));
title('一条多径上的复衰减因子的功率谱(多普勒频谱)')

%%绘图:时变单位脉冲响应
f2=figure;
pulse=abs(at);
index3=1:round(N/1200);
stem3(t(index3),tau,pulse(:,index3),'filled')
ylabel('时域延迟分量')
xlabel('时间')
title('信道冲激响应h(\tau,t)')



%%绘图：输入输出波形
f3=figure;
index=round(N/2):round(N/2+N/6);
subplot 411
stem(tau-min(tau),Ptau,'filled')
title('功率-延迟曲线')

subplot 412

% plot(t(index),real(signal_out(1,index)))
plot(t(index),real(out(index)))
title('输出信号时域')

subplot 413
plot(t(index),real(signal_in(t(index))))
title('输入信号时域')

subplot 414
out=real(out);
[r2,lags2]=xcorr(out,'unbiased');
index414=N-round(40*N/Fs):N+round(40*N/Fs);
F_out=abs(fftshift(fft(r2)));
plot(lags2(index414)*Fs/(2*N),F_out(index414))
title('输出信号频域')


end