test4.m

GPS/GALILEO/BEIDOU兼容接收机自主完好性检测程序

clear all;
close all;
clc;

HAL = 556;
freedom = 8;
tao = 120;
sigma = 10;
P_fa = 0.002/(3600/tao)

T = chi2inv((1-P_fa),freedom)  
T_parity = T*sigma*sigma
Td = sqrt(T_parity/freedom)

% lamada = (556/(sigma*diata_max))^2
lamada = 100;
p_ma=ncx2cdf(T,freedom,lamada);
while(p_ma<0.001)
    p_ma=ncx2cdf(T,freedom,lamada);
    lamada = lamada-0.1;
end
p_ma
lamada
T_min = sqrt(lamada*sigma*sigma/freedom)
diata_ceiling = HAL/(sigma*sqrt(lamada))
arp = sqrt(freedom)*diata_ceiling*Td