I=I1.signals.values;
w=w1.signals.values;
didt=didt1.signals.values;
dwdt=dwdt1.signals.values;
load?U.mat%讀入電壓

%K=[?L?R??K?J?Ff];

R=zeros（55）;
Ry=zeros（51）;
y=U2（2:）;
o=zeros（1297）;
y=[y;o];
K1=zeros（51）;
for?i=20:length（U2）
????
????W=[?didt（i）?I（i）?????w（i）????????0?????????0?????;
??????????0??????0?????????-I（i）???dwdt（i）????w（i）];
????
???R=R+W‘*W;?%R=W‘*W;
????
????Ry=Ry+W‘*y（:i）;%W‘*y（:i）;
????K=R\Ry;%K=inv（R）*Ry
????K1=[K1?K];
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-02-22?16:25??最小二乘法\
?????文件???????32348??2014-02-22?15:59??最小二乘法\test_1.mdl
?????文件?????????482??2014-02-22?16:27??最小二乘法\test_2m.m
?????文件????????1460??2014-02-17?22:02??最小二乘法\U.mat