2

我想模拟代表几何分布的值。我使用下面的代码完成的情节似乎产生了正确的情节。但我希望 x 轴严重偏离位置,并且我希望 x 轴编号为 1、2、3 等,而不是我现在收到的 10、20、30 等。我还想将 Y 轴绘制为对数刻度。我正在尝试获取代码中给出的“X”图。

%Geometric Distribution%
N=100;%Number of simulation
P=0.1;
X=zeros(N,1);%simulation data
Ti=0;%Counter

for Ti=2:N
  U=rand(1);
  a=log10(U);
  b=log10(1-P);
  c=(a/b);
  d=1+round(c);
  X(Ti)=d;
  Ti=Ti+1;
end

t = 0:N-1;
hist(X);
4

2 回答 2

3

hist(X,min(X):max(X))

在此处输入图像描述

于 2011-12-05T09:53:30.700 回答
2

用对数比例绘制条形图很棘手。请参阅:如何使用 log scale 绘制 hist,或使用此:

[n, xout] = hist(X,0:max(X));
bar(xout, n, 'barwidth', 1, 'basevalue', 1);
set(gca,'YScale','log')
于 2011-12-05T10:43:24.117 回答