matlab矩阵操作大全

命令 产生线性等分向量 函数 linspace

格式 y = linspace(a,b) %在(a, b)上产生100个线性等分点 y = linspace(a,b,n) %在(a, b)上产生n个线性等分点 命令 产生对数等分向量 函数 logspace

格式 y = logspace(a,b) %在（ ）之间产生50个对数等分向量 y = logspace(a,b,n) y = logspace(a,pi) 命令 计算矩阵中元素个数

n = numel(a) %返回矩阵A的元素的个数 命令 产生以输入元素为对角线元素的矩阵 函数 blkdiag

格式 out = blkdiag(a,b,c,d,…) %产生以a,b,c,d,…为对角线元素的矩阵 例1-13

>> out = blkdiag(1,2,3,4) out =

1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 命令 友矩阵 函数 compan

格式 A = compan(u) %u为多项式系统向量，A为友矩阵，A的第1行元素为 -u (2:n)/u(1)，其中u (2:n)为u的第2到第n个元素，A

为特征值就是多项式的特征根。 例1-14 求多项式的友矩阵和根 >> u=[1 0 -7 6];

>> A=\%求多项式的友矩阵 A =

0 7 -6 1 0 0 0 1 0

>> eig(A) %A的特征值就是多项式的根 ans = -3.0000 2.0000 1.0000

命令 hadamard矩阵 函数 hadamard

格式 H = hadamard(n) %返回n阶hadamard矩阵 例1-15

>> h=\h =

1 1 1 1 1 -1 1 -1 1 1 -1 -1 1 -1 -1 1 命令 Hankel方阵 函数 hankel

格式 H = hankel(c) %第1列元素为c，反三角以下元素为0。 H = hankel(c,r) %第1列元素为c，最后一行元素为r，如果c的最后一个元素与r的第一个元素不同，交叉位置元素取为c的最后一个元素。 例1-16

>> c=\

c =

1 2 3 r =

7 8 9 10 >> h=\h =

1 2 3 8 2 3 8 9 3 8 9 10 命令 Hilbert矩阵 函数 hilb

格式 H = hilb(n) %返回n阶Hilbert矩阵，其元素为H(i,j)=1/(i+j-1)。 例1-17 产生一个3阶Hilbert矩阵 >> format rat %以有理形式输出 >> H=\H =

1 1/2 1/3 1/2 1/3 1/4 1/3 1/4 1/5 命令 逆Hilbert矩阵 函数 invhilb

格式 H = invhilb(n) %产生n阶逆Hilbert矩阵 命令 Magic(魔方)矩阵 函数 magic

格式 M = magic(n) %产生n 阶魔方矩阵 例1-18

>> M=\M =

8 1 6 3 5 7 4 9 2 命令 Pascal矩阵

函数 pascal

格式 A = pascal(n) %产生n阶Pascal矩阵，它是对称、正定矩阵，它的元素由Pascal三角组成，它的逆矩阵的所有元素都是整数。 A = pascal(n,1) %返回由下三角的Cholesky系数组成的Pascal矩阵

A = pascal(n,2) %返回Pascal(n,1)的转置和交换的形式 例1-19

>> A=\A =

1 1 1 1 1 2 3 4 1 3 6 10 1 4 10 20 >> A=\A =

1 0 0 1 -1 0 1 -2 1 >> A=\A =

1 1 1 -2 -1 0 1 0 0 命令 托普利兹矩阵 函数 toeplitz

格式 T = toeplitz(c,r) %生成一个非对称的托普利兹矩阵，将c作为第1列，将r作为第1 行，其余元素与左上角相邻元素相等。 T = toeplitz(r) %用向量r生成一个对称的托普利兹矩阵 例1-20

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

>> r=[1.5 2.5 3.5 4.5 5.5];

>> T=\T =

1 5/2 7/2 9/2 11/2 2 1 5/2 7/2 9/2 3 2 1 5/2 7/2 4 3 2 1 5/2 5 4 3 2 1 命令 Wilkinson特征值测试阵 函数 wilkinson

格式 W = wilkinson(n) 例1-21

>> W=\W =

3/2 1 0 1 1/2 1 0 1 1/2 0 0 1 >> W=\W =

3 1 0 0 0 1 2 1 0 0 0 1 1 1 0 0 0 1 0 1 0 0 0 1 1 0 0 0 0 1 0 0 0 0 0

%返回n阶Wilkinson特征值测试阵0 0 1 3/2 0 0 0 0 0 0 0 0 1 0 2 1 1 3