cep.m

来自「功能极全的GPS开发工具箱」· M 代码 · 共 24 行

%                                cep.m
%  Scope:  This MATLAB macro computes the value of the circular error probable 
%          (CEP).
%  Usage:  xycep = cep(x,y)
%  Description of parameters:
%          x     - input, error vector containing the errors on x axis
%          y     - input, error vector containing the errors on y axis
%          xycep - output, CEP value
%  Remark: The components of the input vectors and the output CEP value 
%          should have the same unit. 
%  Last update:  06/26/00
%  Copyright (C) 1996-00 by LL Consulting. All Rights Reserved.

function  xycep = cep(x,y)

[nrowx,ncolx] = size(x);
[nrowy,ncoly] = size(y);
if  (nrowx ~= nrowy) | (ncolx ~= 1) | (ncoly ~= 1)
   disp('Error  -  CEP; check the dimensions of the input vectors');
   disp('  ');
   return
end 

xycep = 0.615 * std(x) + 0.562 * std(y);