12 / 2 页下一页

查看: 9207|回复: 19

Newton法和Broyden法求解非线性方程组——求助！！！！！

[复制链接]

字体大小: 正常放大

大象装冰箱

123 主题	303 听众	1万积分

★丶 CUMT丶默默. ≮把心挂起≯

TA的每日心情

	奋斗 2025-1-3 11:37

签到天数: 386 天

[LV.9]以坛为家II

群组: 岩土力学与地下工程

群组: 2014国赛优秀论文解析

群组: 2015国赛优秀论文解析

群组: 2013年美赛优秀论文解

群组: 2017算法基础课程

电梯直达

1^#

发表于 2010-3-4 16:40 |只看该作者 |正序浏览

|招呼Ta 关注Ta

本帖最后由大象装冰箱于 2010-4-7 09:41 编辑

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

请教的问题.png

zan

转播0 淘帖0 分享0 收藏1 支持0 反对0 微信

我一直以为隐了身别人就找不到我了
没有用的
象我这样拉风的男人,无论躲在哪都好象漆黑中的萤火虫一样
那么的鲜明,那么的出众
我那忧郁的眼神、唏嘘的须根,还有牛B的判断力,已经将我深深的出卖了

使用道具举报

wujianjack2

21 主题	97 听众	3110 积分

TA的每日心情

	奋斗 2014-3-2 00:26

签到天数: 243 天

[LV.8]以坛为家I

20^#

发表于 2013-9-22 21:12 |只看该作者 |招呼Ta 关注Ta

尝试用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
在精度允许的情况下，结果还不唯一呢。