function X = EvChanAlgorithm(BSN, MSP, R, Noise)
%EvChanAlgorithm 實現改進Chan算法
%EVCHANALGORITHM 本函數用于實現無線定位中的CHAN算法
%               - BSN  為基站個數，3 < BSN <= 7；
%               - MSP  為移動臺的初始位置, MSx, MSy均為[0,1]之間的數；
%                      特別要注意服務小區與MS之間的關系，MS的位置不能越界。
%               - Noise 測距誤差方差。
%               - R    為小區半徑，單位(meter)；
%               - X    為移動臺經算法處理后的位置.
%See also: EvChanAlgorithm.m

%   李金倫，西南交通大學
%       12 October, 2005, 第一版

%   參數檢查：
if  nargout>1,
    error('Too many output arguments.');
end
if nargin ~= 4,
    error('Wrong number of input arguments.');
end

% 算法開始
MS = R*MSP;
BS = R*NetworkTop(BSN);

% MeaDist
for i = 1: BSN,
    MeaDist(i) = sqrt((MS(1) - BS(1,i))^2 + (MS(2) - BS(2,i))^2) + MeaNoise(Noise);
end

% Chan、Taylor估計位置
EMSC = ChanAlgorithm_d(BSN, MSP, R, Noise,MeaDist);
EMSTC = TaylorAlgorithm_d(BSN, EMSC'/R, R, Noise, MeaDist);
EMSTC = EMSTC';

% 殘差
ResC = Residual(MeaDist,EMSC,BSN,R);
ResTC = Residual(MeaDist,EMSTC,BSN,R);

% 估計位置
EMS = (EMSC/ResC + EMSTC/ResTC)/(1/ResC + 1/ResTC);

% 結果輸出
if nargout == 1,
    X = EMS;
elseif nargout == 0,
    disp(EMS);
end