《数值及分析》的上机实习报告.docVIP

数值分析上机实验报告 （C语言） 姓 名： 班 级：2006级 学 号：2006050349 院 系：机械工程学院 指导教师： 第一题：Househloder变换和线性代数方程组求解 一、程序要求 已知A和b矩阵。 用Household变换，把A化为三对角阵B(并打印B)。 用超松弛法求解Bx=b(取松弛因子ω=1.4,X(0)=0,迭代9次)。 用消去法求解Bx=b。 二、解题算法 1．Household变换 令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+Sign(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．超松驰法 3．消去法 , , , , q1[0]=0 , u1[0]=0 , x[9]=u1[9] 三．源程序 #include math.h #include stdio.h #define N 9 float sign(float x) /*子函数符号函数*/ { if (x0.0) return 1.0; else if (x0.0) return -1.0; else return 0.0; } void housholder(float a[N][N]) /*子函数Household 变换 */ { int i,j,n; float s,a1,k,u[N],y[N],q[N]; for(j=0;jN-2;j++) { s=0.0; for(i=j+1;iN;i++) s=s+pow(a[i][j],2); s=sqrt(s); a1=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*sign(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]/a1; } for(k=0,i=0;iN;i++) k=k+u[i]*y[i]; k=0.5*k/a1; 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(Household变换的结果:); /*输出Household变换的结果*/ for(i=0;iN;++i) for(j=0;jN;++j) { if (j%N==0) printf(\n); printf(%-8.4f,a[i][j]); } printf(\n); } void guss(float a[N][N],float b[N]) /* 子函数高斯消去法*/ { int i,j; float u[N],y[N],q[N],x[80]; u[0]=a[0][0]; y[0]=b[0]; for(i=1;iN;i++) { q[i]=a[i][i-1]/u[i-1]; u[i]=a[i][i]-q[i]*a[i-1][i]; y[i]=b[i]-q[i]*y[i-1]; } x[N-1]=y[N-1]/u[N-1]; for(i=N-1;i=0;i--) x[i]=(y[i]-a[i][i+1]*x[i+1])/u[i]; printf(消去法的结果：\n);