Distribution of the maximum of correlated receiver horizons
ID: distribution-of-the-maximum-of-correlated-receiver-horizons
The distribution of the maximum of correlated receiver horizons must use the joint timing law. Let H=max(H1,H2). Then Pr(H≤h)=Pr(H1≤h,H2≤h), which generally is not the product of the marginal probabilities. The independent special case is derived in channel switching under uncertainty in a combined receiver horizon.
For an explicit correlated model, suppose (H1,H2) has a bivariate normal distribution with common mean μ, common standard deviation s>0, and correlation coefficient −1<ρ<1. With z=(h−μ)/s, the cumulative distribution function and density of the maximum areHere ϕ,Φ are the standard normal distribution density and cdf, and Φ2 is the bivariate cdf. The density follows by differentiating along both boundaries and using the conditional normal law. The official SciPy reference supplies a cdf method parameterized by the mean vector and covariance matrix; it is a numerical evaluation, not an assumption of channel independence.[ref]
FH(h)=Φ2(z,z;ρ),fH(h)=s2ϕ(z)Φ(z1+ρ1−ρ).
SciPy 1.16.2 documentation. scipy.stats.multivariate_normal, Parameters: mean and cov; Methods: cdf, maxpts, abseps, releps
The equal-variance model also gives exact moments. Write H=(S+∣D∣)/2, with S=H1+H2 and D=H1−H2. They are independent normal variables, with Var(D)=2s2(1−ρ). The absolute-normal moment then givesFor an original shared-clock example, let μ=100ns, s=1ns, and ρ=1/2. Then EH≈100.399ns and the standard deviation is about 0.917ns. Positive correlation reduces the expected gain from selecting the larger horizon, while increasing its variance relative to the independent equal-channel case.
EH=μ+sπ1−ρ,Var(H)=s2(1−π1−ρ).
A common-clock contribution to a maximum receiver horizon can be modeled for nonnegative ρ by Hj=μ+s(ρZ0+1−ρZj) with independent standard normal inputs. As ρ tends to one, both channels share the same uncertain horizon: selection does not remove the common error. Actual receiver timing correlations must be justified by calibration and reference reuse; the example neither measures them nor guarantees that the receiver is invertible.
New to topics? Read the docs here!