ls20.m

最小二乘估计

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);