首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 750 毫秒
1.
在地球物理电磁勘探领域有限元数值模拟中,最后都会得到一个大型稀疏的复系数线性方程组,受计算机内存空间的限制,必须根据有限元刚度矩阵的稀疏性对其进行压缩存储。由于电磁场有限元计算的自由度大都在三个以上,因而提出了适合多自由度的块按行压缩稀疏存储方案,并通过存储格式的转换,把块按行压缩方式转换成流行的,大型稀疏矩阵的行压缩存储格式,以便于求解。用求解大型稀疏方程组的Krylov子空间方法中的稳定双共轭梯度(Bicgstab)方法,收敛速度快,精度高,而且稳定性好,结合ilu预处理技术,可以大大提高求解大型稀疏方程组的效率。  相似文献   

2.
介绍了对稀疏矩阵进行压缩存储时,稀疏矩阵相乘运算的基本思想和算法.在此基础上,探讨了避免矩阵中零元素相乘的无效操作,奠定了该算法的矩阵变换基础.  相似文献   

3.
在海量影像数据的管理中 ,数据的压缩是关键技术。基于提升方法构造的整数小波变换 ,采用行索引稀疏存储方案 ,可以对海量影像数据进行高保真的压缩。与现有的地理信息系统基础软件平台结合起来 ,进行二次开发 ,可以建立起正射影像数据库系统框架。通过多级索引指针 ,可以实现影像数据的快速调度与显示 ,实时漫游、局部放大和缩小、多分辨显示等操作。  相似文献   

4.
在三维电阻率的正反演计算中,快速、准确的正演计算是反演的关键。而正演计算往往涉及到求解大型线性方程组Ax=b的问题,通过Lanczos迭代构造出对称三对角阵方程组,并采用正交分解法进行求解,与传统算法相比,此算法占用内存少、收敛速度快、且稳定;针对大型稀疏矩阵的特点,采用简单地记录矩阵的非零元素值及其所在行、列值的方法,来存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。通过理论分析和点电源三维地电场计算实例,阐述该法是地电三维正演计算的有效方法。  相似文献   

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.
赵强  焦玉勇  张秀丽  谢壁婷  王龙  黄刚海 《岩土力学》2019,40(11):4515-4522
非连续变形分析方法(DDA)是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。  相似文献   

11.
PSO-LSSVM模型在位移反分析中的应用   总被引:4,自引:1,他引:3  
邬凯  盛谦  梅松华  李佳 《岩土力学》2009,30(4):1109-1114
提出了一种基于均匀设计原理、最小二乘支持向量机(LSSVM)和粒子群优化算法(PSO)的快速位移反分析方法。该方法利用均匀设计和有限差分法获得学习样本,再用粒子群算法搜索最优的最小二乘支持向量机模型参数。并用最小二乘支持向量机回归模型建立反演参数与监测点位移值之间的非线性映射关系,最后用粒子群算法从全局空间上搜索与实测位移最吻合的反演参数。该反演模型利用了粒子群算法高效简单、均匀设计构造高质量小样本以及最小二乘支持向量机的小样本、泛化性能好的特点。将该模型应用于龙滩水电站左岸地下厂房区岩体地应力场的反演分析中,计算结果与实测的位移值和地应力值均吻合较好,说明了该模型在岩土工程快速反演分析中具有良好的应用价值。  相似文献   

12.
嵌入式空间索引策略   总被引:4,自引:0,他引:4  
谢忠  凤鸣  马常杰 《地球科学》2006,31(5):653-658
嵌入式GIS的移动终端设备被有限的资源所束缚.使得需要使用合适的策略去解决这些问题如有限存储与大数据量之间的矛盾以及适应实时的需求.分别总结了移动终端和GIS空间数据的特点, 如移动终端的效率性和空间数据的复杂性等.将嵌入式GIS应用划分为3类: 电子地图浏览、导航系统和野外数据采集.综合描述已有的索引算法: R树、四叉树、Buddy树, 同时提出了基于主题导向的空间索引策略.实践证明, 该索引策略对已有的索引算法扬长避短, 大大提高嵌入式GIS数据存取、内存使用和CPU性能方面的效率.   相似文献   

13.
跨孔层析成像LSQR算法研究   总被引:2,自引:0,他引:2  
目前广泛应用于工程、环境、水文等领域的探地雷达跨孔层析成像技术具有分辨率高,探测范围广,得到的图像结果直观,便于工作者分析地下地层结构等特点.层析成像的关键在于解系数矩阵为大型稀疏矩阵的方程组,笔者采用LSQR算法.此方法具有速度快,对不适定问题数值稳定以及从迭代过程很容易求得数值分析的数值的优点,用其解大型稀疏矩阵方程组得到的的结果清晰可靠,在各种层析成像中LSQR算法得到了广泛的应用.  相似文献   

14.
有限元水流计算中内存和运行效率初探   总被引:1,自引:0,他引:1       下载免费PDF全文
周杰  汪德爟 《水科学进展》2004,15(5):593-597
有限元在非恒定流计算中耗内存大、耗机时多的问题一直困扰着工程师,限制了有限元的应用.采用分步有限元方法离散二维浅水方程,并初步探索了行指标矩阵压缩存储稀疏矩阵的方法和效率以及预条件双共扼梯度法在求解有限元方法中形成的大型线性方程组中的效率,取得了满意的结果.在一定程度上缓解了有限元在非恒定流计算中存在的耗内存大、耗时多的问题.  相似文献   

15.
利用二次场算法研究了可控源音频大地电磁法二维正演问题。采用有限单元法进行正演模拟,将矩阵压缩存储和共轭梯度解方程方法应用到正演算法中,加快了正演算法的速度,并且将地形因素考虑到正演算法中。通过不同的模型验算,检验了算法的精度。  相似文献   

16.
为实现井—地激发极化数据的快速反演,分析三维有限元正演四面体网格剖分形成的系数矩阵的元素规律,应用MSR非零元素压缩存储和SSOR-PCG方程求解技术,提高了正演执行效率。在正演基础上,利用Jacobian矩阵与正演方程的简单关系,采用共轭梯度法求解最小二乘目标函数的模型修改量,避开Jacobian矩阵的直接计算和存储,实现了三维井-地激电数据的分钟级快速反演。模型正演合成数据的反演结果表明,能基本再现已知的模型。  相似文献   

17.
双列块法用于有限元分析计算中求解大型线代数方程组。目前在有限元计算中最通常使用的求解大型线代数方程组的方法是分块的三角分解法。由于这种方法对方程组系数矩阵的分块受带宽的限制‘以至目前一般计算机的内存贮量限制了大带宽问题的求解,使一些大型的有限元课题,特别是三维问题的解算难以实现。双列块法解决了这个西难。双列块法对方程组系数矩阵实行按计算机许可的容量划分列块,按列分解,可以完全不受带宽的限制。方法是成功的、有效的。用双列块法在PE-3220小型机上仅用686K字节的内存解算了最大半带宽为1746,系数矩阵存贮量为5000K以上字节的三维(124个20节点等参元,8个16节点等参节理元、831个节点)有限元课题,这种方法使小型和高档微机用于大型有限元的解算成为可能。  相似文献   

18.
二维电法有限元正演计算的几个问题   总被引:3,自引:0,他引:3  
本文讨论了利用微型计算机进行二维电法有限元正演计算中节省内存及减少运算量的几项措施,即刚度矩阵存储、强加边界条件处理及有限元方程组解法等问题,并提出了相应的算法。算例表明,此方法是可行的。  相似文献   

19.
系数矩阵存储和线性方程组求解是限制三维电磁积分方程方法发展的主要因素。Zhdanov提出准线性(QL)近似技术,建立了复杂散射场与背景场的线性关系,有效地避免了积分方程中大型线性方程组的求解,但是该算法用于多源问题航空电磁正演模拟时精度不高。因此,本文提出一种基于多重网格准线性(MGQL)近似的算法,并利用系数矩阵的Toeplitz性质存储矩阵和快速傅里叶变换,实现了矩阵与向量的快速乘积、降低了计算复杂度,采用多重网格结合了积分方程方法和准线性近似解法的优点,在保证精度的条件下提高计算速度、减少存储量。针对不同类型网格的模拟实验表明,相比于传统积分方程方法,本文算法在保证计算精度的同时,可以将计算速度极大地提高(>10倍)。  相似文献   

20.
采用一维压缩存储正演计算中的对称稀疏矩阵,进行Cholesky分解,利用分解后二个矩阵的对称性和稀疏性,对占用时间较多的回代过程采用先消去列的方法,实现快速回代.算例表明,采用该方法,对于点源场的求解与传统顺代回代求解法对比可以提高五倍的速度,对于大地电磁的正演问题,提高了二倍的速度.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号