ls20.m

来自「最小二乘估计」· M 代码 · 共 14 行

clear;
clc;
STLS=importdata('ST.txt');  
ALS=[ones(20,1)';STLS(:,2)';STLS(:,3)';STLS(:,4)']';
LLS=STLS(:,1);
XLS=inv(ALS'*ALS)*ALS'*LLS;
rLS=ALS*XLS-LLS;
det=sqrt(rLS'*rLS/16);

ALS13=[ones(13,1)';STLS(1:13,2)';STLS(1:13,3)';STLS(1:13,4)']';
LLS13=STLS(1:13,1);
XLS13=inv(ALS13'*ALS13)*ALS13'*LLS13;
rLS13=ALS13*XLS13-LLS13;
det=sqrt(rLS13'*rLS13/9);