通用性的特殊矩阵和用于专门学科的特殊矩阵
1.通用的特殊矩阵
zeros(m): m*m零矩阵
zeros(m,n): m*n零矩阵
zeros(size(A)): 与矩阵A同样大小的零矩阵
>> A=zeros(2,3)
A =
0 0 0
0 0 0
>> zeros(size(reshape(A,3,2)))
ans =
0 0
0 0
0 0
例1 产生5阶两位随机整数矩阵A,再产生均值为0.6、方差为0.1的5阶正态分布随机矩阵B,验证(A+B)I=IA+BI(I为单位矩阵)。
rand(): (0,1)均匀分布的随机数x
fix(a+(b-a+1)*x): [a,b]均匀分布的随机整数
randn(): 均值为0,方差为1的标准正态分布随机数x
ɑ+ɓx,得到均值为ɑ,方差为ɓ^2的随机数
>> A=fix(10+(99-10+1)*rand(5))
A =
78 73 84 49 54
76 12 72 44 50
45 34 38 78 68
68 14 95 81 73
25 18 13 26 77
>> B=0.6+sqrt(0.1)*randn(5)
B =
0.3269 1.0847 0.2556 0.6272 0.4053
0.6245 0.3566 0.6103 0.1283 0.8366
0.2161 0.7174 0.7747 0.3653 0.5392
0.2479 0.5287 0.9480 0.2643 0.8810
0.5978 0.9533 1.0883 1.3433 0.3581
>> C=eye(5)
C =
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0
0 0 0 1 0
0 0 0 0 1
>> (A+B)*C==C*A+B*C
ans =
5×5 logical 数组
1 1 1 1 1
1 1 1 1 1
1 1 1 1 1
1 1 1 1 1
1 1 1 1 1
2.专门学科的特殊矩阵
>> clear
>> M=magic(3)
M =
8 1 6
3 5 7
4 9 2
>> clear
>> A=vander(1:5)
A =
1 1 1 1 1
16 8 4 2 1
81 27 9 3 1
256 64 16 4 1
625 125 25 5 1
>> clear
>> format rat
>> H=hilb(4)
H =
1 1/2 1/3 1/4
1/2 1/3 1/4 1/5
1/3 1/4 1/5 1/6
1/4 1/5 1/6 1/7
>> p=[1,-2,-5,6];
>> A=compan(p)
A =
2 5 -6
1 0 0
0 1 0
>> format rat
>> clear
>> format rat
>> P=pascal(5)
P =
1 1 1 1 1
1 2 3 4 5
1 3 6 10 15
1 4 10 20 35
1 5 15 35 70
>> inv(P)
ans =
5 -10 10 -5 1
-10 30 -35 19 -4
10 -35 46 -27 6
-5 19 -27 17 -4
1 -4 6