〃Lan_Pst。 发表于 2011-4-18 16:47

四参数拟合方程 的方法

四参数拟合方程 的方法(4 Parameter Logistic)找了很久都未找到.
不知道有没有人知道如何求解.

方程如下:(四参数为A B C D)
Y = (A - D) / (1 + (X / C) ^ B)) + D

给出数据集:
  x y
0.00 1.4735
0.025 1.214
0.05 1.0143
0.1 0.7362
0.25 0.3907
0.75 0.2135

求最合适的四参数ABCD值是多少?

答案是 a = 1.47 b = 1.229 c = 0.085 d = 0.12

forcal 发表于 2011-4-19 08:15

用Forcal求解，有2组解：

!using["fcopt","math","sys"];
init(::Array,max)=
{
  max=6,
  Array=arrayinitns{max,2 :
"
0.00  1.4735
0.025 1.214
0.05  1.0143
0.1   0.7362
0.25  0.3907
0.75  0.2135
"
  }.free()
};
f(A, B, C, D :i,s,X,Y:Array,max)=
{
  s=0,i=0,(i<max).while{
    Array.GA,
    s=s+[  (A - D) / (1 + (X / C) ^ B) + D - Y]^2,
    i++
  },
  sqrt
};
Opt;
2组解（A, B, C, D，误差）：
1.470073083612318         1.229193695234829         8.502453298160233e-002    0.1204919664733336        8.553992238184498e-003

0.1204927998043004        -1.22919582966611         8.50244838406566e-002     1.470073025579681         8.553992239016503e-003

看了楼主的一些帖子，似乎楼主不在于结果，而在于如何获得结果。

一般的优化软件都没有公开其算法，如1stOpt、matlab、lingo等。建议从网上公开发布的文章、算法中获取一些算法思想。


〃Lan_Pst。 发表于 2011-4-19 10:54

forcal 发表于 2011-4-19 08:15 static/image/common/back.gif
用Forcal求解，有2组解：

2组解（A, B, C, D，误差）：


是的，因为需要分析其计算过程。
还是非常感谢。。。

243659583 发表于 2011-4-19 14:03

baozheng{:3_42:}{:3_42:}{:3_42:}{:3_42:}{:3_42:}

plmlor 发表于 2011-10-30 08:22

很好！
顶顶顶！

远东森英—咱北方人的地板

bozsu 发表于 2012-1-10 12:09

祝楼主快乐,大家一齐讨论！


这里有木有买双色球的朋友啊，我有方法能提高中奖率，

想了解或者感兴趣的朋友加扣扣：2590587466

wdgqr 发表于 2012-2-25 16:56

真不知道，还有这么好心的楼主 支持











static/image/common/sigline.gif
女裤品牌
淘宝网女裤
当当购书网
一号店团购
haoping66.com
优尚网
页: [1]
查看完整版本: 四参数拟合方程 的方法