task逆矩阵

逆矩阵
对于n阶方阵A，如果存在AB=BA=I（单位矩阵），称B为A的逆矩阵。逆矩阵存在的充要条件是A的行列式不等于0.
方阵为A，则逆矩阵的函数inv(A).
例1

>> A=[1 2 3 4;2 3 5 8;1 5 3 7;3 4 7 11]
> A =     
 1     2     3     4    
 2     3     5     8   
 1     5     3     7    
 3     4     7    11 
 >> det(A) ans =    -3.0000 
>> inv(A) ans =   
> -0.3333   -7.0000    0.3333    5.0000   
>  0.6667   -2.0000    0.3333    1.0000   
>   1.3333    1.0000   -0.3333   -1.0000  
>    -1.0000    2.0000    0.0000   -1.0000 
>    >> inv(A)*A ans =     
>    1.0000   -0.0000    0.0000   -0.0000 
>      -0.0000    1.0000    -0.0000    0.0000  
>        0.0000    1.0000         0         0      
>        0   0         0    1.0000 
>> A\eye(4)
> ans =  
>  -0.3333   -7.0000    0.3333    5.0000   
>   0.6667   -2.0000    0.3333    1.0000   
>    1.3333    1.0000   -0.3333   -1.0000  
>     -1.0000    2.0000    0.0000   -1.0000
> A*(-1) 
> ans =     -1    -2    -3    -4  
            -2    -3    -5    -8   
            -1    -5    -3    -7    
            -3    -4    -7   -11 

例二2.35矩阵的LU分解
将矩阵A分解为l*U,u为上三角矩阵，l为上三角矩阵
LU常用与求解行列式及线性方程组

矩阵的正交分解
qr分解。q为正交矩阵（该矩阵和它的转置矩阵的的乘积为单位矩阵），r为上三角矩阵
广义逆矩阵
APA=A，PAP=P,则P为广义逆矩阵

S数组与矩阵的幂乘
V=x.^y
j矩阵的水平连接和垂直连接
水平连接：horzcat(A,B)
垂直连接：vertcat(A,B)

J矩阵的复制repmat(A,m,n)
x稀疏矩阵的建立
sparse(A);非零元素直接输入法：S=sparse(i,j,s,m,n);对角线方式A=spdiags(B,d,m,n)
寻找矩阵的非零元素：k=find(X);[i,j]=find(X);[i,j,v]=find(X)