作者：未知 来源：月光软件站 加入时间：2005-5-13　月光软件站

%贝努利概型
%
%用法 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

>>

相关软件：