red superior

función > Descripción del parámetro:

Ruido: Variación del error de rango.

R: Radio de la celda.

Ver también: Chan_3BS.

Detección de parámetros:

if nargout ~=1,

error('¡Demasiados argumentos de salida!');

end

if nargin ~= 3,

error('¡error de argumentos de entrada!');

fin

El algoritmo comienza

MS = R*MSP ;

BS = R*NetworkTop(3);

Una matriz:

X21 = BS(1,2) - BS(1,1);

X31 = BS(1,3) - BS(1,1);

Y21 = BS(2,2) - BS(2, 1);

Y31 = BS(2, 3) - BS(2, 1)

A = inv([X21, Y21; X31, Y31]); p>

Matriz B:

R1 = sqrt((BS(1,1) - MS(1))^2 (BS(2,1) - MS(2))^2) ;

R2 = sqrt((BS(1,2) - MS(1))^2 (BS(2,2) - MS(2))^2);

R3 = sqrt((BS (1,3) - MS(1))^2 (BS(2,3) - MS(2))^2);

R21 = R2 - R1 RuidoMedio (Es necesario agregar ruido

R31 = R3 - R1 MeaNoise(Ruido)

B = [R21;

Matriz C; :

K1 = BS(1,1)^2 BS(2,1)^2;

K2 = BS(1,2)^2 BS(2,2) ^2;

K3 = BS(1,3)^2 BS(2,3)^2;

C = 0,5*[R21^2 - K2 K1; 2 - K3 K1];

Coeficientes de ecuaciones cuadráticas:

a = B'*A'*A*B - 1;

b = B' *A'*A*C C'*A'*A*B;

c = C'*A'*A*C;

Dos raíces de la ecuación:

raíz1 = abs((-b sqrt(b^2 - 4*a*c))/(2*a));

raíz2 = abs((-b - sqrt (b^2 - 4* a*c))/(2*a));

Pruebe las raíces de la ecuación:

if root1 lt,

EMS = -A* (B*raíz1 C);

else

EMS = -A*(B*raíz2 C);

end

Resultado de salida:

if nargout == 1,

X = EMS

else;

disp(EMS);

fin