问题的数学模型与算法.pptVIP

湘潭大学信息工程学院 1.1.4 用遗传算法求解布局问题 (1)??? 基本遗传算法的框架 按照Goldberg提出的框架，基本遗传算法(Simple genetic algorithm,SGA)只包括交叉、变异和选择三种基本遗传算子，可用式(1)所描述的一个8元组表示。 SGA＝｛C,E,P0,M,Φ,Γ,Ψ,Τ｝　 (1) 其中C为个体编码方法，E为个体适应度评价函数，P0为初始群体，M为群体规模大小尺寸，Φ为选择算子，Γ为交叉算子，Ψ为变异算子，Τ遗传运算终止条件。 基本遗传算法的主要步骤如下： Step1　随机产生P0个初始种群，并对每一个种群，求出适应度函数的值。 Step2　判断最优个体是否满足要求，如果满足则输出结果后转Step6，否则执行下一步。 Step3　按交叉概率Pc对个体作交叉运算产生新一定量的个体。 Step4　按变异概率Pm对个体作交叉运算产生新一定量的个体，通过交叉和变异后种群的规模为M。 Step5　对M个种群求出其适应度函数的值，按适应度大小选取其中的P0个种后转Step2。 Step6　算法结束。 将式(3)和式(6)代入式(1)并舍去误差项可得到式(7)。 (7) 式(7)的差分方程的逼近误差为O(τ +h 2)，称此逼近关于τ是一阶的，关于h是2阶的。 对初始条件和边界条件作相应逼近逼近，如式(8)。 (8) 式(7)和式(8)构成边值问题式(1)的显式差分格式。 由式(7)和式(8) 可得到差分方程的解式(9)。 (9) 其中　　　。由式(9)可知：第n+1层任意一个内节点处的值Ujn＋1可以由第n层的3个相邻节点处的值Uj-1n, Ujn, Uj+1n决定。由于式(9)的右边不含Ujn＋1，故称之为显示格式。 隐式差分格式　由泰勒展开公式可得式(10)和式(11)。 (10) (11) 由式(10)和式(11)可得到式(12)。 将式(3)和式(12)代入式(1)并舍去误差项可得到式(13)。 (12) (13) 式(13)和式(8)构成边值问题式(1)的隐式差分格式(14)。 (14) 式(14)是关于第n+1层上内节点处的值U1n＋1, U2n＋1, …,UN-1n＋1的线性方程组。我们可以用追赶法或迭代法解此方程组求第n+1层上的值。由于不能直接解出Ujn＋1，故式(14)称之为隐式格式。 六点对称格式　式(7)两边乘以(1-θ)得到式(15)。 (15) 式(14)两边乘以θ得到式(16)。 (16) 式(15)+式(16)得到式(17) 。 (17) 式(17)用到相邻两层六个点处的函数值，故称为六点格式 。 当θ＝1/2时，式(17)被简化为式(18)，它称为六点对称格式。 它可以看作对点(xj,tn＋? )作中心差商的结果。由于 (18) 因此格式(18)的截断误差为O(τ2+h2) ，即对的逼近阶提高了一次。由后面的证明可知：该格式是无条件稳定的。 　　式(18)和式(8)构成边值问题式(1)的六点对称差分格式(19)。 其中：　　　式(19 )对n=0,1,2,…,J可逐次用追赶法求解。 (19) 李查逊格式　由泰勒展开公式有式(20)和式(21)成立。 (20) (21) 式(20)-式(21)的差 除以2τ可得到式(22)。 (22) 由泰勒展开公式可得式(23)和式(24)。 (23) (24) (25) 式(23)+式(24)的和除以h2可得到式(25)。 将式(22)和式(25)代入(1)得到式(26)，其截断误差为O(τ 2 +h 2) 。 (26) (27) 式(27)称为李查孙格式。 前3种格式只要知道第n层上的已知值Ujn通过解方程(方程组)就可以第n+1层上内节点处的值Ujn＋1。而李查孙格式前二层的值Ujn－1 , Uj-1n , Ujn , Uj+1n才能求第n+1层上的值。 2.2　差分方法的稳定性和收敛性 稳定性问题的提出　　对于第一边值问题 其中 采用如下的显式格式求解 对h=π /20,分别取τ1和τ2使r1=τ1 /h2 =5/11, r2=τ2/h2 =5/9,求得的解如图2.2和图2.3所示。其中线段为边界，曲线为精确解，圆点为差分方程解。 图 2.2 h=π /20, r1=5/11时边界和解曲线图 图 2.3 h=π /20, r1=5/9时边界和解曲线图 稳定性定义　　对于给定的ε＞0，总存在与h,τ无关只依赖于ε的正数δ，使当||U0- U0||= ||V0||δ时，依差分格式得到的解Un,Un满足不等式||Un- Un||＝ ||Vn||ε 对于任意的n(