我的新书《基于股票大数据分析的Python入门实战》于近日上架,在这篇博文向大家介绍我的新书:《基于股票大数据分析的Python入门实战》里,介绍了这本书的内容。这里将摘录出部分内容,用以推广本书,请大家多多支持。
《基于股票大数据分析的Python入门实战 视频教学版》,京东链接:https://item.jd.com/69241653952.html。在这篇博文里,将摘录部分内容。
1 OBV指标的原理以及算法
形象地讲,OBV指标是将成交量与股价的关系数字化,并根据股市的成交量变化情况来衡量股市的上涨或下跌支持力,以此来研判股价的走势。OBV指标的设计是基于如下的原理。
第一,如果投资者对当前股价的看法越有分歧,那么成交量就越大,反之成交量就越小,所以可以用成交量来衡量多空双方的力量。
第二,股价在上升时,尤其是在上升初期,必须要较大的成交量相配合,相反股价在下跌时,无需耗费很大的动量,因此成交量未必放大,甚至下跌阶段成交量会有萎缩趋势。
第三,受关注的股票在一段时间内成交量和股价波动会很大,相反冷门股票成交量和价格波动会比较小。
根据上述原则,OBV的算法如下,主要是以日为单位累积成交量。
当日OBV值 = 本日值 + 前日OBV值
如果本日收盘价高于前一日的收盘价,本日的值为正,反之为负,如果本日收盘价和前一日的收盘价相同,则本日值不参与计算,按照这种规则累积计算成交量。成交量可以选择多种计算单位,OBV用到的是成交手数。参考表12-3,通过范例来了解一下OBV的算法。
表 OBV指标算法的实例表
日期 |
收盘价(元) |
成交量(手) |
当日OBV累计值 |
第1天 |
10 |
10000 |
不计算 |
第2天 |
10.2 |
+11000 |
+11000 |
第3天 |
10.3 |
+12000 |
+23000 |
第4天 |
10.2 |
-10000 |
+13000 |
第5天 |
10.1 |
-5000 |
+8000 |
其中,第一天不计算,第2天的收盘价高于第1天,所以当日OBV是当日成交量(为正数)。第3天收盘价也高于第2天,所以该日的OBV是第2天的值(+11000)加上该日成交量(+12000)。
第4天股票下跌,所以当日的OBV累计值是前日的23000减去当日的成交量,结果是+13000,同理第5天也是下跌,当日的OBV是前日值13000减去当日成交量5000,结果是8000。
之后的OBV值按同理计算,将每日算得的OBV值作为纵坐标,交易的日期作为横坐标,将这些点连接起来就是OBV指标线了。
2 绘制K线、均线和OBV指标图的整合效果图
在绘制K线、均线与OBV指标图时,是从csv文件(其实源于网站爬取的股票交易数据)中的Volume字段获得的成交量,它的单位是“股数”,而计算OBV时成交量的单位是“手”,两者的对应关系是1手等于100股。
在DrawKwithOBV.py范例程序中,将绘制整合的效果图,该范例程序存放在MyDjangoDBProj项目中,与DBUtil.py处于同一目录。为了突出OBV算法,范例程序不导入数据库相关的操作,也不输出日志。
1 # !/usr/bin/env python 2 # coding=utf-8 3 import pandas as pd 4 import matplotlib.pyplot as plt 5 from mpl_finance import candlestick2_ochl 6 # 计算OBV的方法 7 def calOBV(df): 8 # 把成交量换算成万手 9 df['VolByHand'] = df['Volume']/1000000 10 # 创建OBV列,先全填充为0 11 df['OBV'] =0 12 cnt=1 # 索引从1开始,即从第2天算起 13 while cnt<=len(df)-1: 14 if(df.iloc[cnt]['Close']>df.iloc[cnt-1]['Close']): 15 df.ix[cnt,'OBV'] = df.ix[cnt-1,'OBV'] + df.ix[cnt,'VolByHand'] 16 if(df.iloc[cnt]['Close']<df.iloc[cnt-1]['Close']): 17 df.ix[cnt,'OBV'] = df.ix[cnt-1,'OBV'] - df.ix[cnt,'VolByHand'] 18 cnt=cnt+1 19 return df
在第7行的calOBV方法中封装了计算OBV指标的程序逻辑。具体执行步骤是,在第9行中为df对象新增VolByHand列,把成交量转换成“万手”,虽然OBV的计算单位是手,但以此绘制出来的指标图上y轴的OBV数值还是过大,所以这里在除以100的基础上再除以10000,转换成“万手”。
随后在第11行新增OBV列,该列的初始值是0。之后在第13行的while循环中,从第2天开始依次遍历df对象,根据OBV的计算规则给每天的OBV列赋值,比如通过第14行的if语句处理当天收盘价上涨的情况,从第15行的程序代码中可以看到,在上涨情况下,当日的OBV值是前日OBV值加上当日的成交量,在第17行中处理了当日下跌的情况,当日的OBV值是前日值减去当日的成交量。
20 filename='D:\\stockData\ch12\\6004602019-01-012019-05-31.csv' 21 df = pd.read_csv(filename,encoding='gbk') 22 # 调用方法计算OBV 23 df = calOBV(df) 24 # print(df) # 可以去除这段注释以查看结果
在第21行从指定的csv文件中读到600460(士兰微)从20190101到20190531的交易数据,并在第23行调用calOBV方法计算OBV值,在该方法的返回结果存放到df对象中,其中OBV值包含在df['OBV']这一列种。如果要检验计算的OBV结果,可以去掉第24行的注释,使得打印语句生效。
25 figure = plt.figure() 26 # 创建子图 27 (axPrice, axOBV) = figure.subplots(2, sharex=True) 28 # 调用方法,在axPrice子图中绘制K线图 29 candlestick2_ochl(ax = axPrice, 30 opens=df["Open"].values, closes=df["Close"].values, 31 highs=df["High"].values, lows=df["Low"].values, 32 width=0.75, colorup='red', colordown='green') 33 axPrice.set_title("K线图和均线图") # 设置子图标题 34 df['Close'].rolling(window=3).mean().plot(ax=axPrice,color="red",label='3日均线') 35 df['Close'].rolling(window=5).mean().plot(ax=axPrice,color="blue",label='5日均线') 36 df['Close'].rolling(window=10).mean().plot(ax=axPrice,color="green",label='10日均线') 37 axPrice.legend(loc='best') # 绘制图例 38 axPrice.set_ylabel("价格(单位:元)") 39 axPrice.grid(linestyle='-.') # 带网格线 40 # 在axOBV子图中绘制OBV图形 41 df['OBV'].plot(ax=axOBV,color="blue",label='OBV') 42 plt.legend(loc='best') # 绘制图例 43 plt.rcParams['font.sans-serif']=['SimHei'] 44 # 在OBV子图上加上负值效果 45 plt.rcParams['axes.unicode_minus'] = False 46 axOBV.set_ylabel("单位:万手") 47 axOBV.set_title("OBV指标图") # 设置子图的标题 48 axOBV.grid(linestyle='-.') # 带网格线 49 # 设置x轴坐标的标签和旋转角度 50 major_index=df.index[df.index%5==0] 51 major_xtics=df['Date'][df.index%5==0] 52 plt.xticks(major_index,major_xtics) 53 plt.setp(plt.gca().get_xticklabels(), rotation=30) 54 plt.show()
在第27行的程序语句设置了两个子图,其中axPrice用于绘制K线和均线,而axOBV则用于绘制OBV指标图。
从第29行到第39行的程序语句用于绘制K线以及三条均线,这部分代码在之前几章中的范例程序中都讲过,所以不再重复说明。在第41行中通过调用df['OBV'].plot方法绘制OBV指标图。
在绘制OBV子图时请注意两个细节:
第一,在第46行中,在axOBV子图内通过调用set_ylabel方法设置了OBV子图的y坐标标签为“万手”;
第二,通过第45行的程序代码,让OBV子图上的y坐标数字有正有负,如果去掉这行语句,OBV子图上y坐标的数字均为正数。
运行这个范例程序,即可看到如图所示的执行结果。
本书的其它内容,大家可以参考如下的博文,