数学建模社区-数学中国

标题: 四参数拟合方程 的方法 [打印本页]
作者: 〃Lan_Pst。    时间: 2011-4-18 16:47
标题: 四参数拟合方程 的方法
四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.8 K4 w$ H. i; m- s* O% V" K1 @
不知道有没有人知道如何求解.. x! a; I. d, A5 \" k2 x3 z! B

: T9 X9 x) I$ X8 y9 R& V方程如下:(四参数为A B C D)' u7 t3 i4 N' l! {7 m- m! f* O- g
Y = (A - D) / (1 + (X / C) ^ B)) + D( B) {) R6 V* p* q
6 v% O6 F% F0 i: ^% @
给出数据集:' C2 Y3 R" g5 S4 g
  x y
" m' J( a/ n, Y* E+ q0.00 1.4735
7 _% k+ g# v4 h  \  b4 r0.025 1.214
3 |/ B) `+ @; Y# g0.05 1.0143 ( C% ^8 u& }* R- v" P& E
0.1 0.7362 + R6 ?: D2 h+ K! Y$ @
0.25 0.3907
) @" X9 i" z7 Z1 S0.75 0.21357 W$ l% d/ M* A4 Y/ }6 ~+ Q

& ~( C3 }1 C$ ^9 S- Z6 U1 S求最合适的四参数ABCD值是多少?/ m/ b: T! e0 ]2 g- y6 F* x& k

' l6 k" b. r  d3 M7 ~答案是 a = 1.47 b = 1.229 c = 0.085 d = 0.12
作者: forcal    时间: 2011-4-19 08:15
用Forcal求解，有2组解：) |8 d) R8 \! o- N1 i4 q$ D
  1. 6 _. |0 P4 ]7 X
  2. !using["fcopt","math","sys"];: g3 g; O: |/ k. B
  3. init(::Array,max)=
    , _$ ^- i6 y4 R: _
  4. {
    % b1 N$ f1 I" A; E2 Z/ k
  5.   max=6,
    . @# L( c& E, ^  w. N
  6.   Array=arrayinitns{max,2 :. `$ ?9 K" S% S
  7. "
    ) |' n- t1 r$ @2 r/ ~
  8. 0.00  1.4735- `# i6 {/ ^  Q9 o8 C6 i  ^+ g
  9. 0.025 1.214
    0 k7 [- A- ]( T% H6 L
  10. 0.05  1.0143 0 N3 K0 q9 D/ G& }
  11. 0.1   0.7362
    ) H: g, ]! N9 _: R
  12. 0.25  0.3907 1 ^9 ]8 o$ o: U3 l
  13. 0.75  0.2135
    . c! h: q' E1 R( y: t
  14. "
    7 A/ @  j: V+ ^% k
  15.   }.free()
    0 e# a) r) e- P
  16. };) P3 k+ i) z! l
  17. f(A, B, C, D :i,s,X,Y:Array,max)=
    ( l5 Q! |8 U  Y1 a( g. k' q* u
  18. {
    , |% c3 a1 N7 d3 S) ?; d0 T
  19.   s=0,i=0,(i<max).while{
    4 S; c5 M2 c1 ]- F8 E: p
  20.     Array.GA[i*2, &X, &Y],
    % E2 d9 X6 i1 f
  21.     s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,
    % \; K. N: b  u) u! {
  22.     i++
    . h: H* e6 d( i" w
  23.   },2 ^8 Q6 ~/ y0 j7 k! c* s
  24.   sqrt[s/max]
    ( g7 C) _' z; ?: _) p' @) P+ f
  25. };9 b- n$ `6 N2 @; Z( H# K3 C+ C
  26. Opt[HFor("f")];
复制代码
2 }7 e4 f7 ~* g; Q2 r
2组解（A, B, C, D，误差）：
/ f- `9 m/ [4 {/ ~1 ^) ^( O1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-0035 N+ F4 n/ u, X: V$ ~; b) ^

' g" z) b+ C- p/ Q( M" r. I0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-003
7 t1 ]# e4 w0 c0 x( l6 X' S3 n" g/ {+ r5 W, f3 ?  I& U# |+ @8 V
看了楼主的一些帖子，似乎楼主不在于结果，而在于如何获得结果。7 t/ ^; V/ G4 n& I1 P+ H( l
. n# B5 v0 \; E
一般的优化软件都没有公开其算法，如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。
* u. f) _6 k- n3 j' l7 \# ~1 s
/ r: }: A4 L/ l1 l( O2 |
作者: 〃Lan_Pst。    时间: 2011-4-19 10:54
forcal 发表于 2011-4-19 08:15 : M* E4 z2 M% ~3 g- u* |
用Forcal求解，有2组解：
! V, u$ T6 p) T. w# c
; q0 v9 `' T& k3 T- N+ \2组解（A, B, C, D，误差）：

/ [. G0 w  K* e. ^是的，因为需要分析其计算过程。
8 q1 _+ P5 w- h5 Q7 B 还是非常感谢。。。
作者: 243659583    时间: 2011-4-19 14:03
baozheng
作者: plmlor    时间: 2011-10-30 08:22
很好！1 T+ [: |; b( E9 b
顶顶顶！; q) m+ ^1 @$ J1 y( d/ {
0 O1 E' _! i" H9 n
远东森英—咱北方人的地板
作者: bozsu    时间: 2012-1-10 12:09
祝楼主快乐,大家一齐讨论！# e6 e9 \% {1 G

6 l1 e" i) ^2 B4 x7 z0 e: T6 ?% Z' l. F6 q) k# c' I! o( C6 s
这里有木有买双色球的朋友啊，我有方法能提高中奖率，9 I! L4 _4 L' P' z
8 w* a5 n4 H1 i3 Z$ M! l
想了解或者感兴趣的朋友加扣扣：2590587466
作者: wdgqr    时间: 2012-2-25 16:56
真不知道，还有这么好心的楼主 支持
* o" l* E: D' D/ H5 A. F8 a1 ^3 S! b

# _1 p# U% I1 y& l0 [) v! F* ~2 O" r- `) `8 l0 {: G3 }+ m

1 [  ~& q9 _1 ?9 q
1 v% E# ?/ A/ g( M$ N# k! s5 @3 `5 H! \  K. H$ s, i8 X, L

" S' A7 s3 g0 N, o/ C: X% r0 L1 _5 g1 Z; R& Y" q: K

: Z/ F5 Z% R& j0 u$ o% E* p
+ O0 K/ N( A& q) F; ?' i0 X# x" ]6 `! U7 b1 E
& [4 i5 H9 I% a7 J; I; b( _1 s9 Y
女裤品牌) R9 D) v1 ?* s8 s
淘宝网女裤8 c  m0 l- @) t" c8 D; c1 b
当当购书网
# {9 ^4 r9 C( D一号店团购# |; T% v: m0 R
haoping66.com
! z3 }; z' f$ H( ?$ H4 [; B; I优尚网



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