22.01.2014 Views

download searchable PDF of Circuit Design book - IEEE Global ...

download searchable PDF of Circuit Design book - IEEE Global ...

download searchable PDF of Circuit Design book - IEEE Global ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

~----- --------<br />

140 Gradient Optimization<br />

condition for a minimum. By similar triangles,<br />

G5-GI 0,<br />

G5 = a . (5.64)<br />

The increase beyond the 0 step just taken is a- 0" so that<br />

- GI 5<br />

a-o,=o,Z, where Z- G5-GI . (5.6 )<br />

Fletcher limits the extrapolation to be no more than four times the prior step;<br />

Le., Z in (5.65) is limited to 4. The variable names employed correspond to the<br />

program code to follow.<br />

Figure 5.20 shows that a minimum has been bounded in 0 when either<br />

F(o) > F i or when the slope is positive. Suppose that this occurs at 0 = >-.. There<br />

are now four pieces <strong>of</strong> information: the two function values, F(O) = F and<br />

F(>-')=F9; and the two slopes F(0)=G5 and F(>-')=GI. These four items<br />

enable the fit <strong>of</strong> a cubic function, which can interpolate the minimum between<br />

the bounds. The cubic function approximates a flat spring fitted to the known<br />

function values and slopes, provided that the slopes are small. Davidon (1959)<br />

suggested the following formulation, and it has been widely applied since then.<br />

Suppose that the fitting function has the form<br />

Then, at a=;\,<br />

h(a)=ao+a,0+a20 2 +a30 3 . (5.66)<br />

F9= F +G5 ·>-.+a2>-.2+ a3>-.3,<br />

GI =G5+ 2a 2 >-.+ 3a 3 >-.2.<br />

The last two equations can be solved for coefficients a 2 and a 3 :<br />

3(F9-F) ->-.(205 +GI)<br />

a2=<br />

>-.2<br />

2(F-F9)H(G5+Gl)<br />

a3=<br />

>-.3<br />

It is convenient to define the constant z as<br />

3(F- F9)<br />

z >-. +GI+G5.<br />

(5.67)<br />

(5.68)<br />

(5.69)<br />

(5.70)<br />

(5.71 )<br />

The cubic interpolation step in 0 is then obtained by differentiating (5.66) and<br />

equating that to zero. The root <strong>of</strong> the resulting equation that is between 0 = 0<br />

and 0 =>-. is thus obtained after considerable algebra:<br />

A I-(GI+W-z)<br />

0=>-' 2W+GI-G5 '<br />

where an additional defined constant is<br />

(5.72)<br />

W=(z2-G5XGlj'/2. (5.73)<br />

The forms <strong>of</strong> these equations are designed to minimize cancellation by

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!