雅可比行列式【2】Jacobian行列式的意义
2.1 线性变换将面积伸缩
对于一个\(\R^2\to\R^2\)的线性变换:
\[
T(x,y)=
\left[
\begin{array}{c}
4x-2y\\
2x+3y
\end{array}
\right]
\]
设区域\(S_1=\{(x,y)|0\leq x,y\leq1\}\),若想要求\(\iint_{S_1}T(x,y)\ d\sigma\).可以通过基底表示单位正方形:\(e_1=(1,0)',e_2=(0,1)'\),则:
\[
S_1=\{xe_1+ye_2|0\leq x,y\leq1\}
\]
设\(A\)为线性变换\(T\)参考标准基地的表示矩阵,即有:
\[
T(xe_1+ye_2)=A(xe_1+ye_2)=xAe_1+yAe_2=xa_1+ya_2
\]
于是:
\[
T(S_1)=\{xa_1+ya_2|0\leq x,y\leq1\}
\]
这表明\(T(S_1)\)是以\(A=(a_1,a_2)\)表示的平行四边形,二阶行列式的绝对值为平行四边形的面积,因此\(v(T(S_1))=|detA|\)。这个结果表明平行四边形\(S_1\)经过线性变换\(T\),面积伸缩了\(|detA|\)倍。
2.2 Jacobian行列式的意义
if \(F:\R^n\to\R^n\) is derivable, then the Jacobian matrix is in \(n\times n\) form in which we could express a number of it. We set the n is equal to 2, and vector function is: \(F:u\to x\)
\[
det\ J(u,v)=
\left|
\begin{matrix}
\frac{\partial x}{\partial u}&\frac{\partial x}{\partial v}\\
\frac{\partial y}{\partial u}&\frac{\partial y}{\partial v}\\
\end{matrix}
\right|=\frac{\partial x}{\partial u}\frac{\partial y}{\partial v}-\frac{\partial x}{\partial v}\frac{\partial y}{\partial v}
\]
若令\(R=\{r_1,r_2\}\),其中\(r_1=(du,0)',r_2=(0,dv)'\)表示长方形,则\(F(R)=\{F(u)|u\in R\}\)近似如下面向量所表示的平行四边形:
\[
J(u,v)(du,0)'=
\left|
\begin{matrix}
\frac{\partial x}{\partial u}&\frac{\partial x}{\partial v}\\
\frac{\partial y}{\partial u}&\frac{\partial y}{\partial v}\\
\end{matrix}
\right|(du,0)'=
\left|
\begin{matrix}
\frac{\partial x}{\partial u}du\\
\frac{\partial y}{\partial u}du
\end{matrix}
\right|\\
J(u,v)(0,dv)'=\left|\begin{matrix}\frac{\partial x}{\partial u}&\frac{\partial x}{\partial v}\\\frac{\partial y}{\partial u}&\frac{\partial y}{\partial v}\\\end{matrix}\right|(0,dv)'=\left|\begin{matrix}\frac{\partial x}{\partial v}dv\\\frac{\partial y}{\partial v}dv\end{matrix}\right|
\]
若令\(dA\)表示平行四边形\(F(R)\)的面积, 因为二阶行列式的行向量所形成的平行四边形面积等于行列式的绝对值,则:
\[
dA=\left|
det
\left[
\begin{matrix}
\frac{\partial x}{\partial u}du&\frac{\partial x}{\partial v}dv\\
\frac{\partial y}{\partial u}du&\frac{\partial y}{\partial v}dv\\
\end{matrix}
\right]
\right|=
\left|
det
\left[
\begin{matrix}
\frac{\partial x}{\partial u}&\frac{\partial x}{\partial v}\\
\frac{\partial y}{\partial u}&\frac{\partial y}{\partial v}\\
\end{matrix}
\right]
\right|dudv=|det\ J(u,v)|dudv
\]
所以微笑区域\(R\)经过向量函数\(F:R\to F(R)\),其面积伸缩了\(|det\ J(u,v)|\)倍。对于\(f:\R^2\to\R\)我们可以得出变换积分公式:
\[
\int_{F(R)} f(x,y)dxdy=\int_{R} f(x(u,v),y(u,v))\left|J(u,v)\right|dudv\\
|J(u,v)|=\left|\frac{\partial(x,y)}{\partial(u,v)}
\right|
\]
雅可比行列式【2】Jacobian行列式的意义的更多相关文章
- Jacobian矩阵和Hessian矩阵
1.Jacobian矩阵 在矩阵论中,Jacobian矩阵是一阶偏导矩阵,其行列式称为Jacobian行列式.假设 函数 $f:R^n \to R^m$, 输入是向量 $x \in R^n$ ,输出为 ...
- 《Linear Algebra and Its Applications》-chaper3-行列式-行列式初等变换
承接上一篇文章对行列式的引入,这篇文章将进一步记录关于行列式的有关内容,包括如下的几个方面: (1)行列式3个初等变换的证明. (2)转置行列式与原行列式相等的证明. (3)定理det(AB) = d ...
- 【原创】开源Math.NET基础数学类库使用(15)C#计算矩阵行列式
本博客所有文章分类的总目录:[总目录]本博客博文总目录-实时更新 开源Math.NET基础数学类库使用总目录:[目录]开源Math.NET基础数学类库使用总目录 上个月 ...
- 开源Math.NET基础数学类库使用(15)C#计算矩阵行列式
原文:[原创]开源Math.NET基础数学类库使用(15)C#计算矩阵行列式 本博客所有文章分类的总目录:http://www.cnblogs.com/asxinyu/p ...
- 矩阵&行列式
# 代数 排列 对换,对于一个排列操作,对于一个偶排列一次对换之后变为奇排列 反之变为偶排列 行列式 N阶行列式室友N^2个数aij(i,j = 1,2,3,...n) 行列式的数=\(\sum_ { ...
- Eigen求矩阵行列式 及 行列式本质
转置.伴随.行列式.逆矩阵 小矩阵(4 * 4及以下)eigen会自动优化,默认采用LU分解,效率不高 #include <iostream> #include <Eigen/Den ...
- C++中计算矩阵的行列式
使用eigen库: 求行列式: #include <iostream> #include <Eigen/Dense> using namespace std; using na ...
- 方阵行列式并行化计算(OpenMP,MPI),并计算加速比
00][100].在创建方阵时,方阵的阶数N(N<100)由外部输入.然后用两层"for循环"来给方阵 p左上角 N×N个位置赋值.具体实现如下: /* * 定义矩阵阶数N ...
- 洛谷P7112 行列式求值
行列式求值 这是一个让你掉头发的模板题 行列式的定义 行列式 (\(\texttt{Determinant}\)) 是一个函数定义,取值是一个标量. 对一个 \(n\times n\) 的矩阵 \(A ...
随机推荐
- <深入理解redis>读书笔记
chapter2 键管理与数据结构 对大多数redis解决方案而言,键的命名设计至关重要.对于管理来说,内存消耗和redis性能都与数据结构设计相关.所以对开发者而言,最好有数据结构的命名文档规范. ...
- Day 26:Dom4j修改xml
Dom4j:Writing a document to a file import java.io.File; import java.io.FileOutputStream; import org. ...
- Gym - 101158C Distribution Center
题意:n个传送带,传送带i运送编号为i的物品,机器人可以负责把传送带i上的物品放到传送带i + 1上,也可以把传送带i + 1上的物品放到传送带i上,机器人分布在传送带上x轴的不同位置,问每个传送带最 ...
- Dockerfile详解,以及构建自定义镜像
Dockerfile使用 前面的操作我们一直下载下载官方已经构建好的镜像,直接下载后就可以run,如果我们想要在镜像中添加自己的应用,比如在tomcat中添加自己的app,构建一个自定义的镜像,那么我 ...
- 028-PHP常用数学函数abs和acos和asin
<?php print(abs(-));//打印绝对值 // 从 -1 到1打印acos函数的值 print("<TABLE BORDER=\"1\"> ...
- 实验吧-隐写术-刷新 刷新 快刷新(f5-steganography)
题目:刷新 刷新 快刷新 其实这就算是很好的提示了,刷新的时候就是F5嘛,这个题就是F5隐写,用f5-steganography来解题. 题中的图片: f5-steganography下载链接:htt ...
- JavaScript获取当前日期前后7天的日期
function fun_date(aa){ var date1 = new Date(), time1=date1.getFullYear()+"-"+(date1.getMon ...
- C# Process类详解
C# Process类详解 Process[] processes = Process.GetProcessesByName(current.ProcessName); 根据进程名字找到所有进程,返回 ...
- 2016蓝桥杯省赛C/C++A组第七题 剪邮票(暴力+并查集)
题意:有12张连在一起的12生肖的邮票.现在你要从中剪下5张来,要求必须是连着的.(仅仅连接一个角不算相连) 分析:暴力+并查集. 1.记录下每个数字所在位置. 2.先枚举各不相同的5个数的所有可能情 ...
- AVCodec 结构体
typedef struct AVCodec { // 标示Codec 的名字, 比如,"h264" "h263" 等. const char *name; / ...