Info
This question is closed. Reopen it to edit or answer.
what is wrong with this code?
1 view (last 30 days)
Show older comments
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
2 Comments
Walter Roberson
on 31 Mar 2020
What leads you to say there is a problem with the code? Is there an error message?
Answers (1)
Jamal Nasir
on 31 Mar 2020
the variable temp1 is not single value. it is a vector so need to check the equation
1 Comment
This question is closed.
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!