MeteoInfoLab脚本示例：曲线拟合

MeteoInfo

MeteoInfoLab中mipylib.numeric中增加了fitting包，用来进行曲线拟合，目前有多项式、幂函数、指数函数拟合功能。

多项式拟合：

1. from mipylib.numeric import fitting
   
2. 
   
3. x = linspace(0, 4*pi, 10)
   
4. y = sin(x)
   
5. 
   
6. #Plot data points
   
7. plot(x, y, 'ro', fill=False, size=1)
   
8. 
   
9. #Use polyfit to fit a 7th-degree polynomial to the points
   
10. r = fitting.polyfit(x, y, 7)
    
11. 
    
12. #Plot fitting line
    
13. xx = linspace(0, 4*pi, 100)
    
14. p = r[0]
    
15. yy = fitting.polyval(p, xx)
    
16. plot(xx, yy, '-b')
    
17. title('Polynomial fitting example')

from mipylib.numeric import fitting<br /> <br /> x = linspace(0, 4*pi, 10)<br /> y = sin(x)<br /> <br /> #Plot data points<br /> plot(x, y, 'ro', fill=False, size=1)<br /> <br /> #Use polyfit to fit a 7th-degree polynomial to the points<br /> r = fitting.polyfit(x, y, 7)<br /> <br /> #Plot fitting line<br /> xx = linspace(0, 4*pi, 100)<br /> p = r[0]<br /> yy = fitting.polyval(p, xx)<br /> plot(xx, yy, '-b')<br /> title('Polynomial fitting example')

幂函数拟合：

1. from mipylib.numeric import fitting
   
2. 
   
3. fn = 'D:/Temp/ascii/PM&vis-1.txt'
   
4. ncol = numasciicol(fn)
   
5. nrow = numasciirow(fn)
   
6. a = asciiread(fn,shape=(nrow,ncol))
   
7. x=a[:,0]
   
8. y=a[:,1]
   
9. z=a[:,2]
   
10. axes(tickfontsize=17)
    
11. ls=scatter(x,y,s=8,c=z,cmap='NCV_jet',edgecolor=None,cnum=20)
    
12. xlim(0,450)
    
13. ylim(0,30)
    
14. xlabel(r'$\rm{PM_{2.5}} \ (\mu g \ m^{-3})$',fontsize=17)
    
15. ylabel(r'$\rm{Visibility (km)}$',fontname='Arial',fontsize=17)
    
16. colorbar(ls,fontsize=17,label='RH(%)')
    
17. 
    
18. #Pow law fitting
    
19. a,b,r,f = fitting.powerfit(x, y, func=True)
    
20. 
    
21. #Plot fitting line
    
22. xx = linspace(x.min(), x.max(), 100)
    
23. #yy = a*pow(xx, b)
    
24. yy = fitting.predict(f, xx)
    
25. plot(xx, yy, '-b', linewidth=2)
    
26. text(250, 20, r'$y = ' + '%.4f' % a + 'x^{%.4f' % b + '}$', fontsize=16)
    
27. text(250, 18, r'$r^2=%.4f' % r + '$', fontsize=16)

from mipylib.numeric import fitting<br /> <br /> fn = 'D:/Temp/ascii/PM&vis-1.txt'<br /> ncol = numasciicol(fn)<br /> nrow = numasciirow(fn)<br /> a = asciiread(fn,shape=(nrow,ncol))<br /> x=a[:,0]<br /> y=a[:,1]<br /> z=a[:,2]<br /> axes(tickfontsize=17)<br /> ls=scatter(x,y,s=8,c=z,cmap='NCV_jet',edgecolor=None,cnum=20)<br /> xlim(0,450)<br /> ylim(0,30)<br /> xlabel(r'$\rm{PM_{2.5}} \ (\mu g \ m^{-3})$',fontsize=17)<br /> ylabel(r'$\rm{Visibility (km)}$',fontname='Arial',fontsize=17) <br /> colorbar(ls,fontsize=17,label='RH(%)')<br /> <br /> #Pow law fitting<br /> a,b,r,f = fitting.powerfit(x, y, func=True)<br /> <br /> #Plot fitting line<br /> xx = linspace(x.min(), x.max(), 100)<br /> #yy = a*pow(xx, b)<br /> yy = fitting.predict(f, xx)<br /> plot(xx, yy, '-b', linewidth=2)<br /> text(250, 20, r'$y = ' + '%.4f' % a + 'x^{%.4f' % b + '}$', fontsize=16)<br /> text(250, 18, r'$r^2=%.4f' % r + '$', fontsize=16)

zihan6688

MeteoInfo越来越强大了，支持

猫猫2003

老师，你展示下你的PM&vis-1.txt里面的数据格式吗？

MeteoInfo

猫猫2003 发表于 2017-5-31 10:36
老师，你展示下你的PM&vis-1.txt里面的数据格式吗？

182.82         1.1868912        92
193.41         1.207008        83
159.08         2.8565856        78
146.80         3.2790384        71
210.08         2.414016        74
160.63         3.954388114        82
76.85         3.126725486        80
42.48         7.1615808        76
97.75         4.8079152        71
107.90         3.6411408        83
90.60         5.8137552        79
77.82         7.425973029        63
85.05         4.6872144        76
176.53         1.9513296        83
144.23         1.8708624        80
242.38         1.207008        84
262.58         0.9454896        94
...

静言_GRMC

好强大，试验成功~

猫猫2003

MeteoInfo 发表于 2017-5-31 11:03
182.82         1.1868912        92
193.41         1.207008        83
159.08         2.8565856        78

测试成功了，但没图片上没有函数 。 是怎么回事呢？

猫猫2003

MeteoInfo 发表于 2017-5-31 11:03
182.82         1.1868912        92
193.41         1.207008        83
159.08         2.8565856        78

已经明白，谢谢老师。

葱葱葱白

老师您好，我想请问一下如果我想做直线的拟合用什么函数呢

yanpianpian

图画的好棒，感觉好多东西需要学习

MeteoInfo

葱葱葱白 发表于 2019-4-28 22:25
老师您好，我想请问一下如果我想做直线的拟合用什么函数呢

参考这里：http://bbs.06climate.com/forum.p ... &extra=page%3D1