test4.asv

来自「GPS/GALILEO/BEIDOU兼容接收机自主完好性检测程序」· ASV 代码 · 共 26 行

clear all;
close all;
clc;

diata_max = 2.3;
freedom = 4;
tao = 120;
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
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_max = 
arp = sqrt(freedom)*diata_max*Td