MATLAB数学公式的计算

1、矩阵kron计算

克罗内克积   C=kron（A,B）

https://baike.baidu.com/item/%E5%85%8B%E7%BD%97%E5%86%85%E5%85%8B%E7%A7%AF/6282573?fr=aladdin

2、矩阵的常用算法

https://www.cnblogs.com/zhouqing/p/3304945.html

1. 矩阵的加减乘除和(共轭)转置

(1) 矩阵的加法和减法

　　如果矩阵A和B有相同的维度(行数和列数都相等)，则可以定义它们的和A+B以及它们的差A-B，得到一个与A和B同维度的矩阵C，其中C_ij=A_ij+B_ij或A_ij-B_ij.

      另外Matlab还支持任意一个矩阵A与一个标量s相加，结果为矩阵的每一个元素加减标量，得到一个与A同维度的新的矩阵，即A+s的各个元素为A_ij+s.

(2) 矩阵的乘法

　　如果矩阵A的列数等于矩阵B的行数，则可以将A和B相乘，命令为A*B，得到一个新的矩阵C，C的行数等于A的行数，列数等于B的列数. 由于矩阵的乘法不满足交换律，所以一般A*B不等于B*A.

(3) 矩阵的张量积（tensor product）

　　矩阵A和B的张量积A⊗B可以方便地用kron函数计算，即使用命令kron(A,B), 例如

(4) 矩阵的除法

　　在MatLab中，有两个矩阵除法符号，左除\和右除/. 如果A是一个非奇异方阵(nonsingular square, 即满秩方阵)，B的行数与A的行数相等，那么A\B=A^-1B. 如果C的列数与A的列数相等，那么C/A=CA^-1.

　　从另一个角度来看，X=A\B是矩阵方程AX=B的解，X=C/A是矩阵方程XA=C的解. 如果b是一个行数与A的行数相等的列向量，则向量x=A\b是线性方程组 Ax=b的解. 且在矩阵方程AX=B中，A可以是一个m×n的矩阵，如果m=n则有唯一解；如果m<n则有多个解，Matlab会返回一个基础解；如果m>n则会返回一个最小方阵解.

(5) 矩阵的转置和共轭转置

　　在Matlab中，矩阵的共轭转置用撇号’表示，如果不需要对元素进行共轭运算，仅仅只对矩阵进行转置，则在撇号之前输入一个点号，即 .’ . 对于实数矩阵A，A’和A.’是相同的.

2. 矩阵元素操作运算

　　矩阵的运算既可以是如前所述的正常的整体运算，也可以是矩阵对应的元素依次进行标量运算，也叫数组运算，即把矩阵看做是二维数组. 对矩阵进行数组运算后得到的结果是一个与参与运算的矩阵维度相同的新矩阵，.这种元素间的算术运算的前提是参与运算的两个矩阵的维数要相同.

      对于加法和减法，元素操作运算和矩阵运算没有差别，而对于乘、除和幂运算符，相应的数组运算符是在一般的算术运算符前面加上一个点号，如

+       -       .*       ./       .\       .^

其中，A./B 是指A中的元素除以B中相应的元素，即A./B 的第i行第j列的元素(A./B)_ij=A_ij/B_ij，而(A.\B)_ij=B_ij\A_ij. 这些元素运算符的使用例子如下所示：

　　

　　

　　在Matlab中预定义的数学标准函数，如sin(x), abs(x)等都是基于对矩阵元素的运算. 如果函数f(x)是这样的一个函数，A是一个m×n的矩阵，其元素是a_ij ,那么 f(A)也是一个m×n的矩阵，其第i行第j列的元素为f(a_ij)，例如

　　

其中pi是Matlab的预定义变量，值为π，i也是预定义变量，表示复数的单位.

3. 常用的矩阵函数

　　矩阵函数是指参数为矩阵的函数，函数结果可能是一个标量值也可能是一个函数或者向量. Matlab中常用的矩阵函数包括：

　　(1) rank(A): 求矩阵A的秩，即A中线性无关的行数或者线性无关的列数.

　　

　　(2) det(A): 求矩阵A的行列式值.

　　

　　(3) inv(A): 如果A是一个非奇异(nonsingular)矩阵，则inv(A)返回A的逆矩阵. 另外还可以用左除A\eye(n)或右除eye(n)/A来计算A的逆，且在Matlab中用左除或右除来计算逆所花的计算时间比用inv函数要少，也比inv具有更好的容错性(error-detection properties).

　　(4) dot(x,y): 求同维度的向量x和y的内积/点积. 若A和B是两个具有相同维度的矩阵，则dot(A,B)是计算A和B对应列的内积，结果是一个行向量，这个行向量的列数等于A或B的列数. 例如

　　

　　(5) cross(x,y): 计算同维度的向量x和y的叉积，结果是一个向量，其方向由右手定则决定，长度等于|x|*|y|sin<x,y>. 若A和B是两个具有相同维度的矩阵，则cross(A,B)是计算A和B对应列的叉积，结果是一个维度与A和B相等的矩阵.

　　(6) kron(A,B): 得到矩阵A和B的张量积.

　　(7) isequal(A,B): 如果矩阵A和B是相同的，即具有相同的维数和相同的内容，则返回1.

　　(8) isreal(A): 判断A是否是一个实矩阵，如果是则返回1，否则返回0.

　　(9) trace(A): 计算方阵A的迹，即对角线元素之和.

　　(10) eig(A): 计算方阵A的特征值，结果是一个列向量，向量中元素的个数等于特征值的个数，即A的维度(A的行数或列数).

　　(11) [U,D]=eig(A): 计算方阵A的特征值和特征向量，得到两个方阵U和D，其中D的对角线元素为A的特征值，U的列向量为A的归一化的(normalized)特征向量，例如

　　

3、求矩阵的范数

https://jingyan.baidu.com/article/a378c960c07500b3282830a0.html

http://blog.youkuaiyun.com/yihaizhiyan/article/details/6904599

nm1=norm(a,1) ，其中norm就是求矩阵范数的函数，1表示的是1范数

nm2=norm(a,2) ，其中norm就是求矩阵范数的函数，2表示的是2范数

nm3=norm(a,inf) ，其中norm就是求矩阵范数的函数，inf表示的是无穷范数

nm4=norm(a,'for') ，其中norm就是求矩阵范数的函数，for表示的是Frobenius范数

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范数 ，矩阵元p阶范数估计需要自己编程求，

计算公式如下

 举个例子吧a=magic(3)sum(sum(abs(a)^4))^(1/4)a = 8 1 6 3 5 7 4 9 2

ans = 19.7411

转自：http://zhidao.baidu.com/question/245196523.html

4、Matlab如何得出矩阵的前k个最大特征值所对应的特征向量？

http://www.ilovematlab.cn/thread-50763-1-1.html

A=[1 3 7;3 8 9;5 4 6];
[V,D]=eig(A);
D=diag(D);
[ D,I]=sort(D,'descend');
[D(1:2) V(:,I(1:2))']

用sort排序，decend是降序，取前面是特征值比较大的；默认是升序，前面是特征值比较小的

5、matlab迭代求解

最基本的循环

1. a=;
   
2. b=;
   
3. d=;
   
4. c=1;
   
5. n=0;
   
6. errf=1;
   
7. while errf>1e-8&n<100%计算结果精度要求1e-8，如果迭代超过100次还不收敛，退出循环。
   
8. y=a*b*c;
   
9. c1=y*d;
   
10. errf=(c1-c)/c1;
    
11. c=c1;
    
12. n=n+1;
    
13. end