瞬变电磁法一维正演问题

MiniTime

一维正演我用的汉克尔变换和余弦变换，时间从10^-6到10^1(s)，但是例如对均匀地层做正演，在时间晚期，就是接近1s及以后的时间范围，计算得到的感应电动势是负值，数量级为10^-10左右，但是正常情况，负值应该为10^-11以及更小数量级，但是C语言编出来的结果不对（所有变量为float，double我也试过，结果也不对），。1）是不是滤波系数不精确？2）时间后期要处理？

◇巳◆伪亽妻

进来看看了，想学习下。

昆仑剑客

好东东！想学习学习一下。

zp了了

学习了，谢谢楼主分享。

yeb9420

小木虫上看到你了！

MiniTime

yeb9420 发表于 2014-4-12 20:41
小木虫上看到你了！

呵呵，就是啊

canefly

这里面有几个小窍门：
1，均匀大地的瞬变电磁场响应值，你自己好好用MATLAB绘绘图，就明白时间与精度关系了。然后自己想办法解决吧。
2。尽量用余玄变换。或者直接点用A.D.chave的CHAVE数值积分程序，那个精度更高。
3.你自己也知道晚期数据的数量级是什么概念。double是必须的。至于编程语言，FORTRAN、MATLAB、C++都不是问题
是武汉地大的话，我可以给你看看算100s甚至更长的程序和漂亮图，但是那有什么意义呢？

MiniTime

canefly 发表于 2014-5-16 16:03
这里面有几个小窍门：
1，均匀大地的瞬变电磁场响应值，你自己好好用MATLAB绘绘图，就明白时间与精度关系 ...

非常感谢你的解释哈

MiniTime

zp了了 发表于 2014-4-10 10:09
学习了，谢谢楼主分享。

请问，经汉克尔滤波的数据晚期会出现负值吗？

jxl860408

路过学习一下