论文查重v1.0part2.docVIP

论文查重v1.0part2 1、相关定义 1.1、-2-1 目标函数的定义 最优化问题的处理,常借助数学而写成以下的形式: .在 xS 的范围内。 .求取 f(x)为最小值时的解 x*。 在上面的表达式里 x 是一个 n 维的向量,而 f(x)则是用来把”好”、” 坏”给量化的纯量函数(Scalar Function,计算结果不是数组的函数),称为目标函数(Objective Function)或代价函数(Cost Function)。如 果原函数的形式是要求取”最大值”通常只要在函数式的前方加一个负号,即可转换为求最小值问题 。 此外,上式中 S 是所有可能解的集合,我们被限制只能在它的范围内寻找最佳解。由于 x 是一个 n 维的向量,XRn,如果 S=Rn,则可以查找的空间没有受到限制,称为没有约束的最优化问题 (Unconstrained Optimization Problem)。如果 S≠Rn,则称为受到约束的最优化问题(Constrained Optimization Problem)。通常约束条件(Constraints)可以写成包括一组等式及一组不等式的数学式: S={x | p(x) = 0, q(x) ≥0 } 其中等式”p(x) = 0”及不等式”q(x)≥0”用来定义 S 的范围。 由于”没有约束的最优化问题”比”受到约束的最优化问题”容易求解,通常我们会以罚函数法 (Method of Penalty Function)先行将后者转换为前者。它的基本做法是:”扩充目标函数,加上能加上约 束条件的项。当约束条件满足时,这些相为零,否则是大于零的值。”也就是说,对于违反约束条件的 答案给予”惩罚”,使得答案往合于约束条件的方向查找。罚函数法按照约束条件的种类可以分为下列 两种做法: (1)等式条件 P(x)=0 可以定义一个罚函数 P(x)满足以下的特性: 当 p(x)=0 时,P(x)=0 当 p(x)≠0 时,P(x)0 最简单的做法是令 p(x)=0 时,P(x)=[ p(x)]2。 (2)不等式条件 q(x)≥0 需要定义一个罚函数 Q(x)满足以下的特性: Fe-1.0Mn-0.1V-0.2C 及 Fe-1.5Al-1.5Cr-1.5C 合金相图计算研究 当 q(x)≥0 时,Q(x)=0 当 q(x)0 最简单的做法是令 Q(x)=[min{ q(x),0}]2。 其中的最小值函数,其性质是: 当 a≤b 时,min{a,b}=a 当 ab 时,min{a,b}=b 一旦等式条件与不等式条件的罚函数 P(x),Q(x)都建立起来后,目标函数 f(x)就可以修正成 F(x): F(x)= f(x)+μ[ P(x)+Q(x) ] 其中正实数 μ 称为惩罚系数(Penalty Parameter),用来调整惩罚函数在目标函数中的比重。 根据以上理论,要求解的问题正是受到约束的最优化问题,应用罚函数法可将方程式转化为如下形 式: 44 ?2 F( x )= f ( x ) +μ1k4=1∑??? i =41∑ f k(i ) x ( ki ) ? f k ??? 2 + μ2 i 4 =1∑ ??? k 4 =1∑ x ki ?1?? (4.23) 按上述做法,可以将等式条件转换成如下的约束方程。x[1]~x[16]所表示的意义已在上一章中说明, x[17]~x[20]则表示某一温度下四个相的量,其值可以由前面 x[1]~x[16]求出,方法也在以下的等式中给 出。由于平时我们表示合金中元素的含量均用质量分数,这里就需要先转换成摩尔分数,其中 W_Fe 表示 Fe 元素在合金中的摩尔百分含量,W_Fe=0.979542,W_Al 表示 Al 元素在合金中的摩尔百分含量, W_Al=0.010105,W_Cr 表示 Cr 元素在合金中的摩尔百分含量,W_Cr=0.001090,W_C 表示 C 元素在 合金中的摩尔百分含量,W_C=0.009263。 q1=x[17]*x[1]+x[18]*x[2]+x[19]*x[3]+x[20]*x[4]- W_Fe q2=x[17]*x[5]+x[18]*x[6]+x[19]*x[7]+x[20]*x[8]- W_Al q3=x[17]*x[9]+x[18]*x[10]+x[19]*x[11]+x[20]*x[12]- W_Cr q4=x[17]*x[13]+x[18]*x[14]+x[19]*x[15]+x[20]*x[16]- W_C p1=x[1]+x[5]+x[9]+x[13]-1 p2=x[2]+x[6]+x[10]+x[14]-1 p3=x[3]+x[7]+x[11]+x[15]-1 p4=x[4]+x[8]+x[12]+x[16]-1 由 x[1]