MATLAB学习笔记--- 特殊矩阵(2)

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

猜你喜欢

转载自blog.csdn.net/naturly/article/details/82780349