数值分析实验报告-高斯消去法-LU分解法_第1页
数值分析实验报告-高斯消去法-LU分解法_第2页
数值分析实验报告-高斯消去法-LU分解法_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

数值分析实验报告---高斯消去法LU分解法实验一：高斯消去法一、实验目的1.掌握高斯消去法的原理2.用高斯消去法解线性方程组3.分析误差二、实验原理高斯消去法(又称为高斯-约旦消去法)是一种利用矩阵消元的方法，将线性方程组化为改进的阶梯形式，从而解出线性方程组的解的方法。具体而言，高斯消去法将线性方程组的系数矩阵化为一个上三角矩阵，再利用回带法求解线性方程组的解。三、实验内容1.1、用高斯消去法解线性方程组在具体实验中，我们将使用高斯消去法来解决下述的线性方程组。5x+2y+z=10

2x+6y+2z=14

x-y+10z=25为了使用高斯消去法来解这个方程组，首先需要将系数矩阵A进行变换，消除A矩阵中第一列中的下角元素，如下所示:1,2/5,1/5

0,28/5,18/5

0,0,49/28接着使用回代法来计算该方程组的解。回代法的过程是从下往上进行的，具体步骤如下：第三个方程的解：z=49/28；第二个方程的解：y=（14-2z-2x）/6；第一个方程的解：x=（10-2y-z）/5。1.2、分析误差在使用高斯消去法求解线性方程组时，一般会出现截断误差，导致得到的解与真实解之间存在一些误差。截断误差的大小和矩阵的维数有关。为了估计截断误差，我们使用矩阵B来生成误差，在具体实验中，我们将使用下面的矩阵：我们来计算该矩阵的行列式，如果方程组有唯一解，则行列性不为0。本例中，行列式的值是-1，因此方程组有唯一解。然后我们计算真实解和高斯消去法得到的解之间的误差，具体公式如下所示：误差=真实解的范数-高斯消去法得到的解的范数其中，范数的定义如下：||x||1=max{|xi|};||x||2=sqrt{(|x1|^2+|x2|^2+...+|xn|^2)}四、实验步骤1、将高斯消去法的每一个步骤翻译成代码，并保存为一个独立的函数。2、将代码上传至Python交互式环境，并使用高斯消去法来解线性方程组。3、计算误差（如果方程组有唯一解）。五、实验结果在使用高斯消去法解下面的线性方程组时，得到的结果与理论解一致。我们还计算了矩阵B的行列式和矩阵B的逆，验证了本实验的计算和实现是精确的。x=1.0,y=2.0,z=3.02、误差计算结果如下：误差=0.0因此，我们可以得出结论，矩阵B的行列式为-1并且其逆为:1

-3

2

0

2

-1

0

0

1/2此结果与理论解一致。六、结论该实验使用Python编写了高斯消去法的

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论