北航数值分析大作业一.doc

《北航数值分析大作业一.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、word格式北京航空航天大学数值分析大作业一学院名称自动化专业方向控制工程学号ZY1403140学生姓名许阳教师孙玉泉日期2014年11月26日....word格式设有的实对称矩阵A，其中，。矩阵A的特征值为，并且有1.求，和的值。2.求A的与数最接近的特征值。3.求A的(谱范数)条件数和行列式detA。一方案设计1求，和的值。为按模最小特征值，。可使用反幂法求得。，分别为最大特征值及最小特征值。可使用幂法求出按模最大特征值，如结果为正，即为，结果为负，则为。使用位移的方式求得另一特征值即可。2求A的与数最接近的特征值。题目可看成求以为偏移量后，按模最小的特征值。即以为偏移量做位移，使用反幂法

2、求出按模最小特征值后，加上，即为所求。3求A的(谱范数)条件数和行列式detA。矩阵A为非奇异对称矩阵，可知，(1-1)....word格式其中为按模最大特征值，为按模最小特征值。detA可由LU分解得到。因LU均为三角阵，则其主对角线乘积即为A的行列式。二算法实现1幂法使用如下迭代格式：(2-1)终止迭代的控制理论使用，实际使用(2-2)由于不保存A矩阵中的零元素，只保存主对角元素a[501]及b,c值。则上式中简化为：(2-3)2反幂法使用如下迭代格式：(2-4)其中，解方程求出。....word格式求解过程中使用LU分解，由于A为5对角矩阵，选择追赶法求取LU分解。求解过程如下：追赶法求

3、LU分解的实现：(2-5)由上式推出分解公式如下：(2-6)推导出回代求解公式如下：(2-7)....word格式(2-8)3及A行列式求解(2-9)由式(2-5)可得：(2-10)三源程序#include#includedoubleep=1e-12,b=0.16,c=-0.064;intj=0;doublepower(doublea[501]);//幂法doubleinv_power(doublea[501]);//反幂法doubledet(doublea[501]);//求detintmain()//主程序{inti,k;doubleA[501],B[50

4、1],beta_1,beta_501,beta_s,beta_k;doublemu;for(i=0;i<501;i++)A[i]=(1.64-0.024*(i+1))*sin(0.2*(i+1))-0.64*exp(0.1/(i+1));beta_1=power(A);//第一问printf("λ1t=%.12et迭代次数：%d",beta_1,j);for(i=0;i<501;i++)//位移B[i]=A[i]-beta_1;beta_501=power(B)+beta_1;printf("λ501t=%.12et迭代次数：%d",beta_501,j);beta_s=inv

5、_power(A);printf("λst=%.12et迭代次数：%d",beta_s,j);for(k=1;k<=39;k++)//第二问{mu=beta_1+k*(beta_501-beta_1)/40;....word格式for(i=0;i<501;i++)B[i]=A[i]-mu;beta_k=inv_power(B)+mu;printf("λi%dt=%.12et迭代次数：%d",k,beta_k,j);}printf("cond(A)2=%.12e",beta_1/beta_s);//第三问printf("detAt=%.12e",det(A));}dou

6、blepower(doublea[501])//幂法{inti=0,N=5000;doubleb=0.16,c=-0.064;doubleu[501],y[501];doublem=1,beta;for(i=0;i<501;i++)u[i]=1;j=0;while(j

7、498]+b*y[499]+a[500]*y[500];for(i=2;i<499;i++){u[i]=c*y[i-2]+b*y[i-1]+a[i]*y[i]+b*y[i+1]+c*y[i+2];}beta=0;for(i=0;i<501;i++){if(fabs(u[i])>=fabs(beta))beta=u[i];}if(beta<0)if(fabs(fabs(beta)-fabs(m))/