数值分析作业答案(第5章)part2
数值分析第五章答案
数值分析第五章答案【篇一：数值分析第五版计算实习题】第二章2-1程序：clear;clc;x1=[0.2 0.4 0.6 0.8 1.0];y1=[0.98 0.92 0.81 0.64 0.38];n=length(y1);c=y1(:);or j=2:n %求差商for i=n:-1:jc(i)=(c(i)-c(i-1))/(x1(i)-x1(i-j+1));endendsyms x df d;df(1)=1;d(1)=y1(1);for i=2:n %求牛顿差值多项式df(i)=df(i-1)*(x-x1(i-1));d(i)=c(i)*df(i);enddisp(4次牛顿插值多项式);p4=vpa(collect((sum(d))),5) %p4即为4次牛顿插值多项式,并保留小数点后5位数 pp=csape(x1,y1, variational);%调用三次样条函数 q=pp.coefs;disp(三次样条函数);for i=1:4s=q(i,:)*[(x-x1(i))^3;(x-x1(i))^2;(x-x1(i));1];s=vpa(collect(s),5)endx2=0.2:0.08:1.08;dot=[1 2 11 12];figureezplot(p4,[0.2,1.08]);hold ony2=fnval(pp,x2);x=x2(dot);y3=eval(p4);y4=fnval(pp,x2(dot));plot(x2,y2,r,x2(dot),y3,b*,x2(dot),y4,co);title(4次牛顿插值及三次样条);结果如下：4次牛顿插值多项式p4 = - 0.52083*x^4 + 0.83333*x^3 - 1.1042*x^2 + 0.19167*x + 0.98 三次样条函数x∈[0.2,0.4]时， s = - 1.3393*x^3 + 0.80357*x^2 - 0.40714*x + 1.04 x∈[0.4,0.6]时，s = 0.44643*x^3 - 1.3393*x^2 + 0.45*x +0.92571 x∈[0.6,0.8]时，s = - 1.6964*x^3 + 2.5179*x^2 - 1.8643*x + 1.3886 x∈[0.8,1.0]时，s =2.5893*x^3 - 7.7679*x^2 + 6.3643*x - 0.80571 输出图如下2-3（1）程序：clear;clc;x1=[0 1 4 9 16 25 36 49 64];y1=[0 1 2 3 4 5 6 7 8];%插值点n=length(y1);a=ones(n,2);a(:,2)=-x1;c=1;for i=1:nc=conv(c,a(i,:));endq=zeros(n,n);r=zeros(n,n+1);for i=1:n[q(i,:),r(i,:)]=deconv(c,a(i,:));%wn+1/(x-xk)enddw=zeros(1,n);for i=1:ndw(i)=y1(i)/polyval(q(i,:),x1(i));%系数endp=dw*q;syms x l8;for i=1:nl8(i)=p(n-i+1)*x^(i-1);enddisp(8次拉格朗日插值);l8=vpa(collect((sum(l8))),5)xi=0:64;yi=polyval(p,xi);figureplot(xi,yi,x1,y1,r*);hold ontitle(8次拉格朗日插值);结果如下：8次拉格朗日插值l8 =- 3.2806e-10*x^8 + 6.7127e-8*x^7 - 5.4292e-6*x^6 +0.00022297*x^5 - 0.0049807*x^4 + 0.060429*x^3 - 0.38141*x^2 +1.3257*x输出图如下：第五章4-1（3）程序：clc;clear;y= @(x) sqrt(x).*log(x);a=0;b=1;tol=1e-4;p=quad(y,a,b,tol);fprintf(采用自适应辛普森积分结果为： %d \n, p);结果如下：采用自适应辛普森积分结果为： -4.439756e-01第九章9-1(a)程序：clc;clear;a=1;b=2;%定义域h=0.05;%步长n=(b-a)/h;y0=1;%初值f= @(x,y) 1/x^2-y/x;%微分函数xn=linspace(a,b,n+1);%将定义域分为n等份 yn=zeros(1,n);%结果矩阵yn(1)=y0;%赋初值%以下根据改进欧拉公式求解for i=1:nxn=xn(i);xnn=xn(i+1);yn=yn(i);yp=yn+h*f(xn,yn);yc=yn+h*f(xnn,yp);yn=(yp+yc)/2;yn(i+1)=yn;endxn=yn;%以下根据经典四阶r-k法公式求解for i=1:nxn=xn(i);yn=yn(i);k1=f(xn,yn);k2=f(xn+h/2,yn+h/2*k1);k3=f(xn+h/2,yn+h/2*k2);k4=f(xn+h,yn+h*k3);yn=yn+h/6*(k1+2*k2+2*k3+k4);yn(i+1)=yn;enddisp(改进欧拉法四阶经典r-k法); disp([xn yn])结果如下：改进欧拉法四阶经典r-k法 110.998870.998850.99577 0.99780.991140.996940.985320.996340.978570.996030.971110.996060.963110.996450.95470.997230.945980.998410.9370510.92798 1.0020.91883 1.00440.90964 1.00730.90045 1.01060.89129 1.01430.88218 1.01840.87315 1.02290.86421 1.02780.85538 1.03310.84665 1.0388（b）程序：clc;clear;a=0;b=1;%定义域h=[0.1 0.025 0.01];%步长y0=1/3;%初值f= @(x,y) -50*y+50*x^2+2*x;%微分函数 xi=linspace(a,b,11);y=1/3*exp(-50*xi)+xi.^2;%准确解 ym=zeros(1,11);for j=1:3【篇二：数值分析(第五版)计算实习题第五章作业】题：lu分解法：建立m文件function h1=zhijielu(a,b)%h1各阶主子式的行列式值[n n]=size(a);ra=rank(a);if ra~=ndisp(请注意：因为a的n阶行列式h1等于零，所以a不能进行lu 分解。
数值分析课后习题及答案
第一章 绪论（12） 第二章 插值法（40-42）2、当2,1,1-=x 时，4,3,0)(-=x f ，求)(x f 的二次插值多项式。
[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。
3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。
X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ，6.01=x ，则693147.0)5.0()(00-===f x f y ，510826.0)6.0()(11-===f x f y ，则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ，从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。
若取4.00=x ，5.01=x ，6.02=x ，则916291.0)4.0()(00-===f x f y ，693147.0)5.0()(11-===f x f y ，510826.0)6.0()(22-===f x f y ，则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ，从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题：1、令00=x ，11=x ，写出x e x y -=)(的一次插值多项式)(1x L ，并估计插值余项。
数值分析课后习题答案
0 1
0 10 1 1 0 0 0 1
0 0 12 1 1 2 0 0 0
1 2
0 0 0 1 1 0
1 2
1 2
1 2
1
0 0 0 1 0
1 2
1 2
0
1 2
1 2
0
0
0
341 1 1
2-5.对矩阵A进行LDLT分解和GGT分解，并求解方程组
Ax=b，其中
16 4 8
1
A 4 5 4 , b 2
8 4 22
3
解
16 A 4
4 5
84
44 11
2-3(1).对矩阵A进行LU分解,并求解方程组Ax=b,其中
2 1 1 A1 3 2
4 ,b6
1 2 2
5
解
2 A 1
1 3
1 2
2 11
22
1
5 2
1
3 21来自，所以 A12
1
2 1 1
5 3
2-2(1).用列主元Gauss消元法解方程组
3 2 6x1 4 10 7 0x2 7 5 1 5x3 6
解
3 2 6 4 10 7 0 7 10 7 0 7
r1r2
消元
10 7 0 7 3 2 6 4 0 0.1 6 6.1
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623
数值分析第二版(丁丽娟)答案
1.32
1.68
2.08
2.52
3.00
解答下列问题 （1）试列出相应的差分表； （2）写出牛顿向前插值公式； （3）用二次牛顿前插公式计算 f（0.225）；
例3已知当 x=-1，0，2，3时，对应的函数值为
，
，
，
，
，求 的四次 Newton 插值多项式。
例4 设 对 n=1,2,3时
，证明：
例5 设 （1）
6 2730.5000 5051.0000 5051.5000
7 10922.5000 23483.0000 23483.5000
8 43690.5000 80827.0000 80827.5000
21.000000000000000 17.000000000000000 16.238095238095237 16.058823529411764 16.014662756598241 16.003663003663004 16.000915583226515
第一章答案
第二章答案
第三章答案
0 0.5 0.5 1 1 2.5000
5.0000 5.5000
第四章答案
2 10.5000 19.0000 19.5000
3 42.5000 91.0000 91.5000
4 170.5000 315.0000 315.5000
5 682.5000 1467.0000 1467.5000
3、 用规范化幂法求
按模最大的特征值和对应的特征向量，取初值
。当特征值有3位小数稳定时停止。
4、 用反幂法求矩阵
练习五
，迭代7次。
的最接近于6 的特征值和对应的特征向量，取初值
数值分析第五章实习题答案
数值分析第五章实习题答案数值分析第五章实习题答案数值分析是一门研究如何使用计算机来解决数学问题的学科。
在数值分析的学习过程中，实习题是非常重要的一部分，通过实习题的练习，可以帮助我们巩固所学的知识，并且提高我们的解题能力。
本文将为大家提供数值分析第五章实习题的答案，希望对大家的学习有所帮助。
第一题：求下列方程的一个正根，并用二分法和牛顿法分别计算根的近似值。
方程：x^3 - 3x + 1 = 0解答：首先，我们可以通过绘制函数图像来初步估计方程的根的范围。
根据图像，我们可以大致确定根在区间[0, 2]之间。
接下来，我们使用二分法来计算根的近似值。
根据二分法的原理，我们将区间[0, 2]等分为两部分，然后判断根在哪一部分。
不断重复这个过程，直到找到根的近似值。
具体计算过程如下：- 将区间[0, 2]等分为两部分，得到中点x = 1。
- 计算方程在x = 1处的函数值f(1) = -1。
- 根据函数值的正负性，我们可以确定根在区间[1, 2]之间。
- 将区间[1, 2]等分为两部分，得到中点x = 1.5。
- 计算方程在x = 1.5处的函数值f(1.5) = 1.375。
- 根据函数值的正负性，我们可以确定根在区间[1, 1.5]之间。
- 重复以上步骤，直到找到根的近似值。
最终得到根的近似值为x ≈ 1.365。
接下来，我们使用牛顿法来计算根的近似值。
牛顿法是一种迭代法，通过不断逼近根的位置来计算根的近似值。
具体计算过程如下：- 选择初始近似值x0 = 1。
- 计算方程在x = 1处的函数值f(1) = -1。
- 计算方程在x = 1处的导数值f'(1) = 4。
- 利用牛顿法的迭代公式x1 = x0 - f(x0)/f'(x0)，我们可以得到x1 ≈ 1.333。
- 重复以上步骤，直到找到根的近似值。
最终得到根的近似值为x ≈ 1.365。
通过二分法和牛顿法，我们分别得到了方程x^3 - 3x + 1 = 0的一个正根的近似值为x ≈ 1.365。
电子科技大学-数值分析答案-钟尔杰
| x n +1 − 7 |=
而xn具有n位有效数，故
所以
| x n +1 − 7 |≤
由此得xn+1的误差限
1 2 7
| x n − 7 |2 ≤
1 × × 10 2− 2 n 2 7 4
1
| x n +1 − 7 |≤
1 × 10 1− 2 n 2
故，xn+1是 7 的具有 2n位有效数字的近似值。 三、问题 1．假定 a0，b0是非负实数且a0≠b0，按如下递推公式
∑ [ai ∑ b j ]
i =1 j =1
n，仍为( n + 2 ) ( n – 1) / 2。 ，算法输出 11 试构造一个算法，对输入的数据 x0，x1，x2，……，xn，以及x（均为实数） 为 ( x –x0) ( x –x1) ( x –x2)……( x –xn) 的计算结果。 解 算法如下： 第一步：输入x；x0，x1，x2，……，xn，M Å (x – x0 )；k Å 0； 第二步：M Å M×(x – x0 )；k Å k+1； 第三步：判断，若 k ≤ n，则转第二步；否则输出 M，结束。 12 利用级数公式
4
π 1 dx = arctan 1 = 可以计算出无理数π 的值。将定积分表示为积分和 2 4 1+ x
R
H
∫
1
0
xn dx ( n = 1，2，…，20) 的递推 5+ x
关系，并研究递推算法的数值稳定性。 6．计算两个多项式Pn(x)和Qm(x)的乘积多项式Tn+m(x)的方法称为向量的卷积方法。设
第一章 习题解答与问题
一、习题解答 1 设 x>0，x 的相对误差限为 δ，求 ln x 的误差。 解：设 x的准确值为x*，则有 ( | x – x* | /|x*| ) ≤ δ 所以 e(ln x)=| ln x – ln x* | =| x – x* | ×| (ln x)’|x=ξ·≈ ( | x – x* | / | x*| ) ≤ δ 另解： e(ln x)=| ln x – ln x* | =| ln (x / x*) | = | ln (( x – x* + x*)/ x*) | = | ln (( x – x* )/ x* + 1) |≤( | x – x* | /|x*| ) ≤ δ 2 设 x = – 2.18 和 y = 2.1200 都是由准确值经四舍五入而得到的近似值。求绝对误差限 ε( x ) 和 ε( y ) 。 解：| e(x) | = |e(– 2.18)|≤ 0.005，| e(y) | = |e( 2.1200)|≤ 0.00005，所以 ε( x )=0.005， ε( y ) = 0.00005。 3 下近似值的绝对误差限都是 0.005，问各近似值有几位有效数字 x1=1.38，x2= –0.0312，x3= 0.00086 解：根据有效数字定义，绝对误差限不超过末位数半个单位。由题设知，x1，x2， x3有效 数末位数均为小数点后第二位。故x1具有三位有效数字，x2具有一位有效数字，x3具有零位 有效数字。 4 已知近似数 x 有两位有效数字，试求其相对误差限。 解：| er(x) | ≤ 5 × 10– 2 。 5 设 y0 = 28，按递推公式 yn = yn-1 –
数值分析作业答案
第2章 插值法1、当x=1，-1，2时，f(x)=0，-3，4，求f(x)的二次插值多项式。
（1）用单项式基底。
（2）用Lagrange 插值基底。
（3）用Newton 基底。
证明三种方法得到的多项式是相同的。
解：（1）用单项式基底设多项式为:2210)(x a x a a x P ++=，所以：6421111111111222211200-=-==x x x x x x A 37614421111111424113110111)()()(222211200222221112000-=-=---==x x x x x x x x x f x x x f x x x f a 2369421111111441131101111)(1)(1)(12222112002222112001=--=--==x x x x x x x x f x x f x x f a 6565421111111421311011111)(1)(1)(12222112002211002=--=---==x x x x x x x f x x f x x f x a 所以f(x)的二次插值多项式为：2652337)(x x x P ++-= （2）用Lagrange 插值基底)21)(11()2)(1())(())(()(2010210-+-+=----=x x x x x x x x x x x l)21)(11()2)(1())(())(()(2101201------=----=x x x x x x x x x x x l)12)(12()1)(1())(())(()(1202102+-+-=----=x x x x x x x x x x x lLagrange 插值多项式为：372365)1)(1(314)2)(1(61)3(0)()()()()()()(22211002-+=+-⨯+--⨯-+=++=x x x x x x x l x f x l x f x l x f x L所以f(x)的二次插值多项式为：22652337)(x x x L ++-= (3) 用Newton 基底: 均差表如下：Newton 372365)1)(1(65)1(230))(](,,[)](,[)()(21021001002-+=+-+-+=--+-+=x x x x x x x x x x x x f x x x x f x f x N所以f(x)的二次插值多项式为：22652337)(x x x N ++-= 由以上计算可知，三种方法得到的多项式是相同的。
数值分析第五版第5章习题答案
第5章
）矩阵行列式的值很小。
）矩阵的范数小。
）矩阵的范数大。
（7）奇异矩阵的范数一定是零。
答：错误，
∞
•可以不为0。
（8）如果矩阵对称，则|| A||1 = || A||∞。
答：根据范数的定义，正确。
（9）如果线性方程组是良态的，则高斯消去法可以不选主元。
答：错误，不选主元时，可能除数为0。
（10）在求解非奇异性线性方程组时，即使系数矩阵病态，用列主元消去法产生的误差也很小。
答：错误。
对于病态方程组，选主元对误差的降低没有影响。
（11）|| A ||1 = || A T||∞。
答：根据范数的定义，正确。
（12）若A是n n的非奇异矩阵，则
)
(
cond
)
(
cond1-
=A
A。
答：正确。
A是n n的非奇异矩阵，则A存在逆矩阵。
根据条件数的定义有：
1
111111 cond()
cond()()
A A A
A A A A A A A
-
------
=•
=•=•=•
习题
如有侵权请联系告知删除，感谢你们的配合！。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
解 中 ，故不能分解。但由于 ，所以若交换 的第1行与第3行，则可以分解且分解是唯一的。
在 中， ，故不能分解。但 可以分解为
，
其中 ， 为任意常数，且 奇异，故分解不唯一。
对于 ， ，故 可以分解且分解唯一。
。
5.13.求证：(1). ；(2). 。
证明(1).由定义知
故 。
(2).由范数定义，有
证明只需证明 满足向量范数的三个条件。
(1).因 正定对称，故当 ， ；而当 时，
。
(2).对任意 ，有
。
(3).因 正定，故有分解 ，因而
对任意 ，由 的三角不等式有
，故 是 上的向量范数。
又
所以 。
5.14.设 且非奇异，又 设为 上一向量范数，定义
试证明 是 上向量的一种范数。
证明只需证明 满足向量范数的三个条件。
(1).因 非奇异，故对任意 ，有 ，故 ，当且仅当 时，有 。
(2).对任意 ，有
。
(3).对任意 ，有
，故 是 上的向量范数。
5.15.设 为对称正定矩阵，定义 ，试证明 是 上向量的一种范数。
并求出系数矩阵 的行列式(即 )的值。
解
所以解为 ， ， ， 。
5.9.用追赶法解三对角方程组 ，其中
， 。
解设 有分解
，
由公式
其中 ， 分别是系数矩阵的主对角元素及其下边和上边的次对角线元素。
具体计算，可得
， ， ， ， ，
， ， ， 。由，得 ， ， Nhomakorabea ， ；再由
，
得 ， ， ， ， 。
5.11.下述矩阵能否分解为 (其中 为单位下三角矩阵， 为上三角矩阵)若能分解，那么分解是否唯一
数值分析作业答案(第5章) part2
5.6.证明：
(1).如果 是对称正定矩阵，则 也是对称正定矩阵
(2).如果 是对称正定矩阵，则 可以唯一地写成 ，其中 是具有正对角元的下三角矩阵。
证明：
(1).因 是对称正定矩阵，故其特征值 皆大于 ，因此 的特征值 也皆大于 。因此 也皆大于 ，故 是可逆的。又
则 也是对称正定矩阵。
(2).由 是对称正定，故它的所有顺序主子阵均不为零，从而有唯一的杜利特尔分解 。又
其中 为对角矩阵， 为上三角矩阵，于是
由 的对称性，得
由分解的唯一性得
从而
由 的对称正定性，如果设 表示 的各阶顺序主子式，则有
， ，
故
因此
，
其中 为对角元素为正的下三角矩阵。
5.7.用列主元消去法解线性方程组