%%生成高斯白噪聲

clear
clc
%%
%初始化
A=65539;
N=1200;
x0=1;
M=2147483647;
C=0;%混合同余法直接變成乘同余法

%%

%乘法遞推N次
for?k=1:N
????x2=A*x0;
????x1=mod（x2M）;
????v1=x1/2147483647;%將x1中的數除以M得到小于1的隨機數
????v（:k）=v1;
????x0=x1;
end
v2=v;%保存0-1隨機數到v2
ave=mean（v）
var=var（v）
k1=k;
save?v;
%%
%繪圖程序
%k=1:k1;
%plot（kv‘r‘）;
%xlabel（‘k‘）;ylabel（‘v‘）;title（‘（0-1）‘）;
[numval]=hist（v10）;

%%
%生成服從正態分布的隨機數
for?n=1:100
????sum1（:n）=sum（v（:（12*（n-1）+1）:12*n））;%將v中的數分為10組，每組12個求和
????vy（:n）=sum1（:n）-6;%vy近似服從N（01）分布
end
figure
n=1:100;
plot（nvy‘b‘）;
xlabel（‘n‘）;ylabel（‘vy‘）;title（‘服從N（01）的白噪聲‘）





????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????730??2018-01-23?16:53??Gaussian.m

?????文件???????1929??2018-01-23?23:50??Minfit.m

?????文件????????585??2018-01-23?20:05??Mnumber.m

-----------?---------??----------?-----??----

?????????????????3244????????????????????3