Wie man abhängige For-Loops in Matlab vektorisieren

Question

Ich fragte mich, ob jemand mir helfen könnte, diesen Teil meines Codes zu vektorisieren. Hier ist bin_pdf Binomialkoeffizientenfunktion. pb und pd sind Skalarparameter. Danke!

for t=0:min(T,r) 
     for n=0:r-t 
      pp = (bin_pdf(n,pb,r-t) * bin_pdf(t,pd,min(T,r))); 
      pt = pt + t/r * pp; 
      pn = pn + n/r * pp; 
      pc = pc + (r-t-n)/r * pp; 
     end 
    end 

wo

function p = bin_pdf(x,rho,n) 
if (x > n) || (n < 0) 
    p = 0; 
else 
    p = Choosenk(n,x) * rho^x * (1-rho)^(n-x); 
end 

und

function C=Choosenk(n,k) 
if k>n/2 
    k=n-k; 
end; 
C=1; 
for i=0:k-1 
    C=C*(n-i)/(k-i); 
end 
end 

Answer 1

Dies ist die vektorisiert Choosenk ist. Wenn Sie es nicht woanders anrufen, könnten Sie es in das bin_pdf integrieren.

function C=Choosenk(n,k) 
if k>n/2 
    k=n-k; 
end; 
C=prod((n-k+1:n)./(1:k)); 
end