数值分析应用.doc

数值分析应用 数值计算方法是一门工具性、方法性、整合性的新学科，是各种科学与工程计算领域（如：气象、地震、核能技术、石油探勘、航天工程、 密码解译等）中不可缺少的工具。 1、插值和曲线拟合的应用 插值是数值分析领域的一个主要部分,插值理论能解决由已知的表格数值去查找未知的值。例如我们可以结合插值理论建立插值函数进行插值计算,可以得到甘油在某一温度下的粘度。内插和外插在实际预测汽油价格中的比较,得到外插的稳定性、可信性和精度都不如内插。利用实例,通过分段线性插值得到解决画图中的Runge现象的方法。经常由观测或测试可以得到y(x)的一组离散数据。这里y(x)的解析表达式未知，而且这组离散数据并非完全精确，我们需要在给定的函数类上由这组离散数据做出逼近曲线，这就需要用到函数逼近中的曲线拟合。 人口预测：表格中是我国1950年到2005年的人口数（见中国统计年鉴），试预测未来的人口数。这个问题解决需要用到插值或曲线拟合。 年份 人口(万) 1950 55196 1955 61465 1960 66207 1965 72538 1970 82992 1975 92420 1980 98705 1985 105851 1990 11433 1995 121121 2000 126743 2005 130756 2、数值积分与数值微分 函数求导和定积分在微积分中给出了许多有效的方法。但在实际问题中，我们仅给出了函数在一些离散点的值，解析表达式未知或有解析表达式，但难于求得其原函数。这时候，我们需要利用这些离散点上的信息求出函数导数及积分的近似表达，由此需要数值微分和数值积分的概念和方法。 例如：铝制波纹瓦的长度问题 建筑上用的一种铝制波纹瓦是由机器将一块平整的铝板压制而成。假若要求波纹瓦长 4 英尺，每个波纹的高度(从中心线)为 1 英寸，且每个波纹以近似 2p 这个问题就是要求由函数 f (x) = sin x 给定的曲线从 x = 0 到 x = 48 英寸间的弧长 L，即 上述积分为第二类椭圆积分，无法用普通方法来计算，需要利用数值积分的理论。 3、线性方程组的数值解法 在数值计算中，大量存在着线性方程组的求解问题，例如线性问题，非线性问题线性化，微分方程离散后的求解问题等，最后归结为线性方程组的求解问题，插值与曲线拟合等问题也需要求解线性方程组。目前，科学问题常常需要求解高阶线性方程组，必须借助数值计算方法。线性方程组的求解有直接法和迭代法两大类。低阶方程组用直接法比较有效，直接法只需要有限次运算就能得到精确解。而高阶方程组迭代法比较有效，但迭代法一般不能得到精确解。 例：一个古老的数学问题 问：今有 上禾三秉，中禾二秉，下禾一秉，实三十九斗； 上禾二秉，中禾三秉，下禾一秉，实三十四斗； 上禾一秉，中禾二秉，下禾三秉，实二十六斗。 问上、中、下禾实一秉各几何？ ——《九章算术》 4、非线性方程的数值解法 在工程与科学计算中都遇到大量求解非线性方程f(x)=0的问题。非线性科学是当今科学发展的一个重要研究方向，为了理解各种非线性问题的物理机理, 各种非线性模型(非线性方程) 被陆续提了出来。一般而言，很难得到这些非线性方程的准确解析解, 因此, 更多地得到的是应用计算方法得到的数值解.非线性方程组的数值求解是其中不可缺少的内容。求解非线性方程（组）的常用方法是迭代法。 5、矩阵特征值和特征向量计算 在很多工程和科学计算中，需要计算矩阵特征值和特征向量，矩阵A的特征问题是求数和非零向量，是。 例：谷歌搜索是Google公司开发的互联网搜索引擎。主要提供网页搜索，图片搜索，地图搜索，新闻搜索，博客搜索，论坛搜索，学术搜索等。Pagerank算法是当今网络搜索引擎巨匠Google的核心技术。这项技术的一个重要工作是需要计算大矩阵的特征值和特征向量。 例：Google 搜索引擎，1998 年创立，目前市值近2000亿 Gx = x, eTx =1 G: Google Matrix, “the world’s largest matrix computation” x : PageRank vector “The $25,000,000,000 Eigenvector” —— SIAM Review，2006 6、微分方程数值解法 在科学、技术和工程中的大量数学模型都可以用微分方程来描述, 很多近代自然科学的基本方程本身是微分方程。人们用微分方程来解释、预见各种自然现象，不断取得了显著成绩。遗憾的是，绝大多数微分方程定解问题的解不能解析形式来表示。这就产生了理论和应用的