计算数值实验报告(太原理工大学)

本科实验报告

课程名称：计算数值方法

实验地点：综合楼五层506室

专业班级：计科1002 学号： 2010001414 学生姓名： xxx

指导教师：王峥

2012 年 6 月 20

太原理工大学学生实验报告

2,)

内有零点，反之在（a，(a+b)/2

printf("%5d %5f %5f %5f %5f\n",i,a,b,x,fx);

i++;

if(fx==0) break;

if(fx>0) b=x;

else if(fx<0) a=x;

if((b-a)

break; //进行计算并返值

}while(1);

}

六、结果讨论和分析

二分法：

迭代法：

太原理工大学学生实验报告

}

六、结果讨论和分析

Gauss消元法：列主元素消元法：

LU分解法：追赶法：

太原理工大学学生实验报告

分析讨论：

其实，这两个迭代法是之前迭代法的升级，多了几个迭代式子而已，而且两者相差不大比较简单，所以选择了雅克比迭代法进行求解，但是没有与另一种方法高斯赛德尔迭代法进行实质性的比较。

流程图：

实验地点综合楼五层506室指导教师王峥

太原理工大学学生实验报告

六、结果讨论和分析

分析讨论

由于该程序将A矩阵和x0的值编写在程序中，所以要想修改成其他矩阵时比较麻烦，所以也有一定的局限性。幂法是一种求任意矩阵A的按模最大特征值及其对应特征向量的迭代算法。该方法的最大优点是计算简单，容易在计算机上实现，对稀疏矩阵较为适合，但有时收敛速度很慢。

流程图：

实验地点综合楼五层506室指导教师王峥

太原理工大学学生实验报告

for(int i=1;i<=n;i++)

{

M[i][j]=(M[i][j-1]-M[i-1][j-1])/(M[i][0]-M[i-j+1][0]);

}

}

for(int i=1;i<=n;i++)

{

cout<<"其"<

}

cout<<"请输入x的值：x=";

cin>>xx;

for(int i=0;i

{

X*=xx-M[i][0];

N+=M[i+1][i+2]*X;

P=M[0][1]+N;

}

cout<<"其函数值：y="<

}

六、结果讨论和分析

太原理工大学学生实验报告

分析讨论

数据拟合的具体作法是：对给定的数据（x i,y i）（i=0,1,…,m），在取定的函数类中，

求p(x)属于此函数类，使误差r

i =p(x

i

)- y

i

(i=0,1,…,m)的平方和最小，即∑r i2=∑

（∑p(x i)-y i）2=min从几何意义上讲，就是寻求与给定点（x i ,y i）（i=0,1,…,m）的距离平方和为最小的曲线y=p(x)。说明：上述程序中数据已被事先存入，可在源代码中任意修改想要的数据。

流程图：

实验地点综合楼五层506室指导教师王峥

（注：可编辑下载，若有不当之处，请指正，谢谢!）