原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
举例，给定方程组为：
先用matlab自带函数solve解此方程组，确定牛顿迭代时的初值范围，得到根为：
作图验证：
此组值确为方程的根。
通过观察我们可以发现y的取值必须大于0。这在程序中必须说明，如果迭代过程中y小于0，则此迭代法发散。
误差分析：因为范数等价的原因，我们选择2范数。将两次相邻迭代差的2范数作为误差，存储与一个向量或矩阵中，并作出曲线图，观察迭代过程中误差的变化情况。
如选初值为(12，0.3)，得到误差图形为：
选初值为(12，1.2)，误差图为：
我们可以发现误差在前3-5次迭代的时候迅速下降，但是中间会有上升的过程，直至最后误差达到我们设定的误差值。由此猜想迭代过程可能漏掉了一些根，利用作图，得到曲线如下：
可以发现还有两组根，用牛顿迭代法只能得到一组值，可能是因为所给方程比较特殊，它的定义域中x,y均不能为0，导致函数不连续。
另外，也可能因为函数不连续，导致初值只能在根的附件变化时才能得到收敛的结果。
因此我们不妨将初值选的稍微小一些，如:[1,2]，得到根为
选初值为[-1,4]，得到根为
综上所述，此方程组有三组根，不同的初值，会得到不同的解的情况，也会有不同的误差情况。初值选得越接近真值，误差变化程度越小，迭代次数也越少。
程序在另三个附件中。
同理，对于n元函数方程组，我们有：
误差分析类似二元函数方程组。
文档评论（0）