第16讲 数值算.pptVIP

1 逐项递推算法 逐项递推计算方法可以归结为如下关系式： tk = x*tk-1 k=0，1，2，3 ……n uk = uk-1 + ak*tk t0 = 1 u0 = a0 该算法避免了x的k次幂重复计算，从而提高了整个算法的效率 迭代法步骤 将方程组做如下简单变形： a11 x1 = b1 ― ( + a12 x2 + a13x3 + …… + a1n xn ) a22 x2 = b2 ― ( a21 x1 + a23x3 + …… + a2n xn) ………………………………………………………… ann xn = bn ― ( an1 x1 + an2 x2 + a23x3 + …… ) 测试主函数 int main() { double a[n][n]={1,2,1,-2,2,5,3,-2,-2,-2,3,5,1,3,2,3}; double b[n]={4,7,-1,0}; double x[n]; int i,k; // double a[n][n]={20,2,3,1,8,1,2,-3,15}; // double b[n]={24,12,30}; dinv(a); for(i=0;i=n-1;i++) { x[i]=0.0; for(k=0;k=n-1;k++) x[i]+=a[i][k]*b[k]; } for(i=0;i=n-1;i++) cout\tx[i]; coutendl; return 0; } * * 第16讲 数值计算 掌握多项式的计算方法 掌握多元一次代数方程组的求根计算 掌握非线性方程组的求根计算 掌握积分计算的方法 掌握求逆矩阵的方法 数值数据处理（数值运算） ： 主要是数学计算，即求数值解， 如求方程的根、求函数的定积分等 非数值数据处理（又称非数值运算） ： 最常见的是事务管理领域， 如图书馆管理、银行管理、自动化生产线 管理等 计算机数据处理 设n次多项式的通用形式为： Pn(x) = a0 + a1x + a2x2 + …… + anxn (其中n≥1) 给定一个x的值，求多项式函数Pn(x)的值 。 注意： 在使用计算机运行一个算法时，乘除法 要比加减法消耗时间长，因此在编写算法时， 总是试图减少乘除法的次数。 多项式的计算 double Compute_Poly _1( double a[], int n , double x ) { double result , t ; int i ; t = x ; result = a[0] + a[1] * t ; for ( i = 2 ; i = n ; i ++ ) { t = t * x ; result = result + a[i] * t ; } return result ; } 2 秦九韶算法 将多项式按降次幂写成如下形式： Pn(x) = (……((anx+an–1)x +……+a1)x+a0 仔细观察，不难发现一种计算方法，从里 向外一层一层地计算，即有如下递推关系式： u k = a n （k = n –1，……，2，1，0） u k = u k + 1 * x + a k double Compute_Poly _2( double a[], int n , double x ) { double result ; int i ; result = a[n] ; for ( i = n - 1 ; i = 0 ; i–– ) { result = result * x + a[i] ; } return result ; } 秦九韶算法的运算效率较高，秦九韶算法也 称为霍纳（Horner）方法。 含多个未知元的线性方程组的求解问题，即求 a11