QQ登录

只需要一步,快速开始

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

四参数拟合方程 的方法

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

1

主题

0

听众

3

积分

升级  60%

该用户从未签到

跳转到指定楼层
1#
发表于 2011-4-18 16:47 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.
* s' B$ O; r8 u# V7 G不知道有没有人知道如何求解.8 g5 m0 N" o6 \+ R4 c: {) B( T

% U! p, g% y& T2 c/ ^# n方程如下:(四参数为A B C D)
& e: m6 O5 B2 Y3 X0 |Y = (A - D) / (1 + (X / C) ^ B)) + D9 P! d2 u' s4 }7 O6 Y. z$ V
2 ^. J' i1 E  h) N, E
给出数据集:6 A' x$ \' K" q6 S8 k- @5 r4 w9 N
  x y
3 L( w3 S4 z$ k6 N2 y# w! l0.00 1.4735/ d0 ?( V7 h& R2 Q' Q6 y) n
0.025 1.214
3 n3 z% M( G8 j" w' k0.05 1.0143
3 P/ E* q& t4 }6 H. g) ?8 _0.1 0.7362
0 q  q3 m+ |# m" b2 U# B0.25 0.3907 * w( ~. _% ^# }) {3 y+ r) Z
0.75 0.2135/ o% n% r9 z% t2 e% |

8 h4 s: E6 G# v求最合适的四参数ABCD值是多少?9 Y/ X; M: s1 T7 Y7 v; p
& o0 A& b( E: n2 o# X
答案是 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组解:% I. }5 j# L1 q0 w% M

    1. 0 Y: U7 R) I% _+ m* I* g
    2. !using["fcopt","math","sys"];$ s\\" r. A! D6 f6 w  z# S  t
    3. init(::Array,max)=
    4. ' x% s! e7 ^' y2 \' I( z# i/ m
    5. {; M( z' t# W& A7 ?3 T+ u8 ^# O, y
    6.   max=6,. K: h; K: u+ s1 w\\" a
    7.   Array=arrayinitns{max,2 :, t3 e2 Y# T& }: U, I
    8. "+ z6 @' z: r9 r  z2 S& m( Z9 l
    9. 0.00  1.4735, m0 v$ l4 c1 I* w8 g
    10. 0.025 1.214
    11. 1 ~. p' i\\" c, h/ x0 T7 t: Y. v
    12. 0.05  1.0143 - R- _2 _\\" S8 P! f+ B. |: f6 o( ?
    13. 0.1   0.7362
    14. 6 Y* ~5 `  ^, B5 H. q
    15. 0.25  0.3907
    16. + |! V+ ?+ w5 G5 o
    17. 0.75  0.2135
    18. 0 a9 H0 t. \. ?: p0 w. p0 e
    19. "7 r7 i0 D( ]; |- r
    20.   }.free()
    21. # ~' q; _) J4 B, E
    22. };
    23. 9 x6 G# h( d) Y, V( D5 R  m
    24. f(A, B, C, D :i,s,X,Y:Array,max)=: `- v- Y  y* x/ y6 g9 D  b3 n, Y
    25. {6 b; y# r9 F6 U% a6 \- G* |
    26.   s=0,i=0,(i<max).while{
    27. 9 ]! {3 R! ]. d
    28.     Array.GA[i*2, &X, &Y],
    29. 9 [/ o) y\\" c8 Y6 h7 I5 b+ D. |6 h
    30.     s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,. V' ~  n$ s! ^4 K
    31.     i++5 k0 ]; k. y' x; C* S
    32.   },$ q0 j0 a% y* `/ X$ v5 l
    33.   sqrt[s/max]3 @\\" Q6 e' ^7 E+ `\\" N3 J0 P
    34. };+ z6 o' i3 F2 n3 d! q7 X! {
    35. Opt[HFor("f")];

    * k4 i5 s1 A* e' ]2组解(A, B, C, D,误差):: D( ^/ A0 |$ c  ^
    1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-003) N! D8 V5 t, L
    ( Z3 M8 a8 M$ U9 M# l
    0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-003
    / n- l0 o" M, i. S; D" c3 K1 `, t2 ?
    看了楼主的一些帖子,似乎楼主不在于结果,而在于如何获得结果。
      V! v0 ?7 C3 s+ M4 S0 s7 C
    ' v; Q9 |9 i! J* u/ C一般的优化软件都没有公开其算法,如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。
    8 Y  ]7 }; ^; d) P
    ; I8 k6 c; E! p: x" M; y$ T! G$ F- Y. k3 I6 F
    回复

    使用道具 举报

    1

    主题

    0

    听众

    3

    积分

    升级  60%

    该用户从未签到

    forcal 发表于 2011-4-19 08:15
    ! u( W5 j/ J8 B# _) y6 K* l用Forcal求解,有2组解:- [/ g8 J: B) d$ f- Q

    " D' j. J, g8 M) \9 J1 m2组解(A, B, C, D,误差):
    5 e# \" u0 R/ y
    是的,因为需要分析其计算过程。
    " q! F- p1 K, D 还是非常感谢。。。
    回复

    使用道具 举报

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

    很好!( L7 M/ e- l0 r$ L' Q/ U
    顶顶顶!+ F* D" ]4 _2 F6 \$ o6 g. I

    + _* W, r! p7 k远东森英—咱北方人的地板
    回复

    使用道具 举报

    bozsu        

    0

    主题

    4

    听众

    9

    积分

    升级  4.21%

    该用户从未签到

    祝楼主快乐,大家一齐讨论!
    . T1 c# N6 U0 S  |: ~4 ?+ [7 j- J* m9 s9 {' Y. X

    - r3 @4 Z9 }. y$ w! W" Q. B7 `这里有木有买双色球的朋友啊,我有方法能提高中奖率,
    # k2 u) t, d( e' Y  |
    / @; k/ E9 ~$ o9 h想了解或者感兴趣的朋友加扣扣:2590587466
    回复

    使用道具 举报

    wdgqr        

    0

    主题

    0

    听众

    4

    积分

    升级  80%

    该用户从未签到

    真不知道,还有这么好心的楼主 支持
    % }5 @$ I# I5 Y
    ( v. D" o( A0 w1 X9 P* q9 x
    ' T: G. r& G4 a' H# D- y8 ?: V8 c" I9 n9 p" H

    ! g$ h' j& y- P  D1 o7 k- K
    1 n2 ?. q" f7 i# H. M" d  j) T& X$ D: u9 H) Q

    ' A: M6 U+ W$ o( @% g; y6 X) N) m! }, v0 Q8 B( I! _' [0 K

    7 \# _; S' k8 R1 t2 F2 B) E$ X) ]6 ^- s( T; R
    $ r% s) y  z6 ~1 ]; p) ~! ^' z- _

    / ^4 W: s$ Q( y" r* a女裤品牌
    2 D0 Y) [( A2 w: `+ N# s淘宝网女裤
    6 p# }2 P$ @: K' @- a' k! K当当购书网4 U2 S- U2 Q5 H
    一号店团购
    - `1 ]0 A+ }" hhaoping66.com
    ' q8 N3 E* h- ]  n优尚网
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-8-31 18:20 , Processed in 0.400367 second(s), 85 queries .

    回顶部