Propagação de Variância [Ajustamento de Observações] – Código Matlab

% PROPAGAÇÃO DE VARIÂNCIA
% DADOS
A01=60; %azimute
d01=1000; %distância 1
dpA=10; %precisão A
dpd=0.02; %precisão de d
X0=0;
Y0=0;
A01=A01*pi/180;
% CALCULO DAS COORDENADAS
X1=X0+d01*sin(A01);
Y1=Y0+d01*cos(A01);
% propagacao
D=[ sin(A01) d01*cos(A01) ;
cos(A01) -d01*sin(A01) ];

MVCdA=[ dpd^2 0 ;
0 ((dpA/3600)*pi/180)^2];

MVCX1Y1=D*MVCdA*D’;

dpX1=sqrt(MVCX1Y1(1,1));
dpY1=sqrt(MVCX1Y1(2,2));
corX1Y1=MVCX1Y1(1,2)/(dpX1*dpY1);

X1
dpX1
Y1
dpY1
corX1Y1