That could be a challenge.
dt = datetime('now')+days(1:10);
T1 = table(dt(:), rand(10,1)+(0.1:0.1:1).');
B = [day(T1.Var1) ones(size(T1.Var1))] \ T1.Var2;
LinFit = [day(T1.Var1) ones(size(T1.Var1))] * B;
figure
plot(T1.Var1, T1.Var2, 'p')
hold on
plot(T1.Var1, LinFit, '-r')
hold off
grid
Depending on the time series, the results could appear a bit unusual.
.
