数学建模社区-数学中国

标题: 四参数拟合方程 的方法 [打印本页]
作者: 〃Lan_Pst。    时间: 2011-4-18 16:47
标题: 四参数拟合方程 的方法
四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.  D+ T* t: v$ E0 i8 g$ O
不知道有没有人知道如何求解.  Y! ?3 \3 ~' F

0 D% h0 W8 b, Q, t4 y+ p方程如下:(四参数为A B C D)
( Z! b, U8 Z% L1 G, v( a$ uY = (A - D) / (1 + (X / C) ^ B)) + D
- C' ]0 K2 d: L% x2 G( U
0 i0 Y7 N! C6 w6 ]/ y给出数据集:
& r2 G4 ]- p  u% y9 M" G  x y  A( F' z% M9 i6 e6 h0 x
0.00 1.4735
& ?1 m2 I/ J9 C  a" U5 A3 _# l/ U0.025 1.214
6 |8 I& e* s7 ^) y# ~# d2 O0.05 1.0143
9 s3 `, j5 ?8 s0.1 0.7362 ! Z6 c+ B9 X9 G' ~: L
0.25 0.3907
4 R* @. d1 T; p- R/ K, v, ~7 o5 j0.75 0.2135
2 @6 ?4 O- N7 n2 U$ U* h
9 ~& \  m* J. I求最合适的四参数ABCD值是多少?
  u- ^9 V! Y0 _( t' I
9 R4 R' c1 y+ ]/ E答案是 a = 1.47 b = 1.229 c = 0.085 d = 0.12
作者: forcal    时间: 2011-4-19 08:15
用Forcal求解，有2组解：
: R+ M3 o) j$ h+ s* W* L* C! W. @
  1. ) ~/ [) ]3 i( @
  2. !using["fcopt","math","sys"];
    ! S3 q2 g4 Y1 Z9 K' a+ N/ H
  3. init(::Array,max)=
    5 D& S# ^7 r5 n0 F" A  X4 R
  4. {6 T# ]: R9 H! d& J+ @7 n
  5.   max=6,
    5 R$ F0 p" J7 `- Z+ ?6 a2 p
  6.   Array=arrayinitns{max,2 :
    . P8 }1 \2 r3 o7 M% b8 e( X9 m
  7. "
    4 ?5 p" h. `8 v- }' q9 W
  8. 0.00  1.4735$ D7 d. D' i' E' l- y5 b
  9. 0.025 1.214% p9 ^' r* d4 V  }; Y, F4 j) T
  10. 0.05  1.0143
    / F9 U8 V4 w0 s- r- ^* G
  11. 0.1   0.7362
    - f. Y& i4 o* A# @6 a4 t
  12. 0.25  0.3907
    2 e2 f0 U7 b* w/ T
  13. 0.75  0.2135: r, ]4 v0 t1 A6 j# n
  14. "/ M# D8 j" J; q$ M
  15.   }.free()/ C% [7 Q4 h9 s! T) y6 U
  16. };* T: F$ T/ `' W! c# o
  17. f(A, B, C, D :i,s,X,Y:Array,max)=
    . m' ^( O9 L- W7 `6 F4 T
  18. {: p/ i0 @3 n2 d! f# L3 }
  19.   s=0,i=0,(i<max).while{
    ( ~, b2 L2 F9 Q, |6 M' I# c9 r
  20.     Array.GA[i*2, &X, &Y],, e* j6 _$ X5 }. u! i
  21.     s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,
    . f" a  X3 N9 K: U; x
  22.     i++
    ! V0 F1 _5 V5 g+ e& E' l' L
  23.   },
    ! E5 g) ^) E4 X0 ~
  24.   sqrt[s/max]# b1 |9 ?: p; u4 [  F
  25. };7 P1 L; |* S4 k- ~7 A* F2 B& l$ m! r+ Z" R
  26. Opt[HFor("f")];
复制代码
1 s8 Z. R7 ~  f5 Z8 s3 v9 n! g
2组解（A, B, C, D，误差）：
% E) j) v7 _2 Q. A1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-003/ M( r$ |# g! W
0 G9 n: g3 e; V- j: `
0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-0039 q+ Y. u- k! v. I

! r' n3 x+ C- m; q看了楼主的一些帖子，似乎楼主不在于结果，而在于如何获得结果。
! j$ r* U0 T" i$ e  D: f5 B  I! c& J) m0 `, `4 `" s
一般的优化软件都没有公开其算法，如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。; }3 Q+ M: y9 L

& k& `2 _, s2 Y
/ y# h; y: E8 @/ y( g" r
作者: 〃Lan_Pst。    时间: 2011-4-19 10:54
forcal 发表于 2011-4-19 08:15
' ^. T" K" {! U用Forcal求解，有2组解：, `4 {( i. h6 D+ y- p% a
! z7 m7 J5 K  S6 |
2组解（A, B, C, D，误差）：

; m! z4 q: p" e; `是的，因为需要分析其计算过程。/ z$ c4 T8 s2 {. t3 H9 L
还是非常感谢。。。
作者: 243659583    时间: 2011-4-19 14:03
baozheng
作者: plmlor    时间: 2011-10-30 08:22
很好！
0 o1 I; Q3 v2 A* S; G8 @顶顶顶！! a; [7 _! L; E4 M- q" A7 ?! A* v  w

5 l: b+ {' Z) I: k远东森英—咱北方人的地板
作者: bozsu    时间: 2012-1-10 12:09
祝楼主快乐,大家一齐讨论！
0 D7 T/ D+ ?+ [6 q; w3 h' z" C1 O5 S& L% D- \6 ^
5 A0 x3 V# O) T* t8 Z( P/ o
这里有木有买双色球的朋友啊，我有方法能提高中奖率，
' K" Q* j" T, F2 b; D* F
( x0 F7 k" j" L4 W想了解或者感兴趣的朋友加扣扣：2590587466
作者: wdgqr    时间: 2012-2-25 16:56
真不知道，还有这么好心的楼主 支持
+ {" N7 ^% h$ P. u  s+ K$ v8 [
+ y- [# V' O' f/ \
+ A( x/ b6 n& H1 e8 b  [$ U5 c- m8 p/ H3 `0 x
2 |& u* p$ c6 c' s, }. X/ d

! N* _- G6 A6 I) g7 Q' ~" Y) Q3 t; i5 Y' S- H3 ^0 \  p

) {% e8 w% v) Q* t( @9 w- d7 B. P  X
7 \/ U% {& Y  g$ h8 l$ L+ G2 V+ Y5 F  [. t8 U6 a6 s6 ]- A& N
, V3 o& N% W' p% w

5 l- E5 P1 D- `( W7 f: p
' t% B& s, w. a, X1 T2 T* C( C女裤品牌) b4 k5 z0 l  l+ |6 D
淘宝网女裤
& K2 A5 J3 c7 c1 T! h; Z, _当当购书网
# |# S! L( o* {9 m一号店团购
8 {4 N, s! K$ g, w* l. khaoping66.com5 m2 [/ {) H% J$ O# g5 d2 U
优尚网



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