查看: 3850|回复: 5

Forcal全局优化库FcOpt更新

[复制链接]

字体大小: 正常放大

forcal

45 主题	3 听众	282 积分

升级 91%

TA的每日心情

	难过 2012-8-27 18:22

签到天数: 1 天

[LV.1]初来乍到

电梯直达

1^#

发表于 2011-2-8 12:33 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

本帖最后由 forcal 于 2011-2-8 12:37 编辑

Forcal优化库FcOpt最近进行了更新，全局优化能力显著增强，使用更加方便，同时增加了隐函数全部解，非线性方程（组）全部解求解函数，更新了微分方程参数优化（拟合）实例。
参考：http://www.forcal.net/sysm/forcal9/fchtm/fcopt.htm

zan

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

使用道具举报

musuchan

1 主题	6 听众	40 积分

升级 36.84%

该用户从未签到

2^#

发表于 2011-2-8 12:35 |只看该作者 |招呼Ta 关注Ta

谢谢斑竹~~~~~~~~~~

我也说一句

发表

使用道具举报

杨帆

4 主题	4 听众	626 积分

升级 6.5%

TA的每日心情

	奋斗 2012-12-4 13:38

签到天数: 314 天

[LV.8]以坛为家I

群组: 数学趣味、游戏、IQ等

群组: 数学建模

群组: 我们一定会赢

3^#

发表于 2011-2-8 15:59 |只看该作者 |招呼Ta 关注Ta

谢谢楼主，楼主威武

我也说一句

发表

使用道具举报

forcal

45 主题	3 听众	282 积分

升级 91%

TA的每日心情

	难过 2012-8-27 18:22

签到天数: 1 天

[LV.1]初来乍到

4^#

发表于 2011-2-20 19:09 |只看该作者 |招呼Ta 关注Ta

谢谢楼上两位朋友！

使用Forcal优化库FcOpt做优化时，无论简单的优化，还是复杂的优化（例如含积分的优化、共享参数优化、复数优化、微分方程参数优化等等），基本步骤都是一样的：

1、给出拟合数据，一般保存在数组中。如果没有数据，这一步就省了。
2、编写目标函数，做优化时就是要得到目标函数值最小时的最优参数值。自己写目标函数，似乎是麻烦了些，但非常灵活，前面提到的各种各样的优化都可以通过自己写目标函数来实现。正所谓以不变应万变也。
3、选择优化函数进行优化。

在FcOpt中，我最喜欢用的优化函数是Opt，下面给出一个使用Opt做优化的例子：

例子：拟合公式：y = (p1)+(p2*exp(-p3*x/p5)+p4/(1+p4*p5*x))

p1,p2,p3,p4,p5为待求参数

数据（x， y）
0,             0.928
0.0000098, 1.02
0.0000195, 1.12
0.0000293, 1.25
0.0000391, 1.42
0.0000488, 1.7
0.0000586, 2.01
0.0000684, 2.26
0.0000781, 2.46
0.0000879, 2.63
0.0000977, 2.82
0.0001074, 3.01
0.0001172, 3.2
0.000127,  3.41
0.0001367, 3.59
0.0001465, 3.72
0.0001562, 3.85
0.000166,  3.98
0.0001758, 4.08

Forcal代码：

!using["fcopt","math"];9 e& G& X/ x2 U) }
init(::Array,max)= //准备数据
) |\\" Z6 m Y; h& k, F) S* X
{2 Q5 Y$ B1 K$ n9 O1 U
max=19,
2 D# P; W' u/ M, M9 f: m+ O+ d# G
Array=arrayinit{2,max,2 :
+ ~8 H: U- |+ Q- P1 c- K! } k
0, 0.928,
$ w3 g' C; l6 C& H7 ^' a) V
0.0000098, 1.02,
0 t; T* U8 T6 D6 `. k; {6 S1 D. V- F# @
0.0000195, 1.12,( u6 ~7 s3 Y- I' l# y8 B8 F: r
0.0000293, 1.25,& w\\" Q\\" [1 F5 h& B2 R0 X
0.0000391, 1.42,
\\" x; G\\" B1 _6 f$ l
0.0000488, 1.7,
1 F% @$ u1 Q/ t( n& E& i9 z+ ]) z
0.0000586, 2.01,
: X* `0 h8 ]: y/ X* O
0.0000684, 2.26,
{3 ?2 O/ M5 Q. ?# S: P, l
0.0000781, 2.46,
4 F% C8 Q4 S2 T; o0 f K7 G
0.0000879, 2.63,* f, R. Q\\" L% z\\" J% D7 n\\" y+ V
0.0000977, 2.82,
! C! @; N3 o9 O1 X$ J9 [7 `' ?
0.0001074, 3.01,
% v& E7 J. X% t\\" `
0.0001172, 3.2,3 q% K- A* R; A& C1 Y0 }9 m
0.000127, 3.41,4 G3 S: z0 s/ m/ b+ m
0.0001367, 3.59,
0 I0 H' o5 q% Z a3 m
0.0001465, 3.72,8 @9 O# e* A2 k; R
0.0001562, 3.85,
& i/ U; ^0 l: i8 l\\" Y; h2 y$ u
0.000166, 3.98,: \9 ?9 k, H3 q: D' `3 ?' Z0 O
0.0001758, 4.08$ q e! \8 Q4 X2 d) {! R' r* p5 ?
}.free()7 t+ n) }4 }4 W+ e! V) U
};! e3 {! h2 j3 }) P7 |
f(p1,p2,p3,p4,p5 : i,s,x,y : Array,max)= //目标函数( u# y\\" {, t' L8 d+ k+ j# { J+ Z/ c
{
( S5 v' ?2 B, d! V, Y3 y
s=0,i=0,(i<max).while{
/ `+ T7 p4 M! X5 y3 |
x=Array[i,0], y=Array[i,1],
( E' f& ?/ b! ]( c
s=s+[(p1)+(p2*exp(-p3*x/p5)+p4/(1+p4*p5*x))-y]^2,
: |% C, N6 U2 E3 i* r1 }& t
i++
\\" V- i, z$ [+ a p2 r3 S
},
4 T3 T! i\\" R G) b. g4 l4 ?
sqrt[s/max]# B' g% [* J; K1 ~0 L8 x; a
};) u# I' x6 s- `0 a
Opt[HFor("f"), optwaysimdeep, optwayconfra, optdeep,50]; //比较难的优化，我就这样使用Opt函数

!using["fcopt","math"];9 e& G& X/ x2 U) }

init(::Array,max)= //准备数据 
) |\\" Z6 m  Y; h& k, F) S* X
{2 Q5 Y$ B1 K$ n9 O1 U

max=19, 
2 D# P; W' u/ M, M9 f: m+ O+ d# G
    Array=arrayinit{2,max,2 : 
+ ~8 H: U- |+ Q- P1 c- K! }  k
        0,         0.928, 
$ w3 g' C; l6 C& H7 ^' a) V
        0.0000098, 1.02, 
0 t; T* U8 T6 D6 `. k; {6 S1 D. V- F# @
        0.0000195, 1.12,( u6 ~7 s3 Y- I' l# y8 B8 F: r

0.0000293, 1.25,& w\\" Q\\" [1 F5 h& B2 R0 X

0.0000391, 1.42, 
\\" x; G\\" B1 _6 f$ l
        0.0000488, 1.7, 
1 F% @$ u1 Q/ t( n& E& i9 z+ ]) z
        0.0000586, 2.01, 
: X* `0 h8 ]: y/ X* O
        0.0000684, 2.26, 
  {3 ?2 O/ M5 Q. ?# S: P, l
        0.0000781, 2.46, 
4 F% C8 Q4 S2 T; o0 f  K7 G
        0.0000879, 2.63,* f, R. Q\\" L% z\\" J% D7 n\\" y+ V

0.0000977, 2.82, 
! C! @; N3 o9 O1 X$ J9 [7 `' ?
        0.0001074, 3.01, 
% v& E7 J. X% t\\" `
        0.0001172, 3.2,3 q% K- A* R; A& C1 Y0 }9 m

0.000127,  3.41,4 G3 S: z0 s/ m/ b+ m

0.0001367, 3.59, 
0 I0 H' o5 q% Z  a3 m
        0.0001465, 3.72,8 @9 O# e* A2 k; R

0.0001562, 3.85, 
& i/ U; ^0 l: i8 l\\" Y; h2 y$ u
        0.000166,  3.98,: \9 ?9 k, H3 q: D' `3 ?' Z0 O

0.0001758, 4.08$ q  e! \8 Q4 X2 d) {! R' r* p5 ?

}.free()7 t+ n) }4 }4 W+ e! V) U

};! e3 {! h2 j3 }) P7 |

f(p1,p2,p3,p4,p5 : i,s,x,y : Array,max)= //目标函数( u# y\\" {, t' L8 d+ k+ j# {  J+ Z/ c

{ 
( S5 v' ?2 B, d! V, Y3 y
    s=0,i=0,(i<max).while{ 
/ `+ T7 p4 M! X5 y3 |
        x=Array[i,0], y=Array[i,1], 
( E' f& ?/ b! ]( c
        s=s+[(p1)+(p2*exp(-p3*x/p5)+p4/(1+p4*p5*x))-y]^2, 
: |% C, N6 U2 E3 i* r1 }& t
        i++ 
\\" V- i, z$ [+ a  p2 r3 S
    }, 
4 T3 T! i\\" R  G) b. g4 l4 ?
    sqrt[s/max]# B' g% [* J; K1 ~0 L8 x; a

};) u# I' x6 s- `0 a

Opt[HFor("f"), optwaysimdeep, optwayconfra, optdeep,50]; //比较难的优化，我就这样使用Opt函数

结果（需求解几次，前面的数是最优参数，最后一个数是目标函数值）：

6.855486033798995 4.813449005068769 -54298054.6899329 -10.72898361613976 -1516.502649185066 3.337716353170544e-002

感觉做这个优化题，是稍稍有些难度的，换句话说，如果题目难度不超过此题，Forcal可以轻松求解。

喜欢优化的朋友可以尝试一下。

目前我知道1stOpt可以轻松求解此题，还有一位优化高手也可以轻松求解此题，还想知道更多可以求解此题的朋友或软件。

我也说一句

发表

使用道具举报

forcal

45 主题	3 听众	282 积分

升级 91%

TA的每日心情

	难过 2012-8-27 18:22

签到天数: 1 天

[LV.1]初来乍到

5^#

发表于 2011-2-23 21:41 |只看该作者 |招呼Ta 关注Ta

本帖最后由 forcal 于 2011-2-25 19:50 编辑

美国国家标准与技术研究院（NIST: National Institute of Standards and Technology）提供有一套27道非线性拟合测试题，世界上几乎所有著名的数据分析软件包都以能通过该套测试题集为验证标准。Forcal全部能以任意随机初始值求得最优解。

参考：http://www.forcal.net/yyhz/optnist.htm