计算贝努利(bernoulli)概型的MATLAB函数

%贝努利概型
%
%用法 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中计算的结果:

>> c=bernoulli(3,3,0.8)
c =
    0.5120

>>