数值代数实验报告1

暨南大学本科实验报告专用纸

课程名称 数值代数 成绩评定 实验项目名称 (填写实验所属章节名称) 指导教师 刘娟 实验项目编号 07010029901 实验项目类型 验证 实验地点 机房

学生姓名 王伟文 学号 2010051727 学院 信息科学技术学院数学 系 信息与计算科学 专业 2010级 实验时间 2011年 9月 1 日~12月30日 温度24℃

第一章

1.实验选题:(写出你在该项目所选的实验题目)

第一章上机练习1

2.谈谈你对该算法的理解:(简单谈一下你是如何理解该算法的?)

对算法的理解:先将84阶的矩阵A分解为一个下三角矩阵L和上三角矩阵U,先考虑下三角形方程组Ly=b,利用前代法解出y,再考虑上三角形方程组Ux=y;利用回代法解出x。 列主元gauss消去法

Ax b PA LU,Ly Pb,Ux y;

3.实验内容(将实验程序及其实验结果粘贴,最好对程序各部分注释清楚,比如设置了哪

些函数,这些函数的输入输出是什么,具有什么功能?) 实验程序

对矩阵A进行LU分解的程序

function [ L,U ] = LUfac( A ) for k=1:n-1

A(k+1:n,k)=A(k+1:n,k)/A(k,k);

A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n); end

L=tril(A,0); for i=1:n L(i,i)=1; end U=triu(A,0); end

利用前代法解出y值的程序

function [ b ] = TSL( L,b )

你可能喜欢

  • 数值线性代数
  • 线性代数答案
  • 线性代数方程组求解
  • 高斯列主元消去法
  • 数据库期末考试题及答案
  • 考研学校排名
  • 牛顿法解非线性方程组
  • 矩阵特征值

数值代数实验报告1相关文档

最新文档

返回顶部