如何求这组数据的拟合函数 - V2EX
V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
khowarizmi
V2EX    数学

如何求这组数据的拟合函数

  •  
  •   khowarizmi 2015 年 7 月 10 日 5143 次点击
    这是一个创建于 3861 天前的主题,其中的信息可能已经有所发展或是发生改变。
    给了几组点如下:
    (1,21.9059330191461)
    (2,12.3412139792904)
    (5,6.60182614356538)
    (10,4.6877695119712)
    (25,3.54156118026073)
    (50,3.22440644770007)
    (100,3.08252143576504)
    (250,2.99905966403855)

    画出来是这个样子:


    用什么函数去拟合比较好?

    如果我假设 y = a(x + b)^c + d, 应该如何求这些参数?
    7 条回复    2015-10-21 14:38:16 +08:00
    giskard
        1
    giskard  
       2015 年 7 月 10 日   1
    用origin可以很容易做这样的拟合,或者,也可以用python的scipy包进行最小二乘拟合,比如http://andyhuzhill.github.io/python/scipy/2013/04/10/scipy-leastsq-usage/
    chlx
        2
    chlx  
       2015 年 7 月 10 日
    就最小二乘 + 梯度下降吧
    omengye
        3
    omengye  
       2015 年 7 月 24 日
    matlab算了一下,从图像上看显然是反比例函数比较好一些吧.
    假设方程是y=a/x+b这种形式,可以算出来a=17.7722,b=3.1730,当然这个不一定是最佳的.
    lz可以自行画一下看看效果.
    khowarizmi
        4
    khowarizmi  
    OP
       2015 年 7 月 25 日
    @omengye 我最后用的分段,一部分用 log 拟合,一部分用多项式。
    omengye
        5
    omengye  
       2015 年 7 月 25 日 via Android
    @khowarizmi 哦log 的话配合多项式也是可以的,拟合效果并不能完全决定用什么函数去拟合的,这个得跟实际情况有关联。
    miemiekurisu
        6
    miemiekurisu  
       2015 年 8 月 25 日 via Android   1
    有这么几种思路和方法: 广义线性回归,非线性回归,样条。
    分段处理其实就是样条,数据点少的时候用样条精度略差,比较难把握分段数据点临界点,数据多的时候运算速度快,拟合效果好(因为分段处理了),非线性回归拟合和计算都比较复杂。
    广义线性回归求系数比较容易,系数就是最小二乘解
    onerhao
        7
    onerhao  
       2015 年 10 月 21 日
    数据太少了。不过可以考虑基于核(比如高斯核)的线性回归。
    关于     帮助文档     自助推广系统     博客     API     FAQ     Solana     2617 人在线   最高记录 6679       Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 30ms UTC 04:51 PVG 12:51 LAX 20:51 JFK 23:51
    Do have faith in what you're doing.
    ubao msn snddm index pchome yahoo rakuten mypaper meadowduck bidyahoo youbao zxmzxm asda bnvcg cvbfg dfscv mmhjk xxddc yybgb zznbn ccubao uaitu acv GXCV ET GDG YH FG BCVB FJFH CBRE CBC GDG ET54 WRWR RWER WREW WRWER RWER SDG EW SF DSFSF fbbs ubao fhd dfg ewr dg df ewwr ewwr et ruyut utut dfg fgd gdfgt etg dfgt dfgd ert4 gd fgg wr 235 wer3 we vsdf sdf gdf ert xcv sdf rwer hfd dfg cvb rwf afb dfh jgh bmn lgh rty gfds cxv xcv xcs vdas fdf fgd cv sdf tert sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf shasha9178 shasha9178 shasha9178 shasha9178 shasha9178 liflif2 liflif2 liflif2 liflif2 liflif2 liblib3 liblib3 liblib3 liblib3 liblib3 zhazha444 zhazha444 zhazha444 zhazha444 zhazha444 dende5 dende denden denden2 denden21 fenfen9 fenf619 fen619 fenfe9 fe619 sdf sdf sdf sdf sdf zhazh90 zhazh0 zhaa50 zha90 zh590 zho zhoz zhozh zhozho zhozho2 lislis lls95 lili95 lils5 liss9 sdf0ty987 sdft876 sdft9876 sdf09876 sd0t9876 sdf0ty98 sdf0976 sdf0ty986 sdf0ty96 sdf0t76 sdf0876 df0ty98 sf0t876 sd0ty76 sdy76 sdf76 sdf0t76 sdf0ty9 sdf0ty98 sdf0ty987 sdf0ty98 sdf6676 sdf876 sd876 sd876 sdf6 sdf6 sdf9876 sdf0t sdf06 sdf0ty9776 sdf0ty9776 sdf0ty76 sdf8876 sdf0t sd6 sdf06 s688876 sd688 sdf86