数值分析报告1.docVIP

《数 值 分 析》 计 算 报 告 第 一 次 作 业 系 （院）： 宇航学院 专 业： 航空宇航推进理论与工程 学生姓名： 向纪鑫 学 号： SY1415313 2014年11月 一．计算方案的设计 1.初始化并压缩存储矩阵A 定义矩阵，由矩阵性质可知。根据计算实习说明书要求，所有矩阵的零元素都不存储，所以将矩阵压缩存储为矩阵。在矩阵c中检索矩阵的带内元素的方法为：。 2.求和的值 首先，用幂法求出矩阵按模最大的特征值。然后，根据原点平移法，同样根据幂法求出矩阵按模最大的特征值。当时，，；当时，，；当时，或，若，则，若，则。最后，综上所述，只要求出和，通过比较就可以求出和。 求,,的值 用反幂法就可以直接求出。结合上一步求出的，从而很容易求出。其中，用反幂法求的过程中用到了Doolittle分解，即把矩阵写成,是上三角矩阵，所以的值就是所有矩阵主对角线元素的乘积。 求的与最接近的特征值 用带平移的反幂法求出的值然后加上就是，即：求出矩阵按模最小的特征值，。 计算源程序 File1 初始化的函数（shushihua） #includestdio.h #includemath.h void chushihua(double A[5][501]) { int i; 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)); A[0][i]=A[4][i]=-0.064; A[1][i]=A[3][i]=0.16; } A[0][0]=A[0][1]=A[1][0]=A[4][500]=A[4][499]=A[3][500]=0; } File 2 幂法的函数 （mifa） #includestdio.h #includemath.h double mifa(double A[5][501]) { int i,j,t,k=0; double g,B=1,c=2,u[501],y[501]; for(i=0;i501;i++) { u[i]=1; } do { g=0; c=B; for(i=0;i501;i++) { g=u[i]*u[i]+g; } g=sqrt(g); for(i=0;i501;i++) { y[i]=u[i]/g; } for(i=0;i501;i++) { u[i]=0; for(j=0;j501;j++) { t=i-j+2; if(t=0t=4) { u[i]=A[t][j]*y[j]+u[i]; } } } B=0; for(i=0;i501;i++) { B=y[i]*u[i]+B; } k=k++; } while(fabs((B-c)/B)1e-12); return(B); } File 3 Doolittle分解 （doolittle） #includestdio.h #includemath.h void doolittle(double A[5][501]) { int k,i,j,t,ik; double S[501],Max,B=1; for(k=1;k=501;k++) //*用选主元Doolittle分解法解方程组（LU分解）*// { ik=k; S[k-1]=0; for(t=1;t=k-1;t++) if(fabs(k-t)3) S[k-1]=A[k-t+2][t-1]*A[t-k+2][k-1]+S[k-1]; S[k-1]=A[2][k-1]-S[k-1]; Max=S[k-1]; for(i=k;i=501;i++) { S[i-1]=0; for(t=1;t=k-1;t++) if(fabs(i-t)3fabs(t-k)3) S[i-