hi there. can some fix error in this code.

1 次查看（过去 30 天）

rana zain 2020-8-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/587009-hi-there-can-some-fix-error-in-this-code

评论： rana zain 2020-8-31

function G = FI_I(I)
%segma_N=10.^-7;
a = 2;     %Alpha
b = 2;     %Beta
c = 0.2;   %Gamma
d = 0.8;   %Big omega prime
k = 1;
z = 4;      %Propagation Distance
%segma_db = 8;  %Atmospheric attenuation coefficient
segma = (8*0.115129255);
%p_db = 0:1:30; %power
%p = 10.^(p_db/10);
%R = 0.5;    %Optoelectronic conversion factor
s = 0.3;    %Zero boresight 
I1 = exp(-segma*z)*0.11512925; %Atmospheric attenuation
Aperture_radius = 0.1;   %Aperture_radius
segma_s = 0.2; %Jitter standard deviation
omega_z = 2.5;   %Beam Width
v = (sqrt(pi)*Aperture_radius)/(sqrt(2)*omega_z);
omega_zeq = sqrt((omega_z.^2*(sqrt(pi)*erf(v)))/(2*v*exp(-v.^2))); 
g = omega_zeq/(2*segma_s);
A0 = (erf(v)).^2;
A = ((2*a.^(a/2))/(c.^(1+a/2)*gamma(a)))*(((c*b)/(c*b+d)).^(b+a/2));
ak = (nchoosek(b-1,k-1))*(((c*b+d).^(1-k/2))/gamma(k))*((d/c).^(k-1))*((a/b).^(k/2));
%PDF for a non-zero boresight pointing error
part_1 = (((g.^2)*A*exp(-s.^2/(2*segma_s.^2)))./((A0*I1).^(g.^2))).*(I.^(g.^2-1));
Q = integral(@integration,2,10);
G = (part_1)*ak*Q;
end

2 个评论
显示无隐藏无

Walter Roberson 2020-8-31

integration is not defined in your call to integral()

rana zain 2020-8-31

https://drive.google.com/drive/folders/1FjJXWrscCk4JVfjeoAsQ-YwJ0OApYK0P?usp=sharing

I have written three matlab file :

integration.m file for the integration inside equation 14 .
FI_I.m file for the equation 14 and call the integration function (integration.m) inside it.
non_zero_o.m file for finally write a code for plotting the equation 16 where I tried to call the function for equation 14

I am trying to plot a figure similar the given Fig.1 for non-zero boresight pointing error (details is given on the pdf) using the equation 14 and 16. And for Zero-boresight pointing error using the equation 13 and 16.

N.B: In zero-boresight (equation 13) there is a function named Meijer for which I am stuck to write the code.

请先登录，再进行评论。

请先登录，再回答此问题。