共查询到20条相似文献，搜索用时 797 毫秒
1.
在地球物理电磁勘探领域有限元数值模拟中,最后都会得到一个大型稀疏的复系数线性方程组,受计算机内存空间的限制,必须根据有限元刚度矩阵的稀疏性对其进行压缩存储。由于电磁场有限元计算的自由度大都在三个以上,因而提出了适合多自由度的块按行压缩稀疏存储方案,并通过存储格式的转换,把块按行压缩方式转换成流行的,大型稀疏矩阵的行压缩存储格式,以便于求解。用求解大型稀疏方程组的Krylov子空间方法中的稳定双共轭梯度（Bicgstab）方法,收敛速度快,精度高,而且稳定性好,结合ilu预处理技术,可以大大提高求解大型稀疏方程组的效率。 相似文献
2.
针对重力勘探中光滑反演存在的分辨率较低的问题,本文提出一种基于地质体埋深、地层倾向等一定先验信息的局部光滑约束三维反演算法,并提出了一种光滑反演中粗糙度矩阵的存储方式,该存储方式可以将M×N维粗糙度矩阵存储为M×2维,减少了计算机计算内存,且详细介绍了该存储方式下粗糙度矩阵与其他矩阵相乘时,粗糙度矩阵所存储的位置信息的读取方式以及与其他矩阵逐列逐步相乘最终得到计算结果的过程。最后,利用文中提出的算法对理论模型和实测数据进行反演试算,结果表明局部光滑反演算法相比于全局光滑反演结果更加准确,且该算法在一定噪声水平下依然稳定,在实际生产中有效可行。 相似文献
3.
4.
大型对称变带宽方程组的Cholesky分解法 总被引：6，自引：3，他引：3
作者针对地球物理数值模拟中常碰到的大型稀疏变带宽方程组的求解问题,介绍了 大型稀疏变带宽矩阵的存储方法及Ｃｈｏｌｅｓｋｙ分解法,该方法的特点是用二个一维数组,其中一个输入时存储对称稀疏矩阵变带宽内的元素,输出时存储Ｃｈｌｏｅｓｋｙ下三角短阵带宽内的元素;另一个存储对称变带宽矩阵对角线元素在前一维数组中的位置,大大节约了所占的计算内存空间。 相似文献
5.
在许多有限元、有限差分的应用中,待解方程组的系数矩阵是大型稀疏带状阵,方程组的求解一般使用迭代法.与其它存储方式相比,压缩对角存储由于不存矩阵元素的行列索引,对内存的使用最为节省.使用压缩对角存储,沿对角线操作可以完成高效的矩阵—向量乘.由于在以往的文献中,没有提及按行、列操作的算法,压缩对角存储的应用范围受到一定的限制.用行、列、对角版的矩阵—向量乘代表普通意义的行、列、对角方向操作模式,通过等价矩阵推导和伪代码,给出了一个包括按行、列操作算法在内的应用框架.这里运用C语言实现不完全Cholesky分解共轭梯度法解方程,阐述了如何在实际编程中使用这个算法框架.经理论与实验分析表明,对角压缩存储应用于框架中是高效的,因为与使用常用的一维行索引存储格式所编程序相比,同样迭代次数的耗时减少了约25％. 相似文献
6.
井地直流电法三维数值模拟中若干问题研究 总被引：3，自引：0，他引：3
讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全Cholesky分解稳定化的双共轭梯度算法(ICBG),通过与均匀半空间垂直线源解析解的对比,证明了该算法是准确可靠的。 相似文献
7.
所谓近似对称矩阵是指矩阵中仅有极少一部分元素是非对称的,在将对角线以上的非对称元素用其对角线以下的对称元素替代后,该矩阵就变成了一个对称矩阵。在求解非线性连续介质力学问题时常常会遇到近似非对称矩阵。基于Sherman-Morrison定理,给出了一种新的近似非对称矩阵的分解算法。在确保数值稳定性的前提下,无论在求解效率还是在内存开销方面新算法都优于一般的非对称稀疏矩阵的求解器,且仅需对传统的基于LDLT分解的求解器略做修改,即可开发出适应于对称和非对称稀疏矩阵的求解器。最后用一个摩擦接触算例,显示了新算法的优越性。 相似文献
8.
9.
针对传统LSQR反演算法计算时间长、占用内存大的不足,这里提出了一种LSQR的快速算法,通过线性三元组存储稀疏矩阵,实现相应的矩阵运算,使LSQR算法的运算速度及内存占用都有较大地改善。 相似文献
10.
基于压缩感知理论的缺失地震数据重构方法 总被引：1，自引：0，他引：1
压缩感知是一种新的理论,它打破常规尼奎斯特-香农采样定理的制约,利用信号的稀疏特性或可压缩特性,用较少的数据即可重构恢复完整的信号。建立了基于压缩感知理论的缺失地震数据重构模型。首先在与稀疏变换不相关的测量矩阵基础上引入一种约束矩阵,使地震数据的缺失满足或接近高斯随机分布;随机缺失的地震数据变换到稀疏域会产生很多与有效信号不相干的随机噪声,接着通过一种新的自适应阈值迭代算法可以很好地消除稀疏系数中的随机噪声干扰,经过逆稀疏变换即得到重构后的地震数据。Marmousi 2模型测试及实际地震资料处理均验证了该方法的可行性和有效性。重构缺失地震数据取得了较好的效果。 相似文献
11.
Matrix formulation of co-kriging 总被引：11，自引：0，他引：11
Donald E. Myers 《Mathematical Geology》1982,14(3):249-257
The matrix form of the general co-kriging problem is presented. Matrix solutions are given for SRFs with covariance functions, for IRFs of order zero using variograms and for universal co-kriging. General methods for obtaining cross-covariance or cross-variogram models are given. The relationship of the general co-kriging problem to the problem of one under sampled variable is presented. 相似文献
12.
钻头胎体硬度是一项重要技术指标,应与所钻岩石的硬度（可钻性）和研磨性相适应。胎体硬度直接影响钻探技术经济指标的提高和钻进工艺参数的选择。试验研究表明,钻头胎体硬度分布很不均匀,差别很大。首先应从胎体烧结工艺上找原因,烧结压力偏小（5 MPa）可能是首要原因。钻头、工艺参数、操作技术等都很重要,缺一不可。钻进技术经济指标低,不一定都是钻头质量问题,要从多方面进行研究。提倡把经验打钻提高到科学打钻上来。 相似文献
13.
Subhash Lele 《Mathematical Geology》1993,25(5):573-602
Euclidean Distance Matrix Analysis (EDMA) of form is a coordinate free approach to the analysis of form using landmark data. In this paper, the problem of estimation of mean form, variance-covariance matrix, and mean form difference under the Gaussian perturbation model is considered using EDMA. The suggested estimators are based on the method of moments. They are shown to be consistent, that is as the sample size increases these estimators approach the true parameters. They are also shown to be computationally very simple. A method to improve their efficiency is suggested. Estimation in the presence of missing data is studied. In addition, it is shown that the superimposition method of estimation leads to incorrect mean form and variance-covariance structure. 相似文献
14.
15.
基体分离-电感耦合等离子体质谱测定重晶石中超痕量稀土元素 总被引：2，自引：10，他引：2
采用Na2O2、Na2CO3 碱熔重晶石,水提取,氨水沉淀稀土元素,经氢型阳离子交换树脂交换分离富集稀土元素,应用电感耦合等离子体质谱对重晶石中超痕量稀土元素进行测定,用铑作内标元素补偿基体效应和灵敏度漂移。重晶石中钡的分离效率达到99. 999%,检出限为0. 11~6. 9ng/g,相对标准偏差(RSD,n=9)≤13. 7%。 相似文献
16.
17.
Edward R.D. Scott Alan E. Rubin G.Jeffrey Taylor Klaus Keil 《Geochimica et cosmochimica acta》1984,48(9):1741-1757
Matrix material in type 3 chondrites forms rims on chondrules, metal-sulfide aggregates, Ca,Al-rich inclusions and chondritic clasts; it also forms lumps up to a millimeter in size, which may contain coarser silicates. Chondrules of all types were found with internal matrix lumps that appear to have entered the chondrules before the latter had crystallized. Mean concentrations of Mg, Na, Al and Ca in matrix occurrences show up to fivefold variations in a single chondrite. Variations between mean matrix compositions of individual type 3 ordinary chondrites are almost as large and partly reflect systematic differences between H, L and LL matrices. Such variations are probably a result of nebular separation of feldspathic material and ferromagnesian silicates.Compositions of chondrules and their matrix rims are normally unrelated, although rim compositions are correlated with those of matrix lumps inside chondrules. A single chondrule was found with a composition nearly identical to that of its internal matrix lump, suggesting that some chondrules may have formed from matrix material. Matrix lumps are as heterogeneous as chondrules, but mean chondrule and matrix compositions differ, even allowing for possible loss of metallic Fe,Ni during chondrule formation. Since bulk compositions of matrix lumps and rims have probably not changed significantly since their formation except for Fe-Mg exchange, our matrix samples cannot represent typical chondrule precursor materials. 相似文献
18.
19.
本文以鄂尔多斯盆地东部镇川堡含气区的下二叠统下石盒子组为例,提出对等时沉积序列进行最大熵谱分析的方法,以便寻找各地层序列间的历史成因联系和彼此间的亲疏程度,从而较好地解决了诸如古水流方向、河道形态及湖泊岸线位置等问题,研究结果表明,镇川堡含气区下石盒子组沉积物源来自东北付谷、宝德一带,古水流方向北东30-60°,发育两个河流系统,稳定主河道4条.下石盒子组地层序列最大熵谱分析的最佳滤波阶数m=62,谱跨度K的最佳参数为200. 相似文献
20.
Michael W. Davis 《Mathematical Geology》1987,19(2):99-107
An algorithm for producing a nonconditional simulation by multiplying the square root of the covariance matrix by a random vector is described. First, the square root of a matrix (or a function of a matrix in general) is defined. The square root of the matrix can be approximated by a minimax matrix polynomial. The block Toeplitz structure of the covariance matrix is used to minimize storage. Finally, multiplication of the block Toeplitz matrix by the random vector can be evaluated as a convolution using the fast Fourier transform. This results in an algorithm which is not only efficient in terms of storage and computation but also easy to implement. 相似文献