收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

线性方程组预条件技术及在二维三温问题中的应用和实现

张勇  
【摘要】: 线性方程组的求解是当今不少研究领域遇到的一个共同问题,具有重大的理论价值和应用价值。而对于大型线性方程组,其阶数一般上几十万阶,甚至上百万阶,一般都要通过预处理后才能有效求解。目前预条件技术在诸多应用中远没有达到满意的程度。提高线性方程组求解的速度既是理论研究者也是工程技术应用者的追求目标。预条件技术作为线性方程组求解技术的重要组成部分,其重要性在未来还将不断增加。 本文主要研究了一些特殊矩阵和实际问题的预条件技术,同时还设计开发了线性方程组求解技术开发的Matlab软件包,主要工作如下: 1.对M-矩阵、非负矩阵在Jacobi迭代法和Gauss-Seidel迭代法下的预条件技术进行了研究。首先改进了I+S型预条件技术,并提出优化模型进行参数选择;另外利用LU分解方法,提出了两类Gauss型预条件技术,并给出了算法实现及算法比较。通过较充分的理论分析,本文证明了这些预条件技术的合理性和可靠性。数值实验进一步验证了这些预条件技术的明显改进效果和性能。 2.通过构造“好”的预条件子的机理分析,并结合稀疏近似逆技术,提出了I+H_c型预条件技术。根据不同的预条件策略和思路,分别建立了估算I+H_c型预条件子中参数的优化模型。这类预条件技术的明显优点是计算量低,可并行性高的特点。数值实验表明了对于M-矩阵的情形,该类预条件技术联合Krylov子空间迭代法的预处理效果良好。 3.本文对从一类二维PDE问题离散化得到的块三对角矩阵的预条件技术进行了研究。首先通过观察矩阵的逆有关的稀疏模式,提出了块七对角预条件子,并根据基于Frobenius模极小化这类稀疏近似逆技术的一种算法实现:最小残差算法(MR),提出了构造块七对角预条件子的改进的最小残差算法。改进的MR算法的计算量大大低于原有的MR算法,而且改进MR算法预处理效果几乎与MR算法相同。这类预条件子的有效性在Krylov子空间迭代法下进行了测试。数值实验还表明改进的MR算法对于一些非对称情形同样适用。 4.对二维三温问题的预条件技术从两个方向进行了研究。(1)对基于红黑排序下的分块系统的重排技术进行研究。在重排效果较好的hyperplane方法的基础上,提出了一类交替hyperplane(A-hyperplane)重排方法。数值方法表明,相对于RB、hyperplane预条件重排方法,当该重排预条件技术结合ILU(κ)预条件方法时,大大提高了收敛速度。(2)对二维三温问题的全局预处理,应用了一类区域分解方法进行了研究。数值实验表明了在二维三温问题的应用中,区域分解法在增加分解的区域个数以提高并行计算性的同时,迭代次数并没有明显增加,说明了区域分解对二维三温问题具有良好预处理效果。这些工作为今后基于区域分解方法对二维三温问题的并行求解做了一定准备工作。 5.为了提高线性方程组求解方法的开发效率,特别是预条件技术的开发效率,我们开发了LINSOLKIT软件包。本文主要利用了Matlab数学软件的强大功能,并在分析研究了当前一些成熟的基于非Matlab语言开发的线性方程组求解软件包的情况下,借鉴了其中的一些成功设计经验,开发了线性方程组求解的Matlab软件包:LINSOLKIT。该软件包的主要接口包括预条件子构造接口、预条件的迭代法接口、预处理子求解接口等。由于LINSOLKIT是基于Matlab开发,因而LINSOLKIT在新增预条件技术、迭代求解方法等扩充方面都非常方便、高效。本文的大部分数值实验都是直接使用LINSOLKIT提供的图形操作界面完成。另外,本文还基于SparsKit软件包开发了一组Fortran程序用于二维三温问题和二维PDE测试。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 韩凤萍;;对称循环M-矩阵的若干性质[J];大学数学;2011年02期
2 卢飞龙;何希勤;;M矩阵与其逆的Hadamard积的特征值下界[J];辽宁科技大学学报;2010年05期
3 孙仁康;贾振波;;裂纹诊断的数学原理[J];纯粹数学与应用数学;1993年01期
4 单沪军,郑强,吴强;初等r-分块循环矩阵的几个性质[J];山东工业大学学报;1999年03期
5 逄明贤;M-矩阵特征值的若干估计[J];数学杂志;1984年03期
6 谷同祥,迟学斌,刘兴平;稀疏近似逆与多层块ILU预条件技术[J];应用数学和力学;2004年09期
7 石茂忠,张晓东;对角占优矩阵的推广及应用[J];数学研究;1997年02期
8 吴世锦;M-矩阵与广义正定矩阵的关系[J];大学数学;2005年04期
9 林彤,杜宇辉;双对角占优矩阵的性质[J];吉林化工学院学报;2000年02期
10 胡文杰,逄勃;关于M-矩阵等价表征的注记[J];数学研究;2000年01期
11 唐先华;广义正定矩阵与稳定矩阵、M-矩阵的关系[J];衡阳师范学院学报;1994年03期
12 吴池业;黄廷祝;;几种矩阵逆的连续性[J];纯粹数学与应用数学;2009年03期
13 桂曙光;M-矩阵的判定[J];安徽理工大学学报(自然科学版);2004年02期
14 方秀男;汤凤香;;关于r-块对角占优矩阵的对角Schur补[J];渤海大学学报(自然科学版);2010年02期
15 张谋成,黎稳;M-矩阵分裂的迭代矩阵[J];应用数学;1994年04期
16 李庆春,田雪;广义严格对角占优矩阵的判定条件[J];北华大学学报(自然科学版);2000年04期
17 沈光星;对《矩阵对角占优性的推广及应用》一文的注记[J];科技通报;2003年04期
18 裴冀南;;一类非线性积分微分方程的全局指数稳定[J];重庆教育学院学报;2005年06期
19 徐屹;胡乡秋;;广义严格对角占优矩阵判定探究[J];通化师范学院学报;2007年12期
20 斯办更,马万彪;关于时滞微分不等式及其应用的研究[J];内蒙古师范大学学报(自然科学汉文版);1989年03期
中国重要会议论文全文数据库 前10条
1 周洪伟;陈璞;尹家聪;;自动多重子结构方法及应用[A];2007全国结构动力学学术研讨会论文集[C];2007年
2 邓爱明;张明荣;梅华生;;基于特征值的零构件疲劳可靠性统计方法[A];第七届全国MTS材料试验学术会议论文集(一)[C];2007年
3 叶康生;袁驷;;旋转梁振动分析的ODE求解器法[A];第五届全国结构工程学术会议论文集(第一卷)[C];1996年
4 仓艳;袁嗣杰;;通信信号调制识别中的信噪比估计[A];全国第一届嵌入式技术联合学术会议论文集[C];2006年
5 张军英;许进;保铮;;一类离散Hopfield网的特征分析(Ⅰ)—稳定吸引状态[A];1999年中国神经网络与信号处理学术会议论文集[C];1999年
6 周开锡;;基于spss的主成分分析在水库水质分析中的应用[A];四川省首届环境影响评价学术研讨会论文集[C];2009年
7 孙利民;;纵向振动杆特征值摄动优化设计方法[A];第二届全国结构工程学术会议论文集(下)[C];1993年
8 周华莹;金梁;季新生;;一种提高CDMA系统软切换成功率的“特征值”分析法[A];现代通信理论与信号处理进展——2003年通信理论与信号处理年会论文集[C];2003年
9 李宝成;孙峥;;特征值定理向伪谱定理推广的注记[A];江苏省现场统计研究会第八次学术年会论文集[C];2003年
10 杨晓龙;姜慧;林克昌;;锤击预应力高强混凝土管桩承载力特征值计算[A];首届全国建筑结构技术交流会论文集[C];2006年
中国博士学位论文全文数据库 前10条
1 张勇;线性方程组预条件技术及在二维三温问题中的应用和实现[D];电子科技大学;2008年
2 程光辉;大型稀疏线性代数系统迭代解法研究[D];电子科技大学;2008年
3 杭旭登;偏微分方程迭代并行解法与网格优化方法[D];中国工程物理研究院;2004年
4 任志刚;预条件算法及在电磁场数值模拟中的应用[D];电子科技大学;2010年
5 王转德;迭代矩阵的谱分析[D];电子科技大学;2009年
6 田贵贤;图谱理论和几类矩阵的谱与组合特征研究[D];电子科技大学;2009年
7 王瑞瑞;求解某些特殊稀疏线性系统的数值解法[D];厦门大学;2009年
8 颜溯;基于Calderón技术的计算电磁学积分方程方法研究[D];电子科技大学;2011年
9 向华;结构线性方程组的迭代方法与扰动分析[D];复旦大学;2006年
10 吴雅容;关于图的特征值的几个问题的研究[D];华东师范大学;2011年
中国硕士学位论文全文数据库 前10条
1 李瑶;混合类型分裂与高斯类型预条件迭代[D];扬州大学;2011年
2 张传文;新预条件下矩阵的收敛性分析及其比较[D];扬州大学;2010年
3 金一;P2P环境下混合增强型防火墙研究与实现[D];苏州大学;2005年
4 李德凡;基于遗传算法的合作伙伴评价选择研究[D];燕山大学;2006年
5 陈蓓;关于Toeplitz矩阵特征值反问题的研究[D];南京航空航天大学;2007年
6 张俊珏;单圈图依最大特征值的排序[D];解放军信息工程大学;2007年
7 周立广;关于几类微分算子谱的研究[D];内蒙古师范大学;2008年
8 吴春连;收缩定义域上的Sturm-Liouville问题[D];南京理工大学;2008年
9 宋春燕;单圈图的零度[D];新疆大学;2008年
10 周俊强;一类可靠性模型研究中出现的投影算子的表达式及其应用[D];新疆大学;2008年
中国重要报纸全文数据库 前10条
1 大连国资公司 岳阳;一季报有哪些特征值得特别关注[N];上海证券报;2008年
2 萌生;3600点附近将迎来短期反弹[N];上海证券报;2007年
3 顾洪洪;杀毒软件厂商在卖“过期药”?[N];新华每日电讯;2005年
4 本报记者  张伟;解决安全问题需综合手段[N];中国高新技术产业导报;2006年
5 顾洪洪;主动防御技术为反病毒业带来新天地[N];人民公安报;2008年
6 联合证券 王慧;底部夯实阶段显现三大特征[N];中国证券报;2007年
7 楼栋;天量成交导演股指“过山车”[N];中国证券报;2006年
8 北京首放;调整压力日益增大的背后[N];上海证券报;2007年
9 孙颖;“华赛”金奖照获奖资格取消,PS过[N];新华每日电讯;2008年
10 张玉光 方成亮;病毒年增23万何以如此张狂?[N];中国电脑教育报;2007年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978