?? riddling.m
字號:
%產生一個N(0,1)正態分布隨機數
%采用篩選法,精度較高
%其他方法可參閱《現代應用數學手冊--概率統計與隨機過程卷》清華大學出版社馬振華主編
%function y=riddling()
function y=riddling()
sign=0;
while 1
x=rand(1,2);
v1=2*x(1)-1;v2=2*x(2)-1;
s=v1^2+v2^2;
if s<=1
a=sqrt(-2*log(s)/s);
y=v1*a;
sign=1;
end
if sign==1
break;
end
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -