首页 文章

相对频率直方图和概率密度函数

提问于
浏览 次
0

名为DicePlot的函数模拟滚动10个骰子5000次 .

该函数计算每个卷的10个骰子的值的总和，其将是1⇥5000向量，并且绘制相对频率直方图，其中选择区间的边缘，其中直方图中的每个区间表示总和的可能值 . 骰子

将计算1×5000个骰子值总和的平均值和标准偏差，并绘制在相对频率直方图之上的正态分布的概率密度函数（计算平均值和标准偏差） .

下面是我的代码到目前为止 - 我做错了什么？图表显示但顶部没有额外的红线？我看了这样的答案，我认为我不会像高斯函数那样绘制任何东西 .

% function[]= DicePlot()
for roll=1:5000
    diceValues = randi(6,[1, 10]);
    SumDice(roll) = sum(diceValues);
end
distr=zeros(1,6*10);
for i = 10:60
    distr(i)=histc(SumDice,i);
end
bar(distr,1)
Y = normpdf(X)
xlabel('sum of dice values')
ylabel('relative frequency')
title(['NumDice = ',num2str(NumDice),' , NumRolls = ',num2str(NumRolls)]); 
  end

它应该看起来像
enter image description here

但它看起来像

enter image description here

matlab histogram probability dice

2 回答

  • 2

    红线不在那里，因为你没有绘制它 . 查看 normpdf 的文档 . 它计算pdf，它没有绘制它 . 所以问题是如何将这一行添加到图中 . 这个问题的答案是google "matlab hold on" .

    回复于
  • 2

    这里有一些代码可以帮助您朝着正确的方向前进：

    % Normalize your distribution 
normalizedDist = distr/sum(distr);
bar(normalizedDist ,1); 
hold on

% Setup your density function using the mean and std of your sample data
mu   = mean(SumDice);
stdv = std(SumDice);
yy   = normpdf(xx,mu,stdv);
xx   = linspace(0,60);

% Plot pdf
h = plot(xx,yy,'r'); set(h,'linewidth',1.5);
    回复于

相关问题