共查询到20条相似文献,搜索用时 750 毫秒
1.
在地球物理电磁勘探领域有限元数值模拟中,最后都会得到一个大型稀疏的复系数线性方程组,受计算机内存空间的限制,必须根据有限元刚度矩阵的稀疏性对其进行压缩存储。由于电磁场有限元计算的自由度大都在三个以上,因而提出了适合多自由度的块按行压缩稀疏存储方案,并通过存储格式的转换,把块按行压缩方式转换成流行的,大型稀疏矩阵的行压缩存储格式,以便于求解。用求解大型稀疏方程组的Krylov子空间方法中的稳定双共轭梯度(Bicgstab)方法,收敛速度快,精度高,而且稳定性好,结合ilu预处理技术,可以大大提高求解大型稀疏方程组的效率。 相似文献
2.
许志文 《华东地质学院学报》2000,23(1):87-89
介绍了对稀疏矩阵进行压缩存储时,稀疏矩阵相乘运算的基本思想和算法.在此基础上,探讨了避免矩阵中零元素相乘的无效操作,奠定了该算法的矩阵变换基础. 相似文献
3.
在海量影像数据的管理中 ,数据的压缩是关键技术。基于提升方法构造的整数小波变换 ,采用行索引稀疏存储方案 ,可以对海量影像数据进行高保真的压缩。与现有的地理信息系统基础软件平台结合起来 ,进行二次开发 ,可以建立起正射影像数据库系统框架。通过多级索引指针 ,可以实现影像数据的快速调度与显示 ,实时漫游、局部放大和缩小、多分辨显示等操作。 相似文献
4.
5.
井地直流电法三维数值模拟中若干问题研究 总被引:3,自引:0,他引:3
讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全Cholesky分解稳定化的双共轭梯度算法(ICBG),通过与均匀半空间垂直线源解析解的对比,证明了该算法是准确可靠的。 相似文献
6.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术,构造出三对角阵方程组,然后采用正交分解法进行求解,它是Krylov子空间方法中的一种。与传统迭代算法相比,它占用内存少,收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点,采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。理论分析和计算实例显示,此算法是地电三维正演计算的有效方法,为下一步的反演计算打好基础。 相似文献
7.
为提高网格环境下海量空间数据管理与并行化处理效率,将网格环境下的分布并行处理技术与空间索引相融合,提出了一种空间索引框架(grid slot and hash R tree,GSHR-Tree).该索引树结构基于散列hash表和动态空间槽,结合R树结构的范围查询优势和哈希表结构的高效单key查询,分析改进了索引结构的组织和存储.构造了适合于大规模空间数据的网格并行空间计算的索引结构,该索引树算法根据空间数据划分策略,动态分割空间槽,并将它们映射到多个节点机上.每个节点机再将其对应空间槽中的空间对象组织成R树,以大节点R树方式在多个节点上分布索引数据.以空间范围查询并行处理的系统响应时间为性能评估指标,通过模拟实验证明,该GSHR-Tree索引满足了当前网格环境空间索引的需要,并具有设计合理、性能高效的特点. 相似文献
8.
针对重力勘探中光滑反演存在的分辨率较低的问题,本文提出一种基于地质体埋深、地层倾向等一定先验信息的局部光滑约束三维反演算法,并提出了一种光滑反演中粗糙度矩阵的存储方式,该存储方式可以将M×N维粗糙度矩阵存储为M×2维,减少了计算机计算内存,且详细介绍了该存储方式下粗糙度矩阵与其他矩阵相乘时,粗糙度矩阵所存储的位置信息的读取方式以及与其他矩阵逐列逐步相乘最终得到计算结果的过程。最后,利用文中提出的算法对理论模型和实测数据进行反演试算,结果表明局部光滑反演算法相比于全局光滑反演结果更加准确,且该算法在一定噪声水平下依然稳定,在实际生产中有效可行。 相似文献
9.
二维视电阻率断面的快速最小二乘反演 总被引:12,自引:0,他引:12
快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方 相似文献
10.
非连续变形分析方法(DDA)是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。 相似文献
11.
PSO-LSSVM模型在位移反分析中的应用 总被引:4,自引:1,他引:3
提出了一种基于均匀设计原理、最小二乘支持向量机(LSSVM)和粒子群优化算法(PSO)的快速位移反分析方法。该方法利用均匀设计和有限差分法获得学习样本,再用粒子群算法搜索最优的最小二乘支持向量机模型参数。并用最小二乘支持向量机回归模型建立反演参数与监测点位移值之间的非线性映射关系,最后用粒子群算法从全局空间上搜索与实测位移最吻合的反演参数。该反演模型利用了粒子群算法高效简单、均匀设计构造高质量小样本以及最小二乘支持向量机的小样本、泛化性能好的特点。将该模型应用于龙滩水电站左岸地下厂房区岩体地应力场的反演分析中,计算结果与实测的位移值和地应力值均吻合较好,说明了该模型在岩土工程快速反演分析中具有良好的应用价值。 相似文献
12.
13.
14.
15.
利用二次场算法研究了可控源音频大地电磁法二维正演问题。采用有限单元法进行正演模拟,将矩阵压缩存储和共轭梯度解方程方法应用到正演算法中,加快了正演算法的速度,并且将地形因素考虑到正演算法中。通过不同的模型验算,检验了算法的精度。 相似文献
16.
17.
双列块法用于有限元分析计算中求解大型线代数方程组。目前在有限元计算中最通常使用的求解大型线代数方程组的方法是分块的三角分解法。由于这种方法对方程组系数矩阵的分块受带宽的限制‘以至目前一般计算机的内存贮量限制了大带宽问题的求解,使一些大型的有限元课题,特别是三维问题的解算难以实现。双列块法解决了这个西难。双列块法对方程组系数矩阵实行按计算机许可的容量划分列块,按列分解,可以完全不受带宽的限制。方法是成功的、有效的。用双列块法在PE-3220小型机上仅用686K字节的内存解算了最大半带宽为1746,系数矩阵存贮量为5000K以上字节的三维(124个20节点等参元,8个16节点等参节理元、831个节点)有限元课题,这种方法使小型和高档微机用于大型有限元的解算成为可能。 相似文献
18.
二维电法有限元正演计算的几个问题 总被引:3,自引:0,他引:3
本文讨论了利用微型计算机进行二维电法有限元正演计算中节省内存及减少运算量的几项措施,即刚度矩阵存储、强加边界条件处理及有限元方程组解法等问题,并提出了相应的算法。算例表明,此方法是可行的。 相似文献
19.
系数矩阵存储和线性方程组求解是限制三维电磁积分方程方法发展的主要因素。Zhdanov提出准线性(QL)近似技术,建立了复杂散射场与背景场的线性关系,有效地避免了积分方程中大型线性方程组的求解,但是该算法用于多源问题航空电磁正演模拟时精度不高。因此,本文提出一种基于多重网格准线性(MGQL)近似的算法,并利用系数矩阵的Toeplitz性质存储矩阵和快速傅里叶变换,实现了矩阵与向量的快速乘积、降低了计算复杂度,采用多重网格结合了积分方程方法和准线性近似解法的优点,在保证精度的条件下提高计算速度、减少存储量。针对不同类型网格的模拟实验表明,相比于传统积分方程方法,本文算法在保证计算精度的同时,可以将计算速度极大地提高(>10倍)。 相似文献
20.
采用一维压缩存储正演计算中的对称稀疏矩阵,进行Cholesky分解,利用分解后二个矩阵的对称性和稀疏性,对占用时间较多的回代过程采用先消去列的方法,实现快速回代.算例表明,采用该方法,对于点源场的求解与传统顺代回代求解法对比可以提高五倍的速度,对于大地电磁的正演问题,提高了二倍的速度. 相似文献