视频

拓展学习

例1  验证det(A-1)=1/det(A)。

format rat

A=[1,3,2;-3,2,1;4,1,2]

det(inv(A))

1/det(A)

例2  求3~20阶魔方阵的秩。

for n=3:20

    r(n)=rank(magic(n));

end

bar(r)

grid on

axis([2,21,0,20])

[3:20;r(3:20)]

例3  求2~10阶希尔伯特矩阵的条件数。

for n=2:10

    c(n)=cond(hilb(n));

end

format long

c'