close all
clc
clear all
N = 2e4;
NN = 100;

%方差
pn0 = 1;
pn1 = 10;
%pn1 = 1;   

A1 = rand(1,N);        
A2 = rand(1,N);

%高斯分布
ni = pn1*(-2*log(A1)).^(1/2).*cos(2*pi*A2);
nq = pn1*(-2*log(A1)).^(1/2).*sin(2*pi*A2);
n0 = pn0*(-2*log(A1)).^(1/2).*cos(2*pi*A2);
nz1 = ni-j*nq;
%包絡(luò)情況
nz2 = (ni.^2+nq.^2).^(1/2);   

S = randn (1,N);   
S = 10*randn(1,N);
%瑞利分布
nRay = sqrt(-2*pn1.^2*log(1-A1));
%指數(shù)分布
nexp = -log(1-A1);

%統(tǒng)計(jì)直方圖
 x0 = linspace(-5,5,NN);
 x1 = linspace(0,5,NN);


figure(1)
hist(ni,100)
% hold on
figure(2)
hist(S,100)
figure(3)
hist(nz2,100)
%axis([0 5 0 700])
figure(4)
hist(nRay,NN)
figure(5)
hist(nexp,NN)

%高斯分布的自相關(guān)函數(shù)和功率譜
ni_gauss = xcorr(ni,'biased');
ni_Pw = abs(fft(ni_gauss));

figure(6)
plot(ni_gauss),
figure(7)
plot(10*log10(ni_Pw/max(ni_Pw)));

xn = zeros(1,NN);
xn(50) = 1;
figure(8)
plot(abs(fft(xn,128)));

figure(9)
plot(ni)
figure(10)
plot(n0);