典型例题与习题1教材.pptVIP

%% Perform Newton iterations for k=1:maxIter; Z=Z-(f(Z,d)./fprime(Z,d)); end function y=f(x,d); y=(x.^d)-1; end function y=fprime(x,d); y=d*(x.^(d-1)); end 代码片段1: %% Find d roots of unity, and the mask for j=1:d root=exp(2*pi*i/d)^j; % the jth root Mj=abs(Z-root); % distance % Each root gets a unique number in [1,d] mask=(Mj=tol)*j; renderMat=renderMat+mask; end colormap(hsv); % Set the color map imagesc(renderMat) % Render the fractal 代码片段2: 作业 题目1: 您研究领域中的数值分析问题？ * 题目2: 用Newton迭代法法画出最美的图形 标准: 1. 图形美 2. 代码美 要求: 1. m文件 2. 说明文档(word或pdf) */36 《数值分析》典型例题 I 一、二章内容提要 典型例题分析 例题与练习题 实验题介绍 ? ? ? ? 化大为小 化繁为简 化难为易 核心的概念 误差 算法的构造与分析 收敛性 稳定性 复杂度(时间与空间)等 * 有效数字概念 若近似值 x 的绝对误差限是某一位上的半个 单位,该位到 x 的第一位非零数字一共有 n 位,则称近似值 x 有 n 位有效数字。 * 从左向右看第一个非零数 误差限不超过该位的半个单位 n位有效数字 如果x具有n位有效数字, 则相对误差满足: 其绝对误差满足: 如果一个规格化浮点数 则称近似数x具有n位有效数字。 迭代法思想: * 收敛性 收敛速度 Iterate:?To say or do?again?or?again and again 例1.经过四舍五入得出x1=6.1025和x2=80.100,试问它们分别具有几位有效数字? 解： 例2.已知近似数x有两位有效数字,试求其相对误差限。 解：| er(x)|5*10-2 例3.如下近似值的绝对误差限均为0.005,问各近似值有几位有效数值 x1=1.38, x2=-0.0312, x3=0.00086。 例4.二次方程 x2 – 16 x + 1 = 0, 取 求 使具有4位有效数。 解:直接计算 x1≈8 – 7.937 = 0.063 修改算法 4位有效数 计算出的x1 具有两位有效数 例5. 采用迭代法计算 , 取x0 = 7 (k = 0，1，2，……) 若xk具有n位有效数字,求证xk+1具有2n位有效数字。 Ex1：对 是否都有这一性质? 例6. 序列{ yn }满足递推关系 yn = 10yn-1 – 1 (n = 1,2,·····) 若取 y0=21/2 ≈1.41(三位有效数字)。递推计算 y10 时误差有多大？ 思考: 这个计算过程稳定吗？ 例7. 设 y0=28, 按递推公式 yn = yn-1 – (783)1/2/100 (n = 1,2,·····) 计算到y100。若取(783)1/2 ≈27.982(5位有效数字), 试问计算 y100 将有多大的误差？ 例8.设计算球体V允许其相对误差限为 1%,问测量球半径R 的相对误差限最大为多少? 解：由球体计算公式分析误差传播规律 故当球体V 的相对误差限为 1% 时，测量球半径R的相对误差限最大为0.33%。 ? ? 相对误差传播规律 例9. 利用级数 可计算出无理数? 的近似值。由于交错级数的部分和数列Sn在其极限值上下摆动, 试分析为了得到级数的三位有效数字近似值应取多少项求和。 解: 由部分和 只需 ? n 1000时, Sn有三位有效数。 例10.在计算机上对调和级数逐项求和计算 当 n 很大时，Sn 将不随n 的增加而增加。试分析原因。 例11. 证明方程1-x-sinx=0在区间[0,1]上有一根, 使用二分法求误差不大于0.5*10-4的根需要二分多少次？ 提示: f(0)=1,