北航数值分析实习第一题重点分析.docxVIP

数值分析计算实习报告第一题所在班级A1班学生姓名学生学号2015年11月1 算法设计方案1.1矩阵的压缩存储A矩阵是一个带宽为5的带状矩阵，众多的0会占据很多存储空间，因此可以将其压缩存储。存储方法是将矩阵逆时针旋转45度，并调整使列号不变，调整后的矩阵如下：1.2求解最大、最小的特征值及模最小的特征值1.利用幂法计算出矩阵A按模最大的特征值（它是和中的一个）。2.利用幂法计算出矩阵（其中I为单位矩阵）按模最大的特征值，令（它是和中的另一个）。3.比较和的大小，大的是，小的是。4.利用反幂法直接计算出矩阵A按模最小的特征值。1.3求与某值接近的特征值采用原点平移法即可方便求解与某值接近的特征值。对于每个k值，先利用反幂法计算出矩阵按模最小的特征值，令即可。1.4求解条件书和行列式值对于非奇异实对称矩阵A，它的条件数，其中是按模最大的特征值，是按模最小的特征值。行列式的值为上三角矩阵对角线元素的乘积。过在进行反幂法求解中要进行Doolittle分解，Doolittle分解可以得到一个上三角矩阵，其对角线元素的乘积就是行列式值。2 C++程序#includestdio.h#includemath.hvoidcaculateA(doublea[5][501]);//计算出A矩阵，并且将其压缩存储doublemifa2(doublea[5][501]);//用2范数的幂法求按模最大特征值doublefanmifa2(doublea[5][501]);//用2范数的反幂法求按模最小特征值voidludoolittle(doublea[5][501]);//简易存储的矩阵的doolittle分解doublefindlamnear(doublex);//求x处的特征值intmax3(intx,inty,intz);//求三个数的最大值intmin2(intx,inty);//求两个数的最小值doublefanshu2(doublea[]);//向量的2范数doubleinner(doublex[],doubley[]);//向量乘积voidmulti(doublea[5][501],doublex[],doubley[]);//矩阵与向量相乘voidsolve(doublea[5][501],doubley[],doublex[]);//根据doolittle分解的结果，求解方程解向量xvoidmain(){ #defineerr1e-12 printf(数值分析计算实习第一次作业\n); printf(班\n); doublea[5][501]; double lam1,lam501,lamm,lams,temp,cond,det,x; inti; caculateA(a); lam1=mifa2(a); lamm=lam1; for(i=0;i501;i++) a[2][i]=a[2][i]-lam1; lam501=mifa2(a)+lam1; if(lam1lam501) { temp=lam1; lam1=lam501; lam501=temp; } caculateA(a); ludoolittle(a); det=1.0; for(i=0;i501;i++) { det=det*a[2][i]; } lams=fanmifa2(a); cond=fabs(lamm/lams); printf(最小特征值λ1=%.11e\n,lam1); printf(最大特征值λ501=%.11e\n,lam501); printf(按模最小特征值λs=%.11e\n,lams); for(i=0;i39;i++) { temp=lam1+(i+1)*(lam501-lam1)/40; x=findlamnear(temp); printf(λi%d=%.11e\t,i+1,temp); if(i%2==1) printf(\n); } printf(\n); printf(条件数cond(A)=%.11e\n,cond); printf(行列式值det(A)=%.11e\n,det);}voidcaculateA(doublea[5][501]){ doubleb=0.16; doublec=-0.064; inti; for(i=2;i501;i++) a[0][i]=c; for(i=1;i501;i++) a[1][i]=b; for(i=0;i501;i++) a[2][i]=(1.64-0.024*(i+1))*sin(0.2*(i+1))-0.64*exp(0.1/(i+1)); for(i=0;i500;i++) a[3][i]=b; for(i=0;i501-2;i++) a[4][i