数学建模社区-数学中国

标题: 四参数拟合方程 的方法 [打印本页]
作者: 〃Lan_Pst。    时间: 2011-4-18 16:47
标题: 四参数拟合方程 的方法
四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.( g+ z3 K6 T  b0 o* l" Z; S4 @, U
不知道有没有人知道如何求解.5 l: Q3 d. v2 c

  [, m. g$ R  {, a, k, F方程如下:(四参数为A B C D)
2 {/ j$ P) i& ?# ?2 f( r& |Y = (A - D) / (1 + (X / C) ^ B)) + D
" N0 W" a5 E4 l
4 [) f& A9 }! O给出数据集:
, |5 \3 i4 }$ _. Q% x- g/ `  x y* m, i, e; C! w" z+ _' M
0.00 1.4735/ U) N- @- s4 z9 {* f5 a$ O% m
0.025 1.214
$ M, t  q4 p& `* I" ~& |0 o& H0.05 1.0143
0 s% q& X8 _' h; y+ t4 J" K( J0.1 0.7362
8 B" P" N5 q' R" ?) c6 J0.25 0.3907
2 x9 y# R5 D0 V+ y2 z, L4 X- g0.75 0.2135# K7 e4 y3 E1 Z4 t
1 D  m: d5 X, G, ]  L; M
求最合适的四参数ABCD值是多少?: Q* E' i# c; I, ?. ]: f
% K; Y0 d+ @& G/ Q6 k
答案是 a = 1.47 b = 1.229 c = 0.085 d = 0.12
作者: forcal    时间: 2011-4-19 08:15
用Forcal求解，有2组解：7 i4 s8 r! |4 ~. j7 a4 r0 n

  1. & n/ u$ y9 D8 D5 i
  2. !using["fcopt","math","sys"];% O( h. O# K; _& T0 v5 Q3 B
  3. init(::Array,max)=
    $ y7 z1 \& m3 q) S
  4. {. P, U$ n. T1 E1 u  Q) Y
  5.   max=6,
    - ^9 P" f3 K7 b  f
  6.   Array=arrayinitns{max,2 :
    . _/ }, P! S! N
  7. "7 _" t  I+ a" U/ ]
  8. 0.00  1.4735
    3 q4 c: z) ^2 B% F1 A
  9. 0.025 1.214/ q5 R/ L( A) i, }- t* `) s8 O8 ?; v
  10. 0.05  1.0143
    " N  m, A2 Z9 z" b  J2 W, I
  11. 0.1   0.7362
    $ W3 @7 k. ]; N+ {
  12. 0.25  0.3907
    % t' c& `$ U7 e. v3 U& z- j3 x% Z
  13. 0.75  0.2135  Z& |4 K( {  f4 p+ Y, A2 U
  14. ", K# q9 }2 d# h; c% q
  15.   }.free()) ]+ Z, x2 h8 p) c" v
  16. };! ~" ?% [% Y7 u9 _/ l. G1 E
  17. f(A, B, C, D :i,s,X,Y:Array,max)=
    , I4 w' }% [, G) {, a0 ]8 U( o$ q
  18. {
    6 I; H9 ^. m2 g7 H& Q4 W
  19.   s=0,i=0,(i<max).while{
    + \9 p- b* H) A) j5 m
  20.     Array.GA[i*2, &X, &Y],
    7 K# |3 i; p2 j  }& U5 f9 G
  21.     s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,
    & x/ `- b" }. R9 g5 N" I9 u. P
  22.     i++
    ; X& e$ |( E/ ^8 v6 [2 ?" f
  23.   },
    6 D' G  t# t6 J7 X! d
  24.   sqrt[s/max]
    5 f% x% H, l* Y8 Y5 \3 F
  25. };, O6 v# E' a, [4 t
  26. Opt[HFor("f")];
复制代码
; G# _) c5 i1 O
2组解（A, B, C, D，误差）：
( J0 Y, r* z9 {$ z* g2 ]7 d1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-003
% g! e, P! H3 m$ `& s( {; }8 u1 p. R$ |, w0 b% S
0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-003# Z8 D  s) H# {, P
# `1 p! ]* z) l
看了楼主的一些帖子，似乎楼主不在于结果，而在于如何获得结果。9 O; r' W6 ^* p7 {

: g3 r4 |* c+ O+ l# a3 Q# k. o! E5 F一般的优化软件都没有公开其算法，如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。/ j8 G  W2 X! k6 \7 ]4 y' h  p. ?
5 Q# {' z+ Y* h& o; J5 ?
* v. L$ ^6 z4 ^. X
作者: 〃Lan_Pst。    时间: 2011-4-19 10:54
forcal 发表于 2011-4-19 08:15
& a# Q8 [& o" o' X, ?用Forcal求解，有2组解：# F; y8 `& ?; J* e- \* b

8 O8 u* Z7 J7 h+ i- I2组解（A, B, C, D，误差）：
+ Y7 V" s2 ^- \7 E# s! k& h( @/ m
是的，因为需要分析其计算过程。% e" i* a7 A+ C  c4 E8 T8 [. J
还是非常感谢。。。
作者: 243659583    时间: 2011-4-19 14:03
baozheng
作者: plmlor    时间: 2011-10-30 08:22
很好！
8 b1 W$ k2 Z& g& M* c顶顶顶！
1 Y* F) }- i  G1 n8 A5 t3 d) d( f, n3 y1 h( Q
远东森英—咱北方人的地板
作者: bozsu    时间: 2012-1-10 12:09
祝楼主快乐,大家一齐讨论！- `* X+ [) S+ f
) a) J9 |8 @2 I- R% v+ j
# |$ c% h1 i3 y9 C
这里有木有买双色球的朋友啊，我有方法能提高中奖率，) j" C. k9 q% Z. |% b

8 z) c  r# i5 }- f% i想了解或者感兴趣的朋友加扣扣：2590587466
作者: wdgqr    时间: 2012-2-25 16:56
真不知道，还有这么好心的楼主 支持
* g( t; b, U. E. ]6 Y& k0 u' c; s5 m; O" c0 F5 |6 a& K. }8 O3 s

5 r& _1 l2 R9 |4 u# I. _' W3 ^" V* t: y+ U* k4 r

$ r) E; ^  t& Z% F- U6 v% q" D: n8 g8 B% U8 A

$ x# u* v3 L% [1 {4 _0 f. u+ X5 d$ s5 v/ K# ?1 Z* y: y9 T

7 r6 I* t" C" A  R7 [& w1 e' f7 o) r6 O/ P4 o- v7 t

* t% H6 q! s  ^# h6 @4 I+ q( Q% T
) J% p$ b. l# z: b" C! G# |0 A; d! i8 E4 b" ]! T+ Y! g) k
女裤品牌
7 T) {' |5 o' f2 e6 e3 Q淘宝网女裤; N+ p6 @; ?9 A6 Y4 [
当当购书网
+ ^+ ^+ C- W7 {2 L一号店团购
" E. A+ g  m7 q& c4 K) |haoping66.com1 {7 A" W0 t% @
优尚网



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5