反分析原理.pdf

《反分析原理.pdf》由会员分享，可在线阅读，更多相关《反分析原理.pdf（13页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、1 反分析的原理和计算方法1 概述地下工程开挖过程中，岩土体性态、水土压力和支护结构的受力状态都在不断变化，采用确定不变的力学参数分析不断变化的体系的力学状态，显然不可能得到预想的效果。利用软件提供的反分析方法以现场位移或内力增量量测值等为依据，借助优化反分析方法确定地层性态参数值，并将可使以这些参数值为输入量算得的测点 位移计算值与实测值相比误差为最小的量作为优化反分析解，尔后将其用作预测计算分析的依据。位移反分析方法可分为 正反分析法 和逆反分析法 两类。后者为正分析的逆过程，计算过程简单，但须先建立求逆公式和编制相应的程序，适用性差。前者为正分析计算的 优化逼近过程，一般通过不断修正未知

2、数的试算值逼近和求得优化解，计算机运作时间虽长，但可利用原有正算程序进行计算，便于处理各种类型的反分析问题，并可用于各类非线性问题的分析，适用性强。地下结构的施工常采用分步开挖、分步支护的方式，其位移、结构内力及岩土层应力等随着施工阶段的变化呈现出一种动态响应过程。因此，有必要将常规的反演分析法与施工模拟过程结合起来，建立一种施工动态反演分析方法。在相同工程及地层条件下，通过利用当前施工阶段量测到的全量或增量信息，来反求地层性态参数和初始地应力参数，进而达到准确预测相继施工阶段的岩土介质和结构的力学状态响应，为施工监控设计提供指导性依据。2 量测信息的种类及表达式在建立的反演分析计算法中,现场

3、量测信息一般用作建立反演计算方程的输入量,因而通常是进行反演计算的主要依据。岩土体在工程施工过程中受到扰动后发生的现象，主要是继续变形和破坏，如果归诸于力学原理，则是岩土体的应力场、应变场、位移场和稳定状态在受到扰动的过程中发生了变化。鉴于受力物体的变形、内力、应力和荷载之间存在依存关系，可以推理如能取得岩土体在受到扰动的过程中发生的应力、应变、内力或位移变化值的量测信息，则可望通过正演计算的逆过程得出初始地应力的量值和作用方向，以及用于描述岩土介质的2 受力变形性态的特性参数。2.1 位移量测信息围岩地层中位移量测分为洞周表面各点的收敛位移量测如拱顶下沉、洞周收敛变形、地表沉降、盾构管片接头

4、相对位移等和围岩域内各点的位移量测，主要为围岩径向多点位移、地表深层沉降、水平位移等。在软土岩土工程中，位移量测主要有地表沉降、围护结构的水平位移、垂直位移、土体测斜、周围建筑物、道路和官线的沉降及水平位移等。位移量又分为绝对位移（相对于不动点）和相对位移（相对于同一测线上的基准测点）两种。2.2 内力量测信息内力量测信息包括扰动应力即由开挖等引起的岩土体应力的变化量和构件（支撑、围护、锚杆及衬砌结构等）轴力、弯矩。其中扰动应力为将来扩展反演量测信息。2.3 压力量测信息压力量测信息包括岩土体内部土压力和结构（喷射混凝土、衬砌、围护结构）与岩土体之间的接触压力两种，为将来扩展反演量测信息。2.

5、4 应变量测信息有开挖引起的应变可分为在洞室壁面上发生的应变和在岩土体内部发生的应变两类。前者称为表面应变，后者称为域内应变。在应变量测中常用的是电阻应变片和千分表，其中前者对量测表面应变和域内应变都适用，后者仅适用于量测表面应变。3 目标函数和适应性函数3.1 目标函数隧道及地下结构施工动态反演过程的量测信息拟采用结构变形、内力及地层水平和垂直变形等，待求未知参数X 可设定为各地层弹性模量和初始地应力参数。关于待求未知量X 的最小二乘目标函数为KiiiiFFwXF10)(（3 1）文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K

6、1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4

7、K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N

8、4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1

9、N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T

10、1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3

11、T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO

12、3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E83 式中：K 为量测信息种类，包括绝对位移、相对位移、结构轴力、弯矩等；iiKjjiKjjjiFFFFF12012*,（3 2）其中：*,jjFF任意两施工阶段测点处对应绝对位移、相对

13、位移、结构轴力或弯矩等的计算值和实测值增量；iK 第 i 种量测信息种类的测点个数；iw 加权常数，一般取iw=1。3.2 适应性函数对于岩土工程的位移优化反分析，在应用遗传算法时，由于目标函数比较小，采用适应性函数)(1)(XFxfitness来区分不同的个体（关于遗传算法，详见下节）。3.3 优化方法反演分析中，优化方法和初始值的选择十分重要，这关系到反演最终能否获得成功（即获得正确合理的反演结果）。3.3.1 单纯形法单纯形法的思想是通过对n维空间上1n顶点的函数值进行比较，通过反射、收缩、延伸来排除函数值最大的点，找到函数值最小的点，并形成新的单纯形，这样逐步逼近极小值点。单纯形是 n

14、 维空间中 n+1 个点构成的体积不为零的多面体，这 n+1 个点称为该单纯形的顶点。顶点的位置由n 维空间中的坐标给出，目标函数f(X)定义于 n维空间中。给定顶点的初值 X1，X2，Xn+1后，可求得顶点处的目标函数值f(Xi)。单纯形形心处的坐标为1111niiXnX（2 1）令 Xh，Xl分别为目标函数值取最大和最小的顶点，单纯形法就是要寻找一个具有较小目标函数值的点来取代顶点Xh，方法是通过三种运算：反射，收缩和文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文

15、档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8

16、文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E

17、8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2

18、E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G

19、2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10

20、G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X1

21、0G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E84 延伸。在反射运算中，新顶点坐标为)(hXXXX（2 2)式中，称为反射系数。在计算目标函数后，如有)()()(hlXfXfXf则以 X替代 Xh构成新的单纯形。如有)()(lXfXf则可以扩大步长，进一步寻找更好的点X)(XXXX（2 3)式

22、中称为扩张系数。这时，对于X点，如有)()(XfXf则以 X置换 Xh，并构成新的单纯形。但是如果有)()(XfXf则以 X置换 Xh并构成新的单纯形。如果对于反射后得到的点X，有)()(iXfXf，ih则新的 Xh将是相应于目标函数f(Xh)和 f(X)中较低者。设该点为Xh，用收缩算法寻找新点)(XXXXhc（2 4)式中，为收缩系数。如有)()(hcXfXf则以 Xc 置换 Xh构成新的单纯形。若)()(hcXfXf则以下式取代单纯形的全部顶点)(21iliiXXXX，i=0，1，n（2 5)文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:

23、CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码

24、:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编

25、码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档

26、编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文

27、档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8

28、文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E

29、8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E85 得到新的单纯形。上式实际上是缩小原来的单纯形，并使最好点仍为缩小后的单纯形的一个顶点。重复上述单纯形的算法，单纯形的尺寸将会不断缩小，直至缩小到指定的精度范围以内。

30、3.3.2 阻尼最小二乘法阻尼最小二乘法在给定参数初值的领域内，把函数通过泰勒级数展开，通过反复迭代逐渐逼近目标函数的极小值，得到参数的最优解，增加阻尼因子，大大改善了系数矩阵的求逆条件，为了进一步减少初始参数的影响、增加解的稳定性以及收敛速度，具体过程和算法如下：假设原方程为：dGx（2 6）式中：G、x、d分别为系数矩阵、参数矩阵和实测数据阵。目标函数：niniiiiuuxfGxdGxdxF1122)()()()(（2 7）式中：iu、iu 分别为位移实测值和有限元计算值；),(21nixxxuu；n为实测值的个数；TmxxxX,21，m为参数个数。2|)(|min)()(min)(min

31、xfxfxfxFT（2 8）nijiijxxfxfxxF1)()(2)(),2,1(mj（2 9）mxxFxxFxxFxF)()()()(21=2mnmmnnxxfxxfxxfxxfxxfxxfxxfxxfxxf)()()()()()()()()(212222111211)()()(21xfxfxfn=)()(2xfxJT（2 10）文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编

32、码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档

33、编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文

34、档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8

35、文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E

36、8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2

37、E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G

38、2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E86 矩阵)(xJ为)(xf在 x处的 Jacobi矩阵。将)(xfi在点)(kx处 Taylor 展开到一次项：)()()()()()(kkikiixxxfxfxf（2 11）nxknkkkxJxfxfxfxf21)()(2)(1)()()()()()()()()(xJxfk（2 12）)(xF)()()()()()()(kkkTkx

39、JxfxJ=0（2 13）得迭代公式：)()()()()()(1)()()()1(kTkkTkkkxfxJxJxJxx)()(kkpx（2 14）为了保证收敛于最优解，减少初值的影响，对（2-16）式进行了改进，增加步长因子得迭代方程：)()()()1(kkkkpxx（2 15）使得)()()()()()(kkkkxFpxF同时不断地调整以改变搜索步长，增加解的稳定性和收敛速度。在（2-8）式中，要求对称半正定矩阵)()()()(kTkxJxJ是非奇异的，由于)(xfi的复杂的非线性，这一要求并不总能满足，造成)()()()(kTkxJxJ是病态的或接近病态的，导致收敛速度极慢或计算终止，为此

40、，进行了改进，增加阻尼因子，增大矩阵)()()()(kTkxJxJ的主对角线元素，迭代方程为：)()()()()()(1)()()()()1(kTkkTkkkkxfxJIxJxJxx（2 16）Jacobi矩阵元素的求解，用有限差分代替一阶导数：)()(kxjixxf=jkjkijkjkixxfxxf)()()()()()(（2 17）根据对称矩阵)()()()(kTkxJxJ的正交分解，可以分解为：)()()()(kTkxJxJ=TRR（2 18）文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10

41、 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I1

42、0 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I

43、10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1

44、I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U

45、1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3

46、U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q

47、3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E87 为)()()()(kTkxJxJ的特征值构成的对角线矩阵，而R为)()()()(kTkxJxJ的特征向量矩阵，且满足IRRRRTT（2 19）)()()()(kTkxJxJ+I=TRR+I=TRR（2 20）r00

48、21（2 21）i为)()()()(kTkxJxJ的特征值，矩阵)()()()(kTkxJxJ+I的条件数为：cond)()()()(kTkxJxJ+I=minmaxminmax=cond)()()()(kTkxJxJ（2 22）3.3.3 遗传算法遗传算法是模拟自然进化过程搜索全局最优解的方法。遗传算法的优越性主要表现在它在搜索过程中不容易陷入局部最优解，即使在所定义的适应函数是不连续的、非规则的或有噪声的情况下，它也能以很大的概率找到全局最优解。遗传算法象撒网一洋，在参变量空间中进行搜索，由串组成的群体在遗传算子的作用下，同时对空间中不同的区域进行采样计算，从而构成一个不断变化的群体序列。

49、为了避免陷入局部最优，在遗传算法中还引入了变异，一方面可以在当前解附近找到更好的解；另一方而还可以保持群体的多样性，确保群体继续进化。为了寻找最优解，传统方法是用启发式策略，在单个猜测解的邻域探寻，即使算法中允许偶尔地跳到解空间中更远的部分，这些启发式算法也往往趋向于陷入局部最优。通过保持在解空间不同区域中多个点的搜索，遗传算法以很大概率找到全局最优解。遗传算法中，控制参数的不同对遗传算法的性能产生很大的影响，要想得到遗传算法执行的最优性能，必须确定最优的参数设置。（1）群体规模npopsiz群体规模影响到遗传算法的最终性能和效率。当规模太小时，造成群体的样文档编码:CO3T1N4K1D2 H

50、S4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2 HS4A10Q3U1I10 ZS2Y7X10G2E8文档编码:CO3T1N4K1D2