【matlab】高斯-均匀-指数-瑞利分布

% 【例8-7】利用randan函数产生高斯分布序列。

%  【例8-8】利用rand函数产生服从（a-b）均匀分布的随机序列。

% 【例8-9】服从指数分布的热噪声随机序列的实现。

% 【例8-10】产生瑞利分布的热噪声。

% 【例8-7】利用randan函数产生高斯分布序列。

clear all;close all;

y=randn(2000);

subplot(211);plot(y);

title('服从高斯分布的随机序列信号');

subplot(212);hist(y);

title('服从高斯分布的随机序列信号直方图');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%【例8-8】利用rand函数产生服从（a-b）均匀分布的随机序列。

clear all;close all;

a=2;          % (a-b)均匀分布下限

b=3;      % (a-b)均匀分布上限

fs=1e7;      % 采样率，单位：Hz

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('均匀分布信号直方图');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 【例8-9】 服从指数分布的热噪声随机序列的实现。

clear all; close all;

dba=2.5;      % 指数分布参数

fs=1e7;       % 采样率，单位：Hz

t=1e-3;      % 随机序列长度，单位：s

n=t*fs;

rand('state',0);    % 把均匀分布伪随机发生器设置为0状态

u=rand(1,n);      % 产生（0-1）单位均匀信号

x=log2(1-u)/(-dba); %广义均匀分布与单位均匀分布之间的关系

subplot(211);

plot(0:1/fs:t-1/fs,x);   %输出信号图

xlabel('t(单位：s)');

title('指数分布信号');

subplot(212);hist(x,0:0.05:4); %输出信号的直方图

title('指数分布信号直方图');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%【例8-10】产生瑞利分布的热噪声

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('瑞利分布信号直方图');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  张德丰编著 《MATLAB数字信号处理与应用》  北京:清华大学出版社,2010