《matlab数学实验习题全部答案(胡良剑)》由会员分享，可在线阅读，更多相关《matlab数学实验习题全部答案(胡良剑)（33页珍藏版）》请在金锄头文库上搜索。
1、数学实验答案%Page20,ex1 (5) 等于exp(1),exp(2);exp(3),exp(4) (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1=4,false, 2=3,false, 3=2, ture, 4=1,ture (11) 答案表明：编址第2元素满足不等式(30=20)和编址第4元素满足不等式(40=10) (12) 答案表明：编址第2行第1列元素满足不等式(30=20)和编址第2行第2列元素满足不等式(40=10) %Page20, ex2 (1)a, b, c的值尽管都是1， 但数据类型分别为数值，字符， 逻辑， 注意a与c相
2、等， 但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 %Page20,ex3 r=2;p=0.5;n=12; T=log(r)/n/log(1+0.01*p) T = 11.5813 %Page20,ex4 x=-2:0.05:2;f=x.4-2.x; fmin,min_index=min(f) fmin = -1.3907 %最小值 min_index = 54 %最小值点编址 x(min_index) ans = 0.6500 %最小值点 f1,x1_index=min(abs(f) %求近似根-绝对值最小的点 f1 = 0.0328 x1
3、_index = 24 x(x1_index) ans = -0.8500 x(x1_index)=;f=x.4-2.x; %删去绝对值最小的点以求函数绝对值次小的点 f2,x2_index=min(abs(f) %求另一近似根-函数绝对值次小的点 f2 = 0.0630 x2_index = 65 x(x2_index) ans = 1.2500 %Page20,ex5 z=magic(10) z = 92 99 1 8 15 67 74 51 58 40 98 80 7 14 16 73 55 57 64 41 4 81 88 20 22 54 56 63 70 47 85 87 19 2
4、1 3 60 62 69 71 28 86 93 25 2 9 61 68 75 52 34 17 24 76 83 90 42 49 26 33 65 23 5 82 89 91 48 30 32 39 66 79 6 13 95 97 29 31 38 45 72 10 12 94 96 78 35 37 44 46 53 11 18 100 77 84 36 43 50 27 59 sum(z) ans = 505 505 505 505 505 505 505 505 505 505 sum(diag(z) ans = 505 z(:,2)/sqrt(3) ans = 57.1577
5、46.1880 46.7654 50.2295 53.6936 13.8564 2.8868 3.4641 6.9282 10.3923 z(8,:)=z(8,:)+z(3,:) z = 92 99 1 8 15 67 74 51 58 40 98 80 7 14 16 73 55 57 64 41 4 81 88 20 22 54 56 63 70 47 85 87 19 21 3 60 62 69 71 28 86 93 25 2 9 61 68 75 52 34 17 24 76 83 90 42 49 26 33 65 23 5 82 89 91 48 30 32 39 66 83 8
6、7 101 115 119 83 87 101 115 119 10 12 94 96 78 35 37 44 46 53 11 18 100 77 84 36 43 50 27 59 %Page 40 ex1 先在编辑器窗口写下列M函数，保存为eg2_1.m function xbar,s=ex2_1(x) n=length(x); xbar=sum(x)/n; s=sqrt(sum(x.2)-n*xbar2)/(n-1); 例如 x=81 70 65 51 76 66 90 87 61 77; xbar,s=ex2_1(x) xbar = 72.4000 s = 12.1124 %Page
7、 40 ex2 s=log(1);n=0; while se k=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end a,x,k 计算至k=21可满足精度 %Page 40 ex4 clear;tic;s=0; for i=1:1000000 s=s+sqrt(3)/2i; end s,toc tic;s=0;i=1; while i1.1)+x.*(x=-1.1)-1.1*(x1); p=p+b*exp(-y.2-6*x.2).*(x+y-1).*(x+y=1); p=p+a*exp(-0.75*y.2-3.75*x.2+1.5*x).*(x+y A=1 2
8、 3;4 5 6;7 8 0;C=2 -5 -22;-5 -24 -56;-22 -56 -16; X=lyap(A,C) X = 1.0000 -1.0000 -0.0000 -1.0000 2.0000 1.0000 -0.0000 1.0000 7.0000 %Chapter 3 %Exercise 1 a=1,2,3;b=2,4,3;a./b,a.b,a/b,ab ans = 0.5000 0.5000 1.0000 ans = 2 2 1 ans = 0.6552 %一元方程组x2,4,3=1,2,3的近似解 ans = 0 0 0 0 0 0 0.6667 1.3333 1.000
9、0 %矩阵方程1,2,3x11,x12,x13;x21,x22,x23;x31,x32,x33=2,4,3的特解 %Exercise 2(1) A=4 1 -1;3 2 -6;1 -5 3;b=9;-2;1; rank(A), rank(A,b) %A,b为增广矩阵 ans = 3 ans = 3 %可见方程组唯一解 x=Ab x = 2.3830 1.4894 2.0213 %Exercise 2(2) A=4 -3 3;3 2 -6;1 -5 3;b=-1;-2;1; rank(A), rank(A,b) ans = 3 ans = 3 %可见方程组唯一解 x=Ab x = -0.4706 -0.2941 0 %Exercise 2(3) A=4 1;3 2;1 -5;b=1;1;1; rank(A), rank(A,b) ans = 2 ans = 3 %可见方程组无解 x=Ab x = 0.3311 -0.1219 %最小二乘近似解 %Exercise 2(4) a=2,1,-1,1;1,2,1,-1;1,1,2,1;b=1 2 3;%注意b的写法 rank(a),rank(a,b) ans = 3 ans = 3 %rank(a)=rank(a,b) ab ans = 1
