《化工过程解析跟综合》复习题跟习题期末试卷.ppt

《化工过程分析与综合》复习题及习题期末试卷 仅供学习交流，勿用于商业 不独立的列 (f=1)：与 f 值较大的列相比较，若某列中的非零元素与 f 值较大列的非零元素同行，则该列相对于 f 值大的列不独立。 如：S2的f 值较大，与其余小于它的列相比较，会发现S2的非零元素为 C行和A行。而S1列，C行非零。 S3，A行非零。其余列中无与S2同行的非零的元素，则判别出 S1， S3相对于S2不独立。表示为： S1， S3? S2。 S5， S6 ? S4。 据此，可以从某个给定的初始点出发，沿着函数值下降的方向逐步前进（或后退）直至发现函数值开始上升为止。由两边高中间低的三点函数值，就可以确定极小值所在的初始区间[a0，b0] 2、进退算法 (1) 选定初始点a0与步长h (2) 计算并比较y（a0）和y（a0+h），根据比较结果有前进和后退两种可能： ①?前进计算： ②后退运算： 若y（a0+h） ≥ y（a0+3h） ，令a0=a0+h，h=2h，重复上述前进运算。 例： 求函数 的极小所在区间 初始点a0=1，步长h=1 解： h =a0=1 4、黄金分割法的优缺点 （1）、优点：简单、高效，计算次数少 （2）、缺点：对解析性能好的单峰函数，计算量较大。 例：用黄金分割法求下面问题的最优解 min f(x)=e –x+x2 S.t. –2 ≤ x ≤ 3 精度要求:将原始区间缩短5倍. 解：a= -2 , b=3 x1=b- λ(b- a)=3-0.618034*(3+2) = -0.09017 x2=a+λ(b- a)= -2+0.618034*(3+2) =1.09017 f(x1)= 1.10249 f(x2)= 1.5246 f(x1)= 0.9219 f(x2)=1.5246 a= -2, b = 1.09017 x1=b- λ(b- a)=1.09017-0.618034*(1.09017 +2)= -0.81966 x2=a+λ(b- a)= -2+0.618034*(1.09017 +2) = -0.09017 f(x2)= 1.10249 f(x1)=1.1124 f(x2=-0.09017)= 1.10249 f(x1=-0.81966)=1.1124 a= -0.81966 , b = 1.09017 a= -0.81966 , b = 1.09017 例 Hesses矩阵 解这n+m个方程构成的方程组，可得： [x1*， x2* ， x3* …… xn* ， u1* ， u2* …… um*] 其中的[x1*， x2* ， x3* …… xn* ]就是原问题的最优解。 例： 用lagrange乘子法求解 求函数 Y=3+6X1+7X2-7X12+2X1X2-16X22 的最大值，该函数的等值线如图： 解题步骤 （1） 选定初始点x(0)，如图中A 点。A(1.2,-0.2) 并计算该点的函数值y(0)=-2.4 （2） 进行变量轮换。固定一个变量，对另一个寻优。 对此问题求优，可找到最优点为X1=0.3857(B点）， Y(1)=2.0786 再固定X1＝0.3857，对X2寻优。 2 例：用F-R共轭梯度法，求解 解：由于F-R共轭梯度法计算过程中需原函数的梯度信息，所以应先求出梯度函数 第一次迭代，首先确定搜索方向，利用原函数的导数信息 用该方向搜索新点x(1), 将x(1)代入原函数进行一维搜索 k=1 进行第二次迭代，求 代入原函数 解析法可直接解出 可见用此法求解二元二次函数，只要精度足够，两步即达极值点，本题若用梯度法求解，需进行8次迭代。 六．方法评价 优点：1 方法简单，计算及存储量小 2 收敛速度快 缺点：收敛速度依赖于一维搜索的精确性。 作业：用 算法求解 初始点 Min f(x) =60-10 x1-4 x2+ x12 + x22- x1x2 s.t. h (xi)=x1+ x2-8=0 解： lagrange函数L(x,u) L(x,u)= 60-10 x1-4 x2+ x12 + x22- x1x2+ u（x1+ x2-8） 解上面方程组，可得