QQ登录

只需要一步,快速开始

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

四参数拟合方程 的方法

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

1

主题

0

听众

3

积分

升级  60%

该用户从未签到

跳转到指定楼层
1#
发表于 2011-4-18 16:47 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.% _+ r( a% b9 W2 k2 a
不知道有没有人知道如何求解.
2 m, @3 o6 e& o8 E8 c. a6 t0 K) {9 _% ^. d
方程如下:(四参数为A B C D)
3 C# u& F+ s$ Z3 {5 a6 FY = (A - D) / (1 + (X / C) ^ B)) + D
5 O  y+ T$ V& r! H( i6 s* k" e- v" C8 L. d) e
给出数据集:
' r( `4 }# {: d9 Z  x y
6 G4 s3 q. w1 p& z0.00 1.47358 X1 `# |' E. S2 R& s
0.025 1.214
% g! F" y7 \0 ~0.05 1.0143 6 s& V, x6 R: T0 s7 y  Q0 K
0.1 0.7362 4 e$ C& p- M! L( K
0.25 0.3907 5 a1 m3 e. l  |) \( @
0.75 0.21357 m6 E( f2 {4 g: Q4 w, m$ V/ N9 K
1 k% p7 R1 g4 z) E1 U
求最合适的四参数ABCD值是多少?) U! k3 O, |) ?, ^

. J3 w% A4 D. Y: C0 z. n- i4 p答案是 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组解:" \% z5 I; W7 ?7 q+ O4 L

    1. % d* d; e' d% |! H
    2. !using["fcopt","math","sys"];2 v  p! s6 n+ q% K& q' G4 R
    3. init(::Array,max)=
    4. % V8 Q  v6 N0 R, j4 O6 J6 z4 l; {; h
    5. {
    6. ; s/ N\\" m- o4 L$ P
    7.   max=6,
    8. 9 N. A; N1 o# \5 E$ B1 r$ Q9 n! C
    9.   Array=arrayinitns{max,2 :) b( w2 P0 g# w5 F5 n
    10. "2 f# l7 x5 }- u( A
    11. 0.00  1.47351 F\\" p, y, [# ]1 v& W
    12. 0.025 1.2140 H2 ?: v3 u+ L! e
    13. 0.05  1.0143 1 q5 M2 X7 }+ Z/ Z$ v' {
    14. 0.1   0.7362
    15. ) v0 a/ G5 E$ A& \8 ^# b! ^; W0 x( l2 v
    16. 0.25  0.3907
    17. ' b0 G) v\\" x0 C  q8 g
    18. 0.75  0.2135
    19. 9 l! _. Z; A$ L7 F  c. q2 k
    20. "
    21. 6 z1 s# Z+ U) x. C
    22.   }.free()
    23. % J( T1 d; G\\" W: ^1 \( h\\" {: Y\\" @
    24. };4 [' Q! m5 [( y% c
    25. f(A, B, C, D :i,s,X,Y:Array,max)=
    26. 0 h+ B. e- }3 @, V% T% t; f\\" B5 L
    27. {& P: w: G1 [9 y& M. t5 ]2 C
    28.   s=0,i=0,(i<max).while{
    29. 6 @& ~% j; t# U7 i# n% ?( t3 `. l
    30.     Array.GA[i*2, &X, &Y],: z, W. h3 N/ C1 r
    31.     s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,( K+ O! L! T: N! |. e# b/ w
    32.     i++. l: b$ T4 t$ i\\" a. e* U+ V1 D$ \
    33.   },9 F. B, G8 w2 F* v1 f3 M5 g
    34.   sqrt[s/max]
    35. - K\\" S4 y# {: ]/ h$ h
    36. };
    37. * t: Q0 m, r# q: Z1 x
    38. Opt[HFor("f")];
    2 x) i: b1 v  Q1 v8 s
    2组解(A, B, C, D,误差):& C5 g/ T& G9 a
    1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-0032 m8 G- J; d! O5 ?% \# H/ `

    $ v  C; v! B2 m3 T7 \$ Q* v0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-003$ a; E0 \: u- T+ W2 _) T

    8 B( W, \1 h7 T$ y7 S/ w看了楼主的一些帖子,似乎楼主不在于结果,而在于如何获得结果。: n- K. w2 y3 t1 q0 l
    % ]9 m3 P% @. D: [6 {/ W$ A
    一般的优化软件都没有公开其算法,如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。
    % E6 K6 X4 d( ^( k, ^' ]* s; G7 ^. @) j! H+ r- B  h7 R

    , @5 R. b2 J1 l% p! |: d
    回复

    使用道具 举报

    1

    主题

    0

    听众

    3

    积分

    升级  60%

    该用户从未签到

    forcal 发表于 2011-4-19 08:15 6 ~5 q) {4 T0 J
    用Forcal求解,有2组解:+ I# Y7 j7 x) r( t

    ' s! e6 z6 N; I0 Z2组解(A, B, C, D,误差):
    5 `9 M0 Q8 J" x) A
    是的,因为需要分析其计算过程。/ i8 r) M, {/ J# [. c& g  n# x5 ^3 I, c
    还是非常感谢。。。
    回复

    使用道具 举报

    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]初来乍到

    很好!
    % U) g/ T2 Z7 i9 w5 o/ w* Q顶顶顶!  H- Z' Z( e% _" x9 V
    1 h: [. Q* @8 A
    远东森英—咱北方人的地板
    回复

    使用道具 举报

    bozsu        

    0

    主题

    4

    听众

    9

    积分

    升级  4.21%

    该用户从未签到

    祝楼主快乐,大家一齐讨论!
    # L0 b0 G. X! e( B" f7 E2 w/ t- \" }( w1 W8 G- k1 u
    0 o) v3 w( G2 n# x; j
    这里有木有买双色球的朋友啊,我有方法能提高中奖率,
    " _& g3 `6 I7 i" d* i9 a; Z6 g' Z* E8 E; x
    想了解或者感兴趣的朋友加扣扣:2590587466
    回复

    使用道具 举报

    7#
    无效楼层,该帖已经被删除
    wdgqr        

    0

    主题

    0

    听众

    4

    积分

    升级  80%

    该用户从未签到

    真不知道,还有这么好心的楼主 支持9 A$ Y7 f% F7 w- S- b

    1 u/ X8 W; F0 q- ]7 S* u( G- A( r$ u1 @7 E* x# o  E
    - M& P& z4 g/ x8 B6 N" c

    / U3 p1 u0 {& d3 @- S2 @) W5 \% H, w2 e- L! D- o7 j
    & O5 R& l/ I, e6 V" C. }9 m
    % D" p$ M) J' z% R, n
    4 u. T: b2 V" H0 L0 t& T

    ; B; R/ a0 c5 J1 Q. h: G0 u4 `9 U# s5 o0 Z; d/ N+ W+ e

    9 D. n" Y3 e0 ?& @5 x" G- I6 \# J9 J3 M+ o  d
    女裤品牌
    5 s# L9 V# m. K9 r  M3 r淘宝网女裤9 ^$ U! `) H) u3 g/ }* @" n' b
    当当购书网
    & w+ \. F3 x8 h一号店团购7 @9 R* U' [* G9 a; ^9 z; b
    haoping66.com
    6 F* S. G" P2 l0 R7 k. j- z优尚网
    回复

    使用道具 举报

    9#
    无效楼层,该帖已经被删除
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-6-14 11:50 , Processed in 0.703023 second(s), 96 queries .

    回顶部