can anyone check this code?
3 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Atrolita afra
am 2 Apr. 2020
Kommentiert: Ameer Hamza
am 2 Apr. 2020
function l=LCR_CF(k,c,sigma)
%% pdf of process of various values of system parameters
for x=0:.0625:10
temp=0;
LCR=0;
ro=(x.^2/sigma);
ro=20*log10(ro);
for p=0:10
for n=0:p+1
% closed term expansion
b=nchoosek(p+(1/2),n);
% gamma term calculation
f=factorial(p);
g=(gamma(p+1))*(gamma(c));
g_ma=((ro.^(p+0.5))*(k.^c)*(2*sqrt(2*pi)))/(g*f);
% bessel function term calculation
z=2*sqrt(k*(1+ro));
K = besselk(p-n+c,z);
% Exponential term calculation
exponential=exp(-ro);
% Last term
q=(p-n+c)/2;
last_term=((k/(1+ro)).^q);
% Final pdf function
temp1= b*g_ma*K*exponential*last_term;
temp=temp+temp1;
end
LCR=LCR+temp;
end
l(1,x/(.0625)+1)=LCR;
end
THIS IS THE FUNCTION I AM TRYING TO PLOT BUT I AM GATTING AN ERROR ,""Undefined function 'LCR_CF' for input arguments of type 'double'.""
HOW CAN I SOLVE THIS?
I HAVE CALLED THIS FUNCTION BY THE FOLLOWING CODE:
clc
clear all
close all
syms x ;
%% For different values of k and c. Change the value of k and c according to figure.
sig=1;
y1=LCR_CF(0.2,1.5,sig); % Put the values of k,c,sig accordingly
y2=LCR_CF(0.5,1.5,sig);
y3=LCR_CF(1,1.5,sig);
y4=LCR_CF(1.5,1.5,sig);
y5=LCR_CF(.2,2,sig);
y6=LCR_CF(.5,2,sig);
y7=LCR_CF(1,2,sig);
y8=LCR_CF(1.5,2,sig);
%% Ploting data
x=0:0.0625:10;
figure
% plot(y1,x,'b',y2,x,'b',y3,x,'b',y4,x,'b',y5,x,'g--',...
% y6,x,'g--',y7,x,'g--',y8,x,'g--');
plot(x,y1,'b',x,y2,'b',x,y3,'b',x,y4,'b',x,y5,'g--',...
x,y6,'g--',x,y7,'g--',x,y8,'g--');
xlim([0 10]);
ylim([0 1]);
title('Normalized LCR');
xlabel('x');
ylabel('N(x)');
txt = {'Sigma=1','C=1.5, Solid Line', 'C=2, Dotted Line'};
text(6,6,txt);
txt = {'k=1.5'};
text(1,6,txt);
13 Kommentare
Akzeptierte Antwort
Ameer Hamza
am 2 Apr. 2020
I have corrected the syntax issue in your code. The MATLAB function nchoosek does not work for fractional numbers. However, Now it gives infinity at 0 and also gives complex values. I am not sure whether this is correct or wrong.
function l=LCR_CF(k,c,sigma)
%% pdf of process of various values of system parameters
nck = @(n, k) gamma(n+1)./(gamma(k+1).*gamma(n-k+1));
X = 0:.0625:10;
l = zeros(size(X));
for i=1:numel(X)
x = X(i);
temp=0;
LCR=0;
ro=(x.^2/sigma);
ro=20*log10(ro);
for p=0:10
for n=0:p+1
% closed term expansion
b=nck(p+(1/2),n);
% gamma term calculation
f=factorial(p);
g=(gamma(p+1))*(gamma(c));
g_ma=((ro.^(p+0.5))*(k.^c)*(2*sqrt(2*pi)))/(g*f);
% bessel function term calculation
z=2*sqrt(k*(1+ro));
K = besselk(p-n+c,z);
% Exponential term calculation
exponential=exp(-ro);
% Last term
q=(p-n+c)/2;
last_term=((k/(1+ro)).^q);
% Final pdf function
temp1= b*g_ma*K*exponential*last_term;
temp=temp+temp1;
end
LCR=LCR+temp;
end
l(i)=LCR;
end
3 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Logical finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!