数学建模社区-数学中国

标题: 四参数拟合方程 的方法 [打印本页]
作者: 〃Lan_Pst。    时间: 2011-4-18 16:47
标题: 四参数拟合方程 的方法
四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.
: v+ S' E5 l( C# T; v& S不知道有没有人知道如何求解.
+ I( n& V- F+ d1 k( P& _; _5 R4 n9 C4 N5 b
方程如下:(四参数为A B C D)7 K( C6 j; `! ^$ b$ D) m
Y = (A - D) / (1 + (X / C) ^ B)) + D
8 _+ Z. q" k- ]
; X: T# S# s) y/ G& t给出数据集:
8 \5 S" X8 j3 ~: H" d3 Y  x y
9 Y) M0 E( I0 C$ y. o8 H0.00 1.47357 i# b6 Q) U; ?
0.025 1.214! P+ L- A8 Q* ~2 C
0.05 1.0143 6 P( H- ^$ H% q: z+ b
0.1 0.7362 0 Y; X( |2 L7 s
0.25 0.3907
& T# s2 D5 E& G# Y* T+ N0.75 0.2135' P6 f" L* L  _2 v( W+ o
2 J9 j) f3 U$ X
求最合适的四参数ABCD值是多少?
# Z% f" \- n4 K5 S% H2 z0 M. G/ d0 d7 A& V1 ~" h2 k
答案是 a = 1.47 b = 1.229 c = 0.085 d = 0.12
作者: forcal    时间: 2011-4-19 08:15
用Forcal求解，有2组解：
# `. Z$ r$ Q* U; y, V

  1. 1 w5 G' J( g6 p# _  d" [/ m  t* F
  2. !using["fcopt","math","sys"];
    ; d' \4 u* I2 a% s# j  |
  3. init(::Array,max)=# E8 P  i; A3 `/ }
  4. {7 j" A6 s$ ]& B% d  Y+ |
  5.   max=6,. M$ f! f  `5 I* U; l0 R3 N) \% f
  6.   Array=arrayinitns{max,2 :. O8 p$ P3 M% Z& ?- B
  7. "
    : y! }0 P4 P* L
  8. 0.00  1.4735
    / i$ F8 t, W' m- B; u, L8 I0 {
  9. 0.025 1.214+ T; \& d4 y6 @
  10. 0.05  1.0143 & ^- C+ b5 ?3 L4 l+ g
  11. 0.1   0.7362 ; n  L7 j5 D4 ?. M$ K# a& |( ]& S
  12. 0.25  0.3907
    9 F& Z" Q" H1 T7 e0 J+ h
  13. 0.75  0.2135/ L% X) W5 o* G. B" o5 v# W
  14. "+ [+ Y  D& g# S+ c" U0 ^
  15.   }.free()1 x, S) s$ i4 r* @5 G6 W- ^1 _
  16. };8 W2 v, i; B  C$ R/ O) ?, \2 o
  17. f(A, B, C, D :i,s,X,Y:Array,max)=
    ; z7 i" j4 l4 p
  18. {
    ! Y. g$ M/ N; s9 o1 E# j
  19.   s=0,i=0,(i<max).while{
    . `2 t4 [' \, @# n
  20.     Array.GA[i*2, &X, &Y],
    * L1 ], B3 Y) T& T
  21.     s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,( g1 n- _, p9 {% O  u
  22.     i++
    ) z8 N4 a) f3 {4 k; ^8 _
  23.   },
    " F. }$ L$ n% w+ t" N* j
  24.   sqrt[s/max]
    6 k4 Z" k& }0 t3 K+ x# M
  25. };
    0 c7 [5 H' W2 C, Y3 O9 \
  26. Opt[HFor("f")];
复制代码

" B  E& x1 [: `' S2组解（A, B, C, D，误差）：( e9 I3 ?7 h  x' X8 x
1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-003
" k# C4 y* R2 _' J+ v/ t1 ^" C/ M& w6 ~5 z. X1 L
0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-003. m) F( c' r* H! P8 H7 z

0 p! M7 a9 K' f看了楼主的一些帖子，似乎楼主不在于结果，而在于如何获得结果。, P: M+ \7 u0 L. A7 x) @1 m& Y8 d

. `" R% K/ q" C; M7 t! ~: u一般的优化软件都没有公开其算法，如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。8 B- Z% m: v0 l( A" k8 f! Z
" Z! {% E) o2 N" W) S

0 M- f0 \% t/ H
作者: 〃Lan_Pst。    时间: 2011-4-19 10:54
forcal 发表于 2011-4-19 08:15
# A* C3 c4 u. _" _用Forcal求解，有2组解：: ]$ ?, V, }5 Z; X
; v% p2 R0 V; j
2组解（A, B, C, D，误差）：
. W% F& O- p, v1 D4 L, k0 A
是的，因为需要分析其计算过程。$ Z, O# y- ~" W
还是非常感谢。。。
作者: 243659583    时间: 2011-4-19 14:03
baozheng
作者: plmlor    时间: 2011-10-30 08:22
很好！
% p$ P7 M) B+ Y0 n6 m) }顶顶顶！9 [: z7 l) H1 a9 {6 M/ f5 v

% J0 ]8 {; x( l远东森英—咱北方人的地板
作者: bozsu    时间: 2012-1-10 12:09
祝楼主快乐,大家一齐讨论！: t3 S1 e. S* }, B
7 ]' n! o" u0 s7 Z4 K* K% t7 a3 f

: @% I+ [8 j# c: d8 Y8 V9 M这里有木有买双色球的朋友啊，我有方法能提高中奖率，" l1 {8 v! c; K% _4 [0 r
" n% `7 |8 H! S" a! u' D  ]: \
想了解或者感兴趣的朋友加扣扣：2590587466
作者: wdgqr    时间: 2012-2-25 16:56
真不知道，还有这么好心的楼主 支持
5 c- M: Z. r4 j
7 {4 ~7 j9 P% u8 m; s  ?$ h  A# \7 s' a

5 y6 ?& k# b1 E
) R  t  M" X0 f/ m
) ]+ R3 N' E) S! ?" |, c! ~+ O) |( @9 q5 o3 y* ]; M  n
8 |6 Z2 T! b/ m" y" l$ H

- I% a9 r! J3 }4 M# k
' z! K1 L# Z* T' @% y
; Q/ W$ x; |6 o4 r+ W! ]$ `2 {$ N& e4 [# J- G' \- X

. x: |3 T1 c( g5 ^/ z% o0 s" U: _女裤品牌
9 W; p' \, z  @( l) m) K- D淘宝网女裤- a) Q3 q3 J6 A- R" @  l4 ?
当当购书网
* l, Q1 ~3 a* S) n一号店团购
( v' u& P! o. A! [- y: C% shaoping66.com/ V' {% k% b6 _" U' j
优尚网



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