数值分析第一题.docVIP

第一题 Househloder变换与方程组求解 一、程序要求 已知A 1. 用Householder变换，把A化为三对角阵B(并打印B)。 2. 用超松弛法求解Bx=b(取松弛因子ω=1.4(0)=0,迭代9次)。 3. 用消去法求解Bx=b。 二、解题算法 1．Householder变换 参考《数值分析》课的P147页，其将Householeer变化归纳起来，得到对称矩阵三对角化的算法步骤如下： （１） 令A0=A, aij(1)=aij,已知Ar-1　即Ar-1=(aij(r))； （２） Sr=(∑ni=r+1(air (r)) 2)1/2 ; （３） ar=Sr2+|a(r)r+1,r|Sr; ur=[0,…0,a(r)r+1,r+sing(a(r)r+1,r)Sr, a(r)r+2,r, …,a(r)nr]T （４） yr=Ar-1ur/ar ; （５） kr=0.5*urTyr/ar ; （６） qr=yr-krur ; （７） Ar=Ar-1-(qrurT+urqrT) ; r=1,2….n-2 . 2．超松驰法 参考《数值分析》课的P127页，松弛法进行计算。 其基本思想是迭代，即在方法已求出x(m),x(m-1)的基础上，组合新的序列，从而加快收敛速度。其算式是： xi(m)=(1-ω)xi(m-1)+ω( bijxi(m)+ xj(m-1)+gi) 其中ω是超松弛因子，当ω1时，可以加快收敛速度。 消去中间量 ，就得SOR方法的算式 3．列主元素消去法 列主元素消去法其是对Gauss消去法的改进，其主要思想是保持a[i][i]a[i][j]，ji。其保持了作除数的值是最大的，减少了误差。 三．程序设计 #include math.h #include stdio.h #define N 9 float sing(double z) /*子函数sing函数*/ { if (z0.0) return 1.0; else if (z0.0) return -1.0; else return 0.0; } void housholder(double a[N][N]) /*子函数Householder变换 */ { int i,j,n; double s,b1,k,u[N],y[N],q[N]; for(j=0;jN-2;j++) { s=0; for(i=j+1;iN;i++) s=s+pow(a[i][j],2); s=sqrt(s); b1=s*s+s*fabs(a[j+1][j]); for(n=0;nN;n++) { if (n=j) u[n]=0; else if (n==j+1) u[n]=a[j+1][j]+s*sing(a[j+1][j]); else u[n]=a[n][j]; } for(i=0;iN;i++) { y[i]=0; for(n=0;nN;n++) y[i]=y[i]+a[i][n]*u[n]; y[i]=y[i]/b1; } for(k=0,i=0;iN;i++) k=k+u[i]*y[i]; k=0.5*k/b1; for(i=0;iN;i++) q[i]=y[i]-k*u[i]; for(n=0;nN;n++) for(i=0;iN;i++) a[i][n]=a[i][n]-(q[i]*u[n]+u[i]*q[n]); } printf(Householder bian huan de dui jiao zhen B:); /*输出Householder变换的结果*/ for(i=0;iN;++i) for(j=0;jN;++j) { if (j%N==0) printf(\n);