利用快速松弛反演算法(RRI）快速计算Jacobian矩阵

最新推荐文章于 2024-06-12 21:09:41 发布

odesseyy

最新推荐文章于 2024-06-12 21:09:41 发布

阅读量327

点赞数 1

文章标签：矩阵线性代数

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/abner361/article/details/127882389

版权

mark一下个人理解：

在电磁法勘探中，由于通常计算偏导数矩阵是通过差分的方法对一个频率：

$\frac{\partial F_j}{\partial m_i} = \frac{F_j[(m_{i(n)}+\Delta m_i),(m_{1,n},m_{2,n},\dots ,m_{i-1,n},m_{i+1,n},\dots,m_{N,n}]-F_j(m_n)}{\Delta m_i}$

即是对模型参数 $m_i$ 给顶一个小的增量，其他模型参数固定不变，求取响应的差值和模型增量的比值，即 $\frac{\partial F_j}{\partial m_i}$ 。

然而对于需要多次求解Jacobian矩阵的反演问题，由于需要多次正演，耗时巨大，效率很差。RRI求偏导数的原理：

$\rho_E = \frac{i}{\omega \mu }{(\frac{E_x}{H_y})^2|_{z=0}}$

$V = I\omega \mu \frac{H_y}{E_x}$

=> $\rho_E = \frac{-i\omega\mu}{V^2}|_{z=0}$

对上式进行扰动：

$\delta \rho_E = i\omega\mu \frac{2E^3_0}{(i\omega\mu )^3H^3_0}\delta V$

又因为：

$\delta V(0) = \frac{i\omega\mu}{E^2_0(0)} \int E^2_0(z) \delta \sigma_i(z)dz$

$\frac{\partial \rho_E}{\partial \rho_i}=\frac{\partial \rho_E}{\partial \sigma_i}(\frac{-1}{\rho^2_i})$

最终得到视电阻率对电阻率的偏导数：

$\frac{\partial \rho_E}{\partial \rho_i}=\frac{-2E_0}{i\omega\mu H^3_0 \rho^2_i} \int E^2_0(z)dz$

通过一次正演后，即可知道 $E_0$ 和 $H_0$ ,即：只需要通过一次正演，即可求得偏导数。接下来就是代码编写。

continue...

参考文献：J Torquil Smith,John R Nooker.1991.Rapid Inversion of Two-and Three-Dimensional Magnetotelluric Data[J].Journal of Geophysical Research,96(b3):3905-3922.

关注

1
点赞
踩
1

收藏

觉得还不错? 一键收藏
1
评论
利用快速松弛反演算法(RRI）快速计算Jacobian矩阵

然而对于需要多次求解Jacobian矩阵的反演问题，由于需要多次正演，耗时巨大，效率很差。给顶一个小的增量，其他模型参数固定不变，求取响应的差值和模型增量的比值，即。,即：只需要通过一次正演，即可求得偏导数。接下来就是代码编写。通过一次正演后，即可知道。
复制链接

扫一扫

odesseyy CSDN认证博客专家 CSDN认证企业博客

码龄2年

5: 原创

166万+: 周排名

189万+: 总排名

5842: 访问

: 等级

57: 积分

1: 粉丝

5: 获赞

4: 评论

11: 收藏

私信

关注

热门文章

最新评论

wps中的Latex公式与编号之间的空白消除方法
狂发文章的生煎: 博主，可以求个完整的代码不
稳定双共轭梯度法BiCGSTAB流程
CSDN-Ada助手: 非常感谢博主写的关于稳定双共轭梯度法BiCGSTAB的流程介绍，让我们更好地了解这个算法。我觉得博主可以考虑继续写关于优化算法的文章，比如介绍如何使用BiCGSTAB来解决大规模线性方程组等问题，这样的技术文章对其他用户也会很有帮助。相信下一篇文章会吸引更多读者的关注。期待你的新作品！为了方便博主创作，提高生产力，CSDN上线了AI写作助手功能，就在创作编辑器右侧哦～（https://mp.csdn.net/edit?utm_source=blog_comment_recall ）诚邀您来加入测评，到此（https://activity.csdn.net/creatActivity?id=10450&utm_source=blog_comment_recall）发布测评文章即可获得「话题勋章」，同时还有机会拿定制奖牌。
wps中的Latex公式与编号之间的空白消除方法
Dandelionldc: 想问一下如何打出公式右侧的序号？
利用快速松弛反演算法(RRI）快速计算Jacobian矩阵
CSDN-Ada助手: 你的文章很有意思，不过质量分有点低, 仅为: 12分, 我觉得可以考虑这些改进，说不定有很好的效果。具体如下：(1)使用更多的站内链接；(2)增加除了各种控件外，文章正文的字数；(3)增加条理清晰的目录。

您愿意向朋友推荐“博客详情页”吗？

强烈不推荐
不推荐
一般般
推荐
强烈推荐

提交

最新文章

目录

评论 1

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。