原创力文档- 1、本文档共26页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
matlab第四章 数值计算1
4.1 解线性方程 MATLAB解线性方程一般不用先求逆阵，再做矩阵乘法的方法求解。而是采用类似数量除法的方法来描述一般的线性方程组的解。 ● 系数矩阵左乘未知数 Ax=B x=B\A ● 系数矩阵右乘未知数 xA=B x=B/A 系数矩阵为m行n列矩阵 (1) m=n 为方阵系统，可以求出一个精确的解。 (2) mn 过限定系统，可以求出一个最小方阵的解。 (3) mn 限定不足系统，可以求出一个至多含有m个非零元素的基解。 2.6.1 方阵系统 Eg: A= 2 3 4 Eg: A= 2 3 4 5 4 1 5 4 1 1 3 2 1 3 2 U= 1 B= 4 7 5 2 1 4 7 3 5 2 6 Ax=U x=A-1U=A\U Ax=B x=A-1B=A\B x= -0.8148 x= -1.4444 1.2222 -0.4074 1.6667 2.0000 -1.0000 2.3333 -0.5926 0.2222 1.8889 -0.2963 A*x A*x ans = 1 ans= 4 7 5 2 1 4 7 3 5 2 6 2.6.2 过限定系统 各种曲线拟合的过程中经常遇到过限定系统的线性方程组，MATLAB将两组数据保存在列矢量t与列矢量y中。 Eg: t=[0 0.3 0.8 1.1 1.6 2.3] y=[0.82 0.72 0.63 0.60 0.53 0.50] 认为这组数据大概在如下一个衰减的指数函数曲线上y(t)=C1+C2e-t y是两个矢量的线性组合，一个是常数C1，另一个是e-t，将数据代入函数y(t)=C1+C2e-t得到六个含有两个未知数的等式，将这些数据表示成一个六行两列的矩阵E，把未知参数看作一个矢量c的两个元素方程组为Ec=y。 解： t=[0 0.3 0.8 1.1 1.6 2.3]； y=[0.82 0.72 0.63 0.60 0.53 0.50]； E=[ones(size(t))exp(-t)] E =1.0000 1.0000 1.0000 0.7408 1.0000 0.4493 1.0000 0.3329 1.0000 0.2019 1.0000 0.1003 c=E\y c = 0.4760 0.3413 E*c ans= 0.8173 0.7288 0.6293 0.5896 0.5449 0.5102 4.2 曲线拟合 一．回归法曲线拟合 曲线拟合需要从一些离散的数据中推导出两者之间的数学解析关系，而数据
文档评论（0）