北航数值分析作业.docxVIP

数值分析作业一一、题目设有501×501的矩阵A = 其中;b=0.16,c=-0.064。矩阵A的特征值满足 试求：1. ，和的值；2. A的与最接近的特征值。3. A的（谱范数）条件数和行列式。二、算法设计方案1、矩阵A的压缩存储与检索方法由题目知，矩阵A为带状矩阵，为了节省内存空间，矩阵中的大量0元素不存储，计算过程中出现的0元素也不存储。根据带状矩阵的特点，将该矩阵旋转压缩存储，选择一个5行501列的二维数组C[5][501]存储即可，在C[5][501]中检索A[i][j]的方法为A[i][j] C[i-j+2][j]。2、所求值的计算方法最大特征值，最小特征值，按模最小特征值的计算方法?首先使用一次幂法运算计算得到矩阵A的按模最大的特征值，就是矩阵A的最大或最小特征值。对原矩阵A平移后再一次幂法运算，可以得到另一个按模最大特征值。比较和的符号，正数的是?，负数的是。对矩阵A进行一次反幂法运算，即可得到按模最小特征值。?A与最接近的特征值的计算方法?首先计算出k所对应的μk值，对原矩阵A进行一次移项，即（A-μkI），对新矩阵进行一次反幂法运算，即可得到一个按模最小特征值λi。则原矩阵A与μk值最接近的特征值λik=λi+μk。?A的（谱范数）条件数cond（A）2的计算方法A为对称矩阵，根据条件数的定义可知，，其中为矩阵A的按模最大特征值，为矩阵A的按模最小特征值，这两个值已经求出，代入即可。A的行列式det（A）的计算方法矩阵A的行列式可先对矩阵A进行LU分解后，det（A）等于U所有对角线上元素的乘积。三、程序源代码// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。//#include stdafx.h#include stdio.h#include math.h#include float.h#include stdlib.h#include time.h#include vectorusing namespace std;#define E 1.0e-12//全局变量，定义误差限#define N 501//全局变量，由矩阵的大小决定，后续数组定义都从1开始，0元素所在的位置空置int max2(int a,int b);int min2(int a,int b);int max3(int a,int b, int c);void creat(vector vectordouble d);void chuzhir(vectordouble a);void chuzhi1(vectordouble a,int j);double mifa(vectordouble u,vector vectordouble d,double p);double fanmifa(vectordouble u,vector vectordouble d,double p);void LU(vector vectordouble d);int _tmain(int argc, _TCHAR* argv[]){//double u[501]vectordouble u(501,0);double TZZ1,TZZ501,TZZs,TZZk,cond,det,a,b;//int i;//double C[5][501]vector vectordouble C; C.resize(5);for(int i = 0;i !=5;++i){C[i].resize(N);//使二位向量存储空间先全部赋值为0}creat(C); chuzhi1(u,1);a = mifa(u,C,0);//调用一次幂法先计算出一个特征值b=mifa(u,C,a);//调用一次幂法求解按模最小特征值if(a0){TZZ1=a;TZZ501=b;}else{TZZ501=a;TZZ1=b;}printf(******数值分析大作业一******\n);printf(傅海军 ZY1413205\n);printf(矩阵A最小的特征值TZZ1=%14.11e\n,TZZ1);printf(矩阵A最大的特征值TZZ501=%14.11e\n,TZZ501);//用反幂法求解TZZschuzhi1(u,1);TZZs=fanmifa(u,C,0);printf(矩阵A按模最小的特征值TZZs=%14.11e\n,TZZs);//对uk进行移项，反幂法求解最接近的特征值for (int k = 1; k 40; k++){a=TZZ1+(k*(TZZ501-TZZ1))/40;creat(C);chuzhi1(u,1);TZZk=fanmifa(u,C,a);printf(与μ%02d=%14.11e最接近的特征值TZZ%