第5讲实践课续龙格库塔方法.docVIP

第四讲（续） 3.3 单步法的收敛性与绝对稳定性.3.1 单步法的收敛性定义.1 设y(x)是初值问题(1.2.1)的精确解，是单步法(3.2.2)在处产生的近似解，若则称方法(3.2.2)产生的数值解收敛于. 实际上，定义中是一固定点，当h→0时n→∞，n不是固定的.因显然方法收敛，则在固定点处的整体误差，当p≥1时下面定理给出方法(3.2.2)收敛的条件定理.1　设初值问题(.2.1)的单步法(3.2.2)是p阶方法(p≥1)，且函数对y满足Lipschitz条件，即存在常数L＞0使对均有　　　　则方法(3.2.2)收敛，且. 定理证明略..3.2 绝对稳定性用单步法(3.2.2)求数值解，由于原始数据及计算过程舍入误差影响，实际得到的不是而是，其中是误差，再计算下一步得到以Euler法为例，若令,则　　　　(.3.1) 如果，则从计算到误差不增长，它是稳定的但如果条件不满足就不稳定例.4 y′=-100y,y(0)=1,精确解为，用Euler法求解得若取h=0.025，则,当，而，显然计算是不稳定的.如果用后退Euler法(3.1.5)解此例，仍取h=0.025,则　　　　　　,即显然当，计算是稳定的.　　由此看到稳定性与方法有关，也与有关，在此例中.在研究方法的稳定性时，通常不必对一般的f(x,y)进行讨论，而只针对模型方程　　　　　　　　　　　　　　(.3.2) 这里可能为复数.规定是因为时微分方程(.3.2)本身是不稳定的，而讨论数值方法(3.2.2)的稳定性，必须在微分方程本身稳定的前提下进行.另一方面，对初值问题(1.2.1)，若将f(x,y)在处线性展开，可得于是方程(1.21)可近似表示为它表明用模型方程(3.3.2)是合理的，至于模型方程(3.3.2)中所以用复数λ是因为初值问题(.2.1)如果是方程组，即，则是(m×m)阶矩阵，其特征值可能是复数.当然对单个方程，λ就是实数，此时只要规定＜0即可.用单步法(3.2.2)解模型方程(3.3.2)可得到　　　　　　　　　　　　　　　　　(3.3.3)其中依赖所选方法，如用Euler法则　　　　　　　　　　(3.3.4)此时由(3.3.1)看到误差方程也为，与(3.3.4)是一样的.因此对一般单步法(3.2.2)误差方程也与(3.3.3)一致.下面再考虑二阶R-K方法有对四阶R-K方法，可得 定义.2　将单步法(3.2.2)用于解模型方程(3.3.2)，若得到(3.3.3)中的 则称方法是绝对稳定的.在复平面上复变量满足 的区域，称为方法(3.2.2)的绝对稳定域，它与实轴的交点称为绝对稳定区间.例如对Euler法， 在复平面上是以(-1，0)为圆心，以1为半径的单位圆域内部，当为实数时，则得绝对稳定区间为，因＜0，故有.在例.4中 时方法稳定，而例中取h=0.025故不稳定.对后退Euler法(3.1.5)，因＜0，故，其绝对稳定域是以(1，0)为圆心的单位圆外部，绝对稳定区间为，即对任何h＞0方法都是绝对稳定的.二阶R-K方法的绝对稳定区间为. 三阶R-K方法的绝对稳定区间为. 四阶R-K方法的绝对稳定区间为. 例.5 用经典四阶R-K方法计算初值问题　　　　　　　　　　步长取h=0.1及0.2,给出计算误差并分析其稳定性.解 本题直接按R-K方法(3.2.12)的公式计算.因精确解为，其计算误差如表所示. 从计算结果看到，h=0.2时误差很大，这是由于在λ=-20,h=0.2时λh=-4，而四阶R-K方法的绝对稳定区间为［-2.785,0］,故h=0.2时计算不稳定，误差很大.而h=0.1时=-2，其值在绝对稳定区间［-2.785,0］内，计算稳定，故结果是可靠的. 讲解： 由于微分方程初值问题数值解公式求出的解是一个逐次递推的过程，因此原始数据误差及计算过程舍入误差对解的影响就是数值方法绝对稳定性研究的问题，如果由计算误差不增长，方法就是绝对稳定的.为使问题得到简化通常就是将方法用于解模型方程（.3.2），对于单步法得到的差分方程为，由于模型方程的，代入Euler法，得，对二阶R-K方法，例如，用改进Euler法于是　　　　　对三阶R-K方法有对四阶R-K方法有只要方法，就是绝对稳定的，这时的值当n增大式是减少的，故计算稳定.这时舍入误差影响可忽略不计，而当，则增大，方法不稳定，计算结果是不可靠的.因此用显式单步法必须使，也就是步长选择要满足这一要求.对于隐式的梯形公式将模型方程，即代入得于是　　注意，于是有　　　　，对成立.这就表明对任意步长h，梯形法都是绝对稳定的..4 练习题 2. 对于一阶微分方程初值问题，用Euler法，改进Euler法，二阶R-K方法求解，并作图比较。 3. 对于一阶微分方程初值问题，用Euler法，改进Euler法，