下面介紹矩陣的一些基本操作,包括矩陣的特征值,三角陣,對角陣,矩陣的翻轉等,以及矩陣的一些特性,例如矩陣的秩,矩陣的跡.最後介紹了矩陣的超越函數.

1　　方陣的行列式

1 clear all;
2 A=magic(3)
3 B=[1:3;2 5 7;3 8 7]
4 y1=det(A)
5 y2=det(B)

2　　特征值.特征向量和特征多項式

1 clear all;
2 A=magic(3)
3 %E是特征值
4 E=eig(A)
5 %V的每一列都是特征向量,D的對角線上是特征值
6 [V,D]=eig(A)

 1 clear all;
 2 %p是方程的系數
 3 p=[3
 
 5 2 1]
 4 %求伴隨矩陣
 5 A=compan(p)
 6 %求特征值,特征值就是根
 7 x1=eig(A)
 8 
 9 %用roots的方式求根
10 x2=roots(p)

3　　對角陣

1 clear all;
2 A=rand(3,4)
3 %提取矩陣A的主對角線元素,產生一個列向量
4 b1=diag(A)
5 %0是主對角線,1是往上一個的對角線
6 b2=diag(A,1)
7 %提取第2條對角線的元素組成一個列向量
8 b3=diag(A,2)

4　　上三角陣與下三角陣

1 clear all;
2 A=rand(3,3)
3 %上三角矩陣
4 B1=triu(A)
 
5 %下三角矩陣
6 B2=tril(A)
7 %右上有一個0的三角形
8 B3=tril(A,1)

5　　矩陣的逆和偽逆(廣義逆矩陣)

 1 clear all;
 2 A=magic(3)
 3 %兩行兩列
 4 B=[1 3;2 6]
 5 %求矩陣A的逆矩陣
 6 C=inv(A)
 7 %相乘為單位矩陣
 8 C*A
 9 %求逆矩陣(結果不存在)
10 inv(B)
11 %求廣義逆矩陣
12 D=pinv(B)
13 %相乘結果為E
14 B*D*B

6　　矩陣的秩

1 clear all;
2 A=magic(3)
3 B=[1 2 4;2 3 5;2 4 8]
4 %求秩
 
5 r1=rank(A)
6 r2=rank(B)

7　　矩陣的跡

1 clear all;
2 A=magic(3)
3 %求跡
4 t1=trace(A)
5 %矩陣的特征值
6 eig(A)
7 %矩陣特征值的和
8 t2=sum(eig(A))

8　　矩陣的標準正交基

1 clear all;
2 A=[1 2 3;3 5 7;9 5 8]
3 B=magic(3)
4 %求A的標準正交基
5 C=orth(A)
6 D=orth(B)
7 %正交基的轉置乘以正交基是單位向量
8 C‘*C

9　　LU分解

 1 clear all;
 2 A=[ 2 3 4;8 4 9;5 3 1]
 3 %lu分解
 4 [L1,U1]=lu(A)
 5 %驗證
 6 L1*U1
 7 
 8 %L2下三角矩陣 U2上三角矩陣 P置換矩陣
 9 [L2,U2,P]=lu(A)
10 %下三角矩陣和上三角矩陣合並在矩陣Y中
11 Y1=lu(A)
12 Y2=L2+U2-eye(size(A))

10　　QR分解

1 clear all;
2 A=[ 2 3 4;8 4 9;5 3 1]
3 %QR分解Q正交陣,R1上三角矩陣
4 [Q1,R1]=qr(A)
5 
6 B=[1 2 3 4;3 5 6 2;3 6 9 12]
7 %正交分解
8 [Q2,R2]=qr(B)
9 Q2*R2

6.矩陣的特征和線性代數