From 020f82903d960f379a3b127c0f46d6bf1106cb26 Mon Sep 17 00:00:00 2001 From: Benjamin Franzke Date: Sat, 12 May 2012 15:59:20 +0200 Subject: regress_pred: Add degree param, to fold in quad regression --- regress.m | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) (limited to 'regress.m') diff --git a/regress.m b/regress.m index ad06796..4d13c62 100644 --- a/regress.m +++ b/regress.m @@ -1,4 +1,5 @@ -function mn = regress(x, y) - mn = [ x ones(length(x), 1) ] \ y; +function polynom = regress(x, y, degree) + pot = @(x) x.^[degree:-1:0]; + M = reshape(cell2mat(arrayfun(pot, x(:)', 'UniformOutput', false)), degree+1, [])'; + polynom = M \ y(:); endfunction - -- cgit