task:17 第二章习题

题目

解答

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）同