[Wx,Wy] = SphericalLogIrrPathWeight(xi,dhat,Rytov,[type])
Computes weighting of Cn2 (turbulence strength) over propagation path for the logirradiance covariance for a spherical wave. Can return the large scale (Wx), small scale (Wy), both (Wx,Wy), or Wx + Wy (single output with type = 'both'). This function uses the integral formulation from Andrews and Phillips and can be computationally intensive.
Parameters 
Description 
xi [vector] 
Normalized location along the propagation path. 
dhat [vector] 
Normalized (by sqrt(L/k)) separation of points in the observation plane 
Rytov [scalar] 
Spherical wave Rytov number from SphericalRytov. 
type [string] 
(Optional) Scale of covariance to output if the number of outputs is one. Ignored for two outputs. Specify 'large' (default), 'small', or 'both' for combined. 
Return Values 
Description 
Wx [vector] 
Large scale path weighting (or Wx + Wy if type is 'both' and single output) evaluated along the path 
Wy [vector] 
(Optional) Small scale path weighting evaluated along the path 
>> Atm = AtmStruct(0, 2e4, 0, 100, 'Cn2', 'HV57');
>> wvl = 0.5e6;
>> dhat = [1 2 3 4 5].*0.02./sqrt(Atm.L*wvl/(2*pi));
>> Rytov = SphericalRytov(1,wvl,Atm);
>> W = SphericalLogIrrPathWeight(Atm.z/Atm.L,dhat,Rytov,'both');
Copyright (c) 2009. All rights reserved.

What do you think about this topic? Send feedback!
