(13)MATLAB莱斯(Rician)衰落信道仿真3

文章目录


前言

本文使用复高斯随机过程对莱斯衰落信道进行建模,给出MATLAB仿真代码,并和瑞利分布的PDF理论值进行了对比。


一、复高斯随机过程的莱斯衰落信道模型

在《(11)MATLAB莱斯(Rician)衰落信道仿真2》中,给出以下莱斯衰落信道模型:

该模型可以看作一个复高斯随机过程,该过程的实部和虚部都如从高斯分布,高斯的均值为

标准差为

下面使用该复高斯随机过程对莱斯衰落信道进行MATLAB建模仿真。

二、仿真代码

代码如下:

c 复制代码
clc
close all
clear all

%% 莱斯衰落信道模型
Kdb = -50;                                        % 莱斯因子K的dB值
K = 10^(Kdb/10);                                  % 莱斯因子K转换成线性值
mu = sqrt(K/(2*(K+1)));                           % 均值
sigma = sqrt(1/(2*(K+1)));                        % 标准差

N = 1e5;                                          % 样本数
ric = mu + sigma.*randn(1,N) + 1j*(mu + sigma.*randn(1,N));

% 利用直方图计算莱斯随机变量的概率密度函数估计值
bins_number = 50;
[elements_number,x] = hist(abs(ric),bins_number);
px = elements_number/N/(mean(diff(x)));

%% 瑞利分布概率密度函数(pdf)的理论值
s = sqrt(0.5);
r = 0: 0.1 : 4;                                   % rayleigh随机变量
pdf = (r/s^2).*exp(-r.^2/(2*s^2));                % 理论pdf

% 画图对比
figure()
plot(x,px,'ro','LineWidth', 1.5)
hold on
grid on
plot(r,pdf,'b-' ,'LineWidth', 1.5);
grid on;
title('莱斯分布随机变量的概率密度');
legend('莱斯随机变量的PDF估计值(K=-50dB)' ,'瑞利分布PDF的理论值' );
xlabel('随机变量x');
ylabel('概率密度px');

三、仿真结果画图

四、几点补充说明

(1)关于randn函数

MATLAB中randn函数用于正态分布伪随机数。正态分布,也即高斯分布。若高斯分布的均值为1,标准差为2,要生成该高斯分布的随机变量的100个样本,可以用这行代码实现:r = 1 + 2.*randn(100,1)。

(2)关于莱斯因子的形式

由于莱斯因子K是直射路径与反射路径信号的相对功率,所以,莱斯衰落信道仿真时,莱斯因子一般以dB值的形式给出。

(3)仿真中莱斯因子取值

上面给出的代码中,设置莱斯因子Kdb=-50,得到的随机变量的分布近似瑞利分布。另外,当设置莱斯因子Kdb=15时,得到的pdf估计值将近似高斯分布。


相关推荐
心情好的小球藻几秒前
Python应用进阶DAY9--类型注解Type Hinting
开发语言·python
静心问道4 分钟前
TrOCR: 基于Transformer的光学字符识别方法,使用预训练模型
人工智能·深度学习·transformer·多模态
说私域6 分钟前
基于开源AI大模型、AI智能名片与S2B2C商城小程序源码的用户价值引导与核心用户沉淀策略研究
人工智能·开源
亲持红叶7 分钟前
GLU 变种:ReGLU 、 GEGLU 、 SwiGLU
人工智能·深度学习·神经网络·激活函数
说私域7 分钟前
线上协同办公时代:以开源AI大模型等工具培养网感,拥抱职业变革
人工智能·开源
群联云防护小杜9 分钟前
深度隐匿源IP:高防+群联AI云防护防绕过实战
运维·服务器·前端·网络·人工智能·网络协议·tcp/ip
惜.己12 分钟前
使用python读取json数据,简单的处理成元组数组
开发语言·python·测试工具·json
摘星编程14 分钟前
构建智能客服Agent:从需求分析到生产部署
人工智能·需求分析·智能客服·agent开发·生产部署
不爱学习的YY酱17 分钟前
信息检索革命:Perplexica+cpolar打造你的专属智能搜索中枢
人工智能
Y40900119 分钟前
C语言转Java语言,相同与相异之处
java·c语言·开发语言·笔记