Matlab：確率期待値

Question

公正なコインを10,000回フリップします。 Xを頭の数と尾の数の差に対応させる。

N = 10000; 
i =0; 
r=1/2; 
Q=nchoosek(N,(X+N)/2); 
X=(1,N); 

for i=-N:N 
    P=Q*r.^(X+N)/2*(1-r)^(N-(X+N)/2) %pmf 
    E=sum(abs(X).*P); %expected value 

end 

は、任意の助けを大幅に

Answer 1

をいただければ幸いあなたは1で、単一の行列にすべてのテスト結果を置くことができます。これは私がこれまでにしたものであるX.

の期待値を計算するためにMATLABを使用しますXごとに平均を計算してからXで平均を計算します。

clear 

TAIL=0; HEAD=1; 

NumTests=121; 

NumRollsPerTest=10*1000; 

AllTestsRolls= rand(NumTests,NumRollsPerTest)>0.5 ; %head when rand>0.5 

XperTest=sum(AllTestsRolls==HEAD,2)-sum(AllTestsRolls==TAIL,2);%every row is test so calc per test 

ExpectedX=sum(XperTest)/length(XperTest)