作为作业的一部分,我被要求编写一个Octave函数,该函数将使用参数 p =生成几何RV X 的 m 伪随机结果0.55 以两种不同的方式:
对于2,我的想法是:
然而,我对第一部分感到很遗憾,并且正在努力通过“直接使用随机变量的cdf”来理解问题的含义。
修改
这是我到目前为止所做的:
function x = generate_geometric_cdf(p, m)
% generate geometric RV X with parameter p, m times using the cdf
emperical = zeros(1,m); % allocate array for empirical results of m simulations
for i = 1:m % iterate over m simulations
u = rand; % generate a random number between 0-1
emperical(i) = geocdf(u, p); % store the outcome of the geometric distribution at u
end
end
编辑2
geoinv
是我需要的功能而不是geocdf
。 geocdf(u, p)
返回值为u
的cdf的概率,参数为p
,而geoinv(u, p)
执行反转并返回给定概率u
的随机值参数p
。