Info

此问题已关闭。 请重新打开它进行编辑或回答。

Need help with a problem

1 次查看(过去 30 天)
Gaëtan Poirier
Gaëtan Poirier 2017-10-2
关闭: MATLAB Answer Bot 2021-8-20
There is a program I have one part I have to modify, could anyone help with this?
function y = run (lambda, n_mea, n_smpls_per_mea)
global xacc hist bin_range bin_size n_bins bin_size_inv midpt bin
% define parameters and procedures for histogram:
bin_range = 6.0 * lambda %NOTE: prefactor can be changed
bin_size = 0.1 * lambda %NOTE: prefactor can be changed
n_bins = floor(bin_range/bin_size)+1
bin_size_inv = 1.0/bin_size
% bin is a function: bin(x) provides the bin number for x:
bin =@(x)(min(floor(x*bin_size_inv)+1,n_bins))
% hist accumulates the number of entries in each bin:
hist = zeros(1,n_bins)
xacc = [n_mea]
for m=1:n_mea
xacc(m)=0
for n=1:n_smpls_per_mea
% unifrnd(0,1) gives a random uniform on (0,1)
ranxi = unifrnd(0,1)
% ranexp gives exponential distribution
ranexp = ranxi
% *^^^^^^^^^ <-- MODIFY -- fill in correct form**
xacc(m) = xacc(m) + ranexp
bi = bin(ranexp)
hist(bi) = hist(bi) +1
end
xacc(m) = xacc(m)/n_smpls_per_mea
end
midpt = 0.5*bin_size:bin_size:(n_bins-0.5)*bin_size
end
  8 个评论
Gaëtan Poirier
Gaëtan Poirier 2017-10-2
Error in mkhist (line 13) exact_pdf = exp(-x/lambda)/lambda;
Walter Roberson
Walter Roberson 2017-10-2
mkhist is not called from the code that you posted.

回答(0 个)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by