最优化问题程序设计方法.ppt

《最优化问题程序设计方法.ppt》由会员分享，可在线阅读，更多相关《最优化问题程序设计方法.ppt（49页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、最优化问题程序设计方法现在学习的是第1页，共49页10.1 最优化问题建模一般步骤最优化问题建模一般步骤 一、建立最优化问题的数学模型一、建立最优化问题的数学模型 工程优化问题的数学模型，是要把工程设计中的问题用数学关系式准确表达出来为达到这些要求，所建立起来的数学模型往往都是很复杂的由于工程设计问题各有其特点，所以数学模型也是多种多样的因此，在工程设计中正确地建立数学模型，不仅是一项艰巨复杂的工作，而且也是解决优化设计问题的关键与前提在很多情况下，建立优化问题的数学模型工作一直是一项重要的研究课题 优化数学模型包括三个内容：变量、目标函数及约束条件它们的基本概念和意义已在第一章做了介绍现在学

2、习的是第2页，共49页 二、选择合适的优化方法二、选择合适的优化方法 各种优化方法都有各自的特点和一定的适用场合根据具体的最优化问题，适当地选择优化方法才会有较好的效果 选择优化方法时，主要考虑的因素是：目标函数的维数与连续性；它的一阶、二阶偏导数是否存在，是否易于求得；约束条件是等式约束，还是不等式约束或两者兼有等不同情况 一般地，对于维数较低的问题应选用结构简单、易于编程的方法对于维数较高的问题，效率就显得十分重要，应选择收敛速度较快的方法，对于求导困难或导数不存在的优化问题应选用直接法现在学习的是第3页，共49页 三、制订流程图和编写源程序三、制订流程图和编写源程序 为了使所编写源程序有

3、个正确的思路，必须先根据具体最优化问题制定一个较详细的流程图该图应反应优化计算的步骤及各种运算之间的逻辑关系流程图既便于程序的编制，又便于使用者对程序的阅读 编写源程序是一种技巧较高而且很细致的工作即使是一个较为简单的最优化问题，也需要考虑许多方面的因素若某些优化方法已有比较成熟的源程序，应尽量优先采用，以期缩短编程时间和提高计算的可靠性与有效性现在学习的是第4页，共49页 一个新编制的程序，即使在编写过程中已经作过周密的考虑，也很难在计算机上一次通过，总会发生这样或那样的障碍，可能是语法规则方面的错误，也可能是运行错误等等因此，新编程序必须经过调试和试算后才能确认它的正确性 试算是必要一环，

4、所谓试算，是用一个比较简单的、已经作好标准答案的题目用编好的源程序运算，观察结果是否正确，以期检查程序的正确性，试算通过后再作正式计算，其结果就比较可信了现在学习的是第5页，共49页 分析优化结果的目的在于考证优化结果的正确性与实用性尽管最优化方法本身是一种科学方法，是可以信赖的但由于实际工程问题的复杂性和某些算法在研究上的不完善性，或由于设计者在建模中失误与疏忽，都会导致计算结果与实际情况不相符，甚至有时是荒谬的所以对优化结果要进行分析如果经分析，发现计算结果存在问题，则需寻查原因，进行调整，修改，直至获得完全符合实际情况为止 最后还需指出，一般情况下通过优化计算所得的最优解只能保证是一个局

5、部最优解只有凸规划问题的局部最优解才是全局最优解为了得到全局最优解，只要多选几个分布在不同位置的初始点进行优化计算若所得各解都归于同一解上去，可认为所得解为全局最优解，否则应从这些解中择其目标函数最小者做为全局最优解现在学习的是第6页，共49页10.2常用最优化方法的特点及选用标准 前几章给大家介绍了各种常用最优化方法，它们在实际使用中都有一定效果，各有其特点，看法也不太一致到底哪种算法较好?哪种较差?为了比较它们的特性，我们有必要建立合理的评价准则 最优化方法评价准则主要有以下几方面：（一）可靠性 所谓可靠性是指算法在合理的精度要求下，在一定允许时间内能解出各种不同类型的最优化问题的成功率能

6、够解出的问题越多，则算法的可靠性越好简言之，如果算法对于不同的问题有时好，有时差，随机性太大，实际使用时没有把握，则可靠性就差现在学习的是第7页，共49页 （二）有效性 这是指算法的解题效率而言的有效性常用两种衡量标准其一是用同一题目，在相同的精度要求和初始条件下，比较占用机时数多少其二是在相同精度要求下，计算同一题目获得最优解时所需要的计算目标函数值次数及导数值次数 （三）简便性 简便性包括两个方面的含意一方面是指实现这种算法人们所需要的准备工作量的大小例如，编制程序的复杂程度，程序调试出错率的高低，算法中所用调整参数的多少等等另一方面是指算法所占用存贮单元的数量，如果某些算法占用单元数很大

7、，就会对机型提出特殊要求，显然这对使用者是不方便的现在学习的是第8页，共49页 由上面的三个评价准则可以看出，要断然地肯定某算法最好或最坏是不可能的因为各种算法就上面三个准则作评价时一般是各有长短，而且由于目标函数的多样性，各种算法对不同目标函数所体现出来的准则衡量结果也有差异因此算法的评价实际上是一个比较复杂的问题 下面分别就常用无约束方法和常用约束方法作一概略评论，以供参考现在学习的是第9页，共49页一、常用无约束最优化方法评价准则一、常用无约束最优化方法评价准则表10.1是几种常用无约束方法的比较现在学习的是第10页，共49页二、常用约束最优化方法评价标准二、常用约束最优化方法评价标准表

8、10.2列出了几种常用约束方法的比较现在学习的是第11页，共49页10.3 最优化问题编程的一般过程 本节将讨论有关最优化问题编程的一般过程问题，它是优化工作的一个重要方面我们所关心的是如何针对不同的优化问题编写供使用者调用方便、界面友好的优化子程序或实用程序库，而不是编写最优化方法程序本身 最优化问题的编程一般分为以下几个阶段：一、准备阶段一、准备阶段 确定任务，选用计算机及所用语言准备原始资料，包括最优化问题的数学模型与原始数据等现在学习的是第12页，共49页 二、编制程序二、编制程序 (1)先画出代表整个计算流程的总框图如有必要再画出有关局部过程的详细框图对较大的最优化问题，它给出了整个

9、构思及各分支的相互联系，是必不可少的一步，但有些人忽视这点当然框图也不宜过于详细，否则，反而主次不分，起不到应有作用 (2)根据框图，按算法语言规定逐条编写设计程序特别要强调的是对一个程序必须多次反复核对，每一个符号都要认真对待程序结构有两种：整体式或积木式，一般简单的程序采用只有一个主程序的整体式结构比较简便明了对于大型复杂的最优化问题编程一般采用模块式（积木式）程序设计方法这样主控程序就可以较短，基本体现出总框图的思路，主控程序中需要用到有关的子程序时，只要调用它就行了模块式程序设计具有：各模块，子程序便于分头编制，分段上机，调试，子程序调用方便，易于检查维护等优点 现在学习的是第13页，

10、共49页 三、调试程序三、调试程序 在计算机上编译检查输入的程序有无符号或语法错误有的是程编了，但输入时打错了，有的则是编写时就错的，可按相应的语法错误表逐个加以改正要注意计算机检查错误是一气完成的，前面一句有错误常会引起后面一大串语法错误，有时只要把第一个错误改正后就解决了，当计算机打印出“有错”时首先要认真检查程序的错误，不要以为是计算机误动作 只有认真地把计算机指出的每一个错误全部改正后才叫程序的“语法通过”或“编译通过”要完成这一步应做到认真仔细、有耐心，以及逐渐建立起来的上机经验现在学习的是第14页，共49页 四、上机试算四、上机试算 程序的语法通过只表示语法上没有问题了，但程序中的

11、公式在数学上有无意义（如个别公式中分母出现零值，或结果太大而超界，开平方根项中出现负数等）或公式本身是否有错？执行路线对否？输入数据对否等都没有检查过因此必须先把一个已知其结果的（靠手算或借用别人已算出的）一套原始数据作为输入，然后将通过程序计算出输出结果与已知的结果对照逐一检查，直到每一项都与原来的已知结果相符才算合格，否则计算结果的可靠性是无法保证的 如果一个程序具有多种功能，则必须对所有功能逐一进行试算，直到符合为止 查出运行错误需有一定经验，行之有效的方法是在程序中插入必要的抽查中间结果的输出语句，以便迅速判断出错位置，缩小搜索范围完成试算这一步后，所编程序方可以初步投入使用现在学习的

12、是第15页，共49页 五、改进程序五、改进程序 在使用过程中还可能会碰到一些事先估计不到的意外情况，如原来准备的数组大小超界，曲线不够查等同时也会发现一些值得改进的地方这就要改 编原来的程序，并可进一步明确这个程序的使用范围经过一段使用检验才能认为程序可以正式交付使用修改程序时要特别注意程序中上下的联系，不要只看局部不顾整体而任意改动，否则会越改越乱一般，程序通用性越强，结构越复杂，出错机会也越多如果编制程序时贪快，考虑不周和不仔细，那么调试工作量会大出好几倍一定要使一个错误也没有，计算机才能算出正确的答案现在学习的是第16页，共49页六、在计算机上求解最优化问题的一些要点六、在计算机上求解最

13、优化问题的一些要点 计算机在求解优化问题时，第一次运算就不收敛，这时怎么办?优化计算失败的原因可能是下述情况中的一种：(1)起点（初始点）在非可行域内，因而一开始时就不运算；(2)在移动了一段距离后，停止在非可行域内；(3)是在可行域内移动，但停止在一个显然是非优化的位置上这些困境可能是由下述的一个或几个因素造成的：(1)起点在非可行域内，问题本身一开始就处于病态；(2)起点在非可行域内，所用的优化方法不能处理这种情况；现在学习的是第17页，共49页(3)对所用特定的优化方法，问题处于病态；(4)问题处于过约束状态，不可能有可行域或约束过紧，可行域过小 为摆脱这种困境，要求程序能打印出目标函数

14、、约束函数和变量的终值这样就可发现失败是属于三种情况中的哪一种 如果属第一、二种情况，运算停止在非可行域内，这可能是上述四种原因中的任一种造成的建议首先从第四条着手从有关技术规范的约束条件中，挑出一个或几个怀疑的对象，将它们放宽也可根据打印出来的终值，找出满足的约束条件，将它们放宽这可能要反复多次，直至运算移入可行域内如果仍不成功，应更换所用的优化方法和初值再进行试算现在学习的是第18页，共49页 如果属第一种情况，可先用一种简单的最优化方法找出一个可行点，然后再用另一种方法求出最优解 如果发生第三情况，运算是在可行域内进行的，但后来停顿了这种情况可能是第三条原因造成的求解的问题对所用的算法处

15、于病态，至少在一定范围内是这样的应该用其它的方法用调整程序参数的方法，如调整步长等，有时会起点作用，但一般是不会奏效的现在学习的是第19页，共49页 在进行迭代计算时，如果数据给得不合适，无论怎样迭代也达不到精度要求遇到此种情况，计算机会算个没完，也不打印出任何结果，如不及时停机，永不休止，这就是称为“死循环”在最优化问题求解时迭代出现“死循环”的可能性最大最简单的解决方法是在程序中专设一些计 算 开 关 变 量，如 W；每 迭 代 一 次，它 自 动 加1，超过规定次数时，如 次，则自动停止迭代，转向下一步或结束，同时打印出故障信息，以便查对现在学习的是第20页，共49页10.4 优化问题设

16、计实例 本节主要举出几个优化设计实例，阐明优化设计的一般步骤，包括设计问题的分析，数学模型的建立，优化方法的选用优化结果分析，可供初学者参考 实例一实例一 圆柱螺旋弹簧的优化设计圆柱螺旋弹簧的优化设计 弹簧的优化设计就是在保证满足工作能力要求的条件下，确定一组设计参数，使某个（或某几个）指标达到最优弹簧的类型很多，这里主要介绍常用的压缩圆柱螺旋弹簧的优化设计对其它类型的弹簧，亦可仿此进行现在学习的是第21页，共49页 一、优化设计数学模型的建立一、优化设计数学模型的建立 （一）设计变量 在压缩圆柱螺旋弹簧的优化设计中，涉及的参数较多，如弹簧钢丝直径d、弹簧中径、弹簧总圈数、弹簧工作圈数、弹簧旋

17、、细长、自由高度、节距等等但是，在这些参数中，有些并不是独立的，它们可以表示成其它参数的函数多数情况是把弹簧钢丝直径，中径及工作圈数作为设计变量，当然，也可根据具体要求的不同，选取其它参数作为设计变量 （二）目标函数 在圆柱螺旋弹簧的优化设计中，可以作为优化目标的项目较多，如要求弹簧在满足工作能力条件下，质量最小或外廓尺寸最小；或在一定空间的限制下能贮存的能量最大，或要求弹簧的支态性能最好等（如果同时要求个项目达到某种指标，那就要建立几个目标函数，这就属于多目标优化设计问题了）现在学习的是第22页，共49页（三）约束条件 弹簧优化设计的约束条件，一般有强度约束、刚度约束、振动稳定性约束及弹簧尺

18、寸约束等 二、设计实例 试设计一压缩圆柱螺旋弹簧，要求其质量最小弹簧材料为65Mn，最大工作载荷Pmax=40N，最小工作载荷为0，载荷变化频率fr=25Hz，弹簧寿命为104h，弹簧钢丝直径d的取值范围为14mm，中径D的取值范围为1030mm，工作圈数n不应小于4.5圈，弹簧旋绕比C不应小于4，弹簧一端固定，一端自由，工作温度为50C，弹簧变形量不小于10mm现在学习的是第23页，共49页解解（一）确定优化设计的目标函数 取弹簧钢丝的直径 d、工作圈数 n及中径 为变量即 （10.1）本题的优化目标是使弹簧质量最小，圆柱螺旋弹簧的质量可表示为 （10.2）其中，表示弹簧材料的密度，对于钢材

19、 ；n表示工作圈数；表示死圈数，常取 ，现取 ；表示弹簧中径(单位：mm)；d表示弹簧钢丝直径(单位：mm)将 已 知 参 数 值 代 入 式(10.1)和 式(10.2)，经整理后得目标函数为 现在学习的是第24页，共49页（二）确定约束条件 根据弹簧性能和结构上的要求，可写出如下的约束条件：（1）强度条件 静强度条件静强度条件是指在最大工作载荷作用下，弹簧不应产生永久变形，即 （10.3）其中，K表示弹簧曲度系数，表示弹簧的许用剪应力，取 经代入式(10.3)整理后，得强度约束条件为 现在学习的是第25页，共49页 疲劳强度条件，载荷循环总次数 ，因，属于承受类载荷的弹簧，应进行疲劳强度计

20、算，要求 ，(10.4)其中 为计算安全系数，(10.5)表示钢丝的脉动循环疲劳极限，取 ；表示弹簧中的最小剪应力，因 ，故 Mpa；表示设计安全系数,取 1.2 将 值代入式（10.4）和式（10.5）整理后得 与静强度条件比较，该条件已包含了静强度条件，故略去静强度约束条件，改用 现在学习的是第26页，共49页2）刚度条件 弹簧刚度应满足变形要求即 ，式中 为弹簧刚度变形量(单位mm），其中，G表示弹簧钢丝的剪切弹性模量，其值与温度有关对于工作温度在50左右时可取 ；表示弹簧最小变形量，且由题可知 故得刚度约束条件为 （3）稳定性条件避免压缩弹簧产生失稳现象，对于一端固定、另一端自由的弹簧

21、，要求 ，式中 为弹簧的自由高度 。现在学习的是第27页，共49页代入整理后，得稳定性约束条件为 （4）不发生撞圈现象 由的计算公式可知，此条件自然得到满足（5）不发生共振现象 ，其中，f 表示弹簧固有频率，;表示外载荷变化频率，故得约束条件 （7）对 的取值限制显然，且应取标准值，即1.0，1.2，1.6，2.0，2.5，3.0，3.5，4.0mm等；（式中50为估取的偏大值）；现在学习的是第28页，共49页 由上可知，此压缩圆柱螺旋弹簧的优化设计是一个三维的约束优化问题，其数学模型为现在学习的是第29页，共49页 （三）优化方法与优化结果 因为维数等于3，是个小型约束优化问题，故可用直接优

22、化方法求解在此选用约束坐标轮换法求解取初始设计参数为 初始质量为 优化后的参数经过圆整为 对应的质量为 现在学习的是第30页，共49页实例二实例二 三相电力变压器的优化设计模型三相电力变压器的优化设计模型一、优化设计的要求和设计变量的选取一、优化设计的要求和设计变量的选取 本节介绍了J.Appelbaum给出的模型，铁心结构为阶梯形三心式铁心，如图10.1所示现在学习的是第31页，共49页 对电力变压器的主要要求是：铁损、铜损和阻抗百分值不得超过某个标准值，此外还有最大允许温升的规定，并使变压器制造费用最小下面所介绍的模型中未将最大允许温升作为不等式约束考虑，但规定的设计容量作为不等式约束之一

23、 优化设计变量取以下几个：心柱直径d 单位为cm 窗口高度h 单位为cm 窗口宽度c 单位为cm 铁轭高度y 单位为cm 心柱磁密B 单位为T 铁轭磁密By 单位为T 电流密度J 单位为A/mm2 给定容量P=A1 单位为W现在学习的是第32页，共49页铁损PFe=A2 单位为W铜损PCu=A3 单位为W阻抗百分值铁心体积计算：铁轭体积为 ，式中 为叠压系数；为结构修正系数故得铁心体积与设计变量关系式为 绕组铜体积计算公式为 ，式中为常数 ；为填铜系数现在学习的是第33页，共49页二、不等式约束的推导二、不等式约束的推导（一）设计容量约束式中，为给定输出功率设铁心柱磁通与铁轭磁通相等，则可写成

24、 （二）铁损约束式中，为给定铁损值（三）铜损约束，式中，为给定的铜损；为常数现在学习的是第34页，共49页（四）阻抗百分值约束 ，式中，及 分别为电阻及电抗百分值，对于集中参数组，有式中，为常数于是可得 现在学习的是第35页，共49页 三、目标函数推导三、目标函数推导 优化目标按制造费用最小考虑，其中包括绕组铜组材料费用 ，硅钢片材料费用，及固定费用 （与设计变量无关），即现在学习的是第36页，共49页四、优化设计模型一般形式及计算举例四、优化设计模型一般形式及计算举例于是，三相电力变压器的优化设计模型为：其中，下面是一个设计例子已知 P=10kVA，f=50Hz，10kV/0.3kV，W，=

25、300W，=4%表10.3给出优化设计值与初始值的比较现在学习的是第37页，共49页可见，优化设计结果与初值相比，制造成本下降11%实例三实例三 优化仿真技术用于深部储层解释优化仿真技术用于深部储层解释 最优化方法用于数字测井解释是80年代地球物理测井解释发展的一种新技术，它为综合应用测井资料评价复杂岩性油气藏开辟了一条新途径现在学习的是第38页，共49页 国外Schlumberger Company在80年代开始这方面的研究，并推出了一套被认为是全球通用的GLOBAL程序系统近年来，我国也开展了这方面的研究工作，取得了一些成果，但至今尚未全面投入生产使用 本文针对GLOBAL方法的不足之处，

26、应用数字仿真技术，并采用大系统寻优组合算法对储层作逐点解释在用理论测井值逼近实际测井值 过程中，由计算机自动仿真调整区域参数，从而克服了静态逼近过程中不能实时调整区域参数等不足 此外，还利用设计变量的对称变换、尺度变换对数学模型进行预处理其次，还利用灵敏度分析方法来控制最优解的稳定性 本系统在微机上的PE-3230机上全面实现，现场处理了几口深井，并给出了最终成果图限于篇幅，笔者重点介绍复杂岩性储层解释数学模型和储层解释质量控制函数的建立以及应用效果分析现在学习的是第39页，共49页 一、大庆长垣地区多矿物模型和响应方程的建立一、大庆长垣地区多矿物模型和响应方程的建立 所谓岩石体积模型就是根据

27、测井方法的探测特性和岩石的各种物质在物理性质上的差异，把岩石体积分成几个部分，然后研究每一部分对岩石宏观物理量的贡献，并把岩石的宏观物理量看成是各部分贡献的总和而多矿物体积模型就是把一个复杂的实际地层简化成为由骨架如石英、长石、方解石，岩块、泥质、储集空间（包括残余 油 气 占 据 的 空 间 ，可 动 油 占 据 的 空 间 ，原状地层水占据的空间，等构成根据地球物理多矿物的模型，结合大庆长垣深部地质特点，建立了13个测井响应方程式和响应误差函数式现在学习的是第40页，共49页（一）响应方程 中子测井（），自然伽马能谱测井得到的钍（Th）和钾（K），响应方程如下：现在学习的是第41页，共49

28、页 其中，参数 表示泥质中子孔隙度，表示残余油气饱和度，表示油气的中子孔隙度，表示泥浆滤液密度值，表示泥浆滤液钍含量值，表示泥浆滤液的钾含量值，表示第i种矿物骨架密度值，表示第i种矿物骨架密度值，表示第i种矿物钍含量值，表示第i种矿物钾含量值()；未知量 表示孔隙度，表示泥质百分含量，表示冲洗带含水饱和度，表示含水饱和度，表示第i种矿物体积百分含量()现在学习的是第42页，共49页（二）响应误差函数由专家经验知 ，的响应误差函数如下式：其中，现在学习的是第43页，共49页 此外，根据多矿物体积模型不可得到密度测井深、浅电阻率测井（LLD）、（LLS）、自然伽马测井（GR），岩性密度测井（PEF

29、），井眼补偿声波测井（BHC），电磁波传播测井（EPT）、自然伽马能谱测井中的铀（URAN）等八种测井响应方程及误差函数式，在此不再赘述 二二、复复杂杂岩岩性性储储层层优优化化解解释释数数学学模模型型和和质质量量控控制制函函数数的的建建立立 在建立大庆长垣深部地层储层评价数学模型时，考虑到由于我们只知道作用在系统上的若干输入量（仪器所测的13条实际测井曲线）和输出量（响应方程计算的理论测井值），而对其系统内部结构仅由专家先验知识大致知一些参数的变化范围和哪些测井曲线是起主要作用外，对其系统内部结构却不清楚，这样的系统又叫“灰色系统”对于灰色系统我们只能建立统计数学模型，该模型的各方程式中的储层

30、参数是不知道，或无法准确计算，为了确定这些参数，必须对模型进行参数辨识现在学习的是第44页，共49页 用最优化方法逐点计算储层参数，就是将各种测井资料，误差函数某些地质因素综合成一个多维的信息系统，既然是个受控系统，那么，经常遇到这样一个问题：如何才能使控制系统在满足一定的约束条件下，使其目标函数达到最小，这就是控制系统数字仿真最优化问题下面给出受控制系统储层评价流程图（图10.2）：现在学习的是第45页，共49页 其执行过程是对于某个给定的采样点 （输入）通过系统作用得输出量测井响应 ，计算输入量 在满足一定受控系统条件下的偏差 ，用 表示目标函数达到最小严格的数学推导为：设实测数据 ，理论

31、计算数据 都是随机变量，在大样本情况下，有 。从而 将上式求和，由此推得：现在学习的是第46页，共49页令即得深部复杂岩性储层优化解释数学模型式中，表示测井曲线条数，表示约束条件个数，表示未知量的个数根据 分布的性质现在学习的是第47页，共49页 当给定一个置信概率 时，最小 的区间为 取 ，由此我们建立起检验数字测井优化解释成果质量的准则函数减小非相关函数 当减小非相关函数Rinc在 之间时，输出未知量值最接近储层物理参数值，从而表明数字测井优化解释成果质量可靠现在学习的是第48页，共49页 三、效果分析三、效果分析 根据我们建立的储层优化解释程序系统CDGL，在大庆PE-3230机上处理三口井，分别给出了三口井700m，500m，420m的最终成果图下面主要就ZS5井来说明，通过分析岩屑录井资料和取芯资料知：从Q2d2J基岩，主要矿物成分为石英、长石、岩块、方解石并含有少量泥质，结合岩性分析，ZS5井储层条件不好，只有局部的砂砾岩发育，或有裂缝发育的层段，才有含油气的可能试油资料表明ZS5井在3247.593303.67m油气显示较好，薄片鉴定资料指出ZS5井的泥质含量一般在10%15%之间变化，现在从CDGL处理的ZS5井成果图（图10.3）可以看出，处理结果与上述实际资料基本一致，从而表明了CDGL系统是可信的现在学习的是第49页，共49页