line_nihe.m

来自「这是用于工程测量进行大坝变形分析的软件」· M 代码 · 共 12 行

function [a,b,cigema_squre_line]=line_nihe(x,y)
%设直线为y=ax+b,cigema_squre_line为中误差
for k=1:length(x)
    A(k,1)=x(k);
    A(k,2)=1;
    L(k)=y(k);
end
temp=inv(A'*A)*A'*L';
a=temp(1);
b=temp(2);
V=A*temp-L';
cigema_squre_line=sqrt(V'*V/(length(x)))*1000;