北航数值分析实习题目第一题.docVIP

《数值分析B》大作业一 ZY1515105 ?樊雪松? 一．? 算法设计方案：? 1.矩阵A的存储与检索? 将带状线性矩阵A[501][501]转存为一个矩阵MatrixC[5][501]?。在数组MatrixC[5][501]中检索A的带内元素aij的方法是：A的带内元素aij=C中的元素ci-j+2,j? 。 2.求解λ1，λ501，λs? 1、首先分别使用幂法和反幂法迭代求出矩阵按摸最大和最小的特征值λmax和λmin。λmin即为λs；如果λ max0,则λ501=λmax；如果λmax0,则λ1=λmax。? 2、使用带原点平移的幂法（mifa（）函数），令平移量p=λmax，求出 对应的按摸最大的特征值λ’max，? 如果λ max0,则λ1=λ’max+p；如果λmax0,则λ501=λ’max+p。 3、求解A的与数μk=λ1+k（λ501-λ1）/40 的最接近的特征值λik （k=1,2，…，39）。? ????? 使用带原点平移的反幂法，令平移量p=μk，即可求出与μk最接近的特征值λ ik。? 4、求解A的（谱范数）条件数cond（A）2和行列式detA。 cond（A）2=|λ1/λn|，其中λ1和λn分别是矩阵 A的模最大和最 小特征值。求解矩阵A的行列式，可先对矩阵A进行LU分解后，detA等于U所有对角线上元素的乘积。 二．?源程序? #includestdio.h #includemath.h #includeconio.h //定义A中元素 double C[5][501]; double a[501]; double b; double c; //声明所有函数 void YaSuoJZ(double C[5][501],double a[501],double b,double c) ;//压缩矩阵函数 double mifa(double C[5][501]); //幂法函数 void daizhuangLU(double A[5][501]); //带状矩阵的LU分解 double fanmifa(double C[5][501]);//反幂法函数 //最值函数 int max2(int x,int y); int max3(int x,int y,int z); int min(int x,int y); //最值函数 int max2(int x,int y) //求2个数的最大值 { int z; z=xy?x:y; return(z); } int max3(int x,int y,int z) //求3个数的最大值 { int w; w = z max2(x,y)? z:max2(x,y); return(w); } int min(int x,int y) //求2个数的最小值 { int z; z=xy?y:x; return(z); } //将矩阵A压缩存储在矩阵C中 void YaSuoJZ(double C[5][501],double a[501],double b,double c) { int i; for(i=0;i=500;i++) { if(i=2) C[0][i]=c; else C[0][i]=0; if(i=1) C[1][i]=b; else C[1][i]=0; if(i=499) C[3][i]=b; else C[3][i]=0; if(i=498) C[4][i]=c; else C[4][i]=0; C[2][i]=a[i]; } } //幂法函数：用幂法求矩阵模最大的特征值 double mifa(double C[5][501]) { double u[501]; double y[501]={0},η=0; double β,βk=0; double ε=1; // ε为精度 double sumu=0,sumAY=0; int i,j,k=1; //k为循环次数 for (i=0;i=500;i++) //取任一非零向量u0 u[i] = 1.0; while(ε=1e-12) { for(i=0;i=500;i++) //求u(k-1)的2范数η sumu=sumu+u[i]*u[i]; η=sqrt(sumu); sumu=0; for(i=0