数值分析计算实习作业一.docVIP

数 值 分 析 计算实习题一 学号： 姓名： 院系： 2015年11月5日 数值分析计算实习作业一 PAGE \* MERGEFORMAT9 一、分析 1.1算法分析 题目要求求出： 1）特征值从小到大排列的最小特征值和最大特征值。 2）特征值中模最小的特征值。 3）靠近一组数的一组特征值。 4）矩阵A的条件数cond(A)2。 5）行列式detA。 解决方法： 1）若将所有行列式按模的大小排列则模最大的特征值一定是和中的一个，因此利用幂法求出模最大的特征值。然后利用带原点平移的幂法，将系数矩阵变为即将所有特征值都减去，则特征值按大小顺序排列的次序不变，模最大的特征值依然在整个排列的两端，再用一次幂法得到模最大的特征值，其中为带原点平移的幂法求出的特征值，最后两个特征值、比较大小，大的为，小的为。 2）因为为按模最小的特征值，因此用反幂法可求的其特征值。 3）因为靠近数，因此一定是所有的中模最小的，因此可利用带原点平移的反幂法求出特征值，此时的系数矩阵变为。 4）条件数cond(A)2为模最小的特征值与模最大的特征值的比的绝对值，因此利用1和2中求出的和可解出条件数。 5）可对矩阵A进行LU分解，即则，又因为矩阵L对角线元素为1，则=1，所以，U为上三角阵，行列式为对角线元素的乘积，因此可得A的行列式。 1.2程序分析 1.2.1 因为A为拟三角阵，储存时零元素不储存，因此将矩阵A压缩为5*501的矩阵C A的带内元素=C中的元素 程序中A[5][501]即为压缩后的矩阵。 1.1.2 程序中的B[5][501]为过渡矩阵，在幂法迭代、反幂法迭代以及LU分解中均用矩阵B来计算，计算之间对B进行适当的赋值。先令B=A，调用幂法函数求；再令B=A-I，调用幂法函数求出；再令B=A，调用反幂法函数求；在令B=，调用反幂法函数求，最后令B=A，将B进行LU分解，计算A的行列式。 1.2.3幂法求解过程： 1）取初始向量u0=[1,1,…1]; 2）进行k次迭代，k=1,2,3… 3）对每一次迭代，计算上一次迭代中的的模（2范数），将单位化后赋值给，即，计算矩阵B与向量的乘积，，计算与的内积。 4）如果两次相邻两次迭代的满足：（允许误差），则结束迭代，的值可认为是矩阵B对应的模最大的特征值。如果不满足误差条件，则重复3）的迭代直到达到误差允许值。 1.2.4 反幂法求解过程： 1）取初始向量u0=[1,1,…1];对矩阵B进行LU分解。 2）进行k次迭代，k=1,2,3… 3）对每一次迭代，计算上一次迭代中的的模（2范数），将单位化后赋值给，即，利用LU分解后的矩阵B，求线性方程组，计算与的内积。 4）如果两次相邻两次迭代的满足：（允许误差），则结束迭代，的值可认为是矩阵B对应的模最大的特征值。如果不满足误差条件，则重复3）的迭代直到达到误差允许值。 1.2.5 拟三角矩阵的LU分解过程可参照数值分析书26页的算法得到。 二、程序 整个数值分析在c语言中的程序如下： #include stdio.h #include math.h double mifa(); double fanmifa(); void fenjie();//声明三个函数，mifa是用幂法求模最大的特征值，fanmifa是用反幂法求模最小的特征值，fenjie是对系数矩阵做LU分解 double B[5][501]={0}; double A[5][501]={0};//定义两个系数矩阵为全局变量,B作为求解过程中的系数矩阵使用 double sigma=1e-12;//定义误差允许值 /********主程序********/ void main() { int i,j,t; double lamk,lamm,min,max,miu,lam,x,lams; double cond2; double det=1; for(i=0;i=500;i++)//输入矩阵A { j=i+1; x=0.1/j; A[2][i]=(1.64-0.024*j)*sin(0.2*j)-0.64*exp(x); A[3][i]=0.16; A[4][i]=-0.064; } for(i=1;i=500;i++) { A[1][i]=0.16; } for(i=2;i=500;i++) { A[0][i]=-0.064; } A[3][500]=0; A[4][500]=0; A[4][499]=0; for(i=0;i=4;i++)//先使B=A { for(j=0;j=500;j++) { B[i][j]=A[i][j]; } } lam