Task03:第二章习题

在这里插入图片描述


>> x=linspace(20,100,21)

x =

    20    24    28    32    36    40    44    48    52    56    60    64    68    72    76    80    84    88    92    96   100

在这里插入图片描述

>> x=logspace(0,3,10)

x =

   1.0e+03 *

    0.0010    0.0022    0.0046    0.0100    0.0215    0.0464    0.1000    0.2154    0.4642    1.0000

在这里插入图片描述

>> x=[10 20 30 40 50 60 70 80 90];
>> y=(x/180)*pi

y =

    0.1745    0.3491    0.5236    0.6981    0.8727    1.0472    1.2217    1.3963    1.5708

>> z=sin(y)

z =

    0.1736    0.3420    0.5000    0.6428    0.7660    0.8660    0.9397    0.9848    1.0000

>> sum(z)

ans =

    6.2150

在这里插入图片描述


>> x=[1 2 3 4 5 6 7 8 9 10];
>>  y=(x+x.^2+x.^3)./(1+sin(x./2))

y =

   1.0e+04 *

    0.0002    0.0008    0.0020    0.0044    0.0097    0.0226    0.0615    0.2401    3.6449    2.7023

在这里插入图片描述

>> 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] 

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

在这里插入图片描述

>> A=[1 1 1 1;1 2 2 1;3 1 1 3]

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]

B =

     1     2     3
     3     2     1
     1     3     5
     5     3     1

>> C=A*B

C =

    10    10    10
    14    15    16
    22    20    18

在这里插入图片描述

>> A=magic(4)

A =

    16     2     3    13
     5    11    10     8
     9     7     6    12
     4    14    15     1

>> B=[2 3 4 5;3 4 5 6;4 5 6 7;5 6 7 8]

B =

     2     3     4     5
     3     4     5     6
     4     5     6     7
     5     6     7     8

>> 
>> C=A.*B

C =

    32     6    12    65
    15    44    50    48
    36    35    36    84
    20    84   105     8

在这里插入图片描述

>> A=[1 2 3 4 5;1 3 6 10 15;1 4 10  20 25]

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]

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

在这里插入图片描述

>> A=[64 118 155;121 220 290;178 322 425]

A =

    64   118   155
   121   220   290
   178   322   425

>> B=[3 4 5;6 10 15;10 20 25]

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

在这里插入图片描述

>> A=[1 1 2;1 3 4; 2 4 5]

A =

     1     1     2
     1     3     4
     2     4     5

>> B=[14 18;30 36;39 50]

B =

    14    18
    30    36
    39    50

>> A\B

ans =

     1     5
     3     5
     5     4

在这里插入图片描述

>> A=[2 0 17 1;1 5 3 3;3 1 -1 2;-5 4 8 11]

A =

     2     0    17     1
     1     5     3     3
     3     1    -1     2
    -5     4     8    11

>> det(A)

ans =

  -2.9210e+03

>> A=[1 3 5 7 9;2 4 6 8 10;5 6 7 6 5;10 8 6 4 2;9 7 5 3 1]

A =

     1     3     5     7     9
     2     4     6     8    10
     5     6     7     6     5
    10     8     6     4     2
     9     7     5     3     1

>> det(A)

ans =

     0

在这里插入图片描述

>> A=[1 1 1 1;8 4 2 1;27 9 3 1;64 16 4 1]

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

在这里插入图片描述

>> A=[61 92 17 93 41;79 73 40 91 89]

A =

    61    92    17    93    41
    79    73    40    91    89

>> size(A)

ans =

     2     5
>> A=[-43 28 118 17;-166 -114 -3 -18;12 119 32 72]

A =

   -43    28   118    17
  -166  -114    -3   -18
    12   119    32    72

>> size(A)

ans =

     3     4

在这里插入图片描述

>> A=[16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 1]

A =

    16     2     3    13
     5    11    10     8
     9     7     6    12
     4    14    15     1

>> rank(A)

ans =

     3         %非满秩
>> A=[100 50 33 25 20;50 33 25 20 16;33 25 20 16 14;25 20 16 14 12;20 16 14 12 11]

A =

   100    50    33    25    20
    50    33    25    20    16
    33    25    20    16    14
    25    20    16    14    12
    20    16    14    12    11

>> rank(A)

ans =

     5            %满秩

在这里插入图片描述


>> v=[1 3 5 7]

v =

     1     3     5     7

>> norm(v,2)

ans =

    9.1652

>> v=[1 -2 5 8];
>> norm(v,2)

ans =

    9.6954

在这里插入图片描述

>> A=[3 3 6 2 4;3 9 12 6 8;6 12 15 8 10]

A =

     3     3     6     2     4
     3     9    12     6     8
     6    12    15     8    10

>> norm(A,2)

ans =

   31.1679

>> max(svd(A))

ans =

   31.1679
>> A=[1 0.5 0.3333 0.25;0.5 0.3333 0.25 0.2;0.3333 0.25 0.2 0.1667;0.25 0.2 0.1667 0.1429]

A =

    1.0000    0.5000    0.3333    0.2500
    0.5000    0.3333    0.2500    0.2000
    0.3333    0.2500    0.2000    0.1667
    0.2500    0.2000    0.1667    0.1429

>> norm(A,2)

ans =

    1.5002

在这里插入图片描述


>> A=[5 4 0 0;1 5 4 0;0 1 5 4;0 0 1 5]

A =

     5     4     0     0
     1     5     4     0
     0     1     5     4
     0     0     1     5

>> cond(A)

ans =

    7.0296
>> A=[0 4 4 1;1 5 4 0;0 1 5 4;0 0 1 5]

A =

     0     4     4     1
     1     5     4     0
     0     1     5     4
     0     0     1     5

>> cond(A)

ans =

   17.3128

在这里插入图片描述

>> A=[1 2 3;4 5 6;7 9 11;10 13 16;14 17 20]

A =

     1     2     3
     4     5     6
     7     9    11
    10    13    16
    14    17    20

>> rank(A)

ans =

     2

>> A=[16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 1]

A =

    16     2     3    13
     5    11    10     8
     9     7     6    12
     4    14    15     1

>> rank(A)

ans =

     3

在这里插入图片描述

>> A=[8 1 6;3 5 7;4 9 2]

A =

     8     1     6
     3     5     7
     4     9     2

>> d=eig(A)

d =

   15.0000
    4.8990
   -4.8990

>> [V,D]=eig(A)

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
>> A=[-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0]

A =

   -10   -35   -50   -24
     1     0     0     0
     0     1     0     0
     0     0     1     0

>> d=eig(A)

d =

   -4.0000
   -3.0000
   -2.0000
   -1.0000

>> [V,D]=eig(A)

V =

   -0.9683    0.9429    0.8677   -0.5000
    0.2421   -0.3143   -0.4339    0.5000
   -0.0605    0.1048    0.2169   -0.5000
    0.0151   -0.0349   -0.1085    0.5000


D =

   -4.0000         0         0         0
         0   -3.0000         0         0
         0         0   -2.0000         0
         0         0         0   -1.0000

在这里插入图片描述

>> v=ones(7);
>> v(2:6,2:6)=ones(5)*3;v(3:5,3:5)=ones(3)*5;v(4,4)=7

v =

     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
 >> a=[4 4 4 4 4];b=[5 5 5 5];c=[6 6 6];d=[3 3 3 3];e=[2 2 2 ];
>> v=diag(a)+diag(b,1)+diag(c,2)+diag(d,-1)+diag(e,-2)

v =

     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

在这里插入图片描述

>> A=[16 4 8 4;4 10 8 4;8 8 12 10;4 4 10 12]

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

在这里插入图片描述

>> syms A
syms B
syms C 
syms D
>> G1=horzcat(A,C)
 
G1 =
 
[ A, C]
 
>> G2=horzcat(D,B)
 
G2 =
 
[ D, B]
 
>> G=vertcat(G1,G2)
 
G =
 
[ A, C]
[ D, B]
 
>> inv(G)
 
ans =
 
[  B/(A*B - C*D), -C/(A*B - C*D)]
[ -D/(A*B - C*D),  A/(A*B - C*D)]

在这里插入图片描述

>> v=[1 2 8;1 3 5;2 1 5]

v =

     1     2     8
     1     3     5
     2     1     5

>> inv(v)

ans =

   -0.5000    0.1000    0.7000
   -0.2500    0.5500   -0.1500
    0.2500   -0.1500   -0.0500
>> v=[5 8 13 0 0;2 5 8 13 0;1 2 5 8 13;0 1 2 5 8;0 0 1 2 5]

v =

     5     8    13     0     0
     2     5     8    13     0
     1     2     5     8    13
     0     1     2     5     8
     0     0     1     2     5

>> inv(v)

ans =

   -0.2857   -0.3878    3.2041   -2.1837   -4.8367
    0.0714   -0.1888    0.0204    1.3316   -2.1837
    0.1429    0.2653   -1.2449    0.0204    3.2041
   -0.0714    0.0459    0.2653   -0.1888   -0.3878
         0   -0.0714    0.1429    0.0714   -0.2857

在这里插入图片描述

>> v=[1 1 1 1;1 2 3 4;1 3 6 10;1 4 10 20]

v =

     1     1     1     1
     1     2     3     4
     1     3     6    10
     1     4    10    20

>> w=[v,eye(size(v))]

w =

     1     1     1     1     1     0     0     0
     1     2     3     4     0     1     0     0
     1     3     6    10     0     0     1     0
     1     4    10    20     0     0     0     1

>> rref(w)

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

>> ans(:,5:8)

ans =

     4    -6     4    -1
    -6    14   -11     3
     4   -11    10    -3
    -1     3    -3     1

>> inv(v)

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
>> v=[11 13 0 0;8 11 13 0;0 8 11 13;0 0 8 11]

v =

    11    13     0     0
     8    11    13     0
     0     8    11    13
     0     0     8    11

>> w=[v,eye(size(v))]

w =

    11    13     0     0     1     0     0     0
     8    11    13     0     0     1     0     0
     0     8    11    13     0     0     1     0
     0     0     8    11     0     0     0     1

>> rref(w)

ans =

    1.0000         0         0         0    0.0778    0.0180   -0.1512    0.1787
         0    1.0000         0         0    0.0111   -0.0152    0.1279   -0.1512
         0         0    1.0000         0   -0.0573    0.0787   -0.0152    0.0180
         0         0         0    1.0000    0.0416   -0.0573    0.0111    0.0778

>> ans(:,5:8)

ans =

    0.0778    0.0180   -0.1512    0.1787
    0.0111   -0.0152    0.1279   -0.1512
   -0.0573    0.0787   -0.0152    0.0180
    0.0416   -0.0573    0.0111    0.0778

>> inv(v)

ans =

    0.0778    0.0180   -0.1512    0.1787
    0.0111   -0.0152    0.1279   -0.1512
   -0.0573    0.0787   -0.0152    0.0180
    0.0416   -0.0573    0.0111    0.0778

在这里插入图片描述

>> v=[1 2 3;2 5 8;2 4 6;4 10 16;3 6 9]

v =

     1     2     3
     2     5     8
     2     4     6
     4    10    16
     3     6     9

>> rank(v)

ans =

     2

>> pinv(v)

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
>> v=[1 2 3 0 0;0 3 5 6 0;0 0 5 8 9]

v =

     1     2     3     0     0
     0     3     5     6     0
     0     0     5     8     9

>> rank(v)

ans =

     3

 
>> pinv(v)

ans =

    0.1390   -0.0524    0.0102
    0.1210    0.0321   -0.0244
    0.2064   -0.0039    0.0129
   -0.2324    0.1539    0.0015
    0.0920   -0.1346    0.1026
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值