i-l
77
严+i )胡_£a 局z —£知屮
?
冃 o
川+1)
丫伙+1) _ 丫⑹,-
X
i
一 Xj
a
i
⑹ 如果£⑷二兀⑷-兀*,其中X 是方程组的精确解,求证:
2
1
4 1 4= 2’ 1 1 1
—二兀-~x 2 +x 4 =-.
(a)证明
gE) —g ⑷)=—£
(d) 由此推出,如果A 是具有正对角元素的非奇异矩阵,且高斯一塞德尔方法对任意初始向 量X ⑼是收敛的,则A 是正定阵。
13. 设A 与B 为n 阶矩阵,A 为非奇异,考虑解方程组
其中右,Z? , £, d?丘R o
(a ;義岀下列迭代方法收敛的充要条件
Az{',,+1)=勺 _ BzJ , Azf+" =Z?2 - BZ ;"“(777 > 0);
(b)找岀下列迭代方法收敛的充要条件
Az ;m+1) =b x -Bz^n \Az^n+l) =b 2 -Bz{m+1,(rn>0);
比较两个方法的收敛速度。
14.证明矩阵
对于2
是正定的,而雅可比迭代只对2 2是收敛的。
_
5 1 2 3 _
0 2 0 4 3-12-1
15. 设 L 0 3 ° 7」,试说明人为可约矩阵。
16. 给定迭代过程,严"=CxW + g ,其中C w R nxn 伙=0,1,2,…),试证明:如果C 的 特征值儿(C )= °(’ = 1,2,…),则迭代过程最多迭代n 次收敛于方程组的解。
17. 画出SOR 迭代法的框图。
18. 设A 为不可约弱对角优势阵且0 <0<1,求证:解Ax^b 的SOR 方法收敛。 19. 设4兀=",其中A 为非奇异阵。 (a) 求证4丁4为对称正定阵; (b) 求证cond^A), = (cond(A)2)2 o
第九章矩阵的特征值与特征向量计算
1.用幕法计算下列矩阵的主特征值及对应的特征向量:
'7
3 ~2
'3
-4 3_ A =
3 4 -1 A
2 =
-4 6 3 (a) -2
-1
3
, (b)
3
3
1
当特征值有3位小数稳定时迭代终止。
其中
7=1
j=i
(C)设A 是对称的,二次型
2.方阵T 分块形式为
^11 A? * ?? T : T =
^22 * ?' T
2n
T nn_
其中人(,=1,2,…小)为方阵,T 称为块上三角阵,如果对角块的阶数至多不超过2,则 称T 为准三角形形式,用b (T )记矩阵T 的特征值集合,证明
a (r )= |Ja (7;7).
Z=1
3. 利用反幕法求矩阵
~6 2 f 2 3 1 1 1 1
的最接近于6的特征值及对应的特征向量。 4. 求矩阵
4 0 0 0 3 1 0 1 3
与特征值4对应的特征向量。 5. 用雅可比方法计算
10 1.0 0.5
A= 1.0 1.0 0.25 0.5 0.25 2.0
的全部特征值及特征向量,用此计算结果给出例3的关于p
的最优值。
6. (a )设A 是对称矩阵,X 和珂x h=l )是A 的一个特征值及相应的特征向量,又设P 为 一
个正交阵,使
Px = e
x,并计算 B = PAF^ o
7.利用初等反射阵将
1 3 4 A= 3 1
2 4 2 1
正交相似约化为对称三对角阵。
证明B = PAF^的第一行和第一列除了 X 外其余元素均为零。 (b )对于矩阵
入=9是其特征值,
试求一初等反射阵P ,使
&设A e R nxn,且Q"勺i不全为零,出为使a n = 0的平面旋转阵,试推导计算场山第d 行,第j行元素公式及人磅第i列,第j列元素的计算公式。
9.设人-1是由豪斯荷尔德方法得到的矩阵,又设y是人”-1的一个特征向量。
(a)证明矩阵A对应的特征向量是* = P\P2 -P n-iy .
(b)对于给出的y应如何计算x?
10.用带位移的QR方法计算
全部特征值。
11.试用初等反射阵A分解为QR,其中Q为正交阵,R为上三角阵,
1 1 1 _
2 -1 -1
2-4 5
数值分析习题答案
第一章绪论习题参考答案
£ (lnx) ~ \X I
兀:有5位有效数字,兀;有2位有效数字,兀;有4位有效数字,兀:有5位有效 数字,以有2位有效数字。
£(兀:+ 兀;+x*) Q £(兀:)+ £(兀;)+ £(兀;)nOBxlO"4 +O.5xlO~3 +O.5xlO -3
=1.05xl0-3 g (心;兀;)? 卜対 沢彳)+ X 匐 £(兀;)+ |xfx*| £(£) = 0.214790825
* 1 *
£(—^) — — &(兀4)~ 8.855668 xlO -6
兀|> 兀4
兀4
^00) = 100x^x1x10 3 =|xl0-
cc /ZTTT __ ccr 吃 =28-丁783 = ----------------- 】〕 ? —-— ? 0.01786 =28 + 7783 ?55.982 , 28 + ^783 55.982 f-^00 1 , A7 ----- wdx = --- arc tgN
)N 1 + x 2 2 6
£(x) = f(V?) ?-5_W) = 0.005
2
c
"S )?业-巴-耳 £(S )dg”|£(/) = O.lgM, 绝对误差增
加,相对误差减小。
畑“叽)冷"计算过程不稳定。
/ = (V2-1)6? 0.005051,如果令 72=1.4,则 /; =(A /2-1)6 = 0.004096 ; (V2+1)5
,歩=(3 —2血尸=0.008,儿(3 + 2血尸
,
去=99-70血=1,办的结果最好。
/(30) =-4.094622 ,开平方时用六位函数表计算所得的误差为
?(小_心)?"卜T" _竺g
n \n ?
0.02〃
3V
1 £(v\ I
= = 0.003333
,故t 增加时S 的
心)
/
分另U 代入等价公式 — -l),f 2(x) = -ln(x + Vx 2 +1)
计
算
可
得
C —— =(x+JF —1)£ = 60X 丄 xio" =3X 10~ x-4x^l 2
^ = 1x10^
2
中
£(_/;)= ln(l + —r =)
x_px -1
14.
€? 1 1
Q —— = — x — x 1()7 =8.33x10-7
x+4x^i 60 2 1000000000 ’ w 999999998 ’ nnnnnn x } = ----------- q 1.000000,
x 9 = ----------------- q 1.000000 方程组的真解为 999999999 2
999999999
£(Q=血(1 +
x + y/x~ +1
15.
第二章插值法习题参考答案
1.
n —1
v”(X )=口 &—兀)n (兀—“)
z=0
0< j
.
匕-(兀。,兀],%…_!)= n (羽-勺)
0< j
2.
(兀+ 1)(兀_2) |( 3) (兀_ 1)(兀_2) | 4 (兀―1)(兀+1) (1 +1)(1 —2)(一 厂(—1 — 1)(—1 — 2) * (2-1)(2 + !)
3 线性插值.取 = 0.6, y 0 = —0.693147 ,开=-0.510826
In 0.54 ? 厶(0.54) = y 0 + “ 一 儿? (0.54 -x 0) = -0.620219 尢I —勺
二次插值:取
x 0 = 0.4,兀]=0.5, x 2 = 0.6, y 0 = -0.916291, y x = -0.693147 , y 2 = -0.510826 则 In 0.54
(0.54)
(0.54 — Xj )(0.54 — x 2) (0.54 — x 0 )(0.54 — x 2) (0.54 — )(0.54 —州) =y 0 ------------------------------ + % ------------------------------ + y 2 ------------------------------
(兀-X])(兀-X 2) (X]-兀)(X] -X 2) ■ (兀—X (l )(x 2 -X])
= -0.616707 .
4 K(x) = /(x)—厶(x) = £/〃?(x —Xo)(x — xj 其中
I R } (x) \< — max I cos7x) I - max I (x-x n )(x -x) I
所以总误差界
2泸吟 心网
2 /
、2
訂显』丄x 工]“曲”
2
4
8(60 180丿
而无论用方程一还是方程二代入消元均解得召=100,勺=1.00,结果十分可 靠。
Ay
s
数值分析上机作业
数值分析上机实验报告 选题:曲线拟合的最小二乘法 指导老师: 专业: 学号: 姓名:
课题八曲线拟合的最小二乘法 一、问题提出 从随机的数据中找出其规律性,给出其近似表达式的问题,在生产实践和科学实验中大量存在,通常利用数据的最小二乘法求得拟合曲线。 在某冶炼过程中,根据统计数据的含碳量与时间关系,试求含碳量y 与时间t 的拟合曲线。 二、要求 1、用最小二乘法进行曲线拟合; 2、近似解析表达式为()33221t a t a t a t ++=?; 3、打印出拟合函数()t ?,并打印出()j t ?与()j t y 的误差,12,,2,1 =j ; 4、另外选取一个近似表达式,尝试拟合效果的比较; 5、*绘制出曲线拟合图*。 三、目的和意义 1、掌握曲线拟合的最小二乘法; 2、最小二乘法亦可用于解超定线代数方程组; 3、探索拟合函数的选择与拟合精度间的关系。 四、计算公式 对于给定的测量数据(x i ,f i )(i=1,2,…,n ),设函数分布为 ∑==m j j j x a x y 0)()(? 特别的,取)(x j ?为多项式 j j x x =)(? (j=0, 1,…,m )
则根据最小二乘法原理,可以构造泛函 ∑∑==-=n i m j i j j i m x a f a a a H 1 10))((),,,(? 令 0=??k a H (k=0, 1,…,m ) 则可以得到法方程 ???? ??????? ?=????????????????????????),(),(),(),(),(),(),(),(),(),(),(),(1010101111000100m m m m m m m m f f f a a a ????????????????????? 求该解方程组,则可以得到解m a a a ,,,10 ,因此可得到数据的最小二乘解 ∑=≈m j j j x a x f 0)()(? 曲线拟合:实际工作中,变量间未必都有线性关系,如服药后血药浓度与时间的关系;疾病疗效与疗程长短的关系;毒物剂量与致死率的关系等常呈曲线关系。曲线拟合是指选择适当的曲线类型来拟合观测数据,并用拟合的曲线方程分析两变量间的关系。 五、结构程序设计 在程序结构方面主要是按照顺序结构进行设计,在进行曲线的拟合时,为了进行比较,在程序设计中,直接调用了最小二乘法的拟合函数polyfit ,并且依次调用了plot 、figure 、hold on 函数进行图象的绘制,最后调用了一个绝对值函数abs 用于计算拟合函数与原有数据的误差,进行拟合效果的比较。
最新应用数值分析第四版第一章课后作业答案
第一章 1、 在下列各对数中,x 是精确值 a 的近似值。 3 .14,7/100)4(143 .0,7/1)2(0031 .0,1000/)3(1 .3,)1(========x a x a x a x a ππ 试估计x 的绝对误差和相对误差。 解:(1)0132.00416 .01.3≈= ≈-= -=a e e x a e r π (2)0011.00143 .0143.07/1≈= ≈-=-=a e e x a e r (3)0127.000004 .00031.01000/≈= ≈-=-=a e e x a e r π (4)001.00143 .03.147/100≈= ≈-=-=a e e x a e r 2. 已知四个数:x 1=26.3,x 2=0.0250, x 3= 134.25,x 4=0.001。试估计各近似数的有效位数和误差限,并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。 解:x 1=26.3 n=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2 x 2=0.0250 n=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2 x 3= 134.25 n=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10 -4 x 4=0.001 n=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5 由公式:e r (μ)= e (μ)/∣μ∣≦1/∣μ∣Σn i=1∣?f/?x i ∣δx i e r (μ1)≦1/∣μ1∣[x 2 x 3δx 1+ x 1 x 3δx 2 +x 1 x 2δx 3] =0.34468/88.269275 =0.0039049 e r (μ2)≦1/∣μ2∣[x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3 / x 1δx 4] =0.501937 3、设精确数a>0,x 是a的近似值,x 的相对误差限是0.2,求㏑x 的相对误差限。 解:设=()u f x , ()()()()() ()||||||||||()||()|| | |()||()||||r r r x e u df x e x df x e x e u u dx u dx u x df x x df x x e x x dx u dx u δ= ≈==≤ ()||10.2 (())| |()||ln ln ln r r r r df x x x x f x x x dx u x x x x δδδδ==??==
数值分析试题及答案汇总
数值分析试题 一、 填空题(2 0×2′) 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值,则x 有 2 位 有效数字。 2. 若f (x )=x 7-x 3+1,则f [20,21,22,23,24,25,26,27]= 1 , f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设,‖A ‖∞=___5 ____,‖X ‖∞=__ 3_____, ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ,则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时,若所求节点靠近首节点,应该选用等距节点下牛顿差商 公式的 前插公式 ,若所求节点靠近尾节点,应该选用等距节点下牛顿差商公式的 后插公式 ;如果要估计结果的舍入误差,应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是:=∑=n i i x a 0)( 1 ;所以当 系数a i (x )满足 a i (x )>1 ,计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%,至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ,线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的,其中的残差 r i = (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ,(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时,若在迭代区间存在唯一解,且f (x )
数值分析习题集及答案[1].(优选)
数值分析习题集 (适合课程《数值方法A 》和《数值方法B 》) 长沙理工大学 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2%,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出 它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=( n=1,2,…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对 误差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若
数值分析上机作业
昆明理工大学工科研究生《数值分析》上机实验 学院:材料科学与工程学院 专业:材料物理与化学 学号:2011230024 姓名: 郑录 任课教师:胡杰
P277-E1 1.已知矩阵A= 10787 7565 86109 75910 ?? ?? ?? ?? ?? ??,B= 23456 44567 03678 00289 00010 ?? ?? ?? ?? ?? ?? ?? ?? ,错误!未找到引用源。 = 11/21/31/41/51/6 1/21/31/41/51/61/7 1/31/41/51/61/71/8 1/41/51/61/71/81/9 1/51/61/71/81/91/10 1/61/71/81/91/101/11?????????????????? (1)用MA TLAB函数“eig”求矩阵全部特征值。 (2)用基本QR算法求全部特征值(可用MA TLAB函数“qr”实现矩阵的QR分解)。解:MA TLAB程序如下: 求矩阵A的特征值: clear; A=[10 7 8 7;7 5 6 5;8 6 10 9;7 5 9 10]; E=eig(A) 输出结果: 求矩阵B的特征值: clear; B=[2 3 4 5 6;4 4 5 6 7;0 3 6 7 8;0 0 2 8 9;0 0 0 1 0]; E=eig(B) 输出结果:
求矩阵错误!未找到引用源。的特征值: clear; 错误!未找到引用源。=[1 1/2 1/3 1/4 1/5 1/6; 1/2 1/3 1/4 1/5 1/6 1/7; 1/3 1/4 1/5 1/6 1/7 1/8; 1/4 1/5 1/6 1/7 1/8 1/9;1/5 1/6 1/7 1/8 1/9 1/10; 1/6 1/7 1/8 1/9 1/10 1/11]; E=eig(错误!未找到引用源。) 输出结果: (2)A= 10 7877565861097 5 9 10 第一步:A0=hess(A);[Q0,R0]=qr(A0);A1=R0*Q0 返回得到: 第二部:[Q1,R1]=qr(A1);A2=R1*Q1
数值分析课后题答案
数值分析 第二章 2.当1,1,2x =-时,()0,3,4f x =-,求()f x 的二次插值多项式。 解: 0120121200102021101201220211,1,2, ()0,()3,()4;()()1 ()(1)(2)()()2()()1 ()(1)(2) ()()6 ()()1 ()(1)(1) ()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------= =-+-- 则二次拉格朗日插值多项式为 2 20 ()()k k k L x y l x ==∑ 0223()4() 14 (1)(2)(1)(1)23 537623 l x l x x x x x x x =-+=---+ -+= +- 6.设,0,1,,j x j n =L 为互异节点,求证: (1) 0()n k k j j j x l x x =≡∑ (0,1,,);k n =L (2) ()()0n k j j j x x l x =-≡∑ (0,1,,);k n =L 证明 (1) 令()k f x x = 若插值节点为,0,1,,j x j n =L ,则函数()f x 的n 次插值多项式为0 ()()n k n j j j L x x l x == ∑。
插值余项为(1)1() ()()()()(1)! n n n n f R x f x L x x n ξω++=-= + 又,k n ≤Q (1)()0 ()0 n n f R x ξ+∴=∴= 0()n k k j j j x l x x =∴=∑ (0,1,,);k n =L 0 000 (2)()() (())()()(()) n k j j j n n j i k i k j j j i n n i k i i k j j i j x x l x C x x l x C x x l x =-==-==-=-=-∑∑∑∑∑ 0i n ≤≤Q 又 由上题结论可知 ()n k i j j j x l x x ==∑ ()()0 n i k i i k i k C x x x x -=∴=-=-=∑原式 ∴得证。 7设[]2 (),f x C a b ∈且()()0,f a f b ==求证: 21 max ()()max ().8 a x b a x b f x b a f x ≤≤≤≤''≤- 解:令01,x a x b ==,以此为插值节点,则线性插值多项式为 10 101010 ()() ()x x x x L x f x f x x x x x --=+-- =() ()x b x a f a f b a b x a --=+-- 1()()0()0 f a f b L x ==∴=Q 又
数值分析作业思考题汇总
¥ 数值分析思考题1 1、讨论绝对误差(限)、相对误差(限)与有效数字之间的关系。 2、相对误差在什么情况下可以用下式代替 3、查阅何谓问题的“病态性”,并区分与“数值稳定性”的不同点。 4、取 ,计算 ,下列方法中哪种最好为什么(1)(3 3-,(2)(2 7-,(3) ()3 1 3+ ,(4) ()6 1 1 ,(5)99- , 数值实验 数值实验综述:线性代数方程组的解法是一切科学计算的基础与核心问题。求解方法大致可分为直接法和迭代法两大类。直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法,因此也称为精确法。当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时,Gauss消去法是目前最基本和常用的方法。如若系数矩阵具有某种特殊形式,则为了尽可能地减少计算量与存储量,需采用其他专门的方法来求解。 Gauss消去等同于矩阵的三角分解,但它存在潜在的不稳定性,故需要选主元素。对正定对称矩阵,采用平方根方法无需选主元。方程组的性态与方程组的条件数有关,对于病态的方程组必须采用特殊的方法进行求解。 数值计算方法上机题目1 1、实验1. 病态问题 实验目的: 算法有“优”与“劣”之分,问题也有“好”和“坏”之别。所谓坏问题就是问题本身的解对数据变化的比较敏感,反之属于好问题。希望读者通过本实验对此有一个初步的体会。 数值分析的大部分研究课题中,如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决,当然一般要付出一些代价(如耗用更多的机器时间、占用更多的存储空间等)。 $ r e x x e x x ** * ** - == 141 . ≈)61
东南大学数值分析上机作业汇总
东南大学数值分析上机作业 汇总 -标准化文件发布号:(9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII
数值分析上机报告 院系: 学号: 姓名:
目录 作业1、舍入误差与有效数 (1) 1、函数文件cxdd.m (1) 2、函数文件cddx.m (1) 3、两种方法有效位数对比 (1) 4、心得 (2) 作业2、Newton迭代法 (2) 1、通用程序函数文件 (3) 2、局部收敛性 (4) (1)最大δ值文件 (4) (2)验证局部收敛性 (4) 3、心得 (6) 作业3、列主元素Gauss消去法 (7) 1、列主元Gauss消去法的通用程序 (7) 2、解题中线性方程组 (7) 3、心得 (9) 作业4、三次样条插值函数 (10) 1、第一型三次样条插值函数通用程序: (10) 2、数据输入及计算结果 (12)
作业1、舍入误差与有效数 设∑ =-=N j N j S 2 2 11 ,其精确值为?? ? ??---1112321N N . (1)编制按从小到大的顺序1 1 131121222-? ??+-+-=N S N ,计算N S 的通用程序; (2)编制按从大到小的顺序()1 21 11111222-???+--+-=N N S N ,计算N S 的通用程序; (3)按两种顺序分别计算642101010,,S S S ,并指出有效位数; (4)通过本上机你明白了什么? 程序: 1、函数文件cxdd.m function S=cxdd(N) S=0; i=2.0; while (i<=N) S=S+1.0/(i*i-1); i=i+1; end script 运行结果(省略>>): S=cxdd(80) S= 0.737577 2、函数文件cddx.m function S=cddx (N) S=0; for i=N:-1:2 S=S+1/(i*i-1); end script 运行结果(省略>>): S=cddx(80) S= 0.737577 3、两种方法有效位数对比
应用数值分析(第四版)课后习题答案第9章
第九章习题解答 1.已知矩阵????? ???????=??????????=4114114114,30103212321A A 试用格希哥林圆盘确定A 的特征值的界。 解:,24)2(, 33)1(≤-≤-λλ 2.设T x x x x ),...,,(321=是矩阵A 属于特征值λ的特征向量,若i x x =∞, 试证明特征值的估计式∑≠=≤-n i j j ij ii a a 1λ. 解:,x Ax λ = ∞∞∞∞≤==x A x x Ax i λλ 由 i x x =∞ 得 i n in i ii i x x a x a x a λ=++++ 11 j n j i i ij i ii x a x a ∑≠==-1)(λ j n j i i ij j n j i i ij i ii x a x a x a ∑∑≠=≠=≤=-11λ ∑∑≠=≠=≤≤-n j i i ij i j n j i i ij ii a x x a a 11λ 3.用幂法求矩阵 ???? ??????=1634310232A 的强特征值和特征向量,迭代初值取T y )1,1,1()0(=。 解:y=[1,1,1]';z=y;d=0; A=[2,3,2;10,3,4;3,6,1]; for k=1:100 y=A*z; [c,i]=max(abs(y)); if y(i)<0,c=-c;end
z=y/c if abs(c-d)<0.0001,break; end d=c end 11.0000 =c ,0.7500) 1.0000 0.5000(z 10.9999 =c ,0.7500) 1.0000 0.5000(z 11.0003 =c ,0.7500) 1.0000 0.5000(z 10.9989=c ,0.7500) 1.0000 0.5000(z 11.0040 =c ,0.7498) 1.0000 0.5000(z 10.9859=c ,0.7506) 1.0000 0.5001(z 11.04981 =c ,0.7478) 1.0000 0.4995(z 10.8316 =c ,0.7574) 1.0000 0.5020(z 11.5839 =c ,) 0.7260 1.0000 0.4928 (z 9.4706 =c ,0.8261) 1.0000 0.5280(z 17 = c ,0.5882) 1.0000 0.4118(z 11T (11)10T (10)9T (9)8T (8)7T (7)6T (6)5T (5)4T (4)3T (3)2T (2)1T (1)=========== 强特征值为11,特征向量为T 0.7500) 1.0000 0.5000(。 4.用反幂法求矩阵???? ??????=111132126A 最接近6的特征值和特征向量,迭代初值取 T y )1,1,1()0(=。 解:y=[1,1,1]';z=y;d=0; A=[6,2,1;2,3,1;1,1,1]; for k=1:100 AA=A-6*eye(3); y=AA\z; [c,i]=max(abs(y)); if y(i)<0,c=-c;end z=y/c; if abs(c-d)<0.0001,break; end d=c end d=6+1/c
数值分析试题及答案
一、单项选择题(每小题3分,共15分) 1. 3.142和3.141分别作为π的近似数具有( )和( )位有效数字. A .4和3 B .3和2 C .3和4 D .4和4 2. 已知求积公式 ()()2 1 121 1()(2)636f x dx f Af f ≈ ++? ,则A =( ) A . 16 B .13 C .12 D .2 3 3. 通过点 ()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足( ) A . ()00l x =0, ()110l x = B . ()00l x =0, ()111l x = C .() 00l x =1,()111 l x = D . () 00l x =1,()111 l x = 4. 设求方程 ()0 f x =的根的牛顿法收敛,则它具有( )敛速。 A .超线性 B .平方 C .线性 D .三次 5. 用列主元消元法解线性方程组 1231231 220223332 x x x x x x x x ++=?? ++=??--=? 作第一次消元后得到的第3个方程( ). A . 232 x x -+= B .232 1.5 3.5 x x -+= C . 2323 x x -+= D . 230.5 1.5 x x -=- 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得 分 评卷人 二、填空题(每小题3分,共15分)
1. 设T X )4,3,2(-=, 则=1||||X ,2||||X = . 2. 一阶均差 ()01,f x x = 3. 已知3n =时,科茨系数()()() 33301213,88C C C ===,那么 () 33C = 4. 因为方程()420 x f x x =-+=在区间 []1,2上满足 ,所以()0f x =在区间 内有根。 5. 取步长0.1h =,用欧拉法解初值问题 ()211y y y x y ?'=+?? ?=? 的计算公式 . 填空题答案 1. 9和29 2. ()() 0101 f x f x x x -- 3. 1 8 4. ()()120 f f < 5. ()12 00.1 1.1,0,1,210.11k k y y k k y +???? ?=+? ?=+???? =??L 得 分 评卷人 三、计算题(每题15分,共60分) 1. 已知函数 21 1y x = +的一组数据: 求分 段线性插值函数,并计算 () 1.5f 的近似值. 计算题1.答案 1. 解 []0,1x ∈, ()1010.510.50110x x L x x --=?+?=---% []1,2x ∈,()210.50.20.30.81221x x L x x --=?+?=-+--%
Matlab作业3(数值分析)答案
Matlab作业3(数值分析) 机电工程学院(院、系)专业班组 学号姓名实验日期教师评定 1.计算多项式乘法(x2+2x+2)(x2+5x+4)。 答: 2. (1)将(x-6)(x-3)(x-8)展开为系数多项式的形式。(2)求解在x=8时多项 式(x-1)(x-2) (x-3)(x-4)的值。 答:(1) (2)
3. y=sin(x),x从0到2π,?x=0.02π,求y的最大值、最小值、均值和标准差。 4.设x=[0.00.30.8 1.1 1.6 2.3]',y=[0.500.82 1.14 1.25 1.35 1.40]',试求二次多项式拟合系数,并据此计算x1=[0.9 1.2]时对应的y1。解:x=[0.0 0.3 0.8 1.1 1.6 2.3]'; %输入变量数据x y=[0.50 0.82 1.14 1.25 1.35 1.40]'; %输入变量数据y p=polyfit(x,y,2) %对x,y用二次多项式拟合,得到系数p x1=[0.9 1.2]; %输入点x1 y1=polyval(p,x1) %估计x1处对应的y1 p = -0.2387 0.9191 0.5318 y1 = a) 1.2909
5.实验数据处理:已知某压力传感器的测试数据如下表 p为压力值,u为电压值,试用多项式 d cp bp ap p u+ + + =2 3 ) ( 来拟 合其特性函数,求出a,b,c,d,并把拟合曲线和各个测试数据点画在同一幅图上。解: >> p=[0.0,1.1,2.1,2.8,4.2,5.0,6.1,6.9,8.1,9.0,9.9]; u=[10,11,13,14,17,18,22,24,29,34,39]; x=polyfit(p,u,3) %得多项式系数 t=linspace(0,10,100); y=polyval(x,t); %求多项式得值 plot(p,u,'*',t,y,'r') %画拟和曲线 x = 0.0195 -0.0412 1.4469 9.8267
数值分析整理版试题及答案
数值分析整理版试题及答案
例1、 已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解: (1)k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x x x x x x l x x x x x x x ----= ==-------- ()()()()()()()() ()()021******* ()1211122x x x x x x l x x x x x x x --+-= ==-+---+- ()()()()()()()()()()0122021111 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ (2)一阶均差、二阶均差分别为
[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----= = =----=== --- k x ()k f x 一阶 二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+ +++-=+- 例2、 设2 ()32f x x x =++,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{} span 1,x Φ=的最佳平方逼近多项式。 解: 若{}span 1,x Φ=,则0()1x ?=,1()x x ?=,且()1x ρ=,这样,有
数值分析习题集及答案Word版
数值分析习题集 (适合课程《数值方法A 》和《数值方法B 》) 长沙理工大学 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2%,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对 误差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?
数值分析第四版习题及答案
第四版 数值分析习题 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2%,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=…) 计算到100Y .(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字. 8. 当N 充分大时,怎样求 2 11N dx x +∞ +? ? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±秒的误差,证明当t 增加时S 的绝对误 差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 ln(ln(x x =- 计算,求对数时误差有多大?
数值分析模拟试题
1、 方程组中,,则求解方程组的Jacobi 迭代与Gauss-Seidel 迭代均收敛的a 的范围是___________。 2、,则A 的LDL T 分解中,。 3、,则__________,_______________. 4、已 知,则用复合梯形公式计算求 得,用三点式求得____________. 5、,则_________ ,三点高斯求积公式______________. 6设* 2.40315x =是真值 2.40194x =的近似值,则* x 有________位有效数字。 7 3()1,[0,1,2,3]f x x x f =+-=设 则差商(均差)_____________,[0,1,2,3,4]f =________________。 8 求方程()x f x =根的牛顿迭代格式是__________________。 9.梯形求积公式和复化梯形公式都是插值型求积公式_____(对或错)。 10.牛顿—柯特斯求积公式的系数和()0n n k k C ==∑__________________。 11.用二次拉格朗日插值多项式2()sin0.34L x 计算的值。插值节点和相应的函数值是(0,0),(0.30,0.2955),(0.40,0.3894)。 12.用二分法求方程3()10[1.0,1.5]f x x x =--=在 区间内的一个根,误差限 210ε-=。 13.用列主元消去法解线性方程组 1231231 232346,3525,433032.x x x x x x x x x ++=??++=??++=? 14. 确定求积公式
012()()(0)()h h f x dx A f h A f A f h -≈-++? 。 中待定参数i A 的值(0,1,2)i =,使求积公式的代数精度尽量高;并指出此时求积公式的代数精度。 15、 试求使求积公式的代数精度 尽量高,并求其代数精度。 16.证明区间[a,b]上带权()x ρ的正交多项式(),1,2,n P x n = 的n 个根都是单根,且位于区间(a,b)内。 17.设()()[,],max ()n n a x b f x C a b M f x ≤≤∈=,若取 21cos ,1,2,,222k a b a b k x k n n +--=+= 作节点,证明Lagrange 插值余项有估计式21()max ()!2n n n a x b M b a R x n -≤≤-≤ 18用n=10的复化梯形公式计算时, (1)试用余项估计其误差 (2)用n=10的复化梯形公式计算出该积分的近似值。 19已知方程组AX =f,其中 (1)列出Jacobi 迭代法和Gauss-Seidel 迭代法的分量形式。 (2)求出Jacobi 迭代矩阵的谱半径,SOR 迭代法的最佳松弛参数 和SOR 法 的谱半径(可直接用现有结论) 20试确定常数A ,B ,C 和,使得数值积分公式 有尽可能高的代数精度。试问所得的数值积分公式代数精度是多少? 21证明方程=)(x f x 2-x -3=0在区间(2,3)内有且仅有一个根,并用迭代法求方程在区间(2,3)内的根,精确到小数点后4位。 22设f (1)=2,f (3)=4,f (4)=6,用拉格朗日插值法求f (x )的二次插值多项式P 2(x ),并求f (2)的近似值。
(完整版)数值计算方法上机实习题答案
1. 设?+=1 05dx x x I n n , (1) 由递推公式n I I n n 1 51+-=-,从0I 的几个近似值出发,计算20I ; 解:易得:0I =ln6-ln5=0.1823, 程序为: I=0.182; for n=1:20 I=(-5)*I+1/n; end I 输出结果为:20I = -3.0666e+010 (2) 粗糙估计20I ,用n I I n n 51 5111+- =--,计算0I ; 因为 0095.05 6 0079.01020 201 020 ≈<<≈??dx x I dx x 所以取0087.0)0095.00079.0(2 1 20=+= I 程序为:I=0.0087; for n=1:20 I=(-1/5)*I+1/(5*n); end I 0I = 0.0083 (3) 分析结果的可靠性及产生此现象的原因(重点分析原因)。 首先分析两种递推式的误差;设第一递推式中开始时的误差为000I I E '-=,递推过程的舍入误差不计。并记n n n I I E '-=,则有01)5(5E E E n n n -==-=-Λ。因为=20E 20020)5(I E >>-,所此递推式不可靠。而在第二种递推式中n n E E E )5 1(5110-==-=Λ,误差在缩小, 所以此递推式是可靠的。出现以上运行结果的主要原因是在构造递推式过程中,考虑误差是否得到控制, 即算法是否数值稳定。 2. 求方程0210=-+x e x 的近似根,要求4 1105-+?<-k k x x ,并比较计算量。 (1) 在[0,1]上用二分法; 程序:a=0;b=1.0; while abs(b-a)>5*1e-4 c=(b+a)/2;
数值分析作业
第二章 1. 题目:运用MATLAB编程实现牛顿迭代 2. 实验操作 1、打开MATLAB程序软件。 2、在MATLAB中编辑如下的M程序。 function [p1,err,k,y]=newton(f,df,p0,delta,max) %f 是要求根的方程(f(x)=0); %df 是f(x)的导数; %p0是所给初值,位于x*附近; %delta是给定允许误差; %max是迭代的最大次数; %p1是newton法求得的方程的近似解; %err是p0的误差估计; %k是迭代次数; p0 for k=1:max p1=p0-feval('f',p0)/feval('df',p0); err=abs(p1-p0); p0=p1; k p1 err y=feval('f',p1) if (err> newton('f','df',1.2,10^(-6),20) 3.实验结果
p0 = 1.2000 k =1 p1=1.1030 err=0.0970 y=0.0329 k= 2 p1=1.0524 err=0.0507 y=0.0084 k =3 p1=1.0264 err=0.0260 y=0.0021 k =4 p1=1.0133 err=0.0131 y=5.2963e-004 k =5 p1=1.0066 err=0.0066 y=1.3270e-004 k =6 p1=1.0033 err=0.0033 y=3.3211e-005 k =7 p1=1.0017 err=0.0017 y=8.3074e-006 k =8 p1=1.0008 err=8.3157e-004 y = 2.0774e-006 k =9 p1=1.0004 err=4.1596e-004 y =5.1943e-007 k=10 p1=1.0002 err=2.0802e-004 y= 1.2987e-007 k=11 p1=1.0001 err=1.0402e-004 y =3.2468e-008 k=12 p1=1.0001 err=5.2014e-005 y=8.1170e-009 k=13 p1=1.0000 err=2.6008e-005 y= 2.0293e-009 k=14 p1=1.0000 err=1.3004e-005 y=5.0732e-010 k=15 p1 =1.0000 err=6.5020e-006 y=1.2683e-010 k=16 p1 =1.0000 err=3.2510e-006 y=3.1708e-011 k=17 p1 =1.0000 err=1.6255e-006 y =7.9272e-012 k=18 p1 =1.0000 err =8.1279e-007 y= 1.9820e-012 ans = 1.0000 结果说明:经过18次迭代得到精确解为1,误差为8.1279e-007。
东南大学-数值分析上机题作业-MATLAB版
2015.1.9 上机作业题报告 JONMMX 2000
1.Chapter 1 1.1题目 设S N =∑1j 2?1 N j=2 ,其精确值为 )1 1 123(21+--N N 。 (1)编制按从大到小的顺序1 1 131121222-+ ??+-+-=N S N ,计算S N 的通用程序。 (2)编制按从小到大的顺序1 21 1)1(111222-+ ??+--+-= N N S N ,计算S N 的通用程序。 (3)按两种顺序分别计算64210,10,10S S S ,并指出有效位数。(编制程序时用单精度) (4)通过本次上机题,你明白了什么? 1.2程序 1.3运行结果
1.4结果分析 按从大到小的顺序,有效位数分别为:6,4,3。 按从小到大的顺序,有效位数分别为:5,6,6。 可以看出,不同的算法造成的误差限是不同的,好的算法可以让结果更加精确。当采用从大到小的顺序累加的算法时,误差限随着N 的增大而增大,可见在累加的过程中,误差在放大,造成结果的误差较大。因此,采取从小到大的顺序累加得到的结果更加精确。 2.Chapter 2 2.1题目 (1)给定初值0x 及容许误差ε,编制牛顿法解方程f(x)=0的通用程序。 (2)给定方程03 )(3 =-=x x x f ,易知其有三个根3,0,3321= *=*-=*x x x ○1由牛顿方法的局部收敛性可知存在,0>δ当),(0δδ+-∈x 时,Newton 迭代序列收敛于根x2*。试确定尽可能大的δ。 ○2试取若干初始值,观察当),1(),1,(),,(),,1(),1,(0+∞+-----∞∈δδδδx 时Newton 序列的收敛性以及收敛于哪一个根。 (3)通过本上机题,你明白了什么? 2.2程序