MeteoInfoLab脚本示例：曲线拟合

diva211

老师，请教一下多项式的拟合方程和R该如何编写呢？我用老师上面提供的脚本拟合了一条2次曲线，但是拟合方程和R不知该怎么表达，附上代码和图

1. from mipylib.numeric import fitting
   
2. 
   
3. fn = 'D:/MI/data/RH-O3.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=16)
    
11. ls=scatter(x,y,s=8,c=z,cmap='NCV_jet',edgecolor=None,cnum=20)
    
12. #xlim(0,300)
    
13. #ylim(0,20)
    
14. xlabel('RH(%)',fontname='Arial',fontsize=15)
    
15. ylabel(r'$\rm{O_{3}} \ (\mu g \ m^{-3})
    
16. ,fontsize=15)
    
17. colorbar(ls,fontsize=15,label=r'$\rm{T} \ (^{\circ} \ C)
    
18. )
    
19. 
    
20. #Pow law fitting
    
21. #a,b,r,f = fitting.powerfit(x, y, func=True)
    
22. r = fitting.polyfit(x,y,2)
    
23. 
    
24. #Plot fitting line
    
25. xx = linspace(x.min(), x.max(), 100)
    
26. p = r[0]
    
27. #yy = a*pow(xx, b)
    
28. #yy = fitting.predict(f, xx)
    
29. yy = fitting.polyval(p,xx)
    
30. plot(xx, yy, '-b', linewidth=2)
    
31. #text(70, 250, r'$y = ' + '%.4f' % a + 'x^{%.4f' % b + '}
    
32. , fontsize=16)
    
33. text(70, 230, r'$r^2=%.4f' % r + '
    
34. , fontsize=16)

复制代码

AX32

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

老师，曲线拟合是讨论pm2.5和能见度的关系吧，为什么有第三个变量呢

半颗麦芽糖

王老师您好，我这个拟合函数有一个变量是在幂上的，这该怎么画呢，式中V（D）和D已知。