机械创新实验报告1.docVIP

实验一 一维优化程序的设计、调试与运用 一、实验目的与要求 通过本实验使学生了解常用一维优化方法的基本原理和特点,并通过对某种具体方法的编程调试及验证，加深对该方法基本理论的理解，并培养独立编程的能力。 二、实验方案 用自编优化程序求解下列一维优化问题： minF(X)=4X2 -12X+10 搜索区间：[1,3]，迭代精度E1=10-5 2、minF(X)=X4-5X3+4X2-6X+60 搜索区间：[1,10]，迭代精度E1=10-5 3、minF(X)=(X+1)(X-2)2 搜索区间：[1,10]，迭代精度E1=10-5 4、minF(X)=8X3-2X2-7X+3 搜索区间：[1,10]，迭代精度E1=10-5 程序清单如下： READ A, B, E1 40 NF = 0 50 GOSUB 300 60 PRINT optimum results: 70 PRINT Fmin=; F 80 PRINT X=; X 90 PRINT NF=; NF 100 END 110 DATA 1, 10, 0.00001 300 Subprogram of Golden Section Method 310 S = B - A 320 X1 = A + .618 * (B - A): X = X1: GOSUB 2000: F1 = FX 330 X2 = A + .382 * (B - A): X = X2: GOSUB 2000: F2 = FX 340 IF F1 F2 THEN GOTO 380 350 A = X2: X2 = X1: F2 = F1 360 X1 = A + .618 * (B - A): X = X1: GOSUB 2000: F1 = FX 370 GOTO 400 380 B = X1: X1 = X2: F1 = F2 390 X2 = A + .382 * (B - A): X = X2: GOSUB 2000: F2 = FX 400 IF ABS((B - A) / S) E1 THEN GOTO 340 410 X = .5 * (A + B): GOSUB 2000: F = FX 420 RETURN 2000 Objectve Function 2010 NF = NF + 1 2020 FX = X^4-5*X^3+X^2-6*X+60 2030 RETURN 其中：A———初始搜索区间的下界值； B———初始搜索区间的上界值 E1———迭代精度 实验过程中修改110行DATA后方数据，四次修改依次为： 1、1, 3, 0.00001 2、1, 10, 0.00001 3、-3, 6, 0.00001 4、0, 6, 0.00001 修改2020行FX=后方函数，四次修改依次为： 4*X^2 -12*X+10 2、X^4-5*X^3+4*X^2-6*X+60 3、(X+1)*(X-2)^2 4、8*X^3-2*X^2-7*X+3 三、实验结果及数据处理 1、Fmin= 1 X=1.500121 NF=28 2、Fmin=22.65901 X=3.279831 NF=27 3、Fmin= -0.00000000009822562 X=2.000006 NF=27 4、Fmin= -0.2034246 X=0.6297983 NF=27 1、 搜索区间：[1，3]，迭代精度E1=10-5 计算机输出结果：Fmin=1 X=1.500121 NF=28 2、 搜索区间：[1，10]，迭代精度E1=10-5 计算机输出结果：Fmin=22.65901 X=3.279831 NF=27 3、 搜索区间：[-3，6]，迭代精度E1=10-5 计算机输出结果：Fmin=9.822562*10-11 X=2.000006 NF=27 4、 搜索区间：[0，6]，迭代精度E1=10-5 计算机输出结果：Fmin=-0.2034246 X=0.6297983 NF=27 四、结论 黄金分割法只是一种直接法，利用不停的逼近最优值，在满足迭代精度的条件下就可以把结果当成是最优值。 五、问题与讨论 1、常用一维优化方法有哪些？ 答：常用的一维优化方法有：格点法、 黄金分割法、二次插值法、三次插值法。 2、进退法、黄金分割法、二次插值法基本原理是什么？各有什么特点？ 答：进退法基本原理：取一初始点，若，则在其右方取点；若，则