矩阵与数值分析上机实验报告.docxVIP

实验一设，分别编制从小到大和从大到小的顺序程序分别计算并指出两种方法计算结果的有效位数。1.程序代码从小到大（从大到小）顺序求解#include stdio.h#include math.h#define N 10000//S1000000时 #define N 1000000void main(){float j;float Sn=0.0;//double Sn=0.0for(j=2;j=N;j++) //从大到小的顺序时 for(j=N;j=2;j--){Sn=Sn+1000000.0/(j*j-1);}printf(Sn=%f,Sn);}2.实验结果Sn为 double类型时（精度16位）：从小到大输出结果是S10000=749900.004999S1000000=749999.000001从大到小输出结果是S10000=749900.005000S1000000=749999.000001Sn为float类型时（精度为7位）：从小到大输出结果是S10000=749849.13500S1000000=749849.125001从大到小输出结果是S10000=749899.989583S1000000=749998.9895833.实验分析1）当使用精度较高的数据类型计算时，从大到小与从小到大计算没有太大区别2）当使用精度较低时，计算的先后顺序对结果影响很大，当N=10000时，两个计算结果差为50.854538；当N=1000000时，两个计算结果差为149.874582.3）造成是因为从小到大计算过程中，大的数加小的数时，把小数的舍入结果为0，产生了大数吃小数的现象。舍入误差较大。而从大到小计算过程，先安排小的数参加运算，获得了较高的精度。4）两种运算的精度分别是：以double数据类型从大到小的计算结果作为精确值，分别计算float类型的有效位数从小到大计算：S10000的有效位数为：3位S1000000有效位数是：3位从大到小计算：S10000的有效数位为：7位S1000000有效位数是：7位实验二解线性方程组，分别利用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组，其中常向量为维随机生成的列向量，系数矩阵具有如下形式，其中为阶矩阵，为阶单位矩阵，迭代法计算停止的条件为：，给出时的不同迭代步数．（一）用Jacobi迭代法解线性方程组1.程序代码#include stdafx.h#include stdio.h#include stdlib.h#include math.h#include time.h#include stdlib.h//using namespace std;#define N 20void main(){double A[N*N][N*N]={0};double B[N*N]={0};double X1[N*N]={0};double X2[N*N]={0};double temp=0.0;double maxErr=0.0;double err;bool flag=1;inti;int j;//构造矩阵Afor(i=0;iN*N;i=i+N){for (j=i;ji+N;j++){A[j][j]=4;if (ji+N-1){A[j+1][j]=-1;A[j][j+1]=-1;}if(j+NN*N){A[j+N][j]=-1;A[j][j+N]=-1;}}}//构造常数项B矩阵for(i=0;iN*N;i++){B[i]=1+(int)(10.0*rand()/(RAND_MAX+1.0)); //0~10之间的随机数printf(%f\n,B[i]);}//构造迭代常数矩阵Bfor(i=0;iN*N;i++){B[i]=B[i]/A[i][i];}//构造迭代矩阵Afor(i=0;iN*N;i=i+N){for (j=i;ji+N;j++){if (ji+N-1){A[j+1][j]=-A[j+1][j]/A[j][j];A[j][j+1]=-A[j][j+1]/A[j][j];}if(j+NN*N){A[j+N][j]=-A[j+N][j]/A[j][j];A[j][j+N]=-A[j][j+N]/A[j][j];}A[j][j]=0;}}//计算int count=0;//记录迭代次数while(true){count++;//迭代X的值，如果flag为1，则把结果写入X2中；如果flag为0，把结果写入X1中if(flag==1){flag=0;//flag=0表示结果写入X2for(i=0;iN*N;i++){temp=0;for(j=0;jN*N;j++){temp+=A[i][j]*X1[j];}X2[i]=temp+B[i];}}else {flag=1;//表