变带宽存储形式下对称正定稀疏方程组的解法

来源 :数学学习与研究 | 被引量 : 0次 | 上传用户:zcc8541099
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  【摘要】稀疏方程组一般都是大型系数矩阵方程组,如果按照直接法求解,往往需要耗费大量的内存资源和计算时间.本文则给出了一种使用变带宽压缩存储的方式求解对称正定稀疏方程组的高效解法.
  【关键词】变宽带存储;稀疏方程组;解法
  稀疏方程组一般都是大型系数矩阵方程组,如果按照直接解法的一般原理,往往需要耗费大量的内存资源和计算时间.但是对于一些特殊形式的系数矩阵,我们可以根据其自身特点寻找相应的矩阵存储形式,进而得到一种更为节省内存、提高运算速度的新解法.这里我们以对称正定稀疏方程组为例,应用变带宽存储的方法压缩储存系数矩阵,给出相应方程组的高效算法.
  下用变带宽压缩存储法求解对称正定稀疏线性方程组
  AX=B,(1-1)
  其中A∈Rn×n是对称正定稀疏阵,X,B∈Rn×m(m≥1).
  先对A,再对B分别作如下的分解:
  A=LDLT,(1-2)
  A=LDB~.(1-3)
  其中L为单位下三角矩阵(lii=1,i=1,2,…,n,lij=0,i  设矩阵A各行的第一个非零元素是aimi,i=1,2,…,n.记mij=max(mi,mj),则由式(1-2)和(1-3),根据矩阵乘法规则,推出矩阵L,D,B~的元素的下列计算公式.
  lij=(aij-∑j-1k=mijlikdkljk)/dj,j=mi,mi+1,…,i-1,i=2,3,…,n,
  di=aii-∑i-1k=mil2ikdk,i=1,2,…,n,
  b~ik=(bik-∑i-1j=milijdjb~jk)/di,i=1,2,…,n,k=1,2,…,m.
  将式(1-2)和(1-3)代入式(1-1),得到
  LTX=B~.(1-4)
  这样,方程组(1-1)的求解就归结为方程组(1-4)的求解.
  以下给出对称正定稀疏方程组(1-1)的算法,其中A∈Rn×n是对称正定稀疏阵,X,B∈Rn×m(m≥1).
  1.采用变带宽压缩存储方式
  用一维数组a(1:p)和一维整型数组d(0:n)存放对称矩阵A.
  这里a(1:p)=(a11,a2m2,a22,a3m3,…,a33,…,aimi,…,aii,…,anmn,…,ann),其中p=∑ni=1(i-mi+1)是数组a中元素的个数.另外用d(0)=0,d(i),i=1,2,…,n来标记对角元素aii在数组a中的序号.矩阵A的元素aij处于数组a中的第d(i)-i+j个位置上,即
  aij=a(d(i)-i+j),
  aimi中的列下标mi的计算公式为
  mi=i-(d(i)-d(i-1))+1,i=1,2,…,n.
  矩阵B存放在二维数组b(1:n,1:m)的位置.最终解X也存放在b(1:n,1:m)的位置.
  2.计算矩阵L,D和B~
  L的元素存放在数组a(1:p)中A元素的相应位置,L的对角线元素不存.D的元素di存放在a(1:p)中aii相应的位置.B~存放在b(1:n,1:m)的位置.
  对于i=1,2,…,n
  (1)计算I=d(i)-i;MI←mi=i-(d(i)-d(i-1))+1.
  (2)对于j=MI,MI+1,…,i
  J=d(j)-j;MJ←mj=j-(d(j)-d(j-1))+1;
  M←mij=max(MI,MJ)=max(mi,mj)
  a(I+j):=a(I+j)-∑j-1k=Ma(I+k)·a(d(k))·a(J+k)
  如果j=i,则转(3),否则a(I+j):=a(I+j)/a(d(j));
  b(i,k):=b(i,k)-a(I+j)·a(d(j))·b(j,k),k=1,2,…,m.
  (3)如果a(I+i)=0,则计算停止(这时A非正定)
  否则b(i,k):=b(i,k)/a(d(i)),k=1,2,…,m.
  3.求解方程组LTX=B~
  对于i=n,n-1,…,1
  (1)计算I=d(i)-i;
  MI←mi=i-(d(i)-d(i-1))+1.
  (2)对于j=MI,MI+1,…,i-1
  b(j,k):=b(j,k)-a(I+j)·b(i,k),k=1,2,…,m.
其他文献
高职院校的和谐与稳定是社会和谐与发展的重要组成部分,维护校园稳定的途径与措施有多方面,其中教学是构建和谐稳定校园的着力点。
目的观察阿托品和长托宁治疗有机磷中毒的临床疗效。方法将80例有机磷中毒患者按治疗药物不同分为对照组34例和观察组46例,对照组采用阿托品治疗,观察组采用长托宁治疗,比较
高职院校制度文化建设中要坚持就业导向,突出高职教育的特点。制度文化建设要坚持以树立高职特色的校风为核心,以解放教师思想、唤醒广大学生对制度的认同感为基本点,以制度保证
目的探讨、总结老年慢性肺源性心脏病伴心律失常的临床特点,为其治疗提供相关依据。方法对我科自2003年3月至2009年5月收治的60岁以上老年慢性肺源性心脏病伴心率失常80例的
目的探讨液基细胞学检测和阴道镜联合检查对宫颈病变的诊断价值。方法对1013例宫颈疾病患者进行液基薄层细胞检查并同时采用阴道镜检查及活组织病理检查,以宫颈活组织病理诊
本文探讨了近几十年对词汇搭配概念的不同理解与分类.搭配已由狭义的语义搭配扩展到广义的语篇中的词项同现.词汇搭配在语篇层面上的功能至少有两个:语篇衔接和消除歧义.词汇
目的探讨微型钢板治疗桡骨小头骨折的方法与疗效。方法回顾性分析我院采用微型钢板螺钉内固定治疗的桡骨小头骨折患者32例临床资料。结果疗效:优27例;良3例;中2例,骨折不愈合
种子处理是农业技术的重要措施之一,它是防治农作物苗期病虫害的最主要手段.根据这样的要求,过去任何一种种子的处理方法--拌种、浸种、闷种等是不可能达到的.因此种子包衣技
【摘要】 本文主要就题后反思我思考什么,从围绕容易出错之处展开反思,围绕查找缺漏展开反思,围绕解题的方法和规律展开反思,围绕能否对原题加以变换、引申展开反思四方面指导学生进行题后反思.  【关键词】 数学教学;题后反思;习惯培养;反思途径  不少学生感叹:课堂上的学习,经常是“一听就会,一做就错”,这种情况的出现可能有多方面原因,但我认为最主要的是同学只追求解题的结果,没有从整体上思考问题的习惯,
在简述融合企业文化打造高职特色校园文化必要性和可行性的基础上,重点论述融合企业文化,从创新高职精神文化、发展高职制度文化、丰富高职物质文化、规范高职行为文化等四个