2012-13年数值分析试题（参考）.docVIP

福 建 农 林 大 学 研 究 生 课 程 考 试 卷 授课时间 2012—2013 学年度 第 1 学期 学 号： 1121213004 研究生姓名： 吴翠兰 课 程 名 称： 数 值 分 析 考 试 时 间： 2012-12-25 考 生 成 绩： 授课或主考教师： （签章） 2012-13学年(1)福建农林大学机电学院研究生《数值分析》课程试卷 姓名：吴翠兰 班级：2012级 学号：1121213004 成绩： 1、用对分区间法求方程在区间[1,2]上的一个实根。精度要求为。（10分） 【解】/* 头文件包含区*/ #include math.h #include stdio.h #include stdlib.h /* 定义声明区*/ int N=0,n=0; double a=0,b=0,e=0,x=0; void step_1(void); void step_2(void); void step_3(void); void step_4(void); void step_5(void); void step_6(void); void step_7(void); /* f(x)函数 */ double f(double X) { double result; //输入f(x) result=X+log(X)-2.2; return(result); } /* 子函数 */ void step_1(void) { double f(double X); a=1; b=2; e=0.000001; N=1000; } void step_2(void) { n=0; } void step_3(void) { x=(a+b)/2; } void step_4(void) { if(fabs(f(x))0.000001||((b-a)/2)e){printf(求解成功!\n);printf(方程的根为\nx=%.5f ,x);;while(1);} else step_5(); } void step_5(void) { if(nN)step_6(); else{printf(求解失败...\n);while(1);} } void step_6(void) { if((f(a)*f(x))0)b=x; else a=x; } void step_7(void) { n=n+1;step_3(); } /* 主函数 */ void main() { step_1(); step_2(); while(1) { step_3(); step_4(); step_7(); } } 调试结果： 2、通过对系数矩阵的三角分解求解如下线性方程组。（10分） 【解】 /* 头文件包含区*/ #include math.h #include stdio.h #include stdlib.h #define N 3//矩阵阶数：3 /* 定义声明区*/ int i=0,j=0,k=0,r=0; float sum1=0,sum2=0,sum3=0,sum4=0; float A[N+1][N+1]={0}; float L[N+1][N+1]={0}; float U[N+1][N+1]={0}; float b[N+1]={0}; float y[N+1]={0}; float x[N+1]={0}; /* 主函数 */ void main() { /*矩阵A[i][j]*/ /*矩阵b[k]*/ A[1][1]=2 ;A[1][2]=-1 ;A[1][3]=-1;b[1]=4 ; A[2][1]=3 ;A[2][2]=4 ;A[2][3]=-2;b[2]=11; A[3][1]=3 ;A[3][2]=-2 ;A[3][3]=4 ;b[3]=11; //求矩阵U 矩阵L for(j=1;j=N;j++) { U[1][j]=A[1][j]; } for(i=2;i=N;i++) { L[i][1]=A[i][1]/U[1][1]; } for(r=2;r=N;r++) { for(j=r;j=N;j++) { for(k=1;k=(r-1);k++) { sum1=sum1+L[r][k]*U[k][j]; }