1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2008 - INRIA
//
// This file is distributed under the same license as the Scilab package.
// =============================================================================
a=[1,7;
2,8
4 3];
b=[10;11;-1];
function y=f1(x,m),y=a*x+b;endfunction
[xsol,v]=lsqrsolve([100;100],f1,3);
lsqrsolve: the cosine of the angle between fvec and any column of the jacobian is at most gtol in absolute value.
if size(xsol) <> [2 1] then bugmes();quit;end
if size(v) <> [3 1] then bugmes();quit;end
function y=fj1(x,m),y=a;endfunction
[xsol,v]=lsqrsolve([100;100],f1,3,fj1);
lsqrsolve: the cosine of the angle between fvec and any column of the jacobian is at most gtol in absolute value.
// Data fitting problem
// 1 build the data
a=34;b=12;c=14;
deff('y=FF(x)','y=a*(x-b)+c*x.*x');
X=(0:.1:3)';Y=FF(X)+100*(rand()-.5);
//solve
function e=f2(abc,m)
a=abc(1);b=abc(2),c=abc(3),
e=Y-(a*(X-b)+c*X.*X);
endfunction
[abc,v]=lsqrsolve([10;10;10],f2,size(X,1));
lsqrsolve: relative error between two consecutive iterates is at most xtol.
if size(abc,'*') <> 3 then bugmes();quit;end
if norm(v) < %eps then bugmes();quit;end
|
