数学建模社区-数学中国

标题: 会lingo的老师教教我，对求解非线性方程组编程。谢谢 [打印本页]

作者: li404541784 时间: 2012-7-15 22:53
标题: 会lingo的老师教教我，对求解非线性方程组编程。谢谢
老师教教我，用“集”和for语言对解非线性方程组编程，谢谢，如：
44IG$`8_OR)I)EHJGYA4~U2.jpg

谢谢。

作者: li404541784 时间: 2012-7-15 22:59
哪位会的老师，麻烦教教我。谢谢

作者: zhb 时间: 2012-7-16 09:11
坐等高手。

作者: zhb 时间: 2012-7-16 09:34
mathematica 不可以用么？

作者: darker50 时间: 2012-7-16 09:50
有能力的同学解决下呗，我给大家评分。

作者: zhb 时间: 2012-7-16 11:39

darker50 发表于 2012-7-16 09:50
有能力的同学解决下呗，我给大家评分。

表示用mathematica试过了。。但是据说条件满足不了，解不了，不知道lingo会有怎样的效果。。

作者: li404541784 时间: 2012-7-16 19:52
谢谢各位朋友来，我知道怎么编程了，如下：
sets:
r/1..5/:x;
endsets
min=(x(1)-1)+@sum(r(i)|i#ge#2

x(i-1)-x(i))^i);
@sum(r(i)|i#le#3:x(i)^i)=3*2^0.5+2;
x(2)-x(3)^2+x(4)=2*2^0.5-2;
@for(r

free(x));
@for(r

bnd(-5,x,5));
谢谢关心。

作者: li404541784 时间: 2012-7-16 19:54
系统自动反应

是：（

是： @

作者: heshuangping 时间: 2012-7-17 12:50
嗯不错。。

作者: 123yujingwei 时间: 2012-7-17 16:17
貌似解不出来吧，x5=5，呵呵……

作者: li404541784 时间: 2012-7-17 16:32

123yujingwei 发表于 2012-7-17 16:17
貌似解不出来吧，x5=5，呵呵……

能够解出来：

就这样。

OIACY`19N5$EG~RR(D7{)9D.jpg (117.61 KB, 下载次数: 44)

OIACY`19N5$EG~RR(D7{)9D.jpg

作者: asd3983133 时间: 2012-7-17 19:48
......这个碉堡了......好吧我表示有压力

作者: 雪凌寒霜 时间: 2012-7-17 20:46
model:
min=(x1-1)+(x1-x2)^2+(x2-x3)^3+(x3-x4)^4+(x4-x5)^5;
x1+x2^2+x3^3=3*@sqrt(2)+2;
x2-x3^2+x4=2*@sqrt(2)-2;
x1>-5;
x2>-5;
x3>-5;
x4>-5;
x5>-5;
x1<5;
x2<5;
x3<5;
x4<5;
x5<5;
@free(x1);@free(x2);@free(x3);@free(x4);@free(x5);
end
时间匆忙，也不知道是不是写对了，只能帮到你这些了，其实lingo的命令挺简单的，找个教程一学就会。加油，这是结果
  Local optimal solution found.
  Objective value:                         -0.1937640E-02
  Extended solver steps:                            0
  Total solver iterations:                            4

作者: 123yujingwei 时间: 2012-7-19 10:07

li404541784 发表于 2012-7-17 16:32
能够解出来：

额，单纯形法太强了

作者: wssl103050 时间: 2012-8-25 13:28
model:
sets:
r/1..5/:x;
endsets
min=(x(1)-1)+@sum(r(i)|i#ge#2

x(i-1)-x(i))^i);
@sum(r(i)|i#le#3:x(i)^i)=3*2^0.5+2;
x(2)-x(3)^2+x(4)=2*2^0.5-2;
@for(r(i)

bnd(-5,x,5));
end

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