- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
淮 阴 工 学 院
《数值分析》考试
──基于Matlab的方法综合应用报告
班 级: 金融1121
姓 名: 蒋倩
学 号: 1124104119
成 绩:
数 理 学 院
2014年6月7日
《数值分析》课程考试
2014.6.3
一 课程主要教学内容综述(开卷部分)(70分)
综述要求:
1. 综述内容以本课程实际讲授内容和所使用的教材为基准;
2. 按主题论述,问题提出的背景、动因,解决问题的理论基础和方法;
3. 以数值计算的实例说明方法的应用,并进一步对方法本身和应用结果进行简要的分析和评价;
4. 各数值计算的Matlab程序按顺序编号,附在正文的后面。
*正文用五号字,附件编号,用小五号字,至少五页,可以讨论,要用自己的语言
*参考文献:书(白峰杉,数值计算引论(第二版),北京:高等教育出版社,2010)、电子讲义(章节名,电子讲义,淮阴工学院,2014)
*自己再找一到两个参考文献
*格式:正文、参考文献、附录
正 文
一、线性方程组求解的数值方法
1、高斯消去法:
1.1问题提出的背景、动因[1]:
高斯消去法是求解线性代数方程组的直接方法,假设计算中没有舍入误差,经过有限次算数运算能够给出问题的精确解。
1.2解决问题的理论基础和方法[1,2]:
高斯消去法的基本思想就是将一般的线性方程组化成与之等价的上三角或下三角形式来求解。高斯消去法的消元过程,从代数运算的角度看就是用一个下三角矩阵左乘方程组的系数矩阵,且乘积的结果为上三角矩阵,即
LU 分解的MATLAB实现:
或
1.3例:Ax=b,已知A = [1 -2 0.33 0;2 -4.01 -1 5;-1 3 1 4; -3 6 0 2],b = [-1 1 0 3],用高斯消去法求解x
(程序见附录一)
结果:
X =
-3.70
-1.47
-0.72
0.36
L =
1.00 0 0 0
0.33 1.00 0 0
-0.67 -0.01 1.00 0
-0.33 0 -0.33 1.00
U =
-3.00 6.00 0 2.00
0 1.00 1.00 3.33
0 0 -0.99 6.37
0 0 0 2.79
P =
0 0 0 1.00
0 0 1.00 0
0 1.00 0 0
1.00 0 0 0
x =
-3.70
-1.47
-0.72
0.36
1.4分析和评价:
用高斯消去法(LU分解)和用Matlab求解的结果相同,所以这种方法较为适用。
2、迭代法
2.1问题提出的背景、动因[3]:
通常逆矩阵不易求得,特别是对于大型的线性方程组,需要用迭代法求解。
2.2解决问题的理论基础和方法[3]:
用迭代法求解线性方程组,首先要把线性方程组写成等价的形式
(1)
由迭代格式(1)确定如下的迭代算法:
Jacobian 迭代法:
Gauss-Seidel 迭代法:
2.3例:生成一个矩阵A和一个向量b,用Gauss迭代法和Jacobi迭代法求解方程组X0 =
0.86
0.1
文档评论(0)