共查询到20条相似文献,搜索用时 62 毫秒
1.
基于物理模型的成像算法评价与SASART成像算法 总被引:1,自引:1,他引:0
基于物理模型的图像重建算法评价方法,作者研究设计的SASART算法,给出了算法SVD,CG,LSQR,阻尼LSQR,SIRT,SART及SASART的测试结果,测试数据表明:(1)线性成像方程系统的特性及解结构都对解精度有影响,解模型越粗糙,解的精度越低;(2)自激励联合迭代重重算法(SASART)迭代稳定、抗噪音能力强,用于高噪数据反能获得合理的图像;(3)各种求解算法都具有平滑效应,同时也都会 相似文献
2.
相位解缠是干涉雷达(InSAR)技术测量地形高程的重要步骤之一。由噪声和欠采样导致的相应不一致性是相位解缠的难点。本文详细介绍了各种相应解缠算法,并从理论上比较了各种算法的优缺点,分析了相位解缠算法未来的发展趋势。 相似文献
3.
病态线性方程组的简单迭代解法 总被引:9,自引:0,他引:9
针对地球物理反问题中经常碰到的病态线性方程组的求解问题。本文提出了一种简单迭代(SI)算法,从理论上证明了解序列收敛且收敛到方程组的真解,然后给出了几个算例,将计算结果与对付病态问题能力很强的CG类算法的结果进行了对比,结果表明:SI算法具有极强的抗病态能力,计算精度明显高于CG类算法,但计算速度稍低于后者 相似文献
4.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格(aggregation-based algebraic multigrid, AGMG)算法引入大地电磁三维正演模拟中。首先从准静态条件下的麦克斯韦方程出发,利用交错网格有限体积法进行离散,并采用第一类Dirichlet边界条件形成大型稀疏复线性方程组;然后阐述AGMG算法的粗化策略和套迭代技术,并实施3种不同的AGMG求解算法:1)传统的V循环AGMG算法;2)AGMG预处理共轭梯度(AGMG-CG);3)AGMG预处理广义共轭残差法(AGMG-GCR)。最终实现大地电磁法三维正演模拟。对典型地电模型进行正演模拟,并与已有的大地电磁三维正反演程序(ModEM)进行结果对比,以验证本文算法的准确性。另外,不同剖分网格和极化方式正演模拟结果与准残量最小化(QMR)迭代算法的对比表明,AGMG预处理求解算法(AGMG-CG、AGMG-GCR)不仅能够改善算法的稳定性,而且能够快速有效地求解正演问题;其中AGMG-GCR迭代次数更少,求解速度更快,误差衰减曲线更光滑,在144×152×104网格剖分情况下,相对于现有ModEM程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
5.
二维声波方程稳健迭代速度反演 总被引:1,自引:0,他引:1
从二维声波方程初、边值问题出发,通过引进一个广义幂指数误差分布函数,以及对波动方程的Lippman-Schwinger方程的解进行Bon近似,导出一个残差加权迭代最小二乘(稳健迭代)算法,进而实现对地下介质剖面的速度结构以演计算。模拟计算表明,本方法具有较高的精度和较强的抗干扰能力,是求解声波方程反演问题的一种有效方法,同时也为地下介质的速度分布成像提供一种新的技术。 相似文献
6.
本文介绍一种基于γ-循环阵快速分解方法的、新的基数为2的DFT快速算法;该算法与具有最小乘法数的Preuss算法比较,新算法的乘法数接近Preuss算法,但其加法数仅为它的75%。在Sun-SPARCStation2工作站上,采用Fortran77语言编程实现,结果表明:新算法在阶数(N)小于1024时,速度比分裂基(Split—Radix)算法快。如果进一步考虑算法中的固定常数乘法,那么新算法将会更加有效。 相似文献
7.
ρ(L_(r。w))的界和迭代法收敛性刘光清,辛迪远(长春地质学院基础科学系.长春130026)(四平电视大学.四平136000)本文给出了广义AOR法(广义SSOR,广义SOR,广义G-S及广义西法)迭代矩阵谱半径较理想的上界及其收敛准则,推广了文?.. 相似文献
8.
9.
《物探化探计算技术》2016,(5)
经典位场物性反演是建立在数与模型之间的线性关系,利用离散正演公式将每个剖分单元相对观测点的核函数进行计算,以寻求观测数据与核函数有最好的相关性。该类方法在计算过程最大的问题就是在解巨大超定方程组时,导致求解极不稳定。在前人的研究基础上,提出一种基于先验密度模型约束正则化条件下的重力异常视密度求解方法,通过L曲线法选取最佳的参数,使求解过程中有较高的稳定性和求解精度。这里以一理论模型作为观测数据,并引入噪声作为观测误差,通过算法的对比表明,该方法有很高的计算精度和稳定性,对于需要求解大型超定矩阵的三维模型,在给定先验初始模型的条件下,该方法仍然能够取得很好的效果。 相似文献
10.
在一单行本中,我们提出了一种CSAMT正演模拟算法并说明了这种算法在解释野外数据中的应用效果,在本文中,提出了一种用于二维电阻率模型的CSAMT数据反演算法,这种算法包含一个三维源并可以反演(TE)模式及(TM)模式下的数据,用这种算法对假设的数据进行试算,结果表明,不用做平面波校正就能对远区场,过渡区场及近区场进行反演,最终的模型结果显示,原始电阻率模型可被恢复,而不存在任何由上源于效应而产生的 相似文献
11.
在重力资料反演解释中,正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子,求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中,要求对所求的解答预先有一个初步的推断解释。为此,要求对所解释的问题先进行一番研究,并对综合利用各种方法和手段所获得的先验资料进行对比,提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面,与正则化法和其它线性方程解法的结果作对比,说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。 相似文献
12.
线性方程组的正则化解法及其应用 总被引:3,自引:0,他引:3
在重力资料反演解释中,正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子,求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中,要求对所求的解答预先有一个初步的推断解释。为此,要求对所解释的问题先进行一番研究,并对综合利用各种方法和手段所获得的先验资料进行对比,提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面,与正则化法和其它线性方程解法的结果作对比,说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。 相似文献
13.
14.
15.
In mineral deposit evaluation, the estimation of block values can involve the solution of as many as 100,000 systems of linear equations. To minimize computing costs it is important that the program to solve the linear systems be as efficient as possible. Standard library routines tend to contain very general stable methods using pivoting and iteration to ensure high accuracy solutions. Here, a particular method has been developed that gives a substantial improvement in run time at the occasional expense of accuracy and stability. The method is not new—just a simple application of Gaussian elimination in the symmetric case. 相似文献
16.
利用经典土压力理论设定合理土条推力线范围,根据静力平衡及力矩平衡条件建立以条块界面安全系数 为变量的线性超定方程组,应用MATLAB软件基于最小二乘方法对此方程组求解,依据极大值原理及合理性条件: ,最终获得边坡整体安全系数 。经算例验证得出如下结论:当条块界面与条块整体安全系数相等时,边坡抗剪能力发挥最大;安全系数随推力线位置的升高而减小,但影响并不大;由于考虑了极大值原理,该法求得的安全系数略微偏大于其他方法。 相似文献
17.
Conditioning of coefficient matrices of Ordinary Kriging 总被引:1,自引:0,他引:1
R. J. O''Dowd 《Mathematical Geology》1991,23(5):721-739
The solution of a set of linear equations is central to Ordinary Kriging. Computers are commonly applied because of the amount of data and work involved. There has, until recently, been little attention devoted toward the conditioning of kriging matrices. This article considers implications of conditioning upon numerical stability, instead of on robustness which has been the main focus of past work. The effect of properties of the stationary covariance matrix on the conditioning of the kriging matrix is discussed. The relationship between the covariance and autocorrelation functions allows some conclusions about the conditioning of covariance matrices, based on past work in deconvolution. The conditioning of some coefficient matrices of stationary kriging, defined in terms of either the semivariogram or the covariance, is examined. 相似文献
18.
19.
20.
提出了一种根据非线性程度求解非线性方程组的ABS算法 ,该算法根据曲线的曲率建立非线性程度。初步的数值试验表明 ,多数情况下本文建立的ABS算法比原来的非线性ABS算法收敛快或与原ABS算法迭代次数相同 相似文献