系统优化设计大作业.docVIP

系统优化设计 PAGE1 系统优化设计算法 指导老师：孙xx 姓名： 李x 学号：02009xxx 日期：2012/11/29 目录 黄金分割法 1.1编程原理………………………………………………………3 1.2程序源代码……………………………………………………3 1.3程序框图………………………………………………………6 1.4运行结果………………………………………………………7 1.5实例考核………………………………………………………9 1.6结论……………………………………………………………9 二次插值法 2.1编程原理………………………………………………………10 2.2程序源代码……………………………………………………10 2.3程序框图………………………………………………………12 2.4运行结果………………………………………………………13 2.5实例考核………………………………………………………14 2.6结论……………………………………………………………14 鲍威尔法 3.1编程原理………………………………………………………15 3.2程序源代码……………………………………………………16 3.3程序框图………………………………………………………22 3.4运行结果………………………………………………………23 3.5实例考核………………………………………………………24 3.6结论……………………………………………………………26 黄金分割法 1.1编程原理 黄金分割法的搜索过程是： 1）给出初始搜索区间【a，b】及收敛精度，将赋以0.618 2）按照坐标点计算公式计算a1和a2，并计算其对应的f（a1），f（a2）。 3）根据区间消去原理缩小搜索区间。为了能用原来的坐标点计算公式，需进行区间名称的代换，并在保留区间中计算一个新的试验点及其函数值。 4）检查区间是否缩短到足够小和函数值是否收敛到足够近，如果条件不满足则返回到步骤2。 5）如果条件满足，则取最后两试验点的平均值作为作为极小点的数值近似解。 1.2程序源代码 syms x disp(请输入函数) f=input(请输入函数\nf=) if isempty(f) f=x.^2+2*x end a=input(请输入函数搜索区间的左极限\na=) b=input(请输入函数搜索区间的左极限\nb=) error=input(请输入精度[1e-3]\nerror=) if isempty(error) error=1e-3 end disp(正在用黄金分割法很努力地计算……) [X_end,Y_end]=GOD(f,x,a,b,error) function [X_end,Y_end]=GOD(f,x,a,b,error) l=0.618; X_end=0; x1=b-l*(b-a); x2=a+l*(b-a); y1=subs(f,x,x1); y2=subs(f,x,x2); while abs((b-a)./b)=error||abs((y2-y1)./y2)=error if(y1=y2) a=x1 x1=x2 y1=y2 x2=a+l*(b-a) y2=subs(f,x,x2); else b=x2 x2=x1 y2=y1 x1=b-l*(b-a) y1=subs(f,x,x1); end X_end=(a+b)./2 end X_end=(a+b)./2 Y_end=subs(f,x,X_end) 开始1.3程序框图 开始 给定 给定a、b、 是 否 是 否 结束 结束 1.4运行结果 最终结果与理论解X=-1，Y=-1十分接近，可以认为是正确的。 表1.4.1是迭代12次的结果 可以看出，函数的最优解收敛很快，6次以后的迭代过程中结果变化很小 表1.4.1 1.5实例考核 附录一维搜索考核题minf(t)=t2-10t+36 用我所编的搜索程序得到的结果时 这与理论值t=5 ，f(t)=11吻合地很好，证明了该程序的通用性。 1.6结论 该黄金分割法搜索程序很好地完成了一维问题最优解的搜索，并通过实例得到了印证。可以认为，该程序的编写是成功的，可以用来作为解决一维最优解搜索问题的通用方法。 二、二次插值法 2.1编程原理 二次插值法又称抛物线法。它是利