绘制误差椭圆的方法
首先计算各点误差椭圆三参数,然后绘制误差椭圆。
(本文以十个点为例,即已知10个点,每个点都有10个误差测量点)
计算代码:
for U=2:10:92
str1=[‘B’,num2str(U),’:’,‘B’,num2str(U+9)];
str2=[‘C’,num2str(U),’:’,‘C’,num2str(U+9)];
x=xlsread(‘data1.xlsx’,str1);
y=xlsread(‘data1.xlsx’,str2);
xx=0.0;
yy=0.0;
for i=1:10
xx=xx+x(i);
yy=yy+y(i);
end
x0=xx/10.0;%x坐标均值
y0=yy/10.0;%y坐标均值
xx1=0.0;
yy1=0.0;
xy1=0.0;
for i=1:10
xx1=xx1+(x(i)-x0)^2;
yy1=yy1+(y(i)-y0)^2;
end
Cx2=xx1/10.0;%x坐标方差
Cy2=yy1/10.0;%y坐标方差
Cxy0=cov(x,y);
Cxy=Cxy0(1,2);%x、y协方差
%E:长半轴、F:短半轴、ct:长半轴E方位角
E=sqrt(Cx2+Cy2+((Cx2-Cy2)^2+4CxyCxy))/2;
F=sqrt