- 1、本文档共51页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Matlab数值计算练习
引言
本章将花较大的篇幅讨论若干常见数值计算问题:线性分析、一元和多元函数分析、微积分、数据分析、以及常微分方程求解等。但与一般数值计算教科书不同,本章的讨论重点是:如何利用现有的世界顶级数值计算资源MATLAB。至于数学描述,本章将遵循“最低限度自封闭”的原则处理,以最简明的方式阐述理论数学、数值数学和MATLAB计算指令之间的内在联系及区别。
对于那些熟悉其他高级语言(如FORTRAN,Pascal,C++)的读者来说,通过本章,MATLAB卓越的数组处理能力、浩瀚而灵活的M函数指令、丰富而友善的图形显示指令将使他们体验到解题视野的豁然开朗,感受到摆脱烦琐编程后的眉眼舒展。
对于那些经过大学基本数学教程的读者来说,通过本章,MATLAB精良完善的计算指令,自然易读的程序将使他们感悟“教程”数学的基础地位和局限性,看到从“理想化”简单算例通向科学研究和工程设计实际问题的一条途径。
对于那些熟悉MATLAB基本指令的读者来说,通过本章,围绕基本数值问题展开的内容将使他们体会到各别指令的运用场合和内在关系,获得综合运用不同指令解决具体问题的思路和借鉴。
由于MATLAB的基本运算单元是数组,所以本章内容将从矩阵分析、线性代数的数值计算开始。然后再介绍函数零点、极值的求取,数值微积分,数理统计和分析,拟合和插值,Fourier分析,和一般常微分方程初值问题。本章的最后讨论稀疏矩阵的处理,因为这只有在大型问题中,才须特别处理。
从总体上讲,本章各节之间没有依从关系,即读者没有必要从头到尾系统阅读本章内容。读者完全可以根据需要阅读有关节次。除特别说明外,每节中的例题指令是独立完整的,因此读者可以很容易地在自己机器上实践。
LU分解和恰定方程组的解
LU分解、行列式和逆
(1)LU分解
(2)行列式和逆
恰定方程组的解
【*例5.2.2-1】“求逆”法和“左除”法解恰定方程的性能对比
(1)为对比这两种方法的性能,先用以下指令构造一个条件数很大的高阶恰定方程。
rand(state,12); %选定随机种子,目的是可重复产生随机阵A。
A=rand(100,100)+1.e8; %rand(100,100)生成(100×100)均匀分布随机矩阵。
%每个随机阵元素加的目的是使A阵条件数升高。
x=ones(100,1); %令解向量 x 为全1的100元列向量。
b=A*x; %为使 Ax=b 方程一致,用A和 x 生成 b 向量。
cond(A) %求A阵的条件数。
ans =
1.4426e+012
(2)“求逆”法解恰定方程的误差、残差、运算次数和所用时间
flops(0);tic %浮点运算计数器置0 ;启动计时器Stopwatch Timer
xi=inv(A)*b; % xi 是用“求逆”法解恰定方程所得的解。
ti=toc %关闭计时器,并显示解方程所用的时间。
ci=flops %“求逆”法解方程所用的运算次数
eri=norm(x-xi) %解向量 xi 与真解向量 x 的范-2误差。
rei=norm(A*xi-b)/norm(b) %方程的范-2相对残差
ti =
0.9300
ci =
2070322
eri =
3.0708e-004
rei =
6.6280e-007
(3)“左除”法解恰定方程的误差、残差、运算次数和所用时间
flops(0);tic;xd=A\b; %是用“左除”法解恰定方程所得的解。
td=toc,cd=flops,erd=norm(x-xd),red=norm(A*xd-b)/norm(b)
td =
0.2200
cd =
741872
erd =
3.2243e-004
red =
2.0095e-016
范数、条件数和方程解的精度
【*例5.2.3-1】Hilbert矩阵是著名的病态矩阵。MATLAB中有专门的Hilbert矩阵及其准确逆矩阵的生成函数。本例将对方程近似解和准确解进行比较。所谓n阶Hilbert矩阵的形式是: 。
N=[6 8 10 12 14]; %本例计算的矩阵阶数
for k=1:length(N)
n=N(k); %矩阵的阶
H=hilb(n); %产生n阶Hilbert矩阵
Hi=invhilb(n); %产生完全准确的n阶逆Hilbert矩阵
b=ones(n,1); %生成n阶全1向量
x_approx=H\b; %利用左除H求近似解
x_exact=Hi*b;
您可能关注的文档
- 光电报警器技巧.doc
- 灌注桩分项工程质量技术交底卡(灌注桩)技巧.doc
- 灌浆平洞支护施工方案.技巧.doc
- 灌溉,井房,输变电技巧.doc
- 冠词考点(2016广东专版)技巧.doc
- 贯通建辉酒店服务指南技巧.doc
- 管桩支护施工组织设计技巧.doc
- 管桩首件工程施工方案技巧.doc
- Mapmatrix教学学案.doc
- 中药药理学复(共100分)A卷案例.doc
- 《中国通史》文字稿第12集春秋争霸.docx
- java教程--类与对象-讲义课件(演讲稿).ppt
- Vue应用程序开发-(1).pptx
- 东北师大版社劳动实践与评价指导手册一年级上册主题二活动一寻找五彩的树叶课时课件.pptx
- 外研版英语四年级上册 Module 4 Unit 2 How much is it单元教学设计.docx
- 外研版英语四年级上册Module 4 单元整体教学设计.docx
- 6《上课之前》课件 鄂科技版 心理健康教育一年级.pptx
- 《1~5的认识》说课课件(共25张PPT)人教版一年级上册数学.pptx
- 六《解决问题(1)》说课课件 人教版 三年级上册数学.pptx
- 七《解决问题》说课课件 人教版 二年级上册数学.pptx
文档评论(0)