首先感谢知乎大大的回答
行列式的本质是什么? - 童哲的回答 - 知乎
https://www.zhihu.com/question/36966326/answer/70687817
当然,想这么简单就说服我,我是不干的,我要自己动手试一试知乎大大说的是否正确,于是写了一段代码来验证。
首先考虑一个简单情况
A=[2002],det(A)=4
A
=
[
2
0
0
2
]
,
d
e
t
(
A
)
=
4
现在考虑有一个矩形,它的四个点坐标是 (0,0),(0,1),(1,0),(1,1) ( 0 , 0 ) , ( 0 , 1 ) , ( 1 , 0 ) , ( 1 , 1 ) ,那么面积为1,矩形内所有点大概长这样
![这里写图片描述](https://i-blog.csdnimg.cn/blog_migrate/bf14551104db83c3102b70c1517e3729.png)
现在,将矩形内所有点做关于
A
A
的线性变换(就是左乘A啦),得到的结果长这样(红色的)
emmmm,面积变成了4,果然是倍呢!
再来考虑复杂一点的情况
B=[1324],det(B)=−2
B
=
[
1
2
3
4
]
,
d
e
t
(
B
)
=
−
2
这里 det(B)=−2 d e t ( B ) = − 2 ,是矢量,符号表示方向
还是原来的那个矩形,经过了 B B 的线性变换之后,长这样(绿色的)
![这里写图片描述](https://i-blog.csdnimg.cn/blog_migrate/d058822bc009e40b4824f7545a007c22.png)
通过一通复杂的计算,得到这个绿色的四边形的面积刚好为2,正好是 倍,真是神奇!
最后附上代码
clear all;
close all;
clc;
% create dots
[x, y] = meshgrid(0:0.01:1, 0:0.01:1);
% scatter
figure;
scatter(x(:), y(:),1.0, 'filled');
axis([0 8 0 8]);
hold on
A = [2 0; 0 2];
B = [1 2; 3 4];
p = [x(:)'; y(:)'];
A_output = A*p;
scatter(A_output(1,:), A_output(2,:),1.0, 'filled', 'MarkerEdgeColor',[0.8 0. 0.]);
B_output = B*p;
scatter(B_output(1,:), B_output(2,:),1.0, 'filled', 'MarkerEdgeColor',[0.0 0.8 0.]);