MATLAB 矩阵求值 例题

MATLAB 矩阵求值 例题

eg:验证det(A的逆）=1/det(A)。

format rat
A=[1,3,2;-3,2,1;4,1,2]
det(inv(A))
1/det(A)

矩阵的秩：

eg:求3-20阶魔方阵的在秩

for n=3:20
    r(n)=rank(magic(n));
end
bar(r)
grid on
axis([2,21,0,20])

矩阵的迹：

A=[1,3,2;-3,2,1;4,1,2]
b=trace(A)
t=sum(diag(A))

向量和矩阵的范数：

%norm(V,1):计算V的1–范数。
%norm(v)或norm(V,2):计算V的2–范数。
%norm(V,inf):计算A的无穷–范数。

x=[2 0 1;-1 1 0;-3 3 0]
n=norm(x)
n=norm(x,1)

矩阵的条件数：

%cond(A,1):计算A的1–范数下的条件数。
%cond(A)或cond(A,2):计算A的2–范数下的条件数。
%cond(A,inf):计算A的无穷–范数下的条件数。

eg:求2-10阶希尔伯特矩阵的条件数。

for n=2:10
    c(n)=cond(hilb(n));
end
format long
c'