Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
update lecture 08
  • Loading branch information
rcc02007 committed Feb 14, 2017
1 parent 4b1d857 commit 5128cfd
Show file tree
Hide file tree
Showing 17 changed files with 3,557 additions and 63 deletions.
Binary file removed lecture_07/.lecture_07.md.swo
Binary file not shown.
36 changes: 36 additions & 0 deletions lecture_08/goldmin.m
@@ -0,0 +1,36 @@
function [x,fx,ea,iter]=goldmin(f,xl,xu,es,maxit,varargin)
% goldmin: minimization golden section search
% [x,fx,ea,iter]=goldmin(f,xl,xu,es,maxit,p1,p2,...):
% uses golden section search to find the minimum of f
% input:
% f = name of function
% xl, xu = lower and upper guesses
% es = desired relative error (default = 0.0001%)
% maxit = maximum allowable iterations (default = 50)
% p1,p2,... = additional parameters used by f
% output:
% x = location of minimum
% fx = minimum function value
% ea = approximate relative error (%)
% iter = number of iterations
if nargin<3,error('at least 3 input arguments required'),end
if nargin<4|isempty(es), es=0.0001;end
if nargin<5|isempty(maxit), maxit=50;end
phi=(1+sqrt(5))/2;
iter=0;
while(1)
d = (phi-1)*(xu - xl);
x1 = xl + d;
x2 = xu - d;
if f(x1,varargin{:}) < f(x2,varargin{:})
xopt = x1;
xl = x2;
else
xopt = x2;
xu = x1;
end
iter=iter+1;
if xopt~=0, ea = (2 - phi) * abs((xu - xl) / xopt) * 100;end
if ea <= es | iter >= maxit,break,end
end
x=xopt;fx=f(xopt,varargin{:});
1,032 changes: 969 additions & 63 deletions lecture_08/lecture_08.ipynb

Large diffs are not rendered by default.

0 comments on commit 5128cfd

Please sign in to comment.