基本搞懂程序了,但是我可能没表达清楚,仁兄看看这个程序是否对:
function Result=FcJFd(N,M,x,b,c)
if (M>N)||(x>1)||(x<0)
error('输入错误!!!');
end
A=binopdf(M-1,N-1,x);%代表C(M-1,N-1)
d1=N*b/c*A;
g=0;
for i=1:M-1
B=binopdf(M-1,N-1,x);
g=g+B*(i/M-1);
end
d2=-(1-x)^N+1+sum(g);
Result=c/x/N*(d1-d2);