轻松搞懂均匀分布、⾼斯分布、瑞利分布、莱斯分布(含MATLAB代码) ⽂章⽬录
1、均匀分布
在概率论和统计学中,均匀分布也叫矩形分布,它是对称概率分布,在相同长度间隔的分布概率是等可能的。 均匀分布由两个参数a和b定义,它们是数轴上的最⼩值和最⼤值,通常缩写为U(a,b)。
举个例⼦,掷骰⼦就是⼀个均匀分布,概率论中⼀个很常⽤分布。
MATLAB代码
%%利⽤rand函数产⽣服从(a-b)均匀分布的随机序列。
clear
close all
a=2;%(a-b)均匀分布下限
b=3;%(a-b)均匀分布上限
fs=1e7;%采样率,单位:Hz
s11348
t=1e-3;%随机序列长度,单位:s
n=t*fs;
rand('state',0);%把均匀分布伪随机发⽣器置为0状态
u=rand(1,n);%产⽣(0-1)单位均匀信号 x=(b-a)*u+a;%⼴义均匀分布与单位均匀分布之间的关系
subplot(211);plot(x);%输出信号图
title('均匀分布信号');
subplot(212);hist(x,2:0.01:3);%输出信号的直⽅图刷毛辊
title('均匀分布信号直⽅图');
2、⾼斯分布(正态分布)
正态分布(Normal distribution),也称“常态分布”,⼜名⾼斯分布(Gaussian distribution),最早由棣莫弗(Abraham de Moivre)在求⼆项分布的渐近公式中得到。C.F.⾼斯在研究测量误差时从另⼀个⾓度导出了它。P.S.拉普拉斯和⾼斯研究了它的性质。是⼀个在数学、物理及⼯程等领域都⾮常重要的概率分布,在统计学的许多⽅⾯有着重⼤的影响⼒。
22
若随机变量X服从⼀个数学期望为μ、⽅差为σ的正态分布,记为N(μ,σ)。其概率密度函数为正态分布的期望值μ决定了其位置,其标准差σ决定了分布的幅度。当μ = 0,σ = 1时的正态分布是标准正态分布。
⽣活中很多现象均近似符合整体分布,如整个班级同学的⾝⾼、体重情况。
MATLAB代码
%%正态分布
x=-10:0.01:10;
y=normpdf(x,0,1);%正态分布函数。
figure;
d2x说明书axes1=axes('Pos',[0.10.10.850.85]);
plot(x,y);
]);
set(axes1,'YLim',[-0.010.43],'XLim',[-33
3、瑞利分布
军用伪装网当⼀个随机⼆维向量的两个分量呈独⽴的、有着相同的⽅差的正态分布时,这个向量的模呈瑞利分布。瑞利分布是最常见的⽤于描述平坦衰落信号接收包络或独⽴多径分量接受包络统计时变特性的⼀种分布类型。两个正交⾼斯噪声信号之和的包络服从瑞利分布。
是不是感觉好抽象,有⼈问,为什么会产⽣瑞利分布?这是因为瑞利分布是最常见的⽤于描述平坦衰落信号接收包络或独⽴多径分量接受包络统计时变特性的⼀种分布类型,通信中很多信道的噪声满⾜这种分布,因此经常会⽤到。
MATLAB代码
clear all;close all;
sigma=2;%瑞利分布参数sigma
fs=1e7;%采样率,单位:Hz
t=1e-3;%随机序列长度,单位:s
t1=0:1/fs:t-1/fs;
n=length(t1);
rand('state',0);%把均匀分布伪随机发⽣器置为0状态
u=rand(1,n);%产⽣(0-1)单位均匀信号
x=sqrt(2*log2(1./u))*sigma;%⼴义均匀分布与单位均匀分布之间的关系
subplot(211);plot(x);%输出信号图
xlabel('t(单位:s)');
title('瑞利分布信号');%输出信号的直⽅图
subplot(212);hist(x,0:0.2:20);
title('瑞利分布信号直⽅图');
4、莱斯分布
莱斯分布实际上可以理解为主信号与服从瑞利分布的多径信号分量的和。概率密度函数公式中,R即
为正弦(余弦)信号加窄带⾼斯随机信号的包络,参数A是主信号幅度的峰值,σ^2是多径信号分量的功率,I()是修正的0阶第⼀类贝塞尔函数。
是不是感觉这个更抽象了,那有什么⽤呢,在通信中,有⼀个信号占主要成分的噪声中,信道噪声⼀般呈现莱斯分布。
MATLAB代码:
N =100000;
K =0.5;
const=1/(2*(K+1));
x1=randn(1,N);
x2=randn(1,N);
x=sqrt(const*((x1+sqrt(2*K)).^2+x2.^2));
[y,xi]= ksdensity(x);
气门绞刀plot(xi,y,'bo')
%验证
hold on
sig =1;
v =1;
yn = zeros(size(xi));
k=find(b >0& xi >=0);
if any(k),
xk = xi(k);
% Rician分布的概率密度函数
yp(k)=(xk ./ sig^2).* exp((-xk.^2+ v.^2)./(2*sig^2)).* besselj(0,(xk .*v ./ sig^2)); end
铝空气电池plot(xi,yp,'b')