一维搜索法.docxVIP

一维搜索法——黄金分割法【实验目的】掌握一维收搜索中黄金分割法求解一维最优化问题。【实验学时】2学时【算法原理】要介绍黄金分割法有必要回顾一下古老的黄金分割问题．所谓黄金分割就是将一线段分为二段时，要求整段长L与较长段x的比值正好等于较长段x与较短段的比值（如图1所示），即．于是有，解出其正根．由此可见长段的长度应为全长的0.618倍，而短段的长度应为全长的0.382倍．因为古代的人们认为按0.618的比率来分割线段是最协调，胜似黄金，故称之为黄金分割．图 1用黄金分割法进行一维搜索，其基本思想是在单谷区间内适当插入两点，由此把区间分为三段，然后再通过比较这两点函数值大小，就可以确定是删去最左段还是最右段，或者同时删去左右两段保留中间段．如此继续下去可将单谷区间无限缩小。【算法步骤】现在提出一个问题，就在上如何选取二点使得迭代次数最小而区间缩短最快？要解决这个问题，人们想到对区间选二点等价于将区间长度进行黄金分割，也就是将第一个搜索点取在的0.618处，第二个搜索点取成的对称点即的0.382处（如图2所示），亦即，．图 2计算与的值，并根据与的值的大小关系分情况讨论：（1） 若，说明是好点，于是把区间划掉，保留，则内有一保留点，置新的区间；（2）若，说明是好点，于是应将划掉，保留，则内有保留点，置新的区间．（3）若，则应具体分析，看极小点可能在哪一边再决定取舍．在一般情况下，可同时划掉和，仅保留中间的，置新的区间． 接下来是在留下的区间内找好点．重复上面的步骤，直到搜索区间小于给定的允许误差为止．这样就得到黄金分割法迭代算法：已知，常数0.382，终止限．(1)确定的初始搜索区间．(2)计算．(3)计算．(4) 若，则打印，结束；否则转(5)．(5) 判别是否满足：若满足，则置， 然后转(3)；否则，置，然后转(4)．黄金分割法算法流程如图3所示．图 3【程序代码】%黄金分割法%函数原型X=goldenSearch(f,a,b,eps)%返回值为长度为3的向量X，分别表示最小值点、最小值、搜索次数%f:目标函数,a:搜索区间下限,b:搜索区间上限,eps:终止限function X=goldenSearch(f,a,b,eps)p = 1-(sqrt(5)-1)/2;z = 1;t2 = a+p*(b-a);f2 = f(t2);while 1 t1 = a+b-t2; f1 = f(t1); while 1 if abs(t1-t2)eps x = (t1+t2)/2; y = f(x); X=[x y z]; return; end if z1000 disp(搜索次数超过1000次，未找到最小值，结束) x = NaN; y = NaN; return; end z =z+1; if f1f2 a = t2; t2 = t1; f2 = f1; break; else b = t1; t1 = t2; f1 = f2; t2 = a+p*(b-a); f2 = f(t2); end endend