Lec01(survey).pptVIP

科学计算中的基本概念 获取知识的有效途径 笛卡儿说 只有两种方法使人们得到真正的知识 清晰的直觉—体会思想 必要的推理—数学证明 科学计算中的基本概念 收敛性 (or 复杂度) -----误差估计和分析 -----收敛速度 病态性 稳定性 误差 计算地球的表面积 模型误差：地球被看成是一个球 地球的简单理想模型 测量仪器误差和前面的计算误差 地球的半径要经过测量和计算得到 截断误差：公式中的 是无理数 舍入误差：浮点数的计算 浮点数 一般2进制数 L≤s≤U. 而任意的浮点数 其中 浮点数 浮点数 如果取t=4，L=-4，U=3，这时采用对数坐标，则集合Ｆ为 浮点数 复杂度 回忆: 2阶问题, 3阶问题 考虑一般矩阵的行列式 计算需要的乘法次数 复杂度 指数型算法 算法计算量是问题规模的指数函数 只能够处理规模很小的问题 多项式型算法 算法计算量是问题规模的多项式函数 可以处理规模较大的问题 Complexity Descriptor Data Set Size in Bytes Storage Mode Tiny 102 Piece of Paper Small 104 A Few Pieces of Paper Medium 106 A Floppy Disk Large 108 Hard Disk Huge 1010 Multiple Hard Disks Massive 1012 Robotic Magnetic Tape Storage Silos Super-massive 1015 Distributed Data Archives The Huber-Wegman Taxonomy of Data Set Sizes Complexity O( n ) Calculate Means, Variances, O(n log(n)) Calculate Fast Fourier Transforms O( n 2 ) Solve most Clustering Algorithms O( n 3 ) Solve systme of linear equaton O( a n ) TSP Complexity Complexity Complexity 复杂度----对于直接方法的度量标准 Ax=b 的Gauss 消去法 线性规划问题的Simplex方法 组合优化的问题和方法 收敛性----刻划算法的另外一个重要概念 误差 收敛性 病态性-----刻划模型的概念 考虑如下的问题 f(x)=(x-1)(x-2)…….(x-20) 显然方程 f(x)=0 的解是 1 2 3 4 ……… 19 20 请问: 如下方程的解是什么？ Matlab program p=poly(1:20); ep=zeros(1,21); ep(3)=1.0e-7; re=roots(p+ep) plot(re, g+); hold on plot([1:20],zeros(1,20),r*); hold off ?=10e-5 ?=10e-6 ?=10e-8 稳定性-----刻划算法的关键概念 考虑如下的序列 可以证明 两个算法----有什么差别，哪个可以用?? Algorithm 1 Algorithm 2 Program of algorithm 1 clear; close ep(1)=1 for n=2:100 ep(n)=exp(1.0)-n*ep(n-1) end plot(ep,b*); % plot(1:100,log(ep),b*); Algorithm 1with n=15 Algorithm 1 with n=100 Program of algorithm 2 clear ep(100)=0 for n=100:-1:2 ep(n-1)=(exp(1.0)-ep(n))/n; end plot(ep,b*); Algorithm 2with n=100 Algorithm 2 with n=500 科学结论的取得，不能依靠感觉 简单的计算发现，可以使用的算法是-- Algorithm 2！ 计算中误差并不可怕，重要的是误差在算法