检索结果-南通市图书馆

融合点扩散函数的预条件黏声最小二乘逆时偏移

在线全文

学校读者我要写书评

暂无评论

石油地球物理勘探 2019年第1期54卷 73-83,I0003,I0004页

作者：姚振岸孙成禹喻志超马振东华理工大学放射性地质与勘探技术国防重点学科实验室江西南昌330013 中国石油大学(华东)地球科学与技术学院山东青岛266580 北京大学地球与空间科学学院石油与天然气研究中心北京100871

在传统黏声最小二乘逆时偏移(Q-LSRTM)中,用于反传残差数据的伴随Q传播算子也是衰减的,因此导致其成像分辨率较低。为了改善Q-LSRTM的低分辨率问题,基于点扩散函数(PSF)构建黏声去模糊滤波器,然后将其作为Q-LSRTM迭代过程中的预条件化因... 详细信息

在传统黏声最小二乘逆时偏移(Q-LSRTM)中,用于反传残差数据的伴随Q传播算子也是衰减的,因此导致其成像分辨率较低。为了改善Q-LSRTM的低分辨率问题,基于点扩散函数(PSF)构建黏声去模糊滤波器,然后将其作为Q-LSRTM迭代过程中的预条件化因子,最终实现高分辨的衰减补偿偏移,因此称为预条件Q-LSRTM。模型测试表明,预条件Q-LSRTM能获得更高的分辨率、更均衡的振幅以及更快的收敛速度。Q值偏移模型的敏感性测试表明,与Q-LSRTM相似,为了显著改善成像效果,预条件Q-LSRTM也需要相对精确的Q值模型和速度模型。

关键词：最小二乘逆时偏移黏声点扩散函数去模糊滤波器预条件

一种基于亥姆霍兹分解的大地电磁测深有限元正演预条件解法

在线全文

学校读者我要写书评

暂无评论

地球物理学报 2019年第10期62卷 3898-3911页

作者：郭泽秋董浩中国地质大学(北京)地球物理与信息技术学院北京100083 四川大学灾后重建与管理学院成都610207 地下信息探测技术与仪器教育部重点实验室北京100083 深地科学与工程教育部重点实验室四川大学成都610065

本研究针对大地电磁测深法有限元数值模拟中,迭代法求解线性方程组效率较低的问题,利用亥姆霍兹分解原理,将电场矢量双旋度方程的预条件问题转化为基于矢量位的泊松问题和基于标量位的拉普拉斯问题,并在四面体非结构化棱边元离散的情况... 详细信息

本研究针对大地电磁测深法有限元数值模拟中,迭代法求解线性方程组效率较低的问题,利用亥姆霍兹分解原理,将电场矢量双旋度方程的预条件问题转化为基于矢量位的泊松问题和基于标量位的拉普拉斯问题,并在四面体非结构化棱边元离散的情况下,借助节点元辅助网格离散上述预条件问题,进一步利用代数多重网格方法(AMG)实施求解,最终实现预条件算法.利用经典的COMMEMI理论模型进行试算并与前人的积分方程解进行对比,验证了本文数值模拟程序与预条件方法的正确性和可靠性.此外,利用不同自由度规模的实验模型对这一预条件算法的效率进行了测试.结果表明,这一算法可以有效地提升大地电磁测深法棱边有限元数值模拟迭代法的收敛性,计算效率较通用的不完全LU分解预条件算法明显更高;在较大自由度网格(>1000万)数值模拟计算中,其算法效率及内存占用相对直接解法有较大优势,也使小型工作站上利用较大自由度的有限元网格进行大地电磁测深数值模拟计算成为可能.

关键词：有限单元亥姆霍兹分解大地电磁预条件数值模拟

预条件平方Smith法求解连续Lyapunov方程

在线全文

学校读者我要写书评

暂无评论

华东理工大学学报（自然科学版） 2016年第6期42卷 881-886页

作者：蔡兆克鲍亮初鲁华东理工大学数学系上海200237

探讨了如何数值求解连续时间的Lyapunov矩阵方程AX+XA^T+BB^T=0,给出了一种预条件的平方Smith算法,该算法首先利用交替方向隐式法即ADI法处理连续Lyapunov方程,构造出含ADI参数的对称Stein方程;然后利用平方Smith法迭代产生Krylov子空... 详细信息

探讨了如何数值求解连续时间的Lyapunov矩阵方程AX+XA^T+BB^T=0,给出了一种预条件的平方Smith算法,该算法首先利用交替方向隐式法即ADI法处理连续Lyapunov方程,构造出含ADI参数的对称Stein方程;然后利用平方Smith法迭代产生Krylov子空间中的低秩逼近形式。得到一些数值实验,这些例子表明预条件平方Smith法是非常有效的。

关键词：连续Lyapunov方程预条件平方Smith算法 ADI Krylov空间

薄板弯曲问题边界元法分析中预条件GMRES算法

在线全文

学校读者我要写书评

暂无评论

哈尔滨工程大学学报 2018年第11期39卷 1809-1815页

作者：陈娟肖洪天高广运同济大学地下建筑与工程系上海200092 同济大学岩土及地下工程教育部重点实验室上海200092 山东科技大学土木建筑学院山东省土木工程防灾减灾重点实验室山东青岛266590

为探讨预条件GMRES算法在求解弹性力学问题边界元法方程时的数值特性,以不同约束条件的薄板弯曲问题边界元分析为例,讨论了Jacobi、块Jacobi和对称Gauss-Sediel迭代作为预条件的三种预条件形式对GMRES收敛性的影响,并分析了其与约束条... 详细信息

为探讨预条件GMRES算法在求解弹性力学问题边界元法方程时的数值特性,以不同约束条件的薄板弯曲问题边界元分析为例,讨论了Jacobi、块Jacobi和对称Gauss-Sediel迭代作为预条件的三种预条件形式对GMRES收敛性的影响,并分析了其与约束条件的关系。分析表明:右块Jacobi和Jacobi预条件对该问题加速收敛效果更好,且对同类边界条件简支比固支收敛快,这与边界元法系数矩阵计算中的刚体位移原理有关。并且随问题规模的增大,预条件的GMRES算法效率比Gauss消去优势更明显;当问题自由度接近1万时,后者计算时间为前者的60倍。

关键词： GMRES 边界元预条件收敛性薄板弯曲边界条件

预条件Krylov子空间法求解耦合Sylvester矩阵方程

在线全文

学校读者我要写书评

暂无评论

华东理工大学学报（自然科学版） 2015年第6期41卷 871-876页

作者：徐冬梅鲍亮蔡兆克华东理工大学数学系上海200237

基于Krylov子空间方法,求解耦合Sylvester矩阵方程并给出了预条件全局正交化方法以及预条件全局极小残量法两种方法,同时给出这两种方法的一些理论结果。数值实验的结果表明,采用预条件Krylov子空间方法求解该类方程非常有效。

关键词：预条件 Krylov子空间耦合Sylvester矩阵方程

一类预条件下高阶两参数并行Jacobi型方法的敛散性

在线全文

学校读者我要写书评

暂无评论

一类预条件下高阶两参数并行Jacobi型方法的敛散性

作者：黄江玲陕西师范大学

学位级别：硕士

在科学计算和工程计算等应用领域中,许多实际问题的研究最终可归为大型稀疏线性方程组Ax=b的求解.迭代法则是求此类方程组的有效方法,而迭代矩阵的谱半径又可衡量迭代法的收敛性和收敛速度.因此,本文研究了高阶两参数并行Jacobi型方法(... 详细信息

在科学计算和工程计算等应用领域中,许多实际问题的研究最终可归为大型稀疏线性方程组Ax=b的求解.迭代法则是求此类方程组的有效方法,而迭代矩阵的谱半径又可衡量迭代法的收敛性和收敛速度.因此,本文研究了高阶两参数并行Jacobi型方法(又称高阶2PPJ迭代法)与经典Jacobi迭代法之间的敛散性,并给出了一类预条件P=I+S下高阶2PPJ迭代法的敛散性.主要内容如下:第一章,主要介绍了文中所用到的基本定义和引理,如谱半径,L-阵,M-阵,不可约矩阵以及Perron-Frobenius定理.第二章,讨论了当线性方程组Ax=b的系数矩阵A是满足0预条件矩阵的基础上,考虑当系数矩阵A是满足0预条件矩阵P=I+S,其次运用特征向量法讨论了预条件高阶2PPJ迭代法敛散性,然后通过比较预条件作用前后迭代矩阵的谱半径大小,得到了敛散性比较定理.特别地,当系数矩阵A还满足严格对角占优L-阵时,预条件P=I+S后的A也为严格对角占优L-矩阵,且ρ(J)≤ρ(J)<1,得到在参数0<ω≤1,0<τ<1,阶数P相同时,ρ(Gω,τ,P)≤ρ(Gω,τ,P)<1;P不同时,随着P取值的增加,ρ(Gω,τ,P)越小,即收敛性越好.最后,通过数值算例对结论进行了验证.

关键词： L-阵 Jacobi迭代法高阶2PPJ迭代法预条件谱半径

同方学位论文库评论

在线全文

同方学位论文库

学校读者我要写书评

暂无评论

一种新的预条件伴随状态法初至波走时层析

地球物理学报 2017年第10期60卷 3934-3941页

作者：李勇德董良国刘玉柱同济大学海洋地质国家重点实验室上海200092

伴随状态法初至波走时层析是基于最优化理论的一种层析成像方法,该方法不必进行射线追踪,用两次正演的计算量便可以获得梯度,具有计算效率高、内存占用小等优点,但是其一阶方向在初始模型或观测孔径不理想的情况下往往无法获得正确的反... 详细信息

伴随状态法初至波走时层析是基于最优化理论的一种层析成像方法,该方法不必进行射线追踪,用两次正演的计算量便可以获得梯度,具有计算效率高、内存占用小等优点,但是其一阶方向在初始模型或观测孔径不理想的情况下往往无法获得正确的反演结果,而二阶方向的实现又比较困难且费时.在伴随状态法的基础上,将走时差替换为定值,再次进行反演,便可以得到类似于射线密度的矩阵,用该矩阵的逆可以方便地进行预条件.基于该方法,本文提出了一种简单易行的预条件伴随状态法初至波走时层析的实现方法.理论模型和实际资料处理结果都表明,该方法既保留了伴随状态法初至波走时层析的优点,又可以克服一阶方向的局限,获得良好的反演效果.

关键词：地震层析成像伴随状态法初至波走时近地表速度预条件

在线全文

学校读者我要写书评

暂无评论

预条件下含参数的JOR迭代法敛散性分析

东北师大学报（自然科学版） 2015年第1期47卷 43-47页

作者：王慧勤雷刚宝鸡文理学院数学系陕西宝鸡721013

对于JOR迭代法求解线性方程组Ax=b,运用了预条件加速JOR迭代法的收敛性,在预条件后引入参数α,给出更一般的预条件下含参数形式的JOR迭代方法.证明了这类方法能够加速JOR迭代法的收敛性,找到了参数的最佳取值,并且用数值算例加以验证.

关键词： JOR迭代法收敛性预条件谱半径

基于块结构不完全近似稀疏逆的GPU预条件算法

在线全文

学校读者我要写书评

暂无评论

基于块结构不完全近似稀疏逆的GPU预条件算法

作者：胡逸文信阳师范学院

学位级别：硕士

如何构造一个有效的预条件子,用于加速求解大规模稀疏线性系统的收敛一直都是数值计算领域的研究热点之一。随着线性系统规模的增大,使用不精确的预处理方法求解大规模稀疏线性系统也逐渐开始备受关注,使用该方法可以在时间与精确度之... 详细信息

如何构造一个有效的预条件子,用于加速求解大规模稀疏线性系统的收敛一直都是数值计算领域的研究热点之一。随着线性系统规模的增大,使用不精确的预处理方法求解大规模稀疏线性系统也逐渐开始备受关注,使用该方法可以在时间与精确度之间进行权衡。近年来随着GPU的快速发展,基于GPU使用不精确的方法构造预条件子用于求解大规模稀疏线性系统已经取得了丰硕的成果。目前针对预条件子构造的研究均集中在标量矩阵,但在实际的工程应用中,如多物理场问题,系数矩阵往往是以块状矩阵的形式出现的,虽然块结构的线性系统的数值算法大多是由标量线性系统推导而来,但实现策略以及性能调优等方面仍然存在较大的差异。所以本文面向CPU+GPU的异构计算系统,对快速构造一个有效的块不完全近似稀疏逆的预条件算法进行讨论与研究,主要工作如下:(1)在不完全近似稀疏逆(Incomplete Sparse Approximate Inverses,ISAI)算法的基础上,提出一种块结构的不完全近似稀疏逆(Block Incomplete Sparse Approximate Inverses,Block-ISAI)预条件算法和在GPU加速平台上高效的实现策略。在预处理方面,为避免在预条件步骤中使用强数据依赖的backward-substitution和forward-substitution求解块三角线性系统,基于Block-ILU分解的预处理方法,使用Block-ISAI算法将预条件步骤转化为可高度并行的矩阵向量乘操作。在GPU实现方面,充分考虑了块矩阵的整体稀疏且局部稠密的特性,并结合GPU平台上融合访存和共享内存等优化技术,基于CUDA中基本调度单位warp,提出了一种以warp为单位来分配线程任务的策略,该策略能够保证每一个warp以块为单位并行计算上下三角因子的近似逆的每一列,从而充分发挥GPU细粒度高并发的计算特性。(2)在预处理步骤的求解效率方面和基于Krylov子空间方法的GMRES整体求解时间方面,将本文提出的Block-ISAI预条件算法与基于cu SPARSE的块三角线性求解的预条件算法,在由Intel E5-2640 V4@2.40 GHz处理器与NVIDIA的Tesla V100 GPU组成的异构计算平台上进行了详细的实验对比分析。在选取了来自Suite Sparse矩阵集合中几个典型的矩阵进行测试后,结果表明,虽然本文提出的Block-ISAI算法增加了迭代步数,但大大缩短了迭代中的预处理时间,总的求解时间比基于cu SPARSE预条件算法要少,加速比为1.19～6.69。

关键词：不完全近似稀疏逆预条件块线性系统 GPU Warp