遗忘海岸

江湖程序员 -Feiph(LM战士)

导航

< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5

统计

关于雅可比行列式与积分换元

换元前后微元数目相同,然后我们保证每个微元的积分(就是dxdy * f(x,y) 的简单乘积)相同那么最后的结果必定是一样的。

对于二元情况的证明参考同济高数7版 P151

A

考虑线性方程组

u=ax+by

v=cx+dy

------------------------------

如果在xy平面上取 (0,0),(1,0),(0,1),(1,1)4个点构成一个变长为1的正方形,那么经过

[a   b

c    d] 做变换后会是一个平行四边形。在uv平面上是 <a,b>,<c,d> 两个向量

向量的面积  | <a,b>  x  <c,d> |  = ad-cd  这就表示变换后的面积比原面积是ad-cb/1

 

等于方程组的对应得行列式

 B

x=g(u,v) y=h(u,v) ,  x,y 与 u v不是线性的

但是做全微分后,   dx= Gu du + Gv dv ,  dy=Hu du  +Hv dv

可见微元 dxdy 与 dudv 在指定点(u0,v0)  是成线性关系的。 dxdy 、dudv 面积之比

| Gu   Gv

  Hu    Hv|  即雅可比行列式(行列式不能是0)  即 dxdy/dudv=J   所以做积分变换时 dxdy=  J * dudv

 

考虑  f(x,y) dxdy 积分变换后要保证值一致(微元数一样),由于被积函数一样都是f(x,y) =f( g(u,v),h(u,v)) ,所以当 f(x,y) dxdy =  f(g(u,v),h(u,v)) * J * dudv时才能保证一致。

注意 J^-1  是  u=g(x,y)  ,v=(x,y)    即反函数表示雅克比行列式

他们互为倒数

---------------------------------------------------------

C

已知道u=g(x,y) v=h(x,y)时 J=| (u,v)/(x,y)|

dudv= J dxdy

已知道联合密度函数在区域上的积分是1

那么变换后每个微元 dudv都有一个正好对应的 dxdy,变换前后的密度函数在该区域上的概率取值应该相等即 

fuv(u,v) du dv= fxy(x,y) dxdy =>    fuv(u,v)=   fxy(x,y)  ( dxdy/dudv)=  fxy(x,y) J ^-1

 

posted on   遗忘海岸  阅读(6068)  评论(0编辑  收藏  举报

编辑推荐:
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
阅读排行:
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· .NET周刊【3月第1期 2025-03-02】
· 分享 3 个 .NET 开源的文件压缩处理库,助力快速实现文件压缩解压功能!
· [AI/GPT/综述] AI Agent的设计模式综述
点击右上角即可分享
微信分享提示