《数值分析》上机实习报告5则范文.docx

《数值分析》上机实习报告5则范文

第一篇：《数值分析》上机实习报告数值分析上机实验报告姓名：班级：学号：院系：机械工程学院指导教师：C语言）2006级2006050349（for(j=0;j}voidguss(floata[N][N],floatb[N])/*子函数高斯消去法*/{inti,j;{}printf(“Household变换的结果:”);/*输出Household变换的结果*/for(i=0;ifor(j=0;j}if(j%N==0)printf(“n”);s=0.0;for(i=j+1;i}for(i=0;i}for(k=0,i=0;iq[i]=y[i]-k*u[i];for(i=0;i}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”);/*输出消去法的结果*/for(i=0;i}printf(“n”);}voidsor(floata[N][N],floatb[N])/*子函数超松驰法*/{floata1[N][N];inti,j,k,m;floatx[N],temp[N][N];floatw=1.4,h=0,g=0;for(m=0;m{for(j=1;jq[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];if(i==5)printf(“n”);printf(“x%d=%9.6f”,i,x[i]);}g=(h+b[i-1])*w;x[i-1]=(1-w)*x[i-1]+g;h=0.0;g=0.0;}}printf(“超松弛法的结果：n”);/*输出超松弛法的结果*/for(i=0;i四、运行结果Household变换的结果:12.3841-4.89310.00000.00000.00000.00000.00000.00000.0000-4.893125.39846.49410.0000-0.00000.0000-0.0000-0.0000-0.00000.00006.494120.61158.2439-0.0000-0.0000-0.00000.00000.00000.00000.00008.243923.4229-13.8801-0.00000.0000-0.00000.00000.0000-0.0000-0.0000-13.880129.69834.53450.00000.00000.00000.00000.0000-0.0000-0.00004.534516.00614.88140.0000-0.00000.0000-0.0000-0.00000.00000.00004.881426.0133-4.50360.00000.0000-0.00000.0000-0.00000.00000.0000-4.503621.25404.50450.0000-0.00000.00000.00000.0000-0.00000.00004.504514.5341消去法的结果：x0=1.075799x1=2.275744x2=-2.855514x3=2.293097x4=2.112633x5=-6.423833x6=1.357923x7=0.634245x8=-0.587266超松弛法的结果：x0=1.073409x1=2.272580x2=-2.856600x3=2.292511x4=2.112164x5=-6.422582x6=1.357802x7=0.634259x8=-0.587042五、问题讨论：1．算法分析SOR方法的矩阵形式为：X(m)=(E-ωL)-1((1-ω)E+ωR)x(m-1)+(E-ωL)-1ωg若记Lω=(E-ωL)-1((1-ω)E+ωR)，SOR收敛的充要条件是S(Lω)数组作为子函数参数时，实现的是地址传递，会改变原数组中元素的值，因此主函数调用housholder变换子函数后，原数组a[9][9]中的值已经改变为三对角阵。以下的消去法、超松弛法子函数调用的是已经三对角化后的a[9][9]的值。调用时由于没有定义housholder变换子函数中数组a[][]的大小，运行时出现以下错误提示：errorC2087:'':missingsubscripterrorC2133:'a':unknownsize修改大小为a[N][N]时(N为宏定义常数9)，运行正常。Xn1,Zkmkmax(Xn)sa[0]=a[0][0]-x;if(sa[0]==0)sa[1]=a[0][1]*a[0][1];elsesa[1]=a[1][1]-x-a[0][1]*a[0][1]/sa[0