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'
