23矩阵的条件数与病态方程组讲解课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《23矩阵的条件数与病态方程组讲解课件.ppt》由用户(ziliao2023)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 23 矩阵 条件 病态 方程组 讲解 课件
- 资源描述:
-
1、2.3 矩阵的条件数与病态方程组矩阵的条件数与病态方程组一、矩阵的条件数一、矩阵的条件数 二、线性方程组的性态二、线性方程组的性态三、病态线性方程组的求解三、病态线性方程组的求解2.3 矩阵的条件数与病态方程组矩阵的条件数与病态方程组例例1 方程组方程组0001.4410001.31321xx准确解:准确解:TTxx1,1,21若若A及及b作微小变化,考虑扰动后的方程组:作微小变化,考虑扰动后的方程组:0002.4419999.21321xx准确解:准确解:TTxx10,2,21方程组解的几何解释为:平面上两条接近于平行的直线的方程组解的几何解释为:平面上两条接近于平行的直线的交点,当其中一条
2、直线稍有变化时,新的交点与原交点相交点,当其中一条直线稍有变化时,新的交点与原交点相差很远。差很远。例例2 方程组方程组3133233210957910685657787104321xxxx准确解为:准确解为:TTxxxx)1,1,1,1(),(4321(1)对右端)对右端b作微小扰动:作微小扰动:9.301.339.221.3210957910685657787104321xxxxTTxxxx)1.1,5.4,6.12,2.9(),(4321(2)对系数矩阵)对系数矩阵A作微小扰动:作微小扰动:3133233298.9999.499.6989.998.585604.508.72.71.871
3、04321xxxxTTxxxx)22,34,137,81(),(4321AAxxbbxx4488倍倍15111倍倍02.0001.001.0011.002.000004.008.02.01.000A0.30.9%330.10.3 0 3%3 3%136001136%136016.132.3 矩阵的条件数与病态方程组矩阵的条件数与病态方程组一、矩阵的条件数一、矩阵的条件数 矩阵条件数的定义矩阵条件数的定义矩阵条件数的性质矩阵条件数的性质一、矩阵的条件数一、矩阵的条件数)(111bbAAAAAAAAxx改写(改写(2.22)式:)式:Proof 矩阵条件数的定义:矩阵条件数的定义:矩阵条件数的性质
4、:矩阵条件数的性质:(6)Cond(AB)Cond(A)Cond(B)二、线性方程组的性态二、线性方程组的性态答案:答案:TxAcond)0,2(,104)()1(4Txx)1,1()2(%50%,005.0)3(xxbb希尔伯特(Hilbert)阵njijihhHjinnji,2,1,11,)(,1212111111413121131211nnnnnnHn定义:定义:-最著名的病态矩阵最著名的病态矩阵对称正定矩阵对称正定矩阵在在MATLAB中,函数中,函数hilb()提供了()提供了Hilbert矩阵矩阵)4(4hilbH 10876441053.1)(1049.1)(,1055.1)(Hc
5、ondHcondHcond希尔伯特(Hilbert)阵-最著名的病态矩阵最著名的病态矩阵Hilbert矩阵的条件数:18201715131055109084.1)(10488.8)(106025.1)(,10761.4)(HcondHcondHcondHcond求解病态方程组出现的问题:求解病态方程组出现的问题:例:用例:用MATLAB求解线性方程组求解线性方程组bxHnTneeHb)1,1,1(,输入:输入:,;);1,(*);(;5xbHxnonesHbnhilbHn得:得:ans=1.000,1,000,1.000,1.000,1.000输入:输入:,;);1,(*);(;10 xbHx
6、nonesHbnhilbHn得:得:ans=1.000,1,000,1.000,1.000,0.9999 1.0002,0.9996,1.0004,0.9998,1.000输入:输入:,;);1,(*);(;5xbHxnonesHbnhilbHn得:得:ans=1.000,1,000,1.000,1.000,1.000输入:输入:得:得:ans=1.000,1,000,1.000,1.000,0.999 1.000,0.999,1.000,0.999,1.000,;);1,(*);(;10 xbHxnonesHbnhilbHn输入:输入:得:得:ans=1.000,1,000,1.001,0.
展开阅读全文