第10章计算物理与fortran编程-习题课.ppt-西安交通大学课程中心.pptVIP

随机数的生成-Metropolis算法 Metropolis算法的实现 例3. 用弦割法求方程 在 附近的根( ) 解：取 由迭代公式求得下表 0 0.5 1 0.6 2 0.56532 -0.03468 3 0.56709 0.00178 4 0.56714 0.00001 故 ，满足精度要求． 见Matlab程序 chap10_ex_3_bisection_newton_secant.m （1）利用Matlab编程求解 （2）利用Fortran编程求解 见Fortran程序 chap10_ex_3_bisection_newton_secant.f90 常微分方程边值问题的数值解法 （2） 当 关于 和 是线性时，式（1）为线性两点边值问题 在具体求解常微分方程时，必须附加某种定解条件。定解条件通常有两种，一种是初始条件，另一种是边界条件。与边界条件相应的定解问题称为边值问题。本节介绍求解两点边值问题 （1） 的数值解法。 打靶法 （3） 这里的 为 在 处的斜率。令 ，上述二阶方程可降为一阶方程组 打靶法的基本原理是将两点边值问题（1）转化为下列形式的初值问题 （4） 求根函数F(x) 因此，边值问题变成求合适的 ，使上述方程组初值问题的解满足原边值问题的右端边界条件 ，从而得到边值问题的解。这样，把一个两点边值问题的数值解问题转化为一阶方程组初值问题的数值解问题。方程组初值问题 的所有数值方法在这里都可以使用。问题的关键是如何去找合适的初始斜率的试探值 。 对给定的 ，设初值问题（3）的解为 ，它是 的隐函数。假设 随 是连续变化的，记为 ，于是我们要找的 就是方程 的根。可以用第2章的二分法或迭代法求上述方程的根。比如用弦割法有 （5） 这样，可以按下面简单的计算过程进行求解。先给定两个初始斜率 ，分别作为初值问题（4）的初始条件。用一阶方程组的数值方法求解它们，分别得到区间右端点的函数的计算值 和 。如果 或 ，则以 或 作为两点边值问题的解。否则用割线法（5）求 ，同理得到 ， 再判断它是否满足精度要求 。 如此重复，直到某个 满足 ，此时得到 的 和 就是边值问题的解函数值和它的一阶导数值。上述方程好比打靶， 作为斜率为子弹的发射， 为靶心，故称为打靶法。 打靶法的算法实现 1. 给定初始点的斜率猜测值S1； 2. 用常微分方程的初值问题解法（如RK算法）求解y(b,S1); 3. 给出另一个斜率猜测值S2作为弦割法的第二个启动点; 4. 用常微分方程的初值问题解法（如RK算法）求解y(b,S2); 5. IF (abs(y(b,S1)-y(b))精度 ? ), y(x)=y(x,S1); 6. IF (abs(y(b,S2)-y(b))精度 ? ), y(x)=y(x,S2); else 7. 利用弦割法迭代公式求出下一个斜率S3； 8. IF (abs(y(b,S3)-y(b))精度 ? ), y(x)=y(x,S3); else S1=S2; S2=S3；回到第7步 9. 得到边值问题的解y(x)=y(x,S3); 例 4 用打靶法求解非线性边值问题 要求误差不超过 ，其解析解是 。 解： 对方程进行降阶处理，并转化成初值问题 见Matlab程序 chapter10_ex_4_shooting_method.m （1）利用Matlab编程求解 （2）利用Fortran编程求解 见Fortran程序 chapter10_ex_4_shooting_method.f90 Metropolis (1915-1999) — The Paper was cited 14776 times from 1988 to 2013 虽然前面我们讲述了多种按照规定分布 产生