QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 14543|回复: 6
打印 上一主题 下一主题

四参数拟合方程 的方法

[复制链接]
字体大小: 正常 放大

1

主题

0

听众

3

积分

升级  60%

该用户从未签到

跳转到指定楼层
1#
发表于 2011-4-18 16:47 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.
3 D+ S2 p- R/ Y4 t6 l/ \& N1 r1 P! ~不知道有没有人知道如何求解.
) c5 L) e* n6 ^& G0 h$ K- h+ ]! X1 H7 z& l
方程如下:(四参数为A B C D): [; a. a' a& Q: ~3 \/ ?# \
Y = (A - D) / (1 + (X / C) ^ B)) + D
  o# O; S2 X; ?2 @3 X- I
4 H) G$ z0 x9 \: f, a给出数据集:& |% W0 x6 s" M, c6 Q) v
  x y
9 J/ T" _, X; @( I, U- m0.00 1.4735+ d' Q/ q# y. ^" c% D, Q# G$ z7 W& V
0.025 1.214- [8 ~* S+ d! P
0.05 1.0143 3 X9 [1 o" d$ p# _4 `# W
0.1 0.7362
' [2 F' D2 P- x/ E- e0.25 0.3907
  K3 x6 K5 g! l0.75 0.2135
% r6 F6 b! w+ ~8 `- [  O; K6 \9 ~( e' G1 p( X* S3 }7 X( W5 }
求最合适的四参数ABCD值是多少?* S7 N, }6 \4 a, E+ P% D5 ]

5 a. N2 C+ L7 H5 a5 [& A5 H* Z答案是 a = 1.47 b = 1.229 c = 0.085 d = 0.12
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
forcal 实名认证       

45

主题

3

听众

282

积分

升级  91%

  • TA的每日心情
    难过
    2012-8-27 18:22
  • 签到天数: 1 天

    [LV.1]初来乍到

    用Forcal求解,有2组解:& N8 G+ u2 H9 [% Q% B
    1. / [  {) H' R% {6 @1 Y
    2. !using["fcopt","math","sys"];
    3. / G4 o\\" r+ s* a4 Y7 G7 e\\" g1 f
    4. init(::Array,max)=
    5. ( C- x# `9 D% a# I3 U
    6. {
    7. 4 X4 H( J( W: i6 C5 [: B& c6 O
    8.   max=6,
    9. 5 E# a7 ?1 u9 C
    10.   Array=arrayinitns{max,2 :9 [& g- i* z: `3 [! K; Q9 u
    11. "
    12. 0 \, Y: Y- G\\" {/ V* T! U4 S* z
    13. 0.00  1.4735# w1 v* ~' y6 ^! T8 D
    14. 0.025 1.214
    15. ) j7 {/ A& F; m& `
    16. 0.05  1.0143 : W* v0 r0 j0 ~
    17. 0.1   0.7362
    18. 8 p1 Z) M7 s5 q0 b+ M) H* u  p
    19. 0.25  0.3907 $ q$ N- U- i/ f# _  W4 }
    20. 0.75  0.2135& Q+ J, Y2 v( L\\" q+ F- p4 }2 F
    21. ", X# a* ?% g* C! G
    22.   }.free()
    23. ) }, x* J2 O- \& U
    24. };: v: f  v! p, ?/ S) A3 U7 c
    25. f(A, B, C, D :i,s,X,Y:Array,max)=; P6 R\\" J7 l2 H
    26. {4 i1 {/ c( g4 c2 Z0 F( X
    27.   s=0,i=0,(i<max).while{\\" U8 n3 d) q\\" O2 ~
    28.     Array.GA[i*2, &X, &Y],
    29. . Z9 W0 _: _: r
    30.     s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,3 m3 [% n4 P0 Y7 x
    31.     i++
    32. - b' ]* s( H5 B: M# ?6 J
    33.   },
    34. 7 Q% l, y# _( m! n  M& ^( X
    35.   sqrt[s/max]
    36. . d7 Y% U. q& M/ v. T
    37. };) g. \7 U/ z! N3 m( ?* i
    38. Opt[HFor("f")];

    6 ^) U. u5 g' \7 O2组解(A, B, C, D,误差):  k0 ^& {# S/ r  E5 p' q8 p
    1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-003
    " _" k, @2 D: f" g) i( W% f( u! U7 B& x' O* d
    0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-003
    ( w/ ]( ~1 Q* S* H; N' B/ L
    2 Z& S6 @. J9 i' R& \" f看了楼主的一些帖子,似乎楼主不在于结果,而在于如何获得结果。
    . K# N! N# F: ?& Z$ F( n- \) u$ H' [- q" b# G1 G0 F
    一般的优化软件都没有公开其算法,如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。
    : |  y2 r. ?: k( {% g
    # w8 x* e; k! E( f9 F! J7 R
    * p1 r7 ?) F& r+ h
    回复

    使用道具 举报

    1

    主题

    0

    听众

    3

    积分

    升级  60%

    该用户从未签到

    forcal 发表于 2011-4-19 08:15 7 V8 c2 W5 [& U5 u/ d
    用Forcal求解,有2组解:
    & r6 {& w5 B: {1 Y, Y
    ( l5 w  l: M" S( a- Z5 k) J# n. g2组解(A, B, C, D,误差):

    * x  ~: c' L3 Y5 z是的,因为需要分析其计算过程。/ c9 J( i$ n% t0 P, J' \
    还是非常感谢。。。
    回复

    使用道具 举报

    243659583 实名认证       

    0

    主题

    3

    听众

    124

    积分

    升级  12%

  • TA的每日心情
    开心
    2015-5-13 14:42
  • 签到天数: 38 天

    [LV.5]常住居民I

    社区QQ达人

    回复

    使用道具 举报

    plmlor        

    2

    主题

    4

    听众

    17

    积分

    升级  12.63%

  • TA的每日心情
    开心
    2011-10-31 18:59
  • 签到天数: 2 天

    [LV.1]初来乍到

    很好!, T( w  u, L5 }6 e8 u* z$ e
    顶顶顶!8 b% @0 d. s' x# m0 ^6 }8 Z' E5 l
    3 ^- J3 v  F7 c2 E2 S
    远东森英—咱北方人的地板
    回复

    使用道具 举报

    bozsu        

    0

    主题

    4

    听众

    9

    积分

    升级  4.21%

    该用户从未签到

    祝楼主快乐,大家一齐讨论!; i0 H( |! B. I* A- v, {

    $ _1 @# P5 }' T% [6 s+ u9 P$ w& p# Q& x7 Y3 ?8 z5 A
    这里有木有买双色球的朋友啊,我有方法能提高中奖率,/ B9 v7 B& {/ x4 O9 D* e' y

    2 b- o6 w1 [+ x: k. S3 b想了解或者感兴趣的朋友加扣扣:2590587466
    回复

    使用道具 举报

    wdgqr        

    0

    主题

    0

    听众

    4

    积分

    升级  80%

    该用户从未签到

    真不知道,还有这么好心的楼主 支持& a( n* v1 R1 U2 Y) L) |( J

      F7 x* m$ X8 G6 I8 G; u- I3 _' a$ ^7 t

    : `& Y# L" a4 w) a
    # Z; N& L: B/ [
    ' O/ F3 e0 I- Q) Y0 _
    : @" E! S) ^; k$ Y( f
    * r* S+ V) z9 A& F) g6 R. Z9 {- a  C' M7 @3 T% J
    $ h4 G9 s% p% _; H9 }) |

    , \+ E, c3 V0 m+ ^) p- g
    ) b: w* \, {5 Q* V  d7 d* W5 r1 h$ `% r
    女裤品牌( g; L3 h5 y! \
    淘宝网女裤: V( R# U* H! p! z9 J* A- }2 o
    当当购书网; o4 g( }* z6 _1 d0 u8 e6 ]
    一号店团购
    5 e7 J4 s' e' Hhaoping66.com
    3 ?9 `# M. t2 T  _1 `6 D6 O优尚网
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-9-1 23:44 , Processed in 0.432968 second(s), 85 queries .

    回顶部