|
wolffd@0
|
1 <html>
|
|
wolffd@0
|
2 <head>
|
|
wolffd@0
|
3 <title>
|
|
wolffd@0
|
4 Netlab Reference Manual linemin
|
|
wolffd@0
|
5 </title>
|
|
wolffd@0
|
6 </head>
|
|
wolffd@0
|
7 <body>
|
|
wolffd@0
|
8 <H1> linemin
|
|
wolffd@0
|
9 </H1>
|
|
wolffd@0
|
10 <h2>
|
|
wolffd@0
|
11 Purpose
|
|
wolffd@0
|
12 </h2>
|
|
wolffd@0
|
13 One dimensional minimization.
|
|
wolffd@0
|
14
|
|
wolffd@0
|
15 <p><h2>
|
|
wolffd@0
|
16 Description
|
|
wolffd@0
|
17 </h2>
|
|
wolffd@0
|
18 <CODE>[x, options] = linemin(f, pt, dir, fpt, options)</CODE> uses Brent's
|
|
wolffd@0
|
19 algorithm to find the minimum of the function <CODE>f(x)</CODE> along the
|
|
wolffd@0
|
20 line <CODE>dir</CODE> through the point <CODE>pt</CODE>. The function value at the
|
|
wolffd@0
|
21 starting point is <CODE>fpt</CODE>. The point at which <CODE>f</CODE> has a local minimum
|
|
wolffd@0
|
22 is returned as <CODE>x</CODE>. The function value at that point is returned
|
|
wolffd@0
|
23 in <CODE>options(8)</CODE>.
|
|
wolffd@0
|
24
|
|
wolffd@0
|
25 <p><CODE>linemin(f, pt, dir, fpt, options, p1, p2, ...)</CODE> allows
|
|
wolffd@0
|
26 additional arguments to be passed to <CODE>f()</CODE>.
|
|
wolffd@0
|
27
|
|
wolffd@0
|
28 <p>The optional parameters have the following interpretations.
|
|
wolffd@0
|
29
|
|
wolffd@0
|
30 <p><CODE>options(1)</CODE> is set to 1 to display error values.
|
|
wolffd@0
|
31
|
|
wolffd@0
|
32 <p><CODE>options(2)</CODE> is a measure of the absolute precision required for the value
|
|
wolffd@0
|
33 of <CODE>x</CODE> at the solution.
|
|
wolffd@0
|
34
|
|
wolffd@0
|
35 <p><CODE>options(3)</CODE> is a measure of the precision required of the objective
|
|
wolffd@0
|
36 function at the solution. Both this and the previous condition must be
|
|
wolffd@0
|
37 satisfied for termination.
|
|
wolffd@0
|
38
|
|
wolffd@0
|
39 <p><CODE>options(14)</CODE> is the maximum number of iterations; default 100.
|
|
wolffd@0
|
40
|
|
wolffd@0
|
41 <p><h2>
|
|
wolffd@0
|
42 Examples
|
|
wolffd@0
|
43 </h2>
|
|
wolffd@0
|
44 An example of the use of this function to find the minimum of a function
|
|
wolffd@0
|
45 <CODE>f</CODE> in the direction <CODE>sd</CODE> can be found in <CODE>conjgrad</CODE>
|
|
wolffd@0
|
46 <PRE>
|
|
wolffd@0
|
47
|
|
wolffd@0
|
48 x = linemin(f, xold, sd, fold, lineoptions);
|
|
wolffd@0
|
49 </PRE>
|
|
wolffd@0
|
50
|
|
wolffd@0
|
51
|
|
wolffd@0
|
52 <p><h2>
|
|
wolffd@0
|
53 Algorithm
|
|
wolffd@0
|
54 </h2>
|
|
wolffd@0
|
55
|
|
wolffd@0
|
56 Brent's algorithm uses a mixture of quadratic interpolation and golden
|
|
wolffd@0
|
57 section search to find the minimum of a function of a single variable once
|
|
wolffd@0
|
58 it has been bracketed (which is done with <CODE>minbrack</CODE>). This is adapted
|
|
wolffd@0
|
59 to minimize a function along a line.
|
|
wolffd@0
|
60 This implementation
|
|
wolffd@0
|
61 is based on that in Numerical Recipes.
|
|
wolffd@0
|
62
|
|
wolffd@0
|
63 <p><h2>
|
|
wolffd@0
|
64 See Also
|
|
wolffd@0
|
65 </h2>
|
|
wolffd@0
|
66 <CODE><a href="conjgrad.htm">conjgrad</a></CODE>, <CODE><a href="minbrack.htm">minbrack</a></CODE>, <CODE><a href="quasinew.htm">quasinew</a></CODE><hr>
|
|
wolffd@0
|
67 <b>Pages:</b>
|
|
wolffd@0
|
68 <a href="index.htm">Index</a>
|
|
wolffd@0
|
69 <hr>
|
|
wolffd@0
|
70 <p>Copyright (c) Ian T Nabney (1996-9)
|
|
wolffd@0
|
71
|
|
wolffd@0
|
72
|
|
wolffd@0
|
73 </body>
|
|
wolffd@0
|
74 </html> |
