数值分析李庆杨版
数值分析第五版_李庆扬_王能超_易大义主编课后习题答案

第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求nx 的相对误差。
解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=Q , 1||n p x nx C n n-⋅∴==又((*))(*)r p r x n C x εε≈⋅Q 且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯ 解:*1 1.1021x =是五位有效数字;*20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===g g (*)(*)3(*)r p r r V C R R εεε∴≈=g又(*)1r V ε=Q故度量半径R 时允许的相对误差限为1(*)10.333r R ε=⨯≈6.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=-Q10099Y Y ∴=9998Y Y =9897Y Y =-……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =-,27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯ 100Y ∴的误差限为31102-⨯。
《数值分析-李庆杨》第1章 引言

结束
将[0,1]分为4等分,分别计算4个小曲边梯形的面积的 近似值,然后加起来作为积分的近似值(如图1-1).记被积 函数为 f(x) ,即 f ( x ) 1 1 x
1 h , xi ih, i 0,1,2,3 4
f ( xi ) f ( xi 1 ) Ti h 2
结束
1.2.2 误差的基本概念 定义1 绝对误差,简称误差: e* x * x, 其中x为准确值 , x *为准确值x的近似值 .
. 误差限: * | e* | 的一个上界
例如,毫米尺 765 x 0.5
在工程中常记为:x= x*± * .
如
l=10.2±0.05mm
,R=1500±100Ω
在科学研究、工程实践和经济管理等工作中,存在大 量的科学计算、数据处理等问题.应用计算机解决数值计算 问题是理工科研究生应当具备的基本能力.
§1.1
算法
解决某类数学问题的数值方法称为算法.为使算法能 在计算机上实现,它必须将一个数学问题分解为有限 次的+、-、×、÷运算和一些简单的基本函数运算.
结束
§ 1.2 误差
1.2.1 误差的来源 在运用数学方法解决实际问题的过程中,每一步都可 能带来误差. 1、模型误差 在建立数学模型时,往往要忽视很多次要因 素,把模型“简单化”,“理想化”,这时模型就与真实 背景有了差距,即带入了误差.
2、测量误差 数学模型中的已知参数,多数是通过测量得 到 . 而测量过程受工具、方法、观察者的主观因素、不可 预料的随机干扰等影响必然带入误差.
例3 多项式求值。
2 P ( x ) a a x a x n 0 1 2
an xn
1.4
数值分析课程第五版课后习题答案(李庆扬等)

数值分析课程第五版课后习题答案(李庆扬等)数值分析课程第五版课后习题答案(李庆扬等)第一章:数值分析导论1. 解答:数值分析是一门研究如何使用计算机来解决数学问题的学科。
它包括了从数学理论到计算实现的一系列技术。
数值分析的目标是通过近似的方式求解数学问题,其结果可能不是完全精确的,但是能够满足工程或科学应用的要求。
2. 解答:数值分析在实际应用中起着重要的作用。
它可以用于求解复杂的数学方程、计算机模拟及建模、数据的统计分析等等。
数值分析是科学计算和工程计算的基础,对许多领域都有着广泛的应用,如物理学、经济学、生物学等。
3. 解答:数值方法指的是使用数值计算的方式来求解数学问题。
与解析方法相比,数值方法一般更加灵活和高效,可以处理一些复杂的数学问题。
数值方法主要包括了数值逼近、插值、数值积分、数值微分、线性方程组的求解、非线性方程的求根等。
4. 解答:计算误差是指数值计算结果与精确解之间的差异。
在数值计算中,由于计算机的有限精度以及数值计算方法本身的近似性等因素,都会导致计算误差的产生。
计算误差可以分为截断误差和舍入误差两种。
第二章:数值误差分析1. 解答:绝对误差是指实际值与精确值之间的差异。
例如,对于一个计算出的数值近似解x和精确解x_0,其绝对误差为| x - x_0 |。
绝对误差可以衡量数值近似解的精确程度,通常被用作评估数值计算方法的好坏。
2. 解答:相对误差是指绝对误差与精确解之间的比值。
对于一个计算出的数值近似解x和精确解x_0,其相对误差为| (x - x_0) / x_0 |。
相对误差可以衡量数值近似解相对于精确解的精确度,常用于评估数值计算方法的收敛速度。
3. 解答:舍入误差是由于计算机的有限精度而引起的误差。
计算机中使用的浮点数系统只能表示有限的小数位数,因此在进行数值计算过程中,舍入误差不可避免地会产生。
舍入误差会导致计算结果与精确结果之间存在差异。
4. 解答:误差限度是指对于给定的数值计算问题,所能容忍的误差范围。
数值分析程序 李庆扬版

数值分析部分程序Lagrangefunction y=Lagrange(x,y,xstar) n=length(y);s=0;for i=1:nz1=1;z2=1;for j=1:nif j~=iz1=z1.*(xstar-x(j)); z2=z2.*(x(i)-x(j));endendz=z1./z2;s=s+z*y(i);endy=s;x=[0.5610,0.5628,0.56401,0.56521];y=[0.83741,0.82659,0.82577,0.82495]; xstar=[0.5625,0.5635,0.5645];y=Lagrange(x,y,xstar)Runget=-5:0.1:5;ft=1./(1+t.*t);t1=-5:1:5;ft1=1./(1+t1.*t1);y1=Lagrange(t1,ft1,t);plot(t,ft,t,y1);QJSfunction b=QJS(A,x)n=length(A);b(1)=A(1);c(1)=A(1);for k=1:1:n-1;b(1+k)=b(k)*x+A(k+1);ends=b(n)for i=1:1:n-2;c(i+1)=c(i)*x+b(i+1);ends1=c(n-1)A=[2,0,-3,3,-4];b=QJS(A,-2);s =10s1 =65面例5切比雪夫多项式的零点做插值,没有龙格现象的产生t=-5:0.1:5;ft=1./(1+t.*t);k=0:1:10;t1=5*cos((21-2.*k).*pi./22);ft1=1./(1+t1.*t1);y1=Lagrange(t1,ft1,t);t2=-5:1:5;ft2=1./(1+t2.*t2);y2=Lagrange(t2,ft2,t);plot(t,ft,t,y1,t,y2);希尔伯特矩阵的病态性function H=Hilbert(n)H=zeros(n,n);for i=1:1:n;for j=1:1:n;H(i,j)=1/(i+j-1);endenda=rand(n,1);d=H*a;d(1)=d(1)+0.0001;b=inv(H)*d;e=a-b;norm(e,2)norm(e,inf)function [c,s]=leastp(x)p1=1;p2=inline('x','x');p3=inline('(3*x.^2-1)/2','x');p4=inline('(5*x.^3-3.*x)/2','x');fp1=inline('exp(x)','x');fp2=inline('x.*exp(x)','x');fp3=inline('exp(x).*(3*x.^2-1)/2','x');fp4=inline('exp(x).*(5*x.^3-3.*x)/2','x');c(1)=quad(fp1,-1,1)/2c(2)=3.*quad(fp2,-1,1)/2c(3)=5.*quad(fp3,-1,1)/2c(4)=7.*quad(fp4,-1,1)/2s=c*例7M文件function p=leg(x,n)p(1:n+1)=x;p(1)=1;p(2)=x;for i=3:n+1p(i)=(2*i-3)*p(i-1)*x/(i-1)-(i-2)*p(i-2)/(i-1);p=simplify(p);M文件syms xT1=leg_p(x,3);for i=1:4disp(T1(i))end主程序>> clear>> d_main也可主程序(不要第二个M文件时) syms xT1=leg_p(x,3);for i=1:4disp(T1(i))endM文件function pf=pingfang(f,n) syms xf1=leg_p(x,n);for i=1:n+1g(i)=f1(i)*f;temp(i)=int(g(i),-1,1); a(i)=(2*i-1)*temp(i)/2; endpf=0;for i=1:n+1pf=pf+a(i)*f1(i);Endpf=simplify(pf);pf=vpa(pf,8);M文件f=inline('exp(x)','x');px=pingfang(f,3)主程序>>clear>> pingfang_main也可M文件function pf=pingfang(f,n) syms xf1=leg_p(x,n);for i=1:n+1g(i)=f1(i)*f;temp(i)=int(g(i),-1,1); a(i)=(2*i-1)*temp(i)/2; endpf=0;for i=1:n+1pf=pf+a(i)*f1(i);pf=vpa(pf,8);主程序px=pingfang('exp(x)',3)结果px =.99629405+.99795482*x+.53672140*x^2+.17614068*x^3勒让德逼近:function [P,E]=LRD(a,b)A0=[0,0,0,1];A1=[0,0,1,0];A2=[0,3/2,0,-1/2];A3=[5/2,0,-3/2,0];f=inline('exp(x)','x');fp1=inline('exp(x)','x');fp2=inline('x.*exp(x)','x');fp3=inline('exp(x).*(3*x.^2-1)/2','x');fp4=inline('exp(x).*(5*x.^3-3.*x)/2','x');q1=quad(fp1,a,b);q2=quad(fp2,a,b);q3=quad(fp3,a,b);q4=quad(fp4,a,b);a1=q1*1/2;a2=q2*3/2;a3=q3*5/2;a4=q4*7/2;A=a1.*A0+a2.*A1+a3.*A2+a4.*A3;P=poly2str(A,'x');E1=quad(@(x)exp(2*x),a,b);E2=2.*a1.^2+2/3.*a2.^2+2/5.*a3.^2+2/7.*a4.^2;E=sqrt(E1-E2);end切比雪夫逼近:function P=QBXF1(a,b)A0=[0,0,0,1];A1=[0,0,1,0];A2=[0,2,0,-1];A3=[4,0,-3,0];f1=inline('2.*exp(cos(x))./pi','x');f2=inline('2.*exp(cos(x)).*cos(x)./pi','x');f3=inline('2.*exp(cos(x)).*cos(2.*x)./pi','x'); f4=inline('2.*exp(cos(x)).*cos(3.*x)./pi','x'); q1=quad(f1,a,b);q2=quad(f2,a,b);q3=quad(f3,a,b);q4=quad(f4,a,b);A=q1/2.*A0+q2.*A1+q3.*A2+q4.*A3;P=poly2str(A,'x');end运行程序P=QBXF1(0,pi)。
数值分析第五版_李庆扬

数值分析第五版_李庆扬一、课程基本信息课程中文名称:数值分析课程英文名称:Numerical Analysis课程类别:专业基础课开课学期:秋适用专业:信息与计算科学;应用数学总学时:86学时(其中理论课56学时,上机实习30学时)总学分:5(理论课3学分;上机实习2学分)预修课程(编号):数学分析,高等代数,常微分方程课程简介:本课程是大学本科信息与计算科学和应用数学专业的一门基础课,也是工科研究生的必修课。
本课程的主要内容是研究各种数学问题的数值计算方法的设计、计算误差分析以及有关理论和具体实现的一门数学课程。
是应用数学的重要分支之一。
建议教材:《计算方法》(二版)(邓建中、刘之行),西安,西安交通大学出版社,2001 参考书:[1]数值分析学习指导,关治编,出版社:清华大学出版社,出版时间:2008年;[2]数值分析,何汉林,梅家斌,科学出版社,2007年;[3]《数值计算引论》白峰杉高等教育出版社 2005年[4]《数值分析》(第五版)李庆扬易大义等清华大学出版社 2008年[5]Numerical Analysis,R.Kress,世界图书出版公司20036、数值分析学习辅导习题解析,李宏、徐长发编,华中科技大学出版社,2001年。
二、理论课程教育目标通过本课程的教学使学生能了解现代科学计算中常用的数值计算方法及其基本理论,系统掌握数值分析的基本概念和分析问题、解决问题的基本方法,为运用数值分析的理论知识并为掌握更复杂的现代计算方法打好。
三、理论教学内容与要求(含学时)第一章:计算方法的一般概念(4学时)本章教学内容:理解计算方法的意义、研究内容与方法,理解并掌握误差的概念(包括误差的来源、绝对误差、相对误差),掌握有效数字及舍入误差对计算的影响。
第二章:解线性方程组的直接法(8学时)本章教学内容:1、高斯消去法;选主元的高斯消去法;2、矩阵的LR分解;解三对角方程组的追赶法;解方程组的平方根法;矩阵的求逆;3、方程组的数;病态方程组的判断。
李庆扬数值分析第五版习题答案解析清华大学出版社

又
即计算值比准确值大。
故 在 内至少有三个互异零点,
依此类推, 在 内至少有一个零点。
记为 使
又
其中 依赖于
分段三次埃尔米特插值时,若节点为 ,设步长为 ,即
在小区间 上
16.求一个次数不高于4次的多项式P(x),使它满足
解:利用埃米尔特插值可得到次数不高于4的多项式
设
其中,A为待定常数
从而
17.设 ,在 上取 ,按等距节点求分段线性插值函数 ,计算各节点间中点处的 与 值,并估计误差。
19。观测物体的直线运动,得出以下数据:
时间t(s)
0
0.9
1.9
3.0
3.9
5.0
距离s(m)
0
10
30
50
80
110
求运动方程。
解:
被观测物体的运动距离与运动时间大体为线性函数关系,从而选择线性方程
令
则
则法方程组为
从而解得
故物体运动方程为
20。已知实验数据如下:
19
25
31
38
44
19.0
32.3
将 代入得
由此得矩阵开工的方程组为
求解此方程组,得
又 三次样条表达式为
将 代入得
21.若 是三次样条函数,证明:
若 ,式中 为插值节点,且 ,则
证明:
从而有
第三章 函数逼近与曲线拟合
1. ,给出 上的伯恩斯坦多项式 及 。
解:
伯恩斯坦多项式为
其中
当 时,
当 时,
2.当 时,求证
证明:
若 ,则
3.证明函数 线性无关
解:
采用复化梯形公式时,余项为
数值分析第五版_李庆扬

数值分析第五版_李庆扬数值分析第五版_李庆扬一、课程基本信息课程中文名称:数值分析课程英文名称:Numerical Analysis课程类别:专业基础课开课学期:秋适用专业:信息与计算科学;应用数学总学时:86学时(其中理论课56学时,上机实习30学时)总学分:5(理论课3学分;上机实习2学分)预修课程(编号):数学分析,高等代数,常微分方程课程简介:本课程是大学本科信息与计算科学和应用数学专业的一门基础课,也是工科研究生的必修课。
本课程的主要内容是研究各种数学问题的数值计算方法的设计、计算误差分析以及有关理论和具体实现的一门数学课程。
是应用数学的重要分支之一。
建议教材:《计算方法》(二版)(邓建中、刘之行),西安,西安交通大学出版社,2001 参考书:[1]数值分析学习指导,关治编,出版社:清华大学出版社,出版时间:2008年;[2]数值分析,何汉林,梅家斌,科学出版社,2007年;[3]《数值计算引论》白峰杉高等教育出版社 2005年[4]《数值分析》(第五版)李庆扬易大义等清华大学出版社2008年[5]Numerical Analysis,R.Kress,世界图书出版公司20036、数值分析学习辅导习题解析,李宏、徐长发编,华中科技大学出版社,2001年。
二、理论课程教育目标通过本课程的教学使学生能了解现代科学计算中常用的数值计算方法及其基本理论,系统掌握数值分析的基本概念和分析问题、解决问题的基本方法,为运用数值分析的理论知识并为掌握更复杂的现代计算方法打好。
三、理论教学内容与要求(含学时)第一章:计算方法的一般概念(4学时)本章教学内容:理解计算方法的意义、研究内容与方法,理解并掌握误差的概念(包括误差的来源、绝对误差、相对误差),掌握有效数字及舍入误差对计算的影响。
第二章:解线性方程组的直接法(8学时)本章教学内容:1、高斯消去法;选主元的高斯消去法;2、矩阵的LR分解;解三对角方程组的追赶法;解方程组的平方根法;矩阵的求逆;3、方程组的数;病态方程组的判断。
数值分析第五版_李庆扬__课后习题答案

第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x x e x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。
解:设()n f x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=Q , 1||n p x nx C n n-⋅∴==又((*))(*)r p r x n C x εε≈⋅Q 且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*11.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x x εεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少?解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===g g(*)(*)3(*)r p r r V C R R εεε∴≈=g又(*)1r V ε=Q故度量半径R 时允许的相对误差限为1(*)10.333r R ε=⨯≈6.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=Q10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。
数值分析第五版_李庆扬

数值分析第五版_李庆扬一、课程基本信息课程中文名称:数值分析课程英文名称:Numerical Analysis课程类别:专业基础课开课学期:秋适用专业:信息与计算科学;应用数学总学时:86学时(其中理论课56学时,上机实习30学时)总学分:5(理论课3学分;上机实习2学分)预修课程(编号):数学分析,高等代数,常微分方程课程简介:本课程是大学本科信息与计算科学和应用数学专业的一门基础课,也是工科研究生的必修课。
本课程的主要内容是研究各种数学问题的数值计算方法的设计、计算误差分析以及有关理论和具体实现的一门数学课程。
是应用数学的重要分支之一。
建议教材:《计算方法》(二版)(邓建中、刘之行),西安,西安交通大学出版社,2001 参考书:[1]数值分析学习指导,关治编,出版社:清华大学出版社,出版时间:2008年;[2]数值分析,何汉林,梅家斌,科学出版社,2007年;[3]《数值计算引论》白峰杉高等教育出版社 2005年[4]《数值分析》(第五版)李庆扬易大义等清华大学出版社 2008年[5]Numerical Analysis,R.Kress,世界图书出版公司20036、数值分析学习辅导习题解析,李宏、徐长发编,华中科技大学出版社,2001年。
二、理论课程教育目标通过本课程的教学使学生能了解现代科学计算中常用的数值计算方法及其基本理论,系统掌握数值分析的基本概念和分析问题、解决问题的基本方法,为运用数值分析的理论知识并为掌握更复杂的现代计算方法打好。
三、理论教学内容与要求(含学时)第一章:计算方法的一般概念(4学时)本章教学内容:理解计算方法的意义、研究内容与方法,理解并掌握误差的概念(包括误差的来源、绝对误差、相对误差),掌握有效数字及舍入误差对计算的影响。
第二章:解线性方程组的直接法(8学时)本章教学内容:1、高斯消去法;选主元的高斯消去法;2、矩阵的LR分解;解三对角方程组的追赶法;解方程组的平方根法;矩阵的求逆;3、方程组的数;病态方程组的判断。
数值分析课程第五版课后习题答案(李庆扬等)(OCR)

根是x,,2…,x-,且V。x,x…·,x)=V,Cx6,x…·)(x-x)…(x-x)。
V,(xo,x,…x-x)=11】 -x,)用a-x,)
[证明]由
可得求证。
=V,(Cx8,x,…,xX))11(x-x)
2、当x=1-1,2时,f(x)=0,-3.4,求f(x)的二次插值多项式。
L,(x)=y%((xx6--xx,)((xx-2x-x22))
y=f(x)=f0.5)=-0.693147,y2=f(x)=f(0.6)=-0.510826,则
L2(x)=y。 (x-x)(x-x2)
(x6-x)x-x)
(x-x)(x-x)
(x-x)(x-x2)
(x-xo)(x-x) (x2-xo)(x2-x)
=-0.916291×.(0(.x4-0-.05.)5()x(-00..64)-0.6-.
30—+2—9.x9583x31 ̄02'=0.8336×104
14、试用消元法解方程x组1+10"x=100
x+x2=2
,假定只有三位数计算,问结果是否
可靠?
[解]精确解为x1=0100-*1 10"-2 ,当使用三位数运 算时,得到
x =1,x2=1,结果可靠。
15、已知三角形面积s=s去= absinc,其中c为弧度,0<c< 且测量a,b,c
位有效数字;x=56.430有5位有效数字;x=7×10有2位有效数字。 4、利用公式(3.3)求下列各近似值的误差限,其中x,x;,x,x;均为第3题所给
的数。
(1)x+x2+x:
e(x+x写+x)=>
[解]
E(x)=E(x)+E(x)+E(x;)
3+tx10=1.05×103
(2)xxx;
李庆扬 数值分析第五版 习题答案

第2章 复习与思考题01ii i ii kx x x x 的基函数称为主要性质有 0,()1,k i kx i k()1n l x、什么是牛顿基函数?它与单项式基答:牛顿差值基函数为00101x ),(x x )(x x ),...,(x x )(x x )...(x x )}n 牛顿差值基函数中带有常数项01,,...n x x x ,这有单项式基不同。
阶均差?它有何重要性质 01n 2n 01n 2n -11[,,...,,][,,...,,]n n f x x x f x x x x x xk j 0j 0j-1j j+1j -k x x x x x x x ()...()()...()和k 阶均差的性质0101k-10[,,...,][,,...,]k kf x x x f x x x x x (分子前项多xk )[a,b]上存在阶导数,且节点2n ,[a,b]x ,则1()!f n0()nn n ik k kk k i i ki kx x y l x y x x ,(j 1,2,....,n)个点的牛顿插值多项式01[,,...,]k f x x x ,(k 1,2,....,n)两者的主要差异是未知数不一致。
拉格朗日插值多项式是系数知道,但基函数不知道。
牛顿插值多项式是函数知道,但系数不知道。
与一般多项式基本相同。
y ,其中系数矩阵用下列基底作多项式插值时,120001211112222121...1...1 (1)...n n n n n nnx x x x x x x x x x x x ,无非零元素。
)拉格朗日基底为01{(),(),...,()}n l x l x l x ,已知数为未知数为01{(),(),...,()}n l x l x l x ,则系数矩阵为00101x ),(x x )(x x ),...,(x x )(x x )...(x x )}n ,已,未知数为012{,,,...,}n a a a a ,则系数矩阵为102020211010100...010...01()()...0...............1()()...()n nnnnj j x x x x x x x x x x x x x x x x ,为下三角矩阵,矩阵的上三角元0。
第五版李庆扬数值分析习题答案

1 2 gt ,假定 g 是准确的,而对 t 的测量有 0.1 秒的误差,证明当 t 增 2 加时 S 的绝对误差增加,而相对误差却减少。 1 解: S gt 2 , t 0 2
10.设 S
( S * ) g 2t (t * )
当 t * 增加时, S * 的绝对误差增加
解:
1
1 2 1 * ( x2 ) 10 3 2 1 * ( x3 ) 10 1 2 1 * ( x4 ) 10 3 2 1 * ( x5 ) 10 1 2
( x1* ) 104
* * * (1) ( x1 x2 x4 ) * * * ( x1 ) ( x2 ) ( x4 )
1 1 1 104 103 103 2 2 2 3 1.05 10
* * * (2) ( x1 x2 x3 ) * * * * * * * * * x1 x2 ( x3 ) x2 x3 ( x1 ) x1 x3 ( x2 )
1 1 1 1.1021 0.031 101 0.031 385.6 104 1.1021 385.6 103 2 2 2 0.215
第一章 绪论 1.设 x 0 , x 的相对误差为 ,求 ln x 的误差。 e * x * x * 解:近似值 x* 的相对误差为 = er x* x* 1 而 ln x 的误差为 e ln x * ln x * ln x e* x* 进而有 (ln x*) 2.设 x 的相对误差为 2%,求 x n 的相对误差。 解:设 f ( x) x n ,则函数的条件数为 C p |
则 f (30) ln(30 899) 此时,
李庆扬数值分析第五版习题答案清华大学出版社

李庆扬数值分析第五版习题答案清华大学出版社Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。
解:设()n f x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅ 且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字;*20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x x εεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少解:球体体积为343V R π= 则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=故度量半径R 时允许的相对误差限为1(*)10.333r R ε=⨯≈ 6.设028Y =,按递推公式1n n Y Y -=(n=1,2,…) 计算到100Y27.982≈(5位有效数字),试问计算100Y 将有多大误差解:1n n Y Y -=10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =,若取27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。
数值分析课程第五版课后习题答案(李庆扬等)

第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。
[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%,求n x 的相对误差。
[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而n x 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字:1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字;0031.0*2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56*4=x 有5位有效数字;0.17*5⨯=x 有2位有效数字。
4、利用公式(3.3)求下列各近似值的误差限,其中*4*3*2*1,,,x x x x 均为第3题所给的数。
(1)*4*2*1x x x ++; [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε;(2)*3*2*1x x x ;[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε;(3)*4*2/x x 。
数值分析(李庆扬)第七章)

x1 x* x2 x*
x0 x* x1 x*
x*
x2
( x2 x1 )2 x0 2x1 x2
此种加速需用两次迭代值进行加工。 (P273)
如果将一次改进值作为一步,
则计算公式如下:
校正 再校正
改进
x% k 1 ( xk )
xk 1 ( x% k 1)
xk 1
xk 1
( xk 1 x% k 1 )2 xk 1 2 x% k 1 xk
证:因 ' ( x)在O(x*, * )内连续,且 ' ( x) 1, 故存 在正数L 1, *, 使得对x [x* , x* ],有
' (x) L 1 另一方面,由 (x* ) x*, 又有
(x) x* (x) (x*) L x x* 即 (x) [x* , x* ]。由上面定理知,迭代序列 xn1 ( xn )收敛于x*。
1), 若
f
a
b 2
0
输出根
x
ab 2
, 否则:若
f
a
2
b
0,
令
a1
a
2
b
,b1
b
反之
ab b1 2 , a1 a.
2 ),对 [ a1 ,b1 ] 区间重复1)的计算,并产生 [ a2 ,b2 ],
3),
若
f
ai
2
bi
0,则得到根
x ai bi . 2
二分法的收敛性
二分法产生一个有根区间:
f (2) 1 0
[1.5, 2] 为有根区间。
(1)x x 1 1(x) 因 1.5 1.5 1 1(x) 2 1 2
且
1' (x) 2
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
三、避免误差危害的若干原则
除了分清问题是否病态和算法是否数值稳定外,还要 考虑避免误差危害和防止有效数字损失的如下原则. 1.避免‘大数’除以‘小数’ 例6 仿计算机,采用3位十进制,用消元法求解方程组
1.0010 x 1.00 y 1.00 解: x得, 消 (1.00 1.00105 ) y (2.00 1.00105 )
100 100 项 项 123 0.01 0.01 0.01 0.01 0.01 0.01 123
4.减少运算次数 减少运算次数可以不但节省时间,而且减少舍入误差. 例10 计算多项式的值
Pn ( x) an x n an1x n1 a1x a0 .
* *
§3 误差定性分析、避免误差危害
误差分析简介 概率分析法 向后误差分析法
x g (a1,, an ), x fl g (a1 1,, an n ).
区间分析法
x [ , ], y [ , ], xy
* r 0.005/9.80 0.000005/ 0.00980.
定理1设近似数x * 表示为
x* 10m (a1 a2 101 al 10(l 1) )
* r
(2.1)
其中a1 0 . 若x * 具有n位有效数字,则其相对 误差限为 1 10( n1); 2a1
数值分析
第1章
一、什么是数值分析绪论 Nhomakorabea§1 数值分析的研究对象与特点
数值分析是计算数学的一个主要部分,计算数学是数 学科学的一个分支,它研究用计算机求解各种数学问题 的数值计算方法及其理论与软件实现. 实际问题→数学模型→数值计算方法 →程序设计→上机计算求出结果
二、数值分析的基本内容 1、数值逼近 插值法 函数逼近与曲线拟和 数值积分与数值微分 2、数值代数 线性代数问题(方程组和特征值) 非线性方程(组)数值解法
2.避免两个相近数相减
例7 求解x 2 16 x 1 0.
7 例8 计算A 10( cos 2。. 1 )
又如:当x, y接近时, x lg y ? lg x 1 x ? 当x, x * 接近时, ( x) f ( x*) ? f
3.防止‘大数’吃‘小数’ 例9 仿计算机在3位十进制下,
秦九韶算法:
S n an , S k xS k 1 ak , (k n-1,,0) P ( x) S . n 0
(3.4)
作业
P19, 5, 7, 9.
并估计误差.
I 0 0.6321, ( A) I n 1 nI n1, n 1,2,. * 1 1 e1 I 9 0.0684, ( I 9 ( ) 0.0684) ( B) * 2 10 10 1 (1 I * ), n 9,8,,1. I n1 n n 定义3 一个算法若输入数据有 误差, 而在计算过程中舍入
(2) (1) 10
5
5
1.00105 x 1.00y 1.00 1.00x 1.00y 2.00
5
105 x 1105
1.00001 0.9999899
y
2 105 1105
1.00 10 x 1.00 y 1.00 x* 0.00, y* 1.00 y 1.00 错.为什么,怎么办?
* r
反之, x *的相对误差限为 若 至少具有n位有效数字.
1 10( n1),则x * 2(a1 1)
例3要使 20的相对误差限小于0.1%,要取几位有效数字?
三、数值运算的误差估计
* * 四则运算,设x1, x2为准确值, x1 , x2为近似值,则误差限:
* * * * ( x1 x2 ) ( x1 ) ( x2 ), * * * * * * ( x1 x2 ) | x1 | ( x2 ) | x2 | ( x1 ), * * * * | x1 | ( x2 ) | x2 | ( x1 ) * * ( x1 / x2 ) . * 2 | x2 |
Cramer法则 vs Gauss消去法. 1、面向计算机 3、良好的计算复杂性 4、数值实验
1 x3 ( x ), 2 3
2、可靠的理论分析,保证收敛性、稳定性
四、如何学好数值分析 1、注意掌握基本原理、处理技巧,误差分析
2、注重实际问题,练习、作业 3、积极动手上机实践
§2 数值计算的误差
一、误差来源、分类
模型误差
观测误差
截断误差或方法误差
f (0) f (0) 2 f ( n) (0) n f ( x) Pn ( x) f (0) x x x 1! 2! n! f ( n1) ( ) n1 x 截断误差: Rn ( x) (n 1)! 舍入误差 R 3.14159 0.0000026. 数制转换、机器数.
一元函数f ( x),x为准确值, x * 为近似值,由Taylor公式 f ( x) f ( x*) f ( x*)( x x*)
f ( ) ( x x*) 2 , 2
在x, x * 之间,
得f ( x*)的误差限 ( f ( x*)) | f ( x*) | ( x*).
误差为2%,函数值相对误差为24%.
一般C p 10认为是病态.
其他计算问题也要考虑 条件数,考虑是否病态 .
二、算法的数值稳定性
考虑初始数据误差在计算中的传播问题.
例5
1 1 n x 计算I n e x e dx, n 0,1,, 0
I n 1 nI n1 , n 1,2,, I 0 1 e1.
一、病态问题与条件数
考虑计算函数值问题 ,
f ( x*) f ( x ) f ( x) x x
xf ( x ) f ( x)
Cp,
C p 称为计算函数值问题的 条件数.
例如f ( x) x , C p 10, f (1) 1, f (1.02) 1.24,自变量相对
10
其中a1 0 . 并且 1 x x * 10mn1 (2.2) 2 例1 42.195, 0.0375551, 8.00033, 2.71828,按四舍五 入写出上述各数具有四位有效数字的近似数.
例2 考察三位有效数字重力加速度g, 若以m/s2为单位, g≈9.80m/s2, 1 1 * 10mn1 g 9.80 102 , 2 2 * 1 按(2.1), 0, n 3. 绝对误差限 1 102. m 2 若以km/s2为单位, g≈0.00980m/s2, 1 g 0.00980 105 , 2 * 1 按(2.1), 3, n 3. 绝对误差限 2 105. m 2 而相对误差限相同:
3、常微方程数值解法和偏微方程数值解法
三、数值分析的特点
1 1 n 1 1 ln 2 1 (1) , 2 3 n
ln(1 x) x 2 x3 x , 2 3
ln 1 x 1 x
vs
1 1 1 1 1 1 1 ln 2 2( 3 5 7 ). 3 3 3 5 3 7 3
例如,x 10 1, y 1000 5.
| x|
10%,
* y
| y|
0.5%.
x 3.1415926,
* * 取三位 x3 3.14, 3 0.002, * * 取五位 x5 3.1416, 5 0.000008.
定义2 若近似值x *的误差限是某一位数字 的半个单位, 该位 到x *的第一位非零数字共有 n位,就说x * 有n位有效数字. 即 x* 10m (a1 a2 101 an 10( n1) ) (2.1)
* * 多元函数f ( x1 ,, xn ),x1 ,, xn为准确值x1 ,, xn的近似值, * * 同理得f ( x1 ,, xn )的误差限
f * ( xk ). ( f *) x k 1 k
n
*
例4 场地面积:s ld s (l*) s (d *). ( s*) l d
二、误差、有效数字 定义1 绝对误差,简称误差:
e* x * x, 其中x *为准确值x的近似值 .
. 误差限: * | e* | 的一个上界
例如,毫米尺 765 x 0.5
相对误差:
* er
e * 或 e* e * . , r x* x
* * r | er | 的一个上界. 相对误差限: * x