Main Content

height2range

Convert target height to propagated range

Since R2021b

Description

r = height2range(tgtht,anht,el) returns the propagated range to the target, r, as a function of the target height tgtht, the sensor height anht, and the local elevation angle el assuming a Curved Earth Model with a 4/3 effective Earth radius.

example

r = height2range(tgtht,anht,el,Name=Value) specifies additional inputs using name-value arguments. For example, you can specify a flat Earth model, a curved Earth model with a given radius, or a CRPL Exponential Reference Atmosphere Model with custom values.

example

[r,trueSR,trueEL] = height2range(___,Method="CRPL") also returns the true slant range and the true elevation angle when you specify the Earth model as "CRPL".

Examples

collapse all

Compute the range along the propagated path for a target height of 1 km, an antenna height of 10 meters, and an elevation angle of 2 degrees at the radar. Assume a curved Earth model with a 4/3 effective Earth radius.

r = height2range(1e3,10,2)
r = 
2.7125e+04

Compute the range along the propagated path using the CRPL exponential reference atmosphere. Assume a target height of 1 km, an antenna height of 10 meters, and an elevation angle of 2 degrees at the radar. Additionally, compute the true slant range and the true elevation angle to the target.

[R,SRtrue,elTrue] = height2range(1e3,10,2,Method="CRPL")
R = 
2.7171e+04
SRtrue = 
2.7163e+04
elTrue = 
1.9666

Input Arguments

collapse all

Target height in meters, specified as a nonnegative real-valued scalar or vector. If tgtht is a vector, it must have the same size as the other vector input arguments of height2range. Heights are referenced to the ground.

Data Types: double

Sensor height in meters, specified as a nonnegative real-valued scalar or vector. If anht is a vector, it must have the same size as the other vector input arguments of height2range. Heights are referenced to the ground.

Data Types: double

Local elevation angle in degrees, specified as a real-valued scalar or vector. The local elevation angle is the initial elevation angle of the ray leaving the sensor. If el is a vector, it must have the same size as the other vector input arguments of height2range.

Data Types: double

Name-Value Arguments

Specify optional pairs of arguments as Name1=Value1,...,NameN=ValueN, where Name is the argument name and Value is the corresponding value. Name-value arguments must appear after other arguments, but the order of the pairs does not matter.

Example: Method="CRPL",SurfaceRefractivity=300,RefractionExponent=0.15

Earth model used for the computation, specified as "Curved", "Flat", or "CPRL".

  • "Curved" — Assumes a Curved Earth Model with a 4/3 effective Earth radius, which is an approximation used for modeling refraction effects in the troposphere. To specify another value for the effective Earth radius, use the EffectiveEarthRadius name-value argument.

  • "Flat" — Assumes a Flat Earth Model. In this case, the effective Earth radius is infinite.

  • "CRPL" — Assumes a curved Earth model with the atmosphere defined by the CRPL Exponential Reference Atmosphere Model with a refractivity of 313 N-units and a refraction exponent of 0.143859 km–1. To specify other values for the refractivity and the refraction exponent, use the SurfaceRefractivity and RefractionExponent name value arguments. This method requires el to be positive. For more information, see CRPL Model Geometry.

Data Types: char | string

Effective Earth radius in meters, specified as a positive scalar. If this argument is not specified, height2range calculates the effective Earth radius using a refractivity gradient of –39 × 10–9 N-units/meter, which results in approximately 4/3 of the real Earth radius. This argument applies only if Method is specified as "Curved".

Data Types: double

Surface refractivity in N-units, specified as a nonnegative real-valued scalar. The surface refractivity is a parameter of the CRPL Exponential Reference Atmosphere Model used by height2range. This argument applies only if Method is specified as "CRPL".

Data Types: double

Refraction exponent, specified as a nonnegative real-valued scalar. The refraction exponent is a parameter of the CRPL Exponential Reference Atmosphere Model used by height2range. This argument applies only if Method is specified as "CRPL".

Data Types: double

Output Arguments

collapse all

Propagated range between the target and the sensor in meters, returned as a real-valued scalar or row vector. If r is a vector, it has the same size as the vector input arguments of height2range.

True slant range in meters, returned as a real-valued scalar or row vector. If trueSR is a vector, it has the same size as the vector input arguments of height2range. This argument is available only if Method is specified as "CRPL".

True elevation angle in degrees, returned as a real-valued scalar or row vector. If trueEL is a vector, it has the same size as the vector input arguments of height2range. This argument is available only if Method is specified as "CRPL".

More About

collapse all

Flat Earth Model

The flat Earth model assumes that the Earth has infinite radius and that the index of refraction of air is uniform throughout the atmosphere. The flat Earth model is applicable over short distances and is used in applications like communications, automotive radar, and synthetic aperture radar (SAR).

Given the antenna height ha and the initial elevation angle θ0, the model relates the target height hT and the slant range RT by

hT=ha+RTsinθ0RT=(hTha)cscθ0,

so knowing one of those magnitudes enables you to compute the other. The actual range R is equal to the slant range. The true elevation angle θT is equal to the initial elevation angle.

To compute the ground range G, use

G=(hTha)cotθ0.

Curved Earth Model

The fact that the index of refraction of air depends on height can be treated approximately by using an effective Earth's radius larger than the actual value.

Given the effective Earth's radius R0, the antenna height ha, and the initial elevation angle θ0, the model relates the target height hT and the slant range RT by

(R0+hT)2=(R0+ha)2+RT2+2RT(R0+ha)sinθ0,

so knowing one of those magnitudes enables you to compute the other. In particular,

hT=(R0+ha)2+RT2+2RT(R0+ha)sinθ0R0.

The actual range R is equal to the slant range. The true elevation angle θT is equal to the initial elevation angle.

To compute the ground range G, use

G=R0ϕ=R0arcsinRTcosθ0R0+hT.

A standard propagation model uses an effective Earth's radius that is 4/3 times the actual value. This model has two major limitations:

  1. The model implies a value for the index of refraction near the Earth's surface that is valid only for certain areas and at certain times of the year. To mitigate this limitation, use an effective Earth's radius based on the near-surface refractivity value.

  2. The model implies a value for the gradient of the index of refraction that is unrealistically low at heights of around 8 km. To partially mitigate this limitation, use an effective Earth's radius based on the platform altitudes.

For more information, see effearthradius.

CRPL Exponential Reference Atmosphere Model

Atmospheric refraction evidences itself as a deviation in an electromagnetic ray from a straight line due to variation in air density as a function of height. The Central Radio Propagation Laboratory (CRPL) exponential reference atmosphere model treats refraction effects by assuming that the index of refraction n(h) and the refractivity N decay exponentially with height. The model defines

N=(n(h)1)×106=NseRexph,

where Ns is the atmospheric refractivity value (in units of 10–6) at the surface of the earth, Rexp is the decay constant, and h is the height above the surface in kilometers. Thus

n(h)=1+(Ns×106)eRexph.

The default value of Ns is 313 N-units and can be modified using the SurfaceRefractivity name-value argument in functions that accept it. The default value of Rexp is 0.143859 km–1 and can be modified using the RefractionExponent name-value argument in functions that accept it.

CRPL Model Geometry

When the refractivity of air is incorporated into the curved Earth model, the ray paths do not follow a straight line but curve downward. (This statement assumes standard atmospheric propagation and nonnegative elevation angles.) The true elevation angle $\theta_T$ is different from the initial $\theta_0$. The actual range $R$, which is the distance along the curved path $R'$, is different from the slant range $R_T$.

Given the Earth's radius $R_0$, the antenna height $h_a$, the initial elevation angle $\theta_0$, and the height-dependent index of refraction $n(h)$ with value $n_0$ at $h=0$, the modified model relates the target height $h_T$ and the actual range $R$ by

$$ R=\int_0^{h_T-h_a}{n(h)\,dh}\,
 \left({{1-\left(\frac{\textstyle n_0\cos\theta_0}{\textstyle
 n(h)\left(1+\frac{\textstyle h}{\textstyle
 R_0+h_a}\right)}\right)^2}}\right)^{-1/2}.
$$

When Method is specified as "CRPL", the integral is solved using $n(h)$ from CRPL Exponential Reference Atmosphere Model.

To compute the ground range $G$, use

$$ G=\int_0^{h_T-h_a}\frac{dh}{1+\frac{\textstyle h}{\textstyle
R_0+h_a}}\, \left({{\left(\frac{\textstyle n(h)\left(1+\frac{\textstyle
h}{\textstyle
 R_0+h_a}\right)}{\textstyle
 n_0\cos\theta_0}\right)^2}}-1\right)^{-1/2}.
$$

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using MATLAB® Coder™.

Version History

Introduced in R2021b