combi.m

来自「用matlab实现利用统计混沌方法解决非线性系统时间序列预测的问题」· M 代码 · 共 34 行

function combi(n,p);

cmb = zeros(2,p);

for d = 1:p,
        cmb(1,d) = d;
        cmb(2,d) = n - p + d;
end;
cmb(1,p) = cmb(1,p) - 1;

nbr = 0;

cur = p;

while 0<cur,
        cmb(1,cur) = cmb(1,cur) + 1;
        
        if cmb(1,cur)<cmb(2,cur),
                for i = cur+1:p
                        cmb(1,i) = cmb(1,i-1) + 1;
                end;
                cur = p;
        else,
                cur = cur - 1;
        end;
        
        cmb
        pause;

        nbr = nbr + 1;
end;

nbr
nbr = factorial(n) / factorial(p) / factorial(n-p)