数值分析计算机模拟.docVIP

数 值 分 析 作 业 1已知矩阵 ,,, 求：a + b, ab , ac, ca. 解 键入： a=[1 2 3;5 4 1;2 5 9]; b=[3 5;6 1;2 9]; c=[6 1 2;5 2 8;4 7 1]; a1=a+c, a2=a*b, a3=a*c, a4=c*a a1 = 7 3 5 10 6 9 6 12 10 a2 = 21 34 41 38 54 96 a3 = 28 26 21 54 20 43 73 75 53 a4 = 15 26 37 31 58 89 41 41 28 2已知矩阵 , 用数组算法计算：a. * c, a* c, a.^3和 a^3。 解：键入 a=[1 2 3;5 4 1;2 5 9];c=[6 1 2;5 2 8;4 7 1]; a.*c ans = 6 2 6 25 8 8 8 35 9 ac1=a* c ac1 = 28 26 21 54 20 43 73 75 53 a3=a.^3,a4=a^3 a3 = 1 8 27 125 64 1 8 125 729 a4 = 206 294 364 238 318 364 574 826 1032 3. 设计算法中的注意事项（防止大数吃小数）。 设a=1e+308，b=1.1e+308，c=-1.001e+308，用matlab计算(a+b)+c，a+(b+c)。 解：键入 a=1e+308;b=1.1e+308;c=-1.001e+308;(a+b)+c ans = Inf a+(b+c) ans = 1.0990e+308 4.求方程的根。 解：键入 roots([1 -4 6 -4 1]) ans = 1.0002 1.0000 + 0.0002i 1.0000 - 0.0002i 0.9998 5.对矩阵求奇异值分解。 解：键入 a=[1 2 3;2 3i 4;3 2 5]; [u,s,v]=svd(a) u = -0.4100 + 0.0525i -0.1219 + 0.3923i -0.8068 - 0.0974i -0.5657 - 0.1051i -0.0021 - 0.7846i -0.1241 + 0.1948i -0.7038 + 0.0525i 0.2484 + 0.3923i 0.5261 - 0.0974i s = 8.5212 0 0 0 2.8341 0 0 0 0.5972 v = -0.4287 0.2184 0.8767 -0.2984 + 0.2238i -0.7412 + 0.5559i 0.0388 - 0.0291i -0.8229 + 0.0000i 0.3062 - 0.0000i -0.4787 + 0.0000i 6.已知y=sin x的表格为： x 0 π/6 π/4 π/3 y 0 0.5000 0.7010 0.8660 分别用三次多项式插值，三次样条插值，lagrange插值求sin(π/12),sin(2π/9),sin(π/2)。 解：(1)三次多项式插值 键入; x=[0;pi/6;pi/4;pi/3]; y=[0;0.5;0.701;0.866]; X=[pi/12;2*pi/9;pi/2]; Y=interpl(x,y,X,cuhic) Y= 0.2656 0.6383 1.0668 (2)三次样条插值 键入; x=[0;pi/6;pi/4;pi/3]; y=[0;0.5;0.701;0.866]; X=[pi