文章目录

  • 一、瑞利信道简介
  • 二、从表达式看瑞利分布
  • 三、基于MATLAB的瑞利信道仿真

一、瑞利信道简介

在许多无线通信系统的仿真教材中,我们常常能够看到 “瑞利信道” 的身影,下面我们来解释一下什么是瑞利信道:在城区环境下,发射机与接收机之间的通信往往不能是直达的,从发射机发出的电磁波信号会受到高大建筑物或者是树木的影响,会造成多径效应。因此,接收机所受到的信号是会来自四面八方的,这些多径信号的幅度和相位都是随机变化的,其中,信号包络服从瑞利分布;信号相位为 000~2π2\pi2π之间的均匀分布,因此,这一类信道的衰落我们都可以用瑞利分布来表征。

二、从表达式看瑞利分布

我们从表达式上简单看看:
假设我们发射机要发送的是一个正弦波信号(只携带单一符号):s(t)=A(t)cos(ωct)s(t) = A(t)cos(\omega_ct)s(t)=A(t)cos(ωc​t)
这个信号经过城市环境,接收端受到的,就是多路信号的叠加,而且每一路信号的幅度相位都不一定一样,如下:r(t)=a1(t)cosωc[t−τ1(t)]+a2(t)cosωc[t−τ2(t)]+⋯+an(t)cosωc[t−τn(t)]=∑i=1nai(t)cosωc[t−τi(t)]\begin{aligned} r(t) &= a_1(t)cos\omega_c[t - τ_1(t)]+a_2(t)cos\omega_c[t - τ_2(t)]+\cdots +a_n(t)cos\omega_c[t - τ_n(t)]\\ &=\sum_{i=1}^na_i(t)cos\omega_c[t - τ_i(t)] \end{aligned}r(t)​=a1​(t)cosωc​[t−τ1​(t)]+a2​(t)cosωc​[t−τ2​(t)]+⋯+an​(t)cosωc​[t−τn​(t)]=i=1∑n​ai​(t)cosωc​[t−τi​(t)]​

那么其实我们也可以写成下面的形式:r(t)=∑i=1nai(t)cos[ωct+φi(t)]r(t) = \sum_{i=1}^na_i(t)cos[\omega_ct +φ_i(t)] r(t)=i=1∑n​ai​(t)cos[ωc​t+φi​(t)]
也就是幅度和相位的形式;其中,φi(t)=−ωcτi(t)φ_i(t) = - \omega_cτ_i(t)φi​(t)=−ωc​τi​(t)

接下来,我们用cos的展开公式,可以得到:r(t)=∑i=1nai(t)cosφi(t)cos(ωct)−∑i=1nai(t)sinφi(t)sin(ωct)=X(t)cos(ωct)−Y(t)sin(ωct)\begin{aligned} r(t) &= \sum_{i=1}^na_i(t)cosφ_i(t)cos(\omega_ct)-\sum_{i=1}^na_i(t)sinφ_i(t)sin(\omega_ct)\\ &=X(t)cos(\omega_ct) - Y(t)sin(\omega_ct) \end{aligned} r(t)​=i=1∑n​ai​(t)cosφi​(t)cos(ωc​t)−i=1∑n​ai​(t)sinφi​(t)sin(ωc​t)=X(t)cos(ωc​t)−Y(t)sin(ωc​t)​
其中,X(t)=∑i=1nai(t)cosφi(t)X(t) = \sum_{i=1}^na_i(t)cosφ_i(t)X(t)=∑i=1n​ai​(t)cosφi​(t);Y(t)=∑i=1nai(t)sinφi(t)Y(t) =\sum_{i=1}^na_i(t)sinφ_i(t)Y(t)=∑i=1n​ai​(t)sinφi​(t)

由于 X(t)X(t)X(t) 和 Y(t)Y(t)Y(t) 都是相互独立的随机变量之后,根据中心极限定理,大量独立随机变量之和的分布趋于正态分布。因此,当n足够大时,X(t)X(t)X(t)和Y(t)Y(t)Y(t)都趋于正态分布。通常情况下X(t)X(t)X(t)和Y(t)Y(t)Y(t)的均值为0(由于没有直射路径),方差相等。而同时我们也知道,在MATLAB里面,生成满足正态分布的随机数是非常简单的。

下一步,我们再将上面的式子整理一下,得到:
r(t)=R(t)cos[ωct+φ(t)]r(t) = R(t) cos[\omega_ct +φ(t)] r(t)=R(t)cos[ωc​t+φ(t)]
其中,R(t)=X2(t)+Y2(t)R(t) = \sqrt{X^2(t) + Y^2(t)}R(t)=X2(t)+Y2(t)​;φ(t)=arctanY(t)X(t)φ(t) = arctan \frac{Y(t)}{X(t)}φ(t)=arctanX(t)Y(t)​

我们知道:X(t)X(t)X(t)和Y(t)Y(t)Y(t)趋于正态分布;而我们再看看瑞利分布的定义:

瑞利分布(Rayleigh distribution),当一个随机二维向量的两个分量呈独立的、有着相同的方差的正态分布时,这个向量的模呈瑞利分布。例如,当随机复数的实部和虚部独立同分布于0均值,同方差的正态分布时,该复数的绝对值服从瑞利分布。该分布是以瑞利命名的。

因此,我们可以知道 R(t)R(t)R(t) 就是服从瑞利分布的。

三、基于MATLAB的瑞利信道仿真

许多教材上对瑞利信道的仿真只是画出了瑞利信道理论的概率密度函数,可是在实际工程仿真时,这个理论的概率密度函数还真不太好用,因为我们知道通信系统最简化的模型就是:r(t)=s(t)∗h(t)+n(t)r(t) = s(t)*h(t) + n(t)r(t)=s(t)∗h(t)+n(t)
因此,我们常常希望可以得到信道的 hhh,那么,下面我们给出一种合理且简单的方法获得瑞利信道的衰减系数——即通过MATLAB的randn函数产生两个零均值的高斯随机变量,其包络就满足瑞利分布。

%-----------------------Rayleigh信道衰落系数的计算-------------------------%
clc
clear
close allsigma_ = [0.5 1 2 3 4];
figure;
for i = 1:5sigma = sigma_(i);V = 0:0.01:6;pdf_rayleigh = (2*V/sigma^2).*exp(-V.^2/sigma^2);plot(V, pdf_rayleigh,'-');hold on;N = 10^5;                                         %测试样本点数x = sqrt(sigma^2/2)*randn(2,N)+0;  %产生均值为0,方差为sigma^2/2的高斯随机变量xx = x(1,:) + sqrt(-1)*x(2,:);R = abs(xx);                                      %信号包络(理论上服从瑞利分布)k = 0:0.01:6;                                     %统计区间dk = 0.01;                                        %组距p = [];for i = 1:(length(k)-1)num = length(find(R >= k(i) & R < k(i+1)));   %找到R里面介于k(i)和k(i+1)的元素的个数p = [p num/length(R)];                        %得到R中的值介于k(i)和k(i+1)的概率endpdf = p/dk;                                       %概率除以组距即可得到概率密度k_new = 0.005:0.01:(6-0.005);                     %取一组的中值表征这个区间plot(k_new, pdf,'-.');grid on;set(gca,'GridLineStyle',':','GridColor','k','GridAlpha',1);xlabel('R');ylabel('PDF');title('Rayleigh衰落的概率密度函数');legend('理论值sigma=0.5','仿真值sigma=0.5','理论值sigma=1','仿真值sigma=1',...'理论值sigma=2','仿真值sigma=2','理论值sigma=3','仿真值sigma=3',...'理论值sigma=4','仿真值sigma=4');
end

可以看出,所得到的测试样本很好地满足理论的瑞利信道概率密度函数,因此在未来整个通信系统的仿真时,我们可以用这个方法获取瑞利信道的 hhh。

【通信原理 入坑之路】基于MATLAB的移动通信系统仿真 之 瑞利衰落信道的原理与仿真相关推荐

  1. MATLAB多径衰落信道仿真程序,基于Matlab的移动通信中多径衰落信道的仿真

    移动通信中多径衰落信道的仿真 文章编号:1007-3558(2004)05-0078-04 基于Matlab的移动通信中多径衰落信道的仿真 陈岚,,万国春,,冯志彪2 (1.江西科技师范学院,江西南昌 ...

  2. 【通信原理 入坑之路】—— 理解 4G,5G中的天王山:OFDM【学习笔记4】

    今天将要讨论的部分个人感觉非常绕,博主也是在阅读了大量相关书籍和文献之后才有所体悟.如有表述不当之处欢迎大家在评论区斧正! 在上一篇博文中,我们知道串并转换使得OFDM的符号速率降低了,可是总体符号传 ...

  3. 【通信原理 入坑之路】—— 信号的包络分析2 包络分析的意义以及OQPSK调制的引入

    本篇 B l o g Blog Blog 是接上一篇未完的分析而作的,主要是继续深入理解信号包络以及包络的意义.同时,本篇 B l o g Blog Blog 是<深入浅出通信原理>的学习 ...

  4. 【通信原理 入坑之路】—— 详解IQ调制以及星座图原理

    写在前面:本博客是<深入浅出通信原理>的学习笔记,仅供个人学习参考使用 文章目录 一. IQ调制与解调的原理与过程 1.1 利用旋转向量理解IQ调制(正交调制) 1.2 利用旋转向量理解I ...

  5. 【通信原理 入坑之路】——深入、详细地理解通信里面“卷积”概念

    文章目录 一.多项式乘法和"卷积" 二."卷积"的计算公式 2.1 实例:利用卷积计算两个信号相乘 三.卷积的生动理解 3.1 卷积能够用来解决什么问题? 3. ...

  6. 【通信原理 入坑之路】——信息论部分 之 离散无记忆信源的等长编码

    课本上的讲解公式实在是太多啦!看得人头大,那么本文试图从繁杂的公式中梳理一下几种等长编码的方法,并在最后给出一些解题技巧.我们开始吧! 文章目录 离散无记忆信源的等长编码 1.1 单符号等长编码 1. ...

  7. 【通信原理 入坑之路】—— 数字载波传输系统 之 2PSK(二进制相移键控)2ASK(二进制幅度键控)的超详细分析

    在之前的文章里面,我们讨论过数字基带传输系统,他们指的是传输二进制脉冲成型信号.然而,在绝大部分场合,我们需要将数字基带信号装载到高频载波上(也就是我们所说的调制),通过天线等将高频波发射出去.这就是 ...

  8. 【通信原理 入坑之路】—— 信息论部分 离散无记忆信源的不等长编码 之 霍夫曼编码过程详解

    还记得在上一篇博客中,我们讨论的离散无记忆信源的等长编码吗,我们最后得到的结论是 -- 如果使用等长编码,那么对典型序列编码的编码效率是最高的.然而,要满足较高的编码效率,同时还要保证误码率的情况下, ...

  9. 【通信原理 入坑之路】—— 模拟信号的数字编码 之 PCM编码(A律13折线和μ律15折线的编译码剖析)

    引言: 在我们的数字通信系统中,一般要发送的绝大多数都是模拟信号.而所谓数字通信,传递的就应该是数字信号.因此,我们首先要对这个待发送的模拟信号进行采样(这在我之前<信号与系统>的相关博客 ...

  10. qam已调信号matlab相干解调,通信原理课程设计-基于MATLAB的M-QAM调制及相干解调的设计与仿真.doc...

    通 信 原 理 课 程 设 计 报 告 题 目:基于MATLAB 的M-QAM调 制及相干解调的设计与仿真 班 级:通信工程 1411 姓 名: 成 绩: 日 期:2016 年 12 月 21 日 基 ...

最新文章

  1. 推荐阅读的多核编程技术书籍
  2. NLP之TFTS读入数据:TF之TFTS读入时间序列数据的几种方法
  3. 计算机图形学图形旋转_计算机图形学翻译
  4. Linux运维问题解决(1)——Linux 定时任务 crontab 配置及示例
  5. 举例 微积分 拉格朗日方程_理解微积分真谛:微分方程
  6. 修改TIM文件存储路径
  7. 运放搭建电压电流转换电路分析
  8. 可解释机器学习- LIME模型讲解|interpretable machine learning-LIME
  9. 「IT女神勋章」挑战赛#
  10. 关于计算机的英语作文带翻译100字,求30篇100字英语作文带翻译
  11. widget安卓桌面插件初步使用
  12. javax.servlet.ServletException: java.lang.NoSuchMethodError 解决方法解决办法
  13. Android 加密 AES
  14. 8-14-粘包问题、(通过简单报头)解决粘包问题、定制复杂报头
  15. 【安全知识】——常见杀软对应进程名
  16. 微信小程序中的wxml如何自动缩进、对齐
  17. Silverlight游戏研发手记:(一)差集运算在SLG战斗范围设定中的应用
  18. OSGI框架搭建常见问题即错误
  19. win7怎么调整屏幕亮度_调整Mac外部屏幕亮度小工具
  20. 怎么将webm格式转换成mp4,3招轻松学

热门文章

  1. 小米路由r2d论坛_维修一台亮黄灯的R2D小米路由器2
  2. web安全之逻辑漏洞
  3. 详解 YUV 格式(I420/YUV420/NV12/NV12/YUV422)
  4. 傻瓜式c语言教程,app制作教程简单易学:傻瓜简单自建app,步骤如下
  5. 微信 小程序弹框 简单实现
  6. Inno Setup打包基本笔记
  7. 单点登录Cas如何保证客户端的安全性
  8. 如何将MySQL卸载干净?
  9. [不限速百度云]两款百度网盘破解版分享
  10. Activiti7事件监听