hi there. can some fix error in this code.
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
Réponses (0)
Voir également
Catégories
En savoir plus sur Calculus dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!