Matal线性代数.doc

Matal线性代数 一、常见运算 1、d = det(X) 得到A的行列式 2、k = rank(A) 得到A的秩 3、函数 size 格式 d = size(A) 若A为m*n阶的符号矩阵，则输出结果d=[m，n]。 [m,n] = size(A) %分别返回矩阵A的行数于m，列数于n。 d= size(A, n) %返回由标量n指定的A的方向的维数：n=1为行方向，n=2为列方向。 4、[V,D] = eig(A) %计算A的特征值对角阵D和特征向量V，使AV=VD成立。 5、Y = inv(X) 得到X的逆矩阵 6、范数计算 norm(A,p) Returns sum(abs(A).^p)^(1/p), for any 1 = p = ∞. norm(A) Returns norm(A,2). norm(A,inf) Returns max(abs(A)). norm(A,-inf) Returns min(abs(A)). 一般来说：给定向量X，其范数 norm(X)表示向量X的长度 7、orth(A) 对A的列向量组进行单位正交化 8、rref(A) 对A进行初等行变换，得行最简矩阵 9、b = trace(A) 求A的迹 10、ZR = null(A,r) 求A的基础解系 10、代数运算 A+B A-B A*B A.*B A/B A./B A\B A.\B A^B A.^B A A. 11、生成对角阵 diag([1 2 3]) ans = 1 0 0 0 2 0 0 0 3 12、如何用整数商表示有理数，rat() A=[0.2 0.3 0.25];[a,b]=rat(A) a = 1 3 1 b = 5 10 4 13、如何求余数？rem(12,5)，ans =2 3、Matlab有没有求矩阵行数/列数/维数的函数？ndims(A)返回A的维数size(A)返回A各个维的最大元素个数length(A)返回max(size(A))[m,n]=size(A)如果A是二维数组，返回行数和列数nnz(A)返回A中非0元素的个数 norm 格式 n = norm(X) %X为向量，求欧几里德范数，即。 n = norm(X,inf) %求-范数，即。 n = norm(X,1) %求1-范数，即。 n = norm(X,-inf) %求向量X的元素的绝对值的最小值，即。 n = norm(X, p) %求p-范数，即，所以norm(X,2) = norm(X)。 命令 矩阵的范数 函数 norm 格式 n = norm(A) %A为矩阵，求欧几里德范数，等于A的最大奇异值。 n = norm(A,1) %求A的列范数，等于A的列向量的1-范数的最大值。 n = norm(A,2) %求A的欧几里德范数，和norm(A)相同。 n = norm(A,inf) %求行范数，等于A的行向量的1-范数的最大值 即：max(sum(abs(A)))。 n = norm(A, fro ) %求矩阵A的Frobenius范数， 即sqrt(sum(diag(A*A)))，不能用矩阵p-范数的定义来求。 命令 范数的估计值 函数 normest 格式 nrm = normest(A) %矩阵A的2-范数(欧几里德范数)的估计值，相对误差小于106。 nrm = normest(A,tol) %tol为指定相对误差 [nrm,count] = normest(…) %count给出计算估计值的迭代次数 三、矩阵的变向 （1）矩阵旋转 函数 格式 B = rot90 (A) %将矩阵A逆时针方向旋转90° B = rot90 (A,k) %将矩阵A逆时针方向旋转(k×90°)，k可取正负整数。 例1-50 A=[1 2 3;4 5 6;7 8 9] A = 1 2 3 4 5 6 7 8 9 Y1=rot90(A),Y2=rot90(A,-1) Y1 = %逆时针方向旋转 3 6 9 2 5 8 1 4 7 Y2 = %顺时针方向旋转