MATLAB矩阵运算

常用的函数：

eye(size(A)) 产生与A矩阵同阶的单位矩阵

zeros(m,n) 产生0矩阵

ones(m,n) 产生幺矩阵

rand(m,n) 产生随机元素的矩阵

randn(m,n) 产生均值为0，方差为1的标准正态随机矩阵

size(A) 返回包含两个元素的向量

length(x) 返回向量的长度

diag(X) 1、若X是矩阵，则diag(X)是X的主对角线向量

             2、若X是向量，则diag(X)产生以X为主对角线的对角矩阵

tril(A) 提取一个矩阵的下三角部分

triu(A) 提取一个矩阵的上三角部分

magic、hilb、pascal

矩阵左除"\"要求矩阵行数相同，矩阵右除"/"要求矩阵列数相同

inv(A) 求逆矩阵（方阵）

rank(A) 矩阵的秩

eig(A) 矩阵的特征值与特征向量

rref(A) 矩阵的最简形

det(A) 矩阵的行列式

 

A=[2 3 5;3 4 1;0 9 7]，是一个3*3矩阵，想插入一行r=[1 2 3]时，可以这样做：A=[A;r]这样A就变成了4*3的矩阵

A如下：A=[2 3 5;3 4 1;0 9 7;1 2 3]

假如已有矩阵A为：A=[2 3 5;3 4 1;0 9 7;1 2 3]，想插入一列c=[1;2;3;4]时，可以这样做：A=[A c]

x(i):向量x中的第i个元素

A(i,j):矩阵A中的第i行，第j列元素

按特征值从大到小排序

A = magic(5);
[V,D] = eig(A);
[d,ind] = sort(diag(D));
DS = D(ind,ind);
VS = V(:,ind)

默认情况下，eig 并不总是返回已排序的特征值和特征向量。可以使用 sort 函数将特征值按升序排序，并重新排序相应的特征向量。

使用 diag(D) 从 D 的对角线上提取特征值，然后按升序对得到的向量进行排序。sort 的第二个输出返回索引的置换向量。

使用 ind 对 D 的对角线元素进行重新排序。由于 D 中的特征值对应于 V 的各列中的特征向量，因此还必须使用相同的索引对V的列进行重新排序。