%贝努利概型
%
%用法 bernoulli(m, n, p)
% m 试验次数
% n 事件发生次数
% p 事件发生的概率
%
%Author 张晓辉
%2005-03-13
%
function c=bernoulli(m, n, p)
if m<1 | round(m)~=m
error('m应该为正整数!');
end
[a b]=size(n);
if a*b==1
if n<0 | round(n)~=n | n>m
error('n应该为不大于m非负整数!');
end
if p<=0 | p>=1
error('p应该为0到1之间的实数!');
end
c=nchoosek(m,n)*p^n*(1-p)^(m-n);
else
error('n暂不支持向量!');
end
以下是在MATLAB中计算的结果:
>> help bernoulli
贝努利概型
用法 bernoulli(m, n, p)
m 试验次数
n 事件发生次数
p 事件发生的概率
Author 张晓辉
2005-03-13
>> c=bernoulli(3,3,0.8)
c =
0.5120
>>
