题目
解答
2-1
>> linspace(20,100,21)
ans =
1 至 12 列
20 24 28 32 36 40 44 48 52 56 60 64
13 至 21 列
68 72 76 80 84 88 92 96 100
2-2
>> x=logspace(0,3,10)
x =
1.0e+03 *
1 至 7 列
0.0010 0.0022 0.0046 0.0100 0.0215 0.0464 0.1000
8 至 10 列
0.2154 0.4642 1.0000
2-3
>> for i=2:2:18;
a(10-i/2)=sin(pi/i);
end
>> a
a =
1 至 7 列
0.1736 0.1951 0.2225 0.2588 0.3090 0.3827 0.5000
8 至 9 列
0.7071 1.0000
>> b=a';
>> sum(b)
ans =
3.7489
2-4
>> x=1:10;
>> y=(x+x.^2+x.^3)./(1+sin(x)/2)
y =
1.0e+03 *
1 至 7 列
0.0021 0.0096 0.0364 0.1351 0.2978 0.2999 0.3003
8 至 10 列
0.3907 0.6791 1.5247
2-5
>> A(1,:)=[1 0 0 0 1];
>> A(5,:)=A(1,:);
>> A(2,:)=[0 1 1 1 0];
>> A(3,:)=A(2,:);
>> A(4,:)=A(2,:);
>> A
A =
1 0 0 0 1
0 1 1 1 0
0 1 1 1 0
0 1 1 1 0
1 0 0 0 1
>> det(A)
ans =
0
2-6
>> A=[1 1 1 1;1 2 2 1;3 1 1 3];
>> B=[1 2 3;3 2 1;1 3 5;5 3 1];
>> A*B
ans =
10 10 10
14 15 16
22 20 18
>> B*A
ans =
12 8 8 12
8 8 8 8
19 12 12 19
11 12 12 11
2-7
>> B=[2 3 4 5;3 4 5 6;4 5 6 7;5 6 7 8];
>> A=magic(4);
>> A.*B
ans =
32 6 12 65
15 44 50 48
36 35 36 84
20 84 105 8
2-8
>> A=[1 2 3 4 5;1 3 6 10 15;1 4 10 20 25];
>> B=[17 24 15;23 5 16;4 6 22;10 12 3;11 18 9];
>> C=A*B
C =
170 190 170
375 465 360
624 794 584
>> det(C)
ans =
-6.5700e+04
2-9
>> A=[64 118 155;121 220 290; 178 322 425];
>> B=[3 4 5;6 10 15;10 20 25];
>> A/B
ans =
2.0000 3.0000 4.0000
5.0000 6.0000 7.0000
8.0000 9.0000 10.0000
2-10
>> A=[1 1 2;1 3 4;2 4 5];
>> B=[14 18;30 36;39 50];
>> A\B
ans =
1 5
3 5
5 4
2-11
(1)
>> det([2 0 17 1;1 5 3 3;3 1 -1 2;-5 4 8 11])
ans =
-2.9210e+03
( 2 )同
2-12
>> A=[1 1 1 1;8 4 2 1;27 9 3 1;64 16 4 1];
>> A'
ans =
1 8 27 64
1 4 9 16
1 2 3 4
1 1 1 1
2-13
(1)
>> A=[61 92 17 93 41;79 73 40 91 89];
>> size(A)
ans =
2 5
(2)同
2-14
(1)
>> rank([16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 1])
ans =
3
(2)同
2-15
(1)
>> V=[1 3 5 7];
>> norm(V,2)
ans =
9.1652
(2)同
2-16
(1)
>> norm([3 3 6 2 4;3 9 12 6 8;6 12 15 8 10],2)
ans =
31.1679
(2)同
2-17
(1)
>> cond([5 4 0 0;1 5 4 0;0 1 5 4;0 0 1 5])
ans =
7.0296
(2)同
2-18
(1)
>> rank([1 2 3;4 5 6;7 9 11;10 13 16;14 17 20])
ans =
2
(2)同
2-19
(1)
>> [V,D]=eig([8 1 6;3 5 7;4 9 2])
V =
-0.5774 -0.8131 -0.3416
-0.5774 0.4714 -0.4714
-0.5774 0.3416 0.8131
D =
15.0000 0 0
0 4.8990 0
0 0 -4.8990
(2)同
2-20
(1)
>> A=ones(7);
>> A(2:6,2:6)=3*ones(5);
>> A(3:5,3:5)=5*ones(3);
>> A(4,4)=7*ones(1);
>> A
A =
1 1 1 1 1 1 1
1 3 3 3 3 3 1
1 3 5 5 5 3 1
1 3 5 7 5 3 1
1 3 5 5 5 3 1
1 3 3 3 3 3 1
1 1 1 1 1 1 1
(2)
>> A=diag([2 2 2],-2);
>> B=diag([3 3 3 3],-1);
>> C=diag([4 4 4 4 4]);
>> D=diag([5 5 5 5],1);
>> E=diag([6 6 6],2);
>> F=A+B+C+D+E
F =
4 5 6 0 0
3 4 5 6 0
2 3 4 5 6
0 2 3 4 5
0 0 2 3 4
2-21
>> A=[16 4 8 4;4 10 8 4;8 8 12 10;4 4 10 12];
>> [l,u]=lu(A)
l =
1.0000 0 0 0
0.2500 1.0000 0 0
0.5000 0.6667 0.6667 1.0000
0.2500 0.3333 1.0000 0
u =
16.0000 4.0000 8.0000 4.0000
0 9.0000 6.0000 3.0000
0 0 6.0000 10.0000
0 0 0 -0.6667
>> det(A)
ans =
576.0000
2-22
2-23
(1)
>> inv([1 2 8;1 3 5;2 1 5])
ans =
-0.5000 0.1000 0.7000
-0.2500 0.5500 -0.1500
0.2500 -0.1500 -0.0500
(2)同
2-24
(1)
>> A=[1 1 1 1;1 2 3 4;1 3 6 10;1 4 10 20];
>> inv(A)
ans =
4.0000 -6.0000 4.0000 -1.0000
-6.0000 14.0000 -11.0000 3.0000
4.0000 -11.0000 10.0000 -3.0000
-1.0000 3.0000 -3.0000 1.0000
>> rref([A,eye(size(A))])
ans =
1 0 0 0 4 -6 4 -1
0 1 0 0 -6 14 -11 3
0 0 1 0 4 -11 10 -3
0 0 0 1 -1 3 -3 1
(2)同
2-25
(1)
>> pinv([1 2 3;2 5 8;2 4 6;4 10 16;3 6 9])
ans =
0.2500 -0.2667 0.5000 -0.5333 0.7500
0.0714 -0.0667 0.1429 -0.1333 0.2143
-0.1071 0.1333 -0.2143 0.2667 -0.3214
(2)同