计算实习大作业二.docVIP

计算实习大作业算法设计方案： 按照题目给出的矩阵定义对矩阵赋初值：的特殊性，使用特殊数组进行存储(以后的计算中都使用此数组)；对应的函数为_init()； 与： 首先使用幂法迭代求出矩阵的按模最大的特征值：若，则其必为；于是可知，接下来使用幂法迭代求出矩阵的按模最大的特征值：由矩阵理论可知的特征值按从小到大排列应为、、……、。因此的按模最大的特征值应为。又因为的值已求得，由此可求出。 若，则其必为；于是可知，同(1)，可求出。 使用反幂法迭代计算矩阵按模最小的特征值； 计算矩阵的条件数和行列式，其中，由分解后的矩阵求出。 使用带原点平移的反幂法计算与最接近的特征值共39个(程序中所有迭代的最大次数设为5000)； 输出。 /***************************************** * 计算实习大作业第二题源程序 * 最后修改：2008.04.13 *****************************************/ #include math.h #include stdio.h #define N 501 //矩阵的阶数； #define EPSL 1.0e-12 //精度水平； #define L 5000 //迭代最大次数； #define R 2 #define S 2 #define OUTPUTMODE 1 //输出格式：0--输出至屏幕，1--输出至文件； #define K 39 double u0[N] = {0}, c[R+S + 1][N] = {0}, c_lu[R+S +1][N] = {0}, lmd_k[39] = {0}; double val_b = 0.16, val_c = -0.064, lmd_max, lmd_min, lmd_s, cond, det = 1; int main() { int i; double tmp; FILE *p; void val_init(); double max_lambda(); double min_lambda(); double k_lambda(double); void result_out(int); if(OUTPUTMODE) { p = fopen(c:\\Result.txt, w+); fprintf(p, 计算结果：\n); fclose(p); } val_init(); //对矩阵A进行初始化； if((tmp = max_lambda()) = 0) //若按模最大的特征值大于0，则其必为最大特征值； { lmd_max = tmp; for(i = 0; i N; i++) c[S][i]-= lmd_max; lmd_min = max_lambda() + lmd_max; //求出最小特征值； for(i = 0; i N; i++) c[S][i]+= lmd_max; } else //若按模最大的特征值小于0，则其必为最小特征值； { lmd_min = tmp; for(i = 0; i N; i++) c[S][i]-= lmd_min; lmd_max = max_lambda()+ lmd_min; //求出最大特征值； for(i = 0; i N; i++) c[S][i]+= lmd_min; } lmd_s = min_lambda(); //求出按模最小的特征值； cond = fabs(tmp / lmd_min); //求出矩阵的条件数； for(i = 0; i N; i++) det*= c_lu[S][i]; //求出矩阵的行列式； for(i = 0; i K; i++)