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

计算数值方法实验报告

while((n2>=eps)&&(k>eps; cout<<\输入方程解数：\ cin>>size; a=new double[size];/*动态分配二维数组*/ s=new double*[size]; for(l=0;l>s[l][m]; cout<<\输入系数矩阵\ for(l=0;l>a[l]; cout<<\ jac(s,a,size,eps); return 0; }

11

计算数值方法实验报告

五、结果讨论和分析 分析： 雅克比迭代公式简单，每迭代一次只须计算一次矩阵和向量的乘法。在计算机运算时需要两组存储单元，以存放x(k)及x(k+1)。 实验地点

12

综合楼六层606室 指导教师 王峥 计算数值方法实验报告

学院名称 学生姓名 课程名称 计算机科学与技术 某某 数值计算方法 专业班级 实验日期 实验题目 计算机 2011-6-20 学号 成绩 1111111111 实验四 最小二乘法拟合多项式 一、课题名称 （1）了解矩阵特征值与特征向量问题解法，掌握幂法。 （2）加深对矩阵特征值与特征向量问题求解方法的认识，掌握算法。 （3）会进行误差分析。 二、目的和意义 学习使用最小二乘法拟合多项式 三、计算公式 幂法：由已知的非零向量x0和矩阵A的乘幂构造向量序列{xn}以计算矩阵A的按模最大特征值及其特征向量的方法，称为幂法。 迭代公式： ??y?Axkk?1???mk?max(yk),k?1,2,...?yk?xk?mk?? 结果可取 ??1?mk???yk或?1?xk ?1 四、结构程序设计 #include #include const int N(3); using namespace std; void matrixx(double A[N][N],double x[N],double y[N]) { //求y的值 for(int i=0;i

13

计算数值方法实验报告

double m(double y[N]) { //求m的值，即y的最大值 double max; for(int i=0;iy[i+1]?y[i]:y[i+1]; return max; } void main() { cout<<\矩阵:\cout<<\cout<<\double A[][N]={2,-1,0,-1,2,-1,0,-1,2}; //矩阵A（3*3） double x[N]={1,1,1}; //已知非零向量x0 double y[N]={0,0,0},u[N]={0,0,0},p[N]={0,0,0}; double e=1e-2,delta=1; int k=0; //计数器 while(delta>=e) { for(int q=0;q

14

计算数值方法实验报告

五、结果讨论和分析 分析： 幂法是一种求任意矩阵A的按模最大特征值及其对应特征向量的迭代算法。该方法的最大优点是计算简单，容易在计算机上实现，对稀疏矩阵较为适合，但有时收敛速度很慢。 实验地点

综合楼六层606室 指导教师 王峥 15