what is wrong with this code?

Question

Atrolita afra 2020-3-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/514226-what-is-wrong-with-this-code

关闭： MATLAB Answer Bot 2021-8-20

 
functiony=pdf_rician(k,c,sigma)
 
%% pdf of process of various values of system parameters 
 
for x=0:0.0625:10
temp=0; 
pdf_rici=0;
 
for p=0:10
    for n=0:(p+1)
        %    closed term expansion 
        b=nchoosek(p+0.5,n);  
        
        % gamma term calculation 
        f=factorial(p);
        g=(gamma(p+1))*(gamma(c));
        g_ma= 2*sqrt(2)*sqrt(pi)*((x.^(2*p+1))*(k.^c))/((sigma.^(p+0.5))*g*f);   
        
        % bessel function term calculation  
        z=2*sqrt(k*(sigma+x.^2)/sigma);
        K = besselk(p-n+c,z); 
        
        % Exponential term calculation
        exponential=exp(-x.^2/sigma); 
        % Last term 
        q=(p-n+c)/2;
        last_term=((k*sigma/(sigma+(x.^2))).^q);
        
        % Final pdf function 
        
        temp1= b*g_ma*K*exponential*last_term;
        temp=temp+temp1;
    end
    pdf_rici=pdf_rici+temp;
end
 
y(1,x/(.0625)+1)=pdf_rici; 
 
end 