pailie.m

来自「滤波器的设计以及零极点对应频率响应分析」· M 代码 · 共 24 行

function x=pailie(n,k)
%计算排列组合
l=1;
m=1;
j=1;

if n>k & k~=0
      for i=1:n
          l=l*i;
      end
      for i=1:k;
           m=m*i;
       end

      for i=1:(n-k)
          j=j*i;
    
       end

      x=l/(m*j);

else
    x=1;
end