资源描述：
《数值分析计算实习作业一.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、数值分析计算实习题一学号：姓名：院系：2015年11月5日数值分析计算实习作业一一、分析1.1算法分析题目要求求出：1）特征值从小到大排列的最小特征值和最大特征值。2）特征值中模最小的特征值。3）靠近一组数的一组特征值。4）矩阵A的条件数cond(A)2。5）行列式detA。解决方法：1）若将所有行列式按模的大小排列则模最大的特征值一定是和中的一个，因此利用幂法求出模最大的特征值。然后利用带原点平移的幂法，将系数矩阵变为即将所有特征值都减去，则特征值按大小顺序排列的次序不变，模最大的特征值依然在整个排列的两端，再用一次幂法得到模最
2、大的特征值，其中为带原点平移的幂法求出的特征值，最后两个特征值、比较大小，大的为，小的为。2）因为为按模最小的特征值，因此用反幂法可求的其特征值。3）因为靠近数，因此一定是所有的中模最小的，因此可利用带原点平移的反幂法求出特征值，此时的系数矩阵变为。4）条件数cond(A)2为模最小的特征值与模最大的特征值的比的绝对值，因此利用1和2中求出的和可解出条件数。5）可对矩阵A进行LU分解，即则，又因为矩阵L对角线元素为1，则=1，所以，U为上三角阵，行列式为对角线元素的乘积，因此可得A的行列式。1.2程序分析10数值分析计算实习作业一
3、1.2.1因为A为拟三角阵，储存时零元素不储存，因此将矩阵A压缩为5*501的矩阵CA的带内元素=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,
4、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）如果两次相
5、邻两次迭代的满足：（允许误差），则结束迭10数值分析计算实习作业一代，的值可认为是矩阵B对应的模最大的特征值。如果不满足误差条件，则重复3）的迭代直到达到误差允许值。1.2.5拟三角矩阵的LU分解过程可参照数值分析书26页的算法得到。二、程序整个数值分析在c语言中的程序如下：#include#includedoublemifa();doublefanmifa();voidfenjie();//声明三个函数，mifa是用幂法求模最大的特征值，fanmifa是用反幂法求模最小的特征值，fenjie是对
6、系数矩阵做LU分解doubleB[5][501]={0};doubleA[5][501]={0};//定义两个系数矩阵为全局变量,B作为求解过程中的系数矩阵使用doublesigma=1e-12;//定义误差允许值/********主程序********/voidmain(){inti,j,t;doublelamk,lamm,min,max,miu,lam,x,lams;doublecond2;doubledet=1;for(i=0;i<=500;i++)//输入矩阵A{j=i+1;x=0.1/j;A[2][i]=(1.64-0.
7、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;10数值分析计算实习作业一for(i=0;i<=4;i++)//先使B=A{for(j=0;j<=500;j++){B[i][j]=A[i][j];}}lamk=mifa();//幂法求模最大的特征值for
8、(i=0;i<=4;i++)//赋值计算B=（A-lamk*I）{for(j=0;j<=500;j++){if(i==2)B[i][j]=A[i][j]-lamk;elseB[i][j]=A[i][j];}}lamm=mifa()+lamk;//带
