该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
最小二乘法
%clear all
%xx=[0 0 0 100 100 100 200 200 200];
%yy=[0 100 200 0 100 200 0
100 200];
%zz=[0 50 10 -30 80 -30 10
-40 0];
x=xx;y=yy;z=zz;
sumxz=sum (sum(x.*z));
sumz=sum(sum(z));sumx=sum(x);sumx=sumx(1);sumx2=sum(x.^2);sumx2=sumx2(1);
sumyz=sum(sum(y.*z));sumy=sum(y);sumy=sumy(1);sumy2=sum(y.^2);sumy2=sumy2(1);
p=(sumxz-sumz.*sumx./9)./(9.*sumx2-1.*(sumx.^2))
q=(sumyz-sumz.*sumy./9)./(9.*sumy2-1.*(sumy.^2))
c=sumz./9./9-p.*sumx./9-q.*sumy./9
d=z-p.*x-q.*y;
errorzxerc=max(max(d))-min(min(d))
3.1.3 最小区域法