matlab：（1）做作业中遇到的问题

发表于 | 分类于 |
字数统计 | 阅读时长

1.求矩阵的逆用 inv()

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
>> a=1/27*[-15 9 -27 -33;1 -15 -9 -23;9 0 -18 -18;-7 -3 9 -1]
a =
   -0.5556    0.3333   -1.0000   -1.2222
    0.0370   -0.5556   -0.3333   -0.8519
    0.3333         0   -0.6667   -0.6667
   -0.2593   -0.1111    0.3333   -0.0370
>> b = inv(a)
b =
   -1.0000   -1.0000    3.0000    2.0000
   -0.0000   -2.5000    3.0000    3.5000
   -0.7500   -1.5000    3.0000    5.2500
    0.2500    1.0000   -3.0000   -4.2500

2.显示分数用rats(b)

1
2
3
4
5
6
7
8
>> c = rats(b)
c =
      -1            -1             3             2      
       0           -5/2            3            7/2     
     -3/4          -3/2            3           21/4     
      1/4            1            -3          -17/4

3.用matlab求Ax=0的基础解析的解程序为:

1
2
3
A=[1 1 1 1 -3 -1 1;1 0 0 0 1 1 0;-2 0 0 -1 0 -1 -2];
r=rank(A);
y=null(A, ‘r’ )

得到解为:

1
2
3
4
5
6
7
y=[    0      -1      -1       0
      -1       2       1       1
       1       0       0       0
       0       2       1      -2
       0       1       0       0
       0       0       1       0
       0       0       0       1]

其列向量构成Ax=0的一个基础解系。