请教王老师双X轴的问题

qxsw2016

假设双轴，其底部刻度为['1', '3', '5', '7', '9'] 其顶部刻度对应为其倒数乘以20，即['20.0', '6.7', '4.0', '2.9', '2.2']如果 x1= arange(1,11,2.0) xlbl = ['20.0', '6.7', '4.0', '2.9', '2.2']
那么这条语句xticks(x1,xlbl)，可以标注底部轴刻度值分别为xlbl的值。如何使得 xlbl的值标注到顶部。

MeteoInfo

1. table = readtable('D:/Temp/ascii/temp.csv', delimiter=',', format='%s%f')
   
2. xt = table['temp']
   
3. #print xt
   
4. n = len(xt)
   
5. m = n/2
   
6. ind = arange(m+1)
   
7. ak = zeros(m+1)
   
8. bk = zeros(m+1)
   
9. sksq = zeros(m+1)
   
10. ak[0] = sum(xt)/n
    
11. for k in arange(1,m+1):
    
12.     for t in arange(0,n):   
    
13.         ak[k] = ak[k] + xt[t]*cos((2*pi*k*t)/n)
    
14.         bk[k] = bk[k] + xt[t]*sin((2*pi*k*t)/n)
    
15.     ak[k] = ak[k]*2/n
    
16.     bk[k] = bk[k]*2/n
    
17.     sksq[k] = (ak[k]*ak[k]+bk[k]*bk[k])/2
    
18. #print sksq
    
19. width = 0.1
    
20. ax1=axes()
    
21. y = arange(0,2.0,0.5)
    
22. x1 = arange(1,12,2.0)
    
23. x2 = 20/x1
    
24. ind = ind - width * 0.5
    
25. bar(ind,sksq,width,color='r', edgecolor=None)
    
26. yaxis(tickin=False)
    
27. xaxis(tickin=False)
    
28. yticks(y)
    
29. xticks(x1)
    
30. xlabel('k')
    
31. ylabel(r'$\rm(S^2)
    
32. 
    
33. , fontname='Arial')
    
34. xlim(0,12)
    
35. 
    
36. ax2 = twiny(ax1)
    
37. xaxis(location='top',tickin=False)
    
38. xlabel('T(year)',bold=False)
    
39. xlim(0,12)
    
40. xticks(x1, ['%.1f' % xx for xx in x2])
    
41. title('xx test')
    
42. #antialias(True)

table = readtable('D:/Temp/ascii/temp.csv', delimiter=',', format='%s%f')<br /> xt = table['temp']<br /> #print xt<br /> n = len(xt)<br /> m = n/2<br /> ind = arange(m+1)<br /> ak = zeros(m+1)<br /> bk = zeros(m+1)<br /> sksq = zeros(m+1)<br /> ak[0] = sum(xt)/n<br /> for k in arange(1,m+1):<br />     for t in arange(0,n):    <br />         ak[k] = ak[k] + xt[t]*cos((2*pi*k*t)/n)<br />         bk[k] = bk[k] + xt[t]*sin((2*pi*k*t)/n)<br />     ak[k] = ak[k]*2/n<br />     bk[k] = bk[k]*2/n<br />     sksq[k] = (ak[k]*ak[k]+bk[k]*bk[k])/2<br /> #print sksq<br /> width = 0.1<br /> ax1=axes() <br /> y = arange(0,2.0,0.5)<br /> x1 = arange(1,12,2.0)<br /> x2 = 20/x1<br /> ind = ind - width * 0.5<br /> bar(ind,sksq,width,color='r', edgecolor=None)<br /> yaxis(tickin=False)<br /> xaxis(tickin=False)<br /> yticks(y)<br /> xticks(x1)<br /> xlabel('k')<br /> ylabel(r'$\rm(S^2)<br /> <br /> , fontname='Arial')<br /> xlim(0,12)<br /> <br /> ax2 = twiny(ax1)<br /> xaxis(location='top',tickin=False)<br /> xlabel('T(year)',bold=False)<br /> xlim(0,12) <br /> xticks(x1, ['%.1f' % xx for xx in x2])<br /> title('xx test')<br /> #antialias(True)

软件里有个bug，需要在MeteoInfo QQ群下载axes.py文件，替换 MeteoInfo -> pylib -> mipylib -> plotlib目录中的同名文件。