基于Matlab的无线信道仿真

基于Matlab的无线信道仿真近几年，随着无线通信业务和新兴宽带移动互联网接入业务的快速增长，对无线通信系统的优化显得尤为重要。与有线信道静态和可预测的典型特点相反，在实际中，由于无线信道动态变化且不可预测，无线通信系统的性能在很大程度上取决于无线信道环境，所以对无线信道的准确理解和仿真对设计一个高性能和高频谱效率的无线传输技术显得尤其重要。无线信道的一个典型特征是“衰落”，衰落现象大致可分为两种类型：大尺度衰落和小尺度衰落。其中，大尺度衰落主要在移动设备通过一段较长的距离时体现，它是由信号的损耗（长距离传播）和大的障碍物（如建筑物、中间地形和植物）形成的阴影所引起的，一般分为路径损耗和阴影衰落，另一方面，小尺度衰落是指当移动台在较短距离内移动时，由多条路径的相消或相长干涉所引起信号电平的快速波动，主要表现为多径衰落。它们之间的关系如图1所示。报告中分别对这几种衰落的常见模型进行了总结和仿真。一、大尺度衰落大尺度衰落是在一个较大的范围上考察功率的渐变过程，功率的局部中值随距离变化缓慢。大尺度信道模型主要研究电波传播在时间、空间、频率范围内平均特性。路径损耗路径损耗由发射功率的辐射扩散及信道的传播特性造成，反映在宏观长距离

上。理论上认为，对于相同收发距离，路径损耗相同。其定义为有效发射功率和平均接收功率之间的比值。几种常用的描述大尺度衰落的模型有自由空间模型、对数距离路径损耗模型、Hata-Okumura模型。1.1.1自由空间模型所谓自由空间是指天线周围为无限大真空时的电波传播，它是理想传播条件。电波在自由空间传播时，其能量既不会被障碍物所吸收，也不会产生反射或散射，传播路径上没有障碍物阻挡，到达接收天线的地面反射信号场强也可以忽略不计。P (P (4ndf丫1GGtr其中，P为发射功率，P为接收功率，d为发射端与接收端距离，f为载波tr频率，c为光速取3x108，G为发射端天线增益，G为接收端天线增益。转换tr成分贝表示：PL(dB)=10lgJ=32.45+20lgd+20lgf-10lgGG

s p trr发射端与接收端均是全向天线，G=G=1，得图2：tr图2图2路径损耗随距离、频率变化曲线对数距离路径损耗模型与前面提到的自由空间路径损耗一样，在其他所有实际环境中，平均接收信号功率随距d呈对数方式减小。通过引入随着环境而改变的路径损耗指数n可以修正自由空间模型，从而构造出一个更为普遍的路径损耗衰落模型。4」(dB)=i°H10nlg(+)0其中，d0是一个参考距离，在参考距离或者接近参考距离的位置，路径损耗具有自由空间路径损耗的特点。如表1所示，路径函数主要由传播环境决定，对于不同的传播系统必须确定合适的参考距离d0，例如，在大覆盖范围的蜂窝系统（半径大于10km的蜂窝系统）中，通常会设置d0为1km，对于小区半径为1km的宏蜂窝系统或者具有极小半径的微蜂窝系统，可以分别设置参考距离为100m或1m。环境路径损耗（n）环境路径损耗（n）自由空间2建筑物内视距传输1.6~1.8市区蜂窝2.7~3.5建筑物内障碍物阻挡4~6市区蜂窝阴影3~5工厂内障碍物阻挡2~3表1路径损耗指数对数路径损耗随距离、路径损耗指数变化曲线：图3路径损耗随距离、路径损耗指数变化曲线

Hata-Okumura模型(1)Okumura模型Okumura模型的特点是：以大城市地区准平坦地形的场强中值路径损耗作为基准，对于不同的传播环境和地形条件等因素用校正因子加以修正。Okumura模型中大城市地区准平坦地形的中值路径损耗(dB)由下式给出L=L+A(f,d)-H(h,d)-H(h，f)Mbsm bb mm其中，d为发射端与接收端距离，f为载波频率，Hh为基站天线高度增益，bH为移动台天线高度增益。它们与路径损耗的关系如下：mJS20 305Q70UM?2OQJQO500700IOOU底猾无我也效曲鹿加JS20 305Q70UM?2OQJQO500700IOOU底猾无我也效曲鹿加E图4路径损耗与基站天线高度增益的关系秒胡合无纯而应JdnjHl工三：Jr-Hfe用驾整,里图5路径损耗与移动台天线高度增益的关系⑵Hata模型Hata模型仍然保留了Okumura模型的风格，以市区传播损耗为标准，其他地区在此基础上进行修正。中值路径损耗的经验公式为：L:69.55+26.16lg(f)-13.82lg(H)+(44.9-6.55lg(H))义lg(d)-a(h)M b bm其中，a(h)为修正因子，由所在环境决定。Hata模型适用条件：载波频率m(MHz)f=150-1000MHz；基站高度(m)Hb=30-100m；移动台高度(m)Hm=1-10m；收发天线距离(小区半径)(km)d=1-20km。中小城市修正因子:a(h);(1.11lg(f)-0.7)义H-1.56lg(f)-0.8mm大城市且载波频率f<200MHz时，修正因子：a(h)=8.29(lg(1.54xH))2-1.1

mm大城市且载波频率f>400MHz时，修正因子：a(h)=3.2(lg(11.75xH))2-4.97

mm郊区修正因子:a(h)=2(lg(f/28))2+5.4m农村修正因子:a(h)=40.98+4.78(lg(f))2-18.33lg(f)m仿真结果如下：

图6Hata模型中路径损耗随距离变化曲线(Hb=60,Hm=1,f=200MHz),,▲…中小城市郊区图6Hata模型中路径损耗随距离变化曲线(Hb=60,Hm=1,f=200MHz),,▲…中小城市郊区一，一农村T-大城市图7Hata模型中路径损耗随距离变化曲线(Hb=60,Hm=1,d=2Km)1.1.4布灵顿模型描述大尺度衰落的一种模型。假设发射天线和移动台之间的地面是理想平面大地，并且两者之间的距离远大于发射天线的高度h或移动台的高度h，此时的t r路径损耗公式为：L=120+4010gQ)—2010g(h)—2010g(h)P 10 010t 010r系统设计时一般把接收机高度按典型值 h=1.5m处理，此时路径损耗公式r为：L=116.5+4010g(d)-2010g(h)p ~10 ~10t1.1.5EgLi模型Egli认为不平坦地区的场强等于平面大地反射公式算出的场强加上一个修正值，该修正值为：G=2010g[40010f\J7式中f为工作频率，单位为MHz。Egli模型是从大量实测结果中归纳出来的中值预测公式，属于经验模型。E=E+2010gf幽土]+2010gf400 10(cd2) 10(f)Egli公式可用于计算不平坦地区的传播损耗，计算公式如下：

L=88+40logQ)-2010g(h)-2010g(h)+2010g(f)-Ga 修10 °10t 修10r 修10“L为传播损耗，单位dB；d为距离，单位km;h为发射天线高度，单位ma th为接收天线高度，单位m；f为工作频率，单位MHz；G为地形修正因r子，单位dB。1.2阴影衰落阴影衰落是指移动台在运动中通过不同障碍物的阴影时，构成的接收天线处场强中值的变化，从而引起的衰落。实测数据表明阴影衰落服从对数正态分布。对数正态阴影模型是最常用的统计仿真模型，发射和接收功率比值8=P/P，服tr从对数正态分布：1 (1n8-〃)2f(8；〃;5)= .e26285v2n农村(史)陵艰(史)段喉(史)陵眼10080其均值〃取决于路径损耗和所在区域内建筑物属性，随距离增加，障碍物数量会增加，衰减增加，这里〃近似为路径损耗。标准差5范围在4dB-13dB之间，代表不同的环境，密集城区、城区、郊区塞取8dB，农村5取7dB，扇区和高速5取4dB。仿真结果如下：农村(史)陵艰(史)段喉(史)陵眼1008060 I I I I I I I I I 60□ 10 20 30 40 50 60 70 80 90 100图8不同环境下的阴影衰落波形二、小尺度衰落小尺度衰落是指短期内的衰落，具体指当移动台移动一个较小距离时，接收信号在短期内的快速波动。当多径信号以可变相位到达接收天线时会引起干涉（即相位相同的相长干涉，相位不同的相消干涉）。换句话说，来自本地散射体的大量信号的相对相位关系决定了接收信号的电平波动。而且，每一个多径信号都可能发生变化，而这种变化依赖于移动台和周围物体的速度。总之，小尺度衰落由以下因素决定：多径传播、移动台速度、周围物体的速度和信号的传输带宽。根据不同无线环境，接收信号包络一般服从几种典型分布，如瑞利分布、莱斯分布和Nakagami-m分布。这里主要针对服从瑞利分布的多径信道进行模拟仿真。2.1瑞利信道Jakes模型通过对复正弦波的合成，产生服从给定多普勒谱的瑞利衰落信道。为了使产生的信号幅度近似服从瑞利分布，正弦波的数量必须足够大。此外，必须对每个正弦波发生器进行加权，以便产生想要的多普勒谱。这就是Jakes模型的主要思想。

h(t)=hi(t)+hQh(t)=hi(t)+hQ(t)fm2cos3cos(w t)No No72sin3N图9Jakes模型实现2cos3 cos(w2cos3cos(wt) 1 12cos3cos N2sin3cos(wt) 1 12cos3cos(N0 -0图9介绍了如何实现Jakes模型。假设以均匀方向到达的所有散射分量的射线被近似为N个平面波。定义N0=(N/2-1)/2，其中限定N/2为一个奇数。令0表示第n个平面波的到达角度，取值为0=2nn/N,n=1,2,...,N。如图9所示，n n0TOC\o"1-5"\h\z将N个频率为w的复振荡器的输出求和(w=w，cos0,n=1,2,...,N，每个振0 n ndn 0荡器的输出对应不同的多普勒频移)，然后与频率为w,=2nf的复振荡器的输出d m相加。在复振荡器的总和中，实部hI(t)和虚部hQ(t)可以分别表示为：h(t)=2N0(cos3coswt)+*2cos3coswt1、 、 ，nn Ndn=1和h(t)=2N0(sin3coswt)+v'2sin3coswtQ — ，nn 'Ndn=1其中，3为经过多普勒频移的第n个正弦信号的初始相位，3”为经过最大n N多普勒频移f的正弦信号的初始相位。初始相位的设置必须使衰落信道的相位m服从均匀分布。例如可以设置初始相位为：3N=03=nn/(N+1),n=1,2,...,N

n0 0EJakes模型的复输出可以表示为：h(t)=.0 {h(t)+jh(t)}2N+1IQ0 0其中，E.为衰落信道的平均幅度。经过多普勒频移的正弦信号频率{•}00 nn=1为：为：攻=攻cos0=2nfcos(2nn/N),n=1,2,...,N

ndnm 0经过多普勒频移的正弦数N0必须足够大，以便衰落信道的振幅能够近似服从瑞利分布。N0=8就足够大了。可以证明hI()和hQ(t)满