1. 通用的特殊矩阵
zeros:产生全0矩阵
ones:全1矩阵
eye: 单位矩阵
rand : 0~1之间均匀分布的随机矩阵
randn: 均值为0,方差为1的标准正态分布随机矩阵
调用格式
zeros(m):产生 m x m 零矩阵
zeros(m,n): 产生m x n 零矩阵
zeros(size(A)): 产生与矩阵A同样大小的矩阵
2.
(1)若想得到在任意[a , b]区间上均匀分布的随机数,
则 yi = a + (b - a)xi
例:其中xi是通过rand获得的一组随机数
>> x = 20 + (50 - 20)*rand(5)
x =
44.4417 22.9262 24.7284 24.2566 39.6722
47.1738 28.3549 49.1178 32.6528 21.0714
23.8096 36.4064 48.7150 47.4721 45.4739
47.4013 48.7252 34.5613 43.7662 48.0198
38.9708 48.9467 44.0084 48.7848 40.3621
(2)若已知正态分布随机数xi , 想的到均值、方差,则利用
yi = u + qxi (u为均值 , q为方差)
>> y = 0.6 + sqrt(0.1)*randn(5)
y =
0.8809 1.0549 0.5677 0.5905 0.3269
0.2373 0.7028 0.5236 0.5479 0.6245
0.2620 0.3613 0.7009 0.7985 0.2161
0.3440 1.0333 0.6989 0.9457 0.2479
-0.3311 0.0588 0.3265 0.9508 0.5978
3、
(1)魔方矩阵 magic( n )
每行每列及两条线上的元素之间和都相等
如:
>> m = 100 + magic(5)
m =
117 124 101 108 115
123 105 107 114 116
104 106 113 120 122
110 112 119 121 103
111 118 125 102 109
(2)范德蒙德矩阵 vander( V )
最后一列全为1,倒数第二列为一个指定的向量,其他各列是其最后列与倒数
第二列的乘积
>> A= vander([1;2;3;4])
A =
1 1 1 1
8 4 2 1
27 9 3 1
64 16 4 1
(3)希尔伯特矩阵 hilb( n )
>> h=hilb(4)
h =
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
>> 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
(4)求·希尔伯特矩阵的逆 invhilb( n )
>> H= invhilb(4)
H =
16 -120 240 -140
-120 1200 -2700 1680
240 -2700 6480 -4200
-140 1680 -4200 2800
format short %恢复默认输出格式
(5)特普利矩阵 toeplitz( x , y )
除第一行和第一列外,其他每个元素都与作上角的元素相同
>> T = toeplitz(1:6)
T =
1 2 3 4 5 6
2 1 2 3 4 5
3 2 1 2 3 4
4 3 2 1 2 3
5 4 3 2 1 2
6 5 4 3 2 1
(6)伴随矩阵 compan( p)
>> p = [1,0,-7,6];
>> compan(p)
ans =
0 7 -6
1 0 0
0 1 0
(7) 帕斯卡矩阵 pascal( n )
由杨辉三角形组成的矩阵
>> pascal(6)
ans =
1 1 1 1 1 1
1 2 3 4 5 6
1 3 6 10 15 21
1 4 10 20 35 56
1 5 15 35 70 126
1 6 21 56 126 252