(10)MATLAB莱斯(Rician)衰落信道仿真1

文章目录


前言

首先给出莱斯衰落信道模型,引入了莱斯因子K,并给出莱斯分布的概率密度函数公式。然后导出莱斯分布随机变量的仿真表示式,建立MATLAB仿真代码,并根据莱斯衰落变量估计得到其PDF。


一、莱斯分布随机变量

与瑞利衰落不同,当无线信道中存在一个直射路径信号分量时,接收信号的包络将不再服从瑞利分布,而是服从莱斯分布,此时的小尺度衰落称为莱斯衰落。莱斯衰落时刻i的衰落幅度ri可以表示为:

其中β是直射路径分量的幅度,而xi和yi是服从均值为0、方差为σ^2的平稳高斯随机过程的样本。直射路径信号能量与散射路径信号能量的比值定义了所谓的Rician因子K,其表达式为

莱斯信道的概率密度函数为:

其中I0[.]是第一类零阶修正贝塞尔函数。

已知Rician分布的均方值为 2σ^2(K + 1),其中 σ^2是组成莱斯分布的高斯噪声过程的方差。另外,为了使信号功率和信噪比(SNR)一致,通常需要将莱斯分布的均方值设置为1,即 E{r^2} = 1。在满足E{r^2} = 1的条件下,式(1)可以写成以下形式:

式中,xi和yi是具有方差σ^2=1的零均值平稳高斯随机过程的样本。

接下来,将根据式(4)给出生成莱斯衰落随机变量的MATLAB代码,并根据随机变量计算出其PDF的估计值。

高斯分布随机变量仿真可以参考:
(3)MATLAB生成高斯随机变量及其概率密度函数估计

二、仿真代码与结果

生成莱斯衰落随机变量的MATLAB代码,并根据随机变量计算出其PDF的估计值。

1.仿真代码

莱斯分布随机变量MATLAB代码如下:

c 复制代码
clc
close all
clear all
%% 生成莱斯分布随机变量
Kdb = 1;                      % 莱斯因子K,分贝值
N = 100000;

% 生成莱斯分布随机变量
K = 10^(Kdb/10);              % 分贝值转换成线性值
const = 1/(2*(K+1));
x = randn(1,N);               % 高斯分布随机变量
y = randn(1,N);
r = sqrt(const*((x+sqrt(2*K)).^2 + y.^2));

% 莱斯随机变量的分贝值
rdb = 20*log10(r);
figure()
plot(rdb)
xlim([0 1e3]);
title('莱斯分布随机变量的分贝值');
grid on;
xlabel('samples');
ylabel('amplitude/dB');

%% 莱斯分布的pdf
bins_number = 30;
[elements_number,x] = hist(abs(r),bins_number);
pdfx = elements_number./N./(mean(diff(x)));

% 画图
figure()
plot(x,pdfx,'-*','LineWidth',1.5)
title('莱斯分布的概率密度函数')
grid on
xlabel('x');
ylabel('pdfx');

2.仿真结果画图

仿真代码运行结果画图如下:

(1)莱斯分布随机变量

(2)莱斯分布随机变量概率密度估计值


后续

下一篇文章将会给出莱斯衰落信道的建模方法与MATLAB代码。

相关推荐
满怀10154 分钟前
Python入门(5):异常处理
开发语言·python
攀小黑7 分钟前
Java 多线程加锁 synchronized 关键字 字符串当做key
java·开发语言
每次的天空17 分钟前
Kotlin 作用域函数:apply、let、run、with、also
android·开发语言·kotlin
小林熬夜学编程20 分钟前
【高并发内存池】第八弹---脱离new的定长内存池与多线程malloc测试
c语言·开发语言·数据结构·c++·算法·哈希算法
独好紫罗兰32 分钟前
洛谷题单3-P1980 [NOIP 2013 普及组] 计数问题-python-流程图重构
开发语言·python·算法
独好紫罗兰36 分钟前
洛谷题单3-P1009 [NOIP 1998 普及组] 阶乘之和-python-流程图重构
开发语言·python·算法
Taichi呀37 分钟前
PHP语言基础
android·开发语言·php
Aerkui1 小时前
Python数据类型-int
开发语言·python
莲动渔舟1 小时前
Nyquist插件基础:LISP语法-自定义函数
服务器·开发语言·lisp·音频处理·audacity
莲动渔舟1 小时前
Nyquist插件基础:LISP语法-运算符
开发语言·lisp·音频处理·audacity