Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function [ alpha, beta ] = myFitExp( x,y )
- len = length(x);
- xSquared = [];
- xy = [];
- for(i = 1 : len)
- xSquared = [xSquared, x(i)^2];
- xy = [xy,x(i)*y(i)];
- end
- xAll = 0;
- yAll = 0;
- xSquaredAll = 0;
- xyAll = 0;
- for(i = 1 : len)
- xAll = xAll + x(i);
- yAll = yAll + y(i);
- xSquaredAll = xSquaredAll + xSquared(i);
- xyAll = xyAll +xy(i);
- end
- m = ((len*xyAll)-(xAll*yAll))/((len*xSquaredAll)-xAll^2);
- b = ((yAll - (m*xAll))/len);
- %y' = mx + b
- %m = beta
- %alfa' = b
- alpha = exp(b);
- beta = m;
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement