برنامه نویسی در متلب درس تعدیل و سرشکنی
clc
clear all
l=[59.92;60;60.05;60;60;59.95]
ql=[1 0 0 0 0 0;0 1 0 0 0 0;0 0 1 0 0 0;0 0 0 1 0 0;0 0 0 0 1 0;0 0 0 0 0 1]
w=ql^(-1);%or(invql)
a=[1 0 0;1 -1 0;0 1 -1;0 0 -1]
f=[l(1)+100;l(2);l(3);l(4)-100]
N=a'*w*a
T=a'*w*f
delta=inv(N)*T
fcap=a*delta;
qdelta=inv(N);
qfcap=a*inv(N)*a'
V=f-fcap