Mercurial > hg > ltpda
view m-toolbox/test/test_ao_straightLineFit.m @ 44:409a22968d5e default
Add unit tests
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Tue, 06 Dec 2011 18:42:11 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
%mc fs = 1; m = 1.3; c = 4; x = 1:30; yunits = 'Pa'; ca = ao(plist('type', 'cdata', 'yvals', c.*ones(size(x)), 'yunits', yunits)); xa = ao(plist('type', 'cdata', 'yvals', x, 'yunits', yunits)); n = 1*randn(size(x)); sn = 2; nats = sn.* ao(plist('tsfcn','randn(size(t))', 'fs', fs, 'nsecs', length(n), 'yunits', yunits)); ya = m.*xa + ca + nats; % Fit with straightLineFit, pest output yf_p = straightLineFit(ya, plist('output', 'pest')); yf = yf_p.eval(xa, plist('type', 'tsdata')); % Compare fitted data with original iplot(ya,yf) % Fit with straightLineFit, ao output yf_a = straightLineFit(ya, plist('output', 'ao'));