数学建模社区-数学中国

标题: Newton法和Broyden法求解非线性方程组——求助！！！！！ [打印本页]

作者: 大象装冰箱 时间: 2010-3-4 16:40
标题: Newton法和Broyden法求解非线性方程组——求助！！！！！
本帖最后由大象装冰箱于 2010-4-7 09:41 编辑

下面的这个问题是同学问我的，以前上学时做过的，现在都忘干净了，以前的资料都没有了，想请教下matlab高手，帮忙编一下，万分感谢！

请教的问题.png

作者: 大象装冰箱 时间: 2010-3-4 16:44
请高手不吝赐教，哈哈

作者: bond 时间: 2010-3-4 17:52
我也刚学计算方法，还没学到啊。。。不好意思帮不到忙了。。。

作者: clanswer 时间: 2010-3-4 17:53
我是来做等高手的，呵呵

作者: mnpfc 时间: 2010-3-4 18:06
我也来等matlab高手，哈哈

作者: 吴建宏 时间: 2010-3-5 12:38
这个我研究研究

作者: madio 时间: 2010-3-5 16:44
本人水平有限，Matlab程序编写比较费劲，只好从网上找了一些代码，你修改一下应该就可以用！

Broyden.rar (653 Bytes, 下载次数: 47)

Newton.rar (583 Bytes, 下载次数: 46)

作者: 大象装冰箱 时间: 2010-3-9 09:08
回复 7# madio

谢谢了，madio，嘿嘿，

作者: 数学者 时间: 2010-3-11 23:04
呵呵~

作者: huanxinghua2009 时间: 2010-3-18 19:42
回复 7# madio

xiexie

作者: FrancisWu 时间: 2010-4-4 17:41
回复 7# madio

作者: kviaxlw 时间: 2011-8-24 10:15
O(∩_∩)O谢谢

作者: zhongdan 时间: 2011-8-26 16:46
下下来看看

作者: huangfenli 时间: 2011-8-26 17:29
看一下

作者: ;`﹎_D`.K゛ 时间: 2012-4-4 15:46
kankan

作者: 菁彩无限 时间: 2012-4-8 14:11
晕，刚才还是小学生呢，怎么又成学前班了？

作者: dougualian 时间: 2013-6-12 15:53
找到了，谢谢，，

作者: dougualian 时间: 2013-6-12 15:59
找到谢谢非常感谢

作者: yangyuhonglirui 时间: 2013-9-22 09:58
我也刚入门。。。。。。

作者: wujianjack2 时间: 2013-9-22 21:12
尝试用LINGO写了一下，分两种情况：
情况一：
MIN=X1+X2+X3;
PI=3.1415926535;
3*X1-@COS(X2*X3)-0.5=0;
X1^2-81*(X2+0.1)^2+@SIN(X3)+1.06=0;
@EXP(-X1*X2)+20*X3+(10*PI-3)/3=0;
@FREE(X1);@FREE(X2);@FREE(X3);

运行结果为：
  Global optimal solution found.
  Objective value:                         -0.2302872
  Objective bound:                         -0.2302872
  Infeasibilities:                            0.1365605E-08
  Extended solver steps:                         2
  Total solver iterations:                         276

                     Variable          Value       Reduced Cost
                           X1    0.4981447          0.000000
                           X2    -0.1996059          0.000000
                           X3    -0.5288260          0.000000
                           PI       3.141593             0.000000

                        Row Slack or Surplus    Dual Price
                           1    -0.2302872          -1.000000
                           2       0.000000             0.1584432
                           3       0.000000          -0.3083673
                           4    -0.1365605E-08 -0.6465823E-01
                           5       0.000000          -0.4753296E-01

情况二：
MAX=X1+X2+X3;
PI=3.1415926535;
3*X1-@COS(X2*X3)-0.5=0;
X1^2-81*(X2+0.1)^2+@SIN(X3)+1.06=0;
@EXP(-X1*X2)+20*X3+(10*PI-3)/3=0;
@FREE(X1);@FREE(X2);@FREE(X3);

运行结果如下：
  Global optimal solution found.
  Objective value:                         -0.2359878E-01
  Objective bound:                         -0.2359878E-01
  Infeasibilities:                               0.4289902E-12
  Extended solver steps:                         0
  Total solver iterations:                         77

                     Variable          Value       Reduced Cost
                           X1    0.5000000          0.000000
                           X2       0.000000          0.000000
                           X3    -0.5235988          0.000000
                           PI       3.141593             0.000000

                        Row Slack or Surplus    Dual Price
                           1    -0.2359878E-01       1.000000
                           2       0.000000                -0.1758114
                           3       0.000000                0.3544521
                           4       0.000000             -0.6335628E-01
                           5       0.000000                0.5274341E-01
在精度允许的情况下，结果还不唯一呢。

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)