matlab基础之矩阵备忘录

%%对角阵

diag(A)%取A的主对角线（第0条）元素，产生列向量
diag(A,K)%取A的第K条对角线元素，产生列向量
diag(V)%以向量V为主对角线元素构造对角矩阵
diag(V,K)%以向量V为第K条对角线元素构造对角矩阵

%%三角阵

%上三角
triu(A)%提取A主对角线及以上元素
triu(A,K)%提取A的K对角线及以上元素
% 下三角（同上）
tril(A)
tril(A,K)

%%矩阵操作

A.’%A的转置矩阵
A’%A的转置且对每个数取复共轭
rot90(A,K)%将矩阵A逆时针旋转90°的K倍
fliplr(A)%对矩阵A左右翻转
flipud(A)%对矩阵A上下翻转
inv(A)%对A求逆

%%通用矩阵

%参数为（m,n)或（m）或（sizeof（A））
zeros()%产生全0矩阵
ones()%产生全1矩阵
eyes()%产生对角线为1的矩阵
rand()%产生（0，1）随机分布的矩阵（不包括0和1）
randn()%产生均值为0方差为1的标准正态分布随机矩阵

%%特殊矩阵(例n=4)

A=magic(n)%n阶幻方矩阵

16     2     3    13
 5    11    10     8
 9     7     6    12
 4    14    15     1

A=vander(1:5)%范德蒙矩阵

 1     1     1     1     1
16     8     4     2     1
81    27     9     3     1
256   64    16     4     1
625  125    25     5     1

A=hilb(n)%n阶希尔伯特矩阵，注意提前format rat有理数形式输出

   1              1/2            1/3            1/4     
   1/2            1/3            1/4            1/5     
   1/3            1/4            1/5            1/6     
   1/4            1/5            1/6            1/7 

A=compan( P )%P为多项式系数向量（高次幂在前），生成多项式的伴随矩阵
%x^3 - 2x^2 - 5x + 6 = 0 P=[1 -2 -5 6]

   2              5             -6       
   1              0              0       
   0              1              0   

A=pascal(n)%n阶帕斯卡矩阵

   1              1/2            1/3            1/4     
   1/2            1/3            1/4            1/5     
   1/3            1/4            1/5            1/6     
   1/4            1/5            1/6            1/7