[工学]东北大学数值分析实验报告.doc

数值分析设计实验实验报告 课题一 迭代格式的比较 问题提出 设方程f(x)=x- 3x –1=0 有三个实根 x=1.8793 , x=-0.34727 ,x=-1.53209现采用下面三种不同计算格式，求 f(x)=0的根 x 或x x = x = x = 二、要求 1、编制一个程序进行运算，最后打印出每种迭代格式的敛散情况； 2、用事后误差估计来控制迭代次数，并且打印出迭代的次数； 3、初始值的选取对迭代收敛有何影响； 4、分析迭代收敛和发散的原因。 三、目的和意义 1、通过实验进一步了解方程求根的算法； 2、认识选择计算格式的重要性； 3、掌握迭代算法和精度控制； 4、明确迭代收敛性与初值选取的关系。 四、程序设计流程图 五、源程序代码 #includestdio.h #includemath.h void main() { float x1,x2,x3,q,a,z,p,e=0.00001; x1=-1.0000;x2=-1.0000;x3=1.0000; int i,y=3; printf(0 %f %f %f\n,x1,x2,x3); q=x1-p;a=x2-p;z=x3-p; for(i=1;i=60;i++) { if(qeq(0-e)) goto a; else { p=x1; x1=(3*x1+1)/(x1*x1); printf(%d 1 %f\t,i,x1); q=x1-p; } a: if(aea(0-e)) goto z; else { p=x2; x2=(x2*x2*x2-1)/3; printf(%d 2 %f\t,i,x2); a=x2-p; } z: if(zez(0-e)) goto end; else { p=x3; x3=pow((3*x3+1),1.0/y); printf(%d 3 %f\n,i,x3); z=x3-p; } end:; } } 六。程序运行结果 七. 程序运行结果讨论和分析： 对于迭代格式一、二、三对于初值为-1.0000，-1.0000，1.0000分别迭代了37次，8次，10次，由此可知，简单迭代法的收敛性取决于迭代函数，以及初值x的选取，并且对初值的选取要求较高，需谨慎选取。 课题二 线性方程组的直接算法 一、问题提出 给出下列几个不同类型的线性方程组，请用适当算法计算其解。 设线性方程组 = x= ( 1, -1, 0, 1, 2, 0, 3, 1, -1, 2 ) Gsuss列主元消去法 #include math.h #include conio.h #include stdio.h #define MAX 100 typedef struct{ int row,col; float MAT[MAX][MAX]; float Solution[MAX]; }Matrix; void Gauss(Matrix *M); void MBack(Matrix *M); void MSave(Matrix *M); void MInput(Matrix *M); void MOutput(Matrix *M); void Solution(Matrix *M); void MSort(Matrix *M,int n); void main() { Matrix Mat; MInput(Mat); MSave(Mat); Gauss(Mat); MSave(Mat); if(Mat.row==Mat.col-1) { MBack(Mat); Solution(Mat); } printf(Press any key to halt...); getch(); } void MInput(Matrix *M) { int i,j; printf(输入行数:); scanf(%d,M-row); printf(输入列数:); scanf(%d,M-col); for(i=0;iM-row;i++){ printf(第%d行:,i+1); for(j=0;jM-col;j++){ scanf(%f,M-MAT[i][j]); } } for(i=0;iM-row;i++) M-Solution[i] = 0; } void MOutput(Matrix *M) { int i,j; printf(MATRIX:\n); for(i=0;iM-row;i++){ f