数值计算ppt.pptVIP

求解大规模方程组Ax b面临的挑战 直接法 　有限次运算后可得方程组精确解（如高斯消元法、 LU方法） 迭代法 　从解的某个近似出发，通过构造无穷序列去逼近精确 解的方法（如Jacobi迭代、Gauss-Seidel迭代、 SOR迭代） 计算量与存储量分析 计算量（乘除法） 　　高斯消元法 O n^3 　　一般迭代法 O n^2*k 　Note：n为未知量个数，k为迭代步数 存储量 　　高斯消元法 O n^2 　　一般迭代法 O n 稀疏阵 ; O n^2 稠密阵 计算机主要技术指标 运算速度 　理论峰值速度 CPU主频×每时钟周期执行浮点运算次数 普通家用3.4GHz P4电脑理论峰值 　3.4 ×2 6.8Gflops 　计算规模n 10^9的线性方程组 　直接法：O 10^17 秒,大约O 10^9 年， 　迭代法：O 10^8*k 秒,若迭代k 10^8,大约O 3 年. 计算机主要技术指标（续） 存储容量 　在64位计算机中，一个双精度浮点占8字节， 整型占4字节。 计算规模n 10^9的线性方程组需要内存 　直接法：O 10^9 GB 　迭代法：O 10 GB 稀疏 、O 10^9 GB 稠密 问题与挑战 全球气候预报要在24小时内完成48小时天气预报，至少需要计算635万个网格点，内存需求大于1TB，计算性能要求高达25Tflops/s。 目前的芯片设计水平远远达不到这样的需求。 高性能计算目的 　能计算规模变更大，算的速度变更快 数值计算方法是研究并解决数学问题的数值近似解的方法，简称计算方法，也叫数值分析等。 算法：从给定的已知量出发，经过有限次四则运算及规定的运算顺序，最后求出未知量的数值解，这样构成的完整计算步骤称为算法。 运算量： 一个算法所需的乘除运算总次数. 计算量是衡量一个算法好坏的重要指标。 算法优劣的重要性 例1 已知 a0, a1, a2 ,…, an, x, 计算多项式： 直接计算：运算量（乘） 秦九韶算法（1247年）： 运算量： 例 2 解线性方程组 其中， ? 克兰姆 Cramer 法则： 运算量 乘除 ： ?高斯消元法 Gauss ： 运算量 乘除 Gauss: 3060次； Cramer: 本课程的任务： ● 建立各种数学问题的数值计算算法的方法和理论。 通俗地讲，就是为各种实际问题提供有效的数值近似解方法。 ● 提供在计算机上实际可行的、理论可靠的、计算复杂 性好的各种常用算法。 计算量：一个算法所需的乘除运算总次数，单位是flop. 计算量是衡量一个算法好坏的重要标准。 可行：能够算出结果。可靠：结果是否正确。复杂性：时间复杂性指节省时间，空间复杂性指节省存储空间。  引 言 2000年8月原武汉大学、武汉水利电力大学、武汉测绘科技大学、湖北医科大学合并组建新武汉大学。新武汉大学汇聚了3所被列入国家“211 工程”重点建设的大学、是名副其实的强强联合、优势互补。 引 言 2000年8月原武汉大学、武汉水利电力大学、武汉测绘科技大学、湖北医科大学合并组建新武汉大学。新武汉大学汇聚了3所被列入国家“211 工程”重点建设的大学、是名副其实的强强联合、优势互补。 * * 学习和了解科学计算的桥梁 科学计算突破了实验与 理论方法的局限性 科学计算已和理论、实验 并列为三大科学方法 There are three great branches of science：theory，experiment， and computation. 我们面对的现象和工程问题 越来越复杂越来越困难 例3 数值天气预报 天气预报方程包括：运动方程（牛顿第二运动定律）;热力学方程（热力学第一定律）;连续方程;水汽守恒方程;气体状态方程和静力方程等等一大串非线性偏微分方程。但只有在计算机出现以后，数值天气预报才成为可能。而在使用并行计算机系统之前，由于受处理能力的限制，只能做到24小时天气预报。高性能计算是解决数值预报中大规模科学计算必要手段。采用高性能计算技术，可以从提高分辨率来提高预报精度。 1992年我国的银河机，10亿/秒，作中期天气预报（5-7天）需要3小时。中期天气预报可以减少灾害的百分之40。作全球大气研究计划，需解50万个方程的方程组。 数值方法和数值软件过去50年的主要进展附录 Before 1940 Newton’s method; Gaussian elimination; least squares fitting; Adams and Runge-Kutta formulas;Richardson extrapolation 1940-1970 floa