计算机在化学中应用3.pptVIP

第八章 方程求根 §8.2 二分法 §8.3 迭代法 * 10 READ N, W 20 DIM A(N), B(N), C(N), D(N), E(N), F(N), G(N), H(N), I(N), J(N) 30 FOR I = 1 TO N: READ A(I), B(I), C(I) 32 D(I) = (B(I) + C(I)) / 2 35 NEXT I 40 FOR I = 1 TO N 50 M = D(1): K = 1 60 FOR J = 2 TO N 70 IF D(J) = M THEN 90 80 M = D(J): K = J 90 NEXT J 100 D(K) = W: E(I) = M: F(I) = A(K) 120 NEXT I 125 PRINT XUEHAO, PINGJUN, MINGCI 130 FOR I = N TO 1 STEP -1 132 G(I) = 11 - I 135 PRINT F(I), E(I), G(I) 137 NEXT I D(N)排序前平均成绩； E(N)排序后平均成绩；F(N)学号； G(N)名次。 140 FOR I = 1 TO N 150 M = F(1): K = 1 160 FOR J = 2 TO N 170 IF F(J) = M THEN 190 180 M = F(J): K = J 190 NEXT J 200 F(K) = W: H(I) = M: 210 I(I) = E(K): J(I) = G(K) 220 NEXT I 225 PRINT XUEHAO, PINGJUN, MINGCI 230 FOR I = 1 TO N 235 PRINT H(I), I(I), J(I) 237 NEXT I 240 END 250 DATA 10,10000,8530,92,82,8532,84,100,8542, 79,94,8519,87,92,8503,60,72 260 DATA 8509,67,97,8548,73,96,8533,52,55,8429, 60,62,8442,40,90 计算机程序解方程，只能给出近似解，不能给出解析解。一般分两步进行，第一步求根的初值或存在范围，第二步求精确解。 §8.1 根的初值和存在范围 一般可靠虑以下几种方法。 8.1.1 根据方程的数学性质进行判断。 比如方程： 二次根下的量大于等于0，分母不为0， x-a≥0 b-x＞0 a≤x＜b. 8.1.2根据方程的物理意义进行估计。 可根据方程的物理意义来估计根的值或范围，比如，若方程中变量x代表系统中某种化学物质的摩尔分数，那么必有0≤x≤1。 又如实际气体状态方程 F（P，V，T）=0 是一个复杂的PVT关系（化工热力学？），指定两个量，可用理想气体状态方程 PV=nRT 求第三个量的初值。 8.1.3图解法。 对f（x）=0，直角坐标系中做f（x）→x曲线,与x轴交点就是满足方程的根. 如图： f(x) x f(x) 也可把f(x)=0改成f1(x)=f2(x), f1(x)→x与f2(x)→x分别作图，两曲线交点就是方程的根.图 f(x) x f1(x) f2(x) 8.1.4. 迈步法. f(x) x f(x) a b h x=a，a+h，a+2h……时的函数值 f（a），f（a+h），f（a+2h）……， 直至两个相邻函数不具有相同的符号: f（x）＊f（x+h）≤0 函数f（x）在a，b间连续，且f（a）与f（b）符号不同， 注意h不能选择太大。 如图，若函数f（x）在a，b间连续，且f（a）与f（b）符号不同，则曲线在x=a和x=b之间至少要经过一次横坐标轴，即在a与b间至少有一个实根。 f(x) x f(x) a b h 直至两个相邻函数不具有相同的符号。于是所求的根必定在相应的两个x之间。这种方法实际上就是从左端点x=a出发，按步长h迈步。每迈一步，检查起点x和终点x+h的函数值是否同号，如不同号，即f（x）＊f（x+h）≤0，则所求根必在此x～x+h之间。这种方法称为迈步法，程序在第二节介绍。 应用此法时应注意h不能选择太大。以免在x与x+h间含有两个或两个以上的根。若h过大，在x与x+h间有偶数个根，则f（x）与f（x+h）不异号，这些根会被漏掉。 我们可以选定一个量h，称做步长，然后计算x=a，a+h，a+2h…