最优化上机作业1.docVIP

课程名称： 最优化计算方法 指导教师： 付云姗 班级： 信软2013 姓名： 孙裕道 学号： 1311710119 实验项目名称： 一维优化方法练习 实验目的及要求： 编写程序实现进退法，利用进退法求解极值区间实例； 编写程序实现0.618法，利用0.618法求解极值实例； 编写程序实现牛顿法，来求解一维函数的局部极小值点。 实验原理： 在Matlab环境下，按照要求编写函数和程序，求解实例，直至取得正确的运行结果。 算法1：进退法 S1 给出,h00, 令h:=h0,,,k:=0; S2 令x4=x1+h, 置k:=k+1； S3 若f(x4)f(x1), 则转S4, 否则转S5; S4 令x2=x1,x1=x4,f(x2)=f(x1),f(x1)=f(x4),令 h=2h，转S2； S5 若k=1,则转S6，否则转S7; S6 令h=-h，x2=x4,f(x2)=f(x4)，转S2; S7 令x3=x2,x2=x1,x1=x4，停止计算极小值点包含于[x1,x3]或[x3,x1]. 算法2：0.618法 S1 选定区间[a1,b1]及精度,计算试探点 令 k=1; S2 若,则停止计算。 否则，当时转S3；当时转S4； S3 令, 转 S5; S4令, 转 S5; S5 令k:=k+1，转S2. 算法3：牛顿法 S1 给出精度,令k:=0; S2 若, 停止，极小值点为； S3 令; S4 令k:=k+1, 转S2. 实验内容（方法和步骤）： 题目1 编写程序实现进退法。 利用Matlab编写函数 [xmin,xmax]=Brackeing(f,x0,h0).区间[xmin,xmax]包含函数f的一个局部极值点。 Brackeing.m文件 function [xmin,xmax] = Brackeing(f, x0, h0) h = h0; x1 = x0; k = 0; while 1 k = k + 1; x4 = x1 + h; if f(x4)f(x1) x2 = x1; x1 = x4; h = 2*h; continue; elseif k==1 h = -h; x2 = x4; continue; else x3 = x2; x2 = x1; x1 = x4; break; end end if x1x3 xmin =x1; xmax = x3; else xmin =x3; xmax = x1; end 题目2 利用进退法求解极值区间实例。取初始点x0=0，步长h0=0.1，用进退法求函数 的极值区间。 f.m文件 function y =f(t) y = (t^2 - 1)^2 + (t - 1)^2 +2; end shiyan1_2.m文件 clear x0 = 0; h0 = 0.1; f1 = @f; [xmin, xmax] = Brackeing(f1, x0, h0) 题目3 编写程序实现0.618法。 利用Matlab编写函数 [x,miny]=MinHJ(f,a,b,eps).调用函数得到的x为函数f在区间[a,b]上的局部极小值点。 MinHJ.m文件 function [x, miny] = MinHJ(f, a, b, eps) h = b-a; a0 = a + 0.382 * h; b0 = a + 0.618 * h; while heps if f(a0)f(b0) b = b0;h = b-a; b0 = a0; a0 = a + 0.382 * h; else a = a0; h = b-a; a0 = b0; b0 = a + 0.618 * h; end end x = a0; miny = h; 题目4 利用0.618法求解极值实例。 利用0.618法求下面函数的极小值；，. Shiyan1_4.m文件 clear a = -10; b = 10; eps = 0.1; f1 = @f; [x, miny] = MinHJ(f1, a, b, eps) 题目5 编写程序实现牛顿法。 利用Matlab编写函数 [x,minf]=Newton(f,x0,eps).其中f为目标函数，x0为初始