hi there. can some fix error in this code.
1 次查看(过去 30 天)
显示 更早的评论
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 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!