カーブフィティング

連立方程式を解くような記述で,最小自乗法をつかえる。すごい。

% 非多項式項をもつ線形モデル
% y=Xa -> X\y = X\Xa = a
clc; close all;
t = [0 .3 .8 1.1 1.6 2.3]';
y = [0.5 0.82 1.14 1.25 1.35 1.40]'

X = [ones(size(t)) exp(-t) t.*exp(-t)];
a= X\y

T = [0:0.1:2.5]';
Y = [ones(size(T)) exp(-T) T.*exp(-T)]*a;
plot(T,Y,'-',t,y,'o'), grid on