test2.m

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

clear all;
close all;
clc;

diata_max = 2.163;
freedom = 4;
tao = 10;
sigma = 33;
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
T_min = sqrt(lamada*sigma*sigma/freedom)
p_ma=ncx2cdf(T,freedom,lamada)
arp = sqrt(freedom)*diata_max*Td