文档介绍：实验报告
课程名称 数值分析实验
实验项目 误差分析
专业班级 姓名—学号—
指导教师 成绩
一、 实验目的
了解误差分析对数值计算的重要性。
掌握避免或减小误差的基本方法。
二、 实验设备
安装有C、C++或MATLi]>>yi[i];
cin>>value_x;
void init() ，
int i,j;
for(i=0;i<n+1;i++)
} Si”"
for(i=0;i<n+1;i++)
for(j=0;j<n+1;j++)
fenmu[i][j]=xi[i]-xi[j] float sum=1,result=0;
for(int k=0;k<n+1;k++)
sum=1;
for(int j=0;j<n+1;j++)
( if(j==k) continue;
sum*=fenzi[j]/fenmu[k][j]
lumF
result+=sum;
//printf("result=%lf\n",result) cout<<"result="<<result<<endl;
'main()
getdata();
init();
result();
system("pause");
return 0;
佛山科学技术学院
实验报告
课程名称 数值分析实验
实验项目 高斯消去法解方程
专业班级 09计算机科学与技术2班 姓 名 张国烽 学 号 2009314215
指导教师 成绩 日期
一、实验目的
"一列选主元在计算机上的算法实现。
，程序运行时先输入一个数n表示方程含有的未知数个数，然后输入每个线性方 程的系数和常数，求出线性方程组的解。
二、实验设备
安装有C、C++或MATLAB的计算机。
三、实验原理
设有方程组Ax = b,设A是可逆矩阵。高斯消去法的基本思想就是将矩阵的初等行变换作用于方程组的增 广矩阵B = [A力]，将其中的A变换成一个上三角矩阵，然后求解这个三角形方程组。
将方程组用增广矩阵B= Ab]=（^ ） 表示。
j nx（ n+1）
步骤1：消元过程，对k = 1,2,...,n -1
选主元，找ik^{k, k +1,...,〃}使得
（2）
（3）
（4）
=max a
k <i < n lk
如果a, k= 0，则矩阵A奇异，程序结束；否则执行（3）。
如果广。k，则交换第k行与第i行对应元素位置，a — a
k k kj ikj
消元，对i = k，…,n，计算l没=a* /a^,对j = k + 1,・・・,n +1，
a = a -"a .
a
ik,k
步骤2：回代过程：
（1） 若a = 0,则矩阵奇异，程序结束；否则执行（2）。
nn r
a工
i ,n+1
V j=i+1
(2) x = a /a ;对i = n-1,...,2,1，
计算% =
i
j = k ,•,,,〃 +1。
计算
\
a x
ij j
/ a
ii
（1）
'2x1 + x2 - 3x3 = 1
4 x + x - 2 x = 4
3 x + x