使用验收抑制的正态分布

28次观看(过去30天)
tarek mahjoub.
tarek mahjoub. 于2021年7月30日在6:21
编辑: 亚山于2021年7月30日在18:55
嗨,我想实现以下算法,并在最后绘制X的PDF图。
1.Y1 = -LN(U1),U1是均匀的随机数
2.Y2 = -LN(U2);U2是均匀的随机数
如果Y2≥(y1-1)^ 2/2,则设定| z | = y1;否则返回步骤1。
3.生成U. SET Z = | Z |如果U≤0.5和Z = - | Z |如果你> 0.5
.4。设置x = 0.5 z -2
有没有人有建议如何做到这一点?
1条评论
tarek mahjoub.
tarek mahjoub. 于2021年7月30日在8:57
这是我的尝试。有人可以帮我解决这个代码吗?
n = 100;
z =零(1,n);
absz = zeros(1,n);
x =零(1,n);
功能exp.
Y1 = -log(兰特(1,n));
y2 = -log(rand(1,n));
如果(y1 -1)^ 2/2 <= y2
absz = y1;
别的
exp.
结尾
结尾
u = rand(1,n);
如果你<= 0.5
z = absz;
别的
z = -absz;
结尾
x = z * 0.5 - 2;

登录评论。

答案(1)

亚山
亚山 于2021年7月30日在9:26
编辑:亚山 于2021年7月30日在18:55
这是一个简单的脚本来估计x的pdf。
CLC,清楚
%迭代的数量,这个值越高,pdf越好
%估计
ITT = 1000;
x = nan(ITT,1);
为了n = 1:ITT
尽管1
Y1 = -Log(兰特(1));
y2 = -log(兰特(1));
如果Y2> =(Y1-1)。^ 2/2
absz = y1;
休息
结尾
结尾
u =兰特(1);
如果U> 0.5
z = -absz;
别的
z = absz;
结尾
x(n)= 0.5 * z -2;
结尾
%plot pdf.
直方图(x,'正常化''可能性');
斧头= GCA;
ax.title.string =.'概率分布';
ax.xlabel.string =.'X';
ax.ylabel.string =.'可能性';

社区宝藏狩猎

找到Matlab Central中的宝藏,并发现社区如何帮助您!

开始狩猎!