用matplotlib和pandas绘制股票MACD指标图，并验证化交易策略

我的新书《基于股票大数据分析的Python入门实战》于近日上架，在这篇博文向大家介绍我的新书：《基于股票大数据分析的Python入门实战》里，介绍了这本书的内容。这里将摘录出部分内容，用以推广本书，请大家多多支持。

1 MACD指标的计算方式

从数学角度来分析，MACD指标是根据均线的构造原理，对股票收盘价进行平滑处理，计算出算术平均值以后再进行二次计算，它是属于趋向类指标。

MACD指标是由三部分构成的，分别是：DIF（离差值，也叫差离值）、DEA（离差值平均）和BAR（柱状线）。

具体的计算过程是，首先算出快速移动平均线（EMA1）和慢速移动平均线（EMA2），用这两个数值来测量两者间的差离值（DIF），在此基础上再计算差离值（DIF）N周期的平滑移动平均线DEA（也叫MACD、DEM）线。

如前文所述，EMA1周期参数一般取12日，EMA2一般取26日，而DIF一般取9日，在此基础上，MACD指标的计算步骤如下所示。

第一步：计算移动平均值（即EMA）。

12日EMA1的计算方式是：EMA（12）= 前一日EMA（12）× 11/13 ＋今日收盘价 × 2/13

26日EMA2的计算方式是：EMA（26）= 前一日EMA（26）× 25/27 ＋今日收盘价 ×2 /27

第二步：计算MACD指标中的差离值（即DIF）。

DIF = 今日EMA（12）－今日EMA（26）

第三步：计算差离值的9日EMA（即MACD指标中的DEA）。用差离值计算它的9日EMA，这个值就是差离平均值（DEA）。

今日DEA（MACD）= 前一日DEA × 8/10 ＋今日DIF × 2/10

第四步：计算BAR柱状线。

BAR = 2 × (DIF － DEA)

这里乘以2的原因是，在不影响趋势的情况下，从数值上扩大DIF和DEA差值，这样观察效果就更加明显。

最后，把各点（即每个交易日）的DIF值和DEA值连接起来，就能得到在x轴上下移动的两条线，分别表示短期（即快速，EMA1，周期是12天）和长期（即慢速，EMA2，周期是26天）。而且，DIF和DEA的离差值能构成红、绿两种颜色的柱状线，在x轴之上是红色，而x轴之下是绿色。

2 遍历数据表数据，绘制MACD指标

同K线指标一样，根据不同的计算周期， MACD指标也可以分为日指标、周指标、月指标乃至年指标。在下面的DrawMACD.py范例程序中将绘制日MACD指标，在这个范例程序中可以看到关于数据结构、图形绘制和数据库相关的操作，由于程序代码比较长，下面分段讲解。

1    # !/usr/bin/env python
2    # coding=utf-8
3    import pandas as pd
4    import matplotlib.pyplot as plt
5    import pymysql
6    import sys
7    # 第一个参数是数据，第二个参数是周期
8    def calEMA(df, term):
9        for i in range(len(df)):
10           if i==0: # 第一天
11               df.ix[i,\'EMA\']=df.ix[i,\'close\']
12           if i>0:
13              df.ix[i,\'EMA\']=(term-1)/(term+1)*df.ix[i-1,\'EMA\']+2/(term+1) * df.ix[i,\'close\']
14       EMAList=list(df[\'EMA\'])
15       return EMAList

在第8行到第15行的calEMA方法中，根据第二个参数term，计算快速（周期是12天）和慢速（周期是26天）的EMA值。

具体步骤是，通过第9行的for循环，遍历由第一个参数指定的DataFrame类型的df对象，根据第10行的if条件中，如果是第一天，则EMA值用当天的收盘价，如果满足第12行的条件，即不是第一天，则在第13行中根据8.3.1节的算法，计算当天的EMA值。

请注意，在第11行和第13行中是通过df.ix的形式访问索引行（比如第i行）和指定标签列（比如EMA列）的数值，ix方法与之前loc以及iloc方法不同的是，ix方法可以通过索引值和标签值访问，而loc以及iloc方法只能通过索引值来访问。计算完成后，在第14行把df的EMA列转换成列表类型的对象并在第15行返回。

16   # 定义计算MACD的方法
17   def calMACD(df, shortTerm=12, longTerm=26, DIFTerm=9):
18       shortEMA = calEMA(df, shortTerm)
19       longEMA = calEMA(df, longTerm)
20       df[\'DIF\'] = pd.Series(shortEMA) - pd.Series(longEMA)
21       for i in range(len(df)):
22           if i==0:         # 第一天
23               df.ix[i,\'DEA\'] = df.ix[i,\'DIF\']  # ix可以通过标签名和索引来获取数据
24           if i>0:
25                 df.ix[i,\'DEA\'] = (DIFTerm-1)/(DIFTerm+1)*df.ix[i-1,\'DEA\'] + 2/(DIFTerm+1)*df.ix[i,\'DIF\']
26       df[\'MACD\'] = 2*(df[\'DIF\'] - df[\'DEA\'])
27       return df[[\'date\',\'DIF\',\'DEA\',\'MACD\']]
28       # return df

在第15行到第27行定义的calMACD方法中，将调用第8行定义的calEMA方法来计算MACD的值。具体步骤是，在第18行和第19行通过调用calEMA方法，分别得到了快速和慢速的EMA值，在第20行，用这两个值计算DIF值。请注意，shortEMA和longEMA都是列表类型，所以可以像第20行那样，通过调用pd.Series方法把它们转换成Series类对象后再直接计算差值。

从第21行到第25行的程序语句，也是根据8.3.1节给出的公式计算DEA值，同样要用两条if语句区分“第一天”和“以后几天”这两种情况，在第26行根据计算公式算出MACD的值。

第27行返回指定的列，在后面的代码中还要用到df对象的其他列，此时则可以用如第28行所示的代码返回df的全部列。

29   try:
30       # 打开数据库连接
31       db = pymysql.connect("localhost","root","123456","pythonStock" )
32   except:
33       print(\'Error when Connecting to DB.\')
34       sys.exit()
35   cursor = db.cursor()
36   cursor.execute("select * from stock_600895")
37   cols = cursor.description        # 返回列名
38   heads = []
39   # 依次把每个cols元素中的第一个值放入col数组
40   for index in cols:
41       heads.append(index[0])
42   result = cursor.fetchall()
43   df = pd.DataFrame(list(result))
44   df.columns=heads
45   # print(calMACD(df, 12, 26, 9)) # 输出结果
46   stockDataFrame = calMACD(df, 12, 26, 9)

从第29行到第35行的程序语句，建立了MySQL数据库的连接和获得游标cursor对象，在第36行中，通过select类型的SQL语句，来获取stock_600895表中的所有数据，如8.2节所述，这个数据表中的数据源自雅虎网站。

在第37行中，得到了stock_600895数据表的字段列表。在第40行和第41行的for循环中，把字段列表中的第0行索引元素放入了heads。在第42行和第43行，把从stock_600895数据表中获取的数据放入到df对象。在第44行的程序语句，把包含数据表字段列表的heads对象赋值给df对象的字段。

执行到这里，如果去掉第45行打印语句的注解，就能看到第一列输出的是字段名列表，之后会按天输出与MACD有关的股票指标数据。

在第46行调用了calMACD方法，并把结果赋值给stockDataFrame对象，之后就可以根据stockDataFrame对象中的值开始绘图。

47   # 开始绘图
48   plt.figure()
49   stockDataFrame[\'DEA\'].plot(color="red",label=\'DEA\')
50   stockDataFrame[\'DIF\'].plot(color="blue",label=\'DIF\')
51   plt.legend(loc=\'best\')    # 绘制图例
52   # 设置MACD柱状图
53   for index, row in stockDataFrame.iterrows():
54       if(row[\'MACD\'] >0):   # 大于0则用红色
55           plt.bar(row[\'date\'], row[\'MACD\'],width=0.5, color=\'red\')
56       else:                               # 小于等于0则用绿色
57           plt.bar(row[\'date\'], row[\'MACD\'],width=0.5, color=\'green\')
58   # 设置x轴坐标的标签和旋转角度
59   major_index=stockDataFrame.index[stockDataFrame.index%10==0]
60   major_xtics=stockDataFrame[\'date\'][stockDataFrame.index%10==0]
61   plt.xticks(major_index,major_xtics)
62   plt.setp(plt.gca().get_xticklabels(), rotation=30)
63   # 带网格线，且设置了网格样式
64   plt.grid(linestyle=\'-.\')
65   plt.title("600895张江高科的MACD图")
66   plt.rcParams[\'axes.unicode_minus\'] = False
67   plt.rcParams[\'font.sans-serif\']=[\'SimHei\']
68   plt.show()

在第49和第50行中通过调用plot方法，以折线的形式绘制出DEA和DIF两根线，在第51行中设置了图例。

在第53行到第57行的for循环中，以柱状图的形式依次绘制了每天的MACD值的柱状线，这里用第54行和第56行的if…else语句进行区分，如果row[\’MACD\’]大于0，则MACD柱是红色，反之是绿色。

从第59行到第61行的程序语句设置了x轴的标签，如果显示每天的日期，那么x轴上的文字会过于密集，所以在第59行和第60行进行相应的处理，只显示stockDataFrame.index%10==0（即索引值是10的倍数）的日期。

在第62行设置了x轴文字的旋转角度，在第64行设置了网格的式样，在第65行设置了标题文字，最后在第68行通过调用show方法绘制了整个图形。

请注意，如果不编写第66行的程序语句，那么y轴标签值里的负号就不会显示，这样显示结果就不正确了。读者可以把这条语句注释掉后，再运行一下，看看结果如何。

至此，我们实现了计算并绘制MACD指标线的功能，读者应该掌握了如何获得指定股票在指定时间段内的交易数据，而后可以稍微改写上述的范例程序，绘制出其他股票在指定时间范围内的MACD走势图。

3 MACD与K线均线的整合效果图

MACD是趋势类指标，如果把它与K线和均线整合到一起的话，就能更好地看出股票走势的“趋势性”。在下面的DrawKwithMACD.py范例程序中示范了整合它们的效果，由于程序代码比较长，因而在下面的分析中略去了一些之前分析过的重复代码，读者可以从本书提供下载的范例程序中看到完整的代码。

1    # !/usr/bin/env python
2    # coding=utf-8
3    import pandas as pd
4    import matplotlib.pyplot as plt
5    import pymysql
6    import sys
7    from mpl_finance import candlestick2_ochl
8    from matplotlib.ticker import MultipleLocator
9    # 计算EMA的方法，第一个参数是数据，第二个参数是周期
10   def calEMA(df, term):
11       # 省略具体实现，请参考本书提供下载的完整范例程序
12   # 定义计算MACD的方法
13   def calMACD(df, shortTerm=12, longTerm=26, DIFTerm=9):
14       # 省略中间的计算过程，请参考本书提供下载的完整范例程序
15       return df

从第3行到第8行的程序语句通过import语句导入了必要的依赖包，第10行定义的calEMA方法和DrawMACD.py范例程序中的完全一致，所以就省略了该方法内部的代码。第13行定义计算MACD的calMACD方法和DrawMACD.py范例程序中的同名方法也完全一致，但在最后的第15行，是通过return语句返回整个df对象，而不是返回仅仅包含MACD指标的相关列，这是因为，在后文中需要股票的开盘价等数值来绘制K线图。

16   try:
17       # 打开数据库连接
18       db = pymysql.connect("localhost","root","123456","pythonStock" )
19   except:
20       print(\'Error when Connecting to DB.\')
21       sys.exit()
22   cursor = db.cursor()
23   cursor.execute("select * from stock_600895")
24   cols = cursor.description        # 返回列名
25   heads = []
26   # 依次把每个cols元素中的第一个值放入col数组
27   for index in cols:
28       heads.append(index[0])
29   result = cursor.fetchall()
30   df = pd.DataFrame(list(result))
31   df.columns=heads
32   # print(calMACD(df, 12, 26, 9)) # 输出结果
33   stockDataFrame = calMACD(df, 12, 26, 9)

从第16行到第33行的程序语句把需要的数据放入了stockDataFrame这个DataFrame类型的对象中，之后就可以根据其中的数据画图了，这段程序代码之前分析过，就不再重复讲述了。

34   # 开始绘图，设置大小，共享x坐标轴
35   figure,(axPrice, axMACD) = plt.subplots(2, sharex=True, figsize=(15,8))
36   # 调用方法绘制K线图
37   candlestick2_ochl(ax = axPrice, opens=stockDataFrame["open"].values, closes = stockDataFrame["close"].values, highs=stockDataFrame["high"].values,  lows = stockDataFrame["low"].values, width=0.75, colorup=\'red\', colordown=\'green\')
38   axPrice.set_title("600895张江高科K线图和均线图")   # 设置子图的标题
39   stockDataFrame[\'close\'].rolling(window=3).mean().plot(ax=axPrice,color="red",label=\'3日均线\')
40   stockDataFrame[\'close\'].rolling(window=5).mean().plot(ax=axPrice,color="blue",label=\'5日均线\')
41   stockDataFrame[\'close\'].rolling(window=10).mean().plot(ax=axPrice,color="green",label=\'10日均线\')
42   axPrice.legend(loc=\'best\')       # 绘制图例
43   axPrice.set_ylabel("价格（单位：元）")
44   axPrice.grid(linestyle=\'-.\')     # 带网格线

从第34行到第44行的程序语句绘制了指定时间范围内“张江高科”股票的K线图和均线，这部分代码和第7章drawKMAAndVol.py范例程序中实现同类功能的代码很相似，有差别的是在第35行，第二个子图的名字设置为“axMACD”，在第44行中通过linestyle设置了网格线的样式。

45   # 开始绘制第二个子图
46   stockDataFrame[\'DEA\'].plot(ax=axMACD,color="red",label=\'DEA\')
47   stockDataFrame[\'DIF\'].plot(ax=axMACD,color="blue",label=\'DIF\')
48   plt.legend(loc=\'best\')    # 绘制图例
49   # 设置第二个子图中的MACD柱状图
50   for index, row in stockDataFrame.iterrows():
51       if(row[\'MACD\'] >0):   # 大于0则用红色
52           axMACD.bar(row[\'date\'], row[\'MACD\'],width=0.5, color=\'red\')
53       else:                               # 小于等于0则用绿色
54           axMACD.bar(row[\'date\'], row[\'MACD\'],width=0.5, color=\'green\')
55   axMACD.set_title("600895张江高科MACD")   # 设置子图的标题
56   axMACD.grid(linestyle=\'-.\')                # 带网格线
57   # xmajorLocator = MultipleLocator(10)    # 将x轴的主刻度设置为10的倍数
58   # axMACD.xaxis.set_major_locator(xmajorLocator)
59   major_xtics=stockDataFrame[\'date\'][stockDataFrame.index%10==0]
60   axMACD.set_xticks(major_xtics)
61   # 旋转x轴显示文字的角度
62   for xtick in axMACD.get_xticklabels():
63       xtick.set_rotation(30)
64   plt.rcParams[\'font.sans-serif\']=[\'SimHei\']
65   plt.rcParams[\'axes.unicode_minus\'] = False
66   plt.show()

在上述程序代码中，在axMACD子图内绘制了MACD线，由于是在子图内绘制，因此在第46行和第47行绘制DEA和DIF折线的时候，需要在参数里通过“ax=axMACD”的形式指定所在的子图。

在第59行和第60行中设置了axMACD子图中的x轴标签，由于在第35行中设置了axPrice和axMACD两子图是共享x轴，因此K线和均线所在子图的x轴刻度会和MACD子图中的一样。因为是在子图中，所以需要通过第62行和第63行的for循环依次旋转x轴坐标的标签文字。

在这段代码中其实给出了两种设置x轴标签的方式。如果注释掉第59行和第60行的代码，并去掉第57行和第58行的注释，会发现效果是相同的。

需要说明的是，虽然在第57行和第59行的代码中并没有指定标签文字，但在第37行调用candlestick2_ochl方法绘制K线图时，会设置x轴的标签文字，所以依然能看到x轴上日期的标签。运行这个范例程序后，结果如图所示。

4 MACD指标的指导意义

根据MACD各项指标的含义，可以通过DIF和DEA两者的值、DIF和DEA指标的交叉情况（比如金叉或死叉）以及BAR柱状图的长短与收缩的情况来判断当前股票的趋势。

如下两点是根据DIF和DEA的数值情况以及它们在x轴上下的位置来确定股票的买卖策略。

第一，当DIF和DEA两者的值均大于0（在x轴之上）并向上移动时，一般表示当前处于多头行情中，建议可以买入。反之，当两者的值均小于0且向下移动时，一般表示处于空头行情中，建议卖出或观望。

第二，当DIF和DEA的值均大于0但都在向下移动时，一般表示为上涨趋势即将结束，建议可以卖出股票或观望。同理，当两者的值均小于0，但在向上移动时，一般表示股票将上涨，建议可以持续关注或买进。

如下四点是根据DIF和DEA的交叉情况来决定买卖策略。

第一，DIF与DEA都大于0而且DIF向上突破DEA时，说明当前处于强势阶段，股价再次上涨的可能性比较大，建议可以买进，这就是所谓MACD指标黄金交叉，也叫金叉。

第二，DIF与DEA都小于0，但此时DIF向上突破DEA时，表明股市虽然当前可能仍然处于跌势，但即将转强，建议可以开始买进股票或者重点关注，这也是MACD金叉的一种形式。

第三，DIF与DEA虽然都大于0，但而DIF却向下突破DEA时，这说明当前有可能从强势转变成弱势，股价有可能会跌，此时建议看机会就卖出，这就是所谓MACD指标的死亡交叉，也叫死叉。

第四，DIF和DEA都小于0，在这种情况下又发生了DIF向下突破DEA的情况，这说明可能进入下一阶段的弱势中，股价有可能继续下跌，此时建议卖出股票或观望，这也是MACD死叉的一种形式。

如下两点是根据MACD中BAR柱状图的情况来决定买卖策略。

第一，红柱持续放大，这说明当前处于多头行情中，此时建议买入股票，直到红柱无法再进一步放大时才考虑卖出。相反，如果绿柱持续放大，这说明当前处于空头行情中，股价有可能继续下跌，此时观望或卖出，直到绿柱开始缩小时才能考虑买入。

第二，当红柱逐渐消失而绿柱逐渐出现时，这表明当前的上涨趋势即将结束，有可能开始加速下跌，这时建议可以卖出股票或者观望。反之，当绿柱逐渐消失而红柱开始出现时，这说明下跌行情即将或已经结束，有可能开始加速上涨，此时可以开始买入。

虽然说MACD指标对趋势的分析有一定的指导意义，但它同时也存在一定的盲点。

比如，当没有形成明显的上涨或下跌趋势时（即在盘整阶段），DIF和DEA这两个指标会频繁地出现金叉和死叉的情况，这时由于没有形成趋势，因此金叉和死叉的指导意义并不明显。

又如，MACD指标是对趋势而言的，从中无法看出未来时间段内价格上涨和下跌的幅度。比如在图8-11中，股票“张江高科”在价格高位时，DIF的指标在2左右，但有些股票在高位时，DIF的指标甚至会超过5。

也就是说，无法根据DIF和DEA数值的大小来判断股价会不会进一步涨或进一步跌。有时看似DIF和DEA到达一个高位，但如果当前上涨趋势强劲，股价会继续上涨，同时这两个指标会进一步上升，反之亦然。

因此，在实际使用中，投资者可以用MACD指标结合其他技术指标，比如之前提到的均线，从而能对买卖信号进行多重确认。

5 验证基于柱状图和金叉的买点

在CalBuyPointByMACD.py范例程序中将根据如下原则来验证买点：DIF向上突破DEA（出现金叉），且柱状图在x轴上方（即当前是红柱状态）。

在这个范例程序中，用的是股票“金石资源（代码为603505）从2018年9月到2019年5月的交易数据，程序代码如下。

1    # !/usr/bin/env python
2    # coding=utf-8
3    import pandas as pd
4    import pymysql
5    import sys
6    # 第一个参数是数据，第二个参数是周期
7    def calEMA(df, term):
8      # 省略方法内的代码，请参考本书提供下载的完整范例程序
9    # 定义计算MACD的方法
10   def calMACD(df, shortTerm=12, longTerm=26, DIFTerm=9):
11       # 省略中间计算过程的代码，最后返回的是df，请参考本书提供下载的完整范例程序
12       return df

上述代码的calEMA和calMACD方法和8.3.4节的范例程序中的代码完全一致，所以就不再重复讲述了。

13	def getMACDByCode(code):
14	    try:
15	        # 打开数据库连接
16	        db = pymysql.connect("localhost","root","123456","pythonStock" )
17	    except:
18	        print(\'Error when Connecting to DB.\')
19	        sys.exit()  
20	    cursor = db.cursor()
21	    cursor.execute(\'select * from stock_\'+code)
22	    cols = cursor.description  	# 返回列名
23	    heads = []
24	    # 依次把每个cols元素中的第一个值放入col数组
25	    for index in cols:
26	        heads.append(index[0])
27	    result = cursor.fetchall()
28	    df = pd.DataFrame(list(result))
29	    df.columns=heads
30	    stockDataFrame = calMACD(df, 12, 26, 9)
31	    return stockDataFrame

第13行开始的getMACDByCode方法中包含了从数据表中获取的股票交易数据并返回MACD指标的代码，这部分程序代码与之前DrawKwithMACD.py范例程序中的程序也非常相似，只不过在第21行中是根据股票代码来动态地拼接select语句。该方法在第31行中返回包含MACD指标的stockDataFrame对象。

32	# print(getMACDByCode(\'603505\')) 	# 可去除这条语句的注解以确认数据
33	stockDf = getMACDByCode(\'603505\')
34	cnt=0
35	while cnt<=len(stockDf)-1:
36	    if(cnt>=30):		# 前几天有误差，从第30天算起
37	        try:
38	            # 规则1：这天DIF值上穿DEA
39	            if stockDf.iloc[cnt][\'DIF\']>stockDf.iloc[cnt][\'DEA\'] and stockDf.iloc[cnt-1][\'DIF\']<stockDf.iloc[cnt-1][\'DEA\']:
40	                #规则2：出现红柱，即MACD值大于0
41	                if stockDf.iloc[cnt][\'MACD\']>0:
42	                    print("Buy Point by MACD on:" + stockDf.iloc[cnt][\'date\'])
43	        except:
44	            pass
45	    cnt=cnt+1

如果去掉第32行打印语句的注释，执行后就能确认数据。在第35行到第45行的while循环中，依次遍历了每个交易日的数据。有数据计算的误差，所以在这个范例程序中通过第36行的if语句排除了刚开始29天的数据，从第30天算起。

在第39行的if条件语句中制定了第一个规则，前一个交易日的DIF小于DEA，而且当天DIF大于DEA，即出现上穿金叉的现象。在第41行的if条件语句中制定了第二个规则，即出现金叉的当日，MACD指标需要大于0，即当前BAR柱是红柱状态。运行这个范例程序之后，就能看到如下输出的买点。

Buy Point by MACD on:2018-10-31

Buy Point by MACD on:2019-01-09

Buy Point by MACD on:2019-03-18

Buy Point by MACD on:2019-04-04

Buy Point by MACD on:2019-04-19

下面改写一下之前的范例程序，把股票代码改成603505，把股票名称改为“金石资源”，运行后即可看到如图所示的结果图。

根据图中的价格走势，在表中列出了各买点的确认情况。

表基于MACD得到的买点情况确认表

买点	对买点的分析	正确性
2018-10-31	该日出现DIF金叉，且Bar已经在红柱状态，后市有涨。	正确
2019-01-09	该日出现DIF金叉，且Bar柱开始逐渐变红，后市有涨。	正确
2019-03-18	该日虽然出现金叉，Bar柱也开始变红，但之后几天Bar交替出现红柱和绿柱情况，后市在下跌后，出现上涨情况。	不明确
2019-04-04	该日在出现金叉的同时，Bar柱由绿转红。但之后若干交易日后出现死叉，且Bar柱又转绿，后市下跌。	不正确
2019-04-19	出现金叉，且Bar柱由绿柱一下子变很长，后市有涨。	正确

根据这个范例程序的运行结果，可以得到的结论是：通过MACD指标的确能算出买点，但之前也说过，MACD有盲点，在盘整阶段，趋势没有形成时，此时金叉的指导意义就不是很明显，甚至是错误的。

6 验证基于柱状图和死叉的卖点

参考MACD指标，与8.4.2节描述的情况相反，如果出现如下情况，则可以卖出股票：DIF向下突破DEA（出现死叉），且柱状图向下运动（红柱缩小或绿柱变长）。下面通过股票“士兰微”（代码为600460）从2018年9月到2019年5月的交易数据来验证卖点。

先来做如下的准备工作：在MySQL的pythonStock数据库中创建stock_600460数据表，在之前介绍的InsertDataFromYahoo.py范例程序中，把股票代码改为600460，运行后即可在stock_600460数据表中看到指定时间范围内的交易数据。

验证MACD指标卖点的CalSellPointByMACD.py范例程序与之前CalBuyPointByMACD.py范例程序很相似，下面只分析不同的程序代码部分。

1    # !/usr/bin/env python
2    # coding=utf-8
3    import pandas as pd
4    import pymysql
5    import sys
6    # calEMA方法中的代码没有变
7    def calEMA(df, term):
8        # 省略方法内的程序代码，请参考本书提供下载的完整范例程序
9    # 定义计算MACD的方法内的程序代码也没有变
10   def calMACD(df, shortTerm=12, longTerm=26, DIFTerm=9):
11       # 省略方法内的程序代码，请参考本书提供下载的完整范例程序
12   def getMACDByCode(code):
13       # 和CalBuyPointByMACD.py范例程序中的程序代码一致
14   stockDf = getMACDByCode(\'600460\')
15   cnt=0
16   while cnt<=len(stockDf)-1:
17       if(cnt>=30):         # 前几天有误差，从第30天算起
18           try:
19               # 规则1：这天DIF值下穿DEA
20               if stockDf.iloc[cnt][\'DIF\']<stockDf.iloc[cnt][\'DEA\'] and stockDf.iloc[cnt-1][\'DIF\']>stockDf.iloc[cnt-1][\'DEA\']:
21                   # 规则2：Bar柱是否向下运动
22                   if stockDf.iloc[cnt][\'MACD\']<stockDf.iloc[cnt-1][\'MACD\']:
23                       print("Sell Point by MACD on:" + stockDf.iloc[cnt][\'date\'])
24           except:
25               pass
26       cnt=cnt+1

上述代码中的calEMA、calMACD和getMACDByCode三个方法和CalBuyPointByMACD.py范例程序中的代码完全一致，所以本节仅仅是给出了这些方法的定义，不再重复讲述了。

在第14行中通过调用getMACDByCode方法，获取了600460（士兰微）的交易数据，其中包含了MACD指标数据。在第16行到第26行的while循环中通过遍历stockDf对象，计算卖点。

具体的步骤是，通过第17行的if条件语句排除了误差比较大的数据，随后通过第20行的if语句判断当天是否出现了DIF死叉的情况，即前一个交易日的DIF比DEA大，但当前交易日DIF比DEA小。

当满足这个条件时，再通过第22行的if语句判断当天的Bar柱数值是否小于前一天的，即判断Bar柱是否在向下运动。当满足这两个条件时，通过第23行的代码输出建议卖出股票的日期。运行这个范例程序代码后，可看到如下输出的卖点。

Sell Point by MACD on:2018-10-11

Sell Point by MACD on:2018-11-29

Sell Point by MACD on:2018-12-06

Sell Point by MACD on:2019-02-28

Sell Point by MACD on:2019-04-04

再修改前文提到的DrawKwithMACD.py范例程序，把股票代码改为600460，把股票名称改成“士兰微”，运行后即可看到如图所示的结果图。

图股票“士兰微”的K线、均线整合MACD的走势图

再根据价格走势，在表中列出了各卖点的确认情况。

表基于MACD得到的卖点情况确认表

卖点	对卖点的分析	正确性
2018-10-11	1. 该日出现DIF死叉，且DIF和DEA均在x轴下方，Bar由红转绿，且绿柱持续扩大。 2. 虽然能验证该点附近处于弱势，但由于此点已经处于弱势，所以后市价位跌幅不大。	不明确
2018-11-29	1. 在DIF和DEA上行过程中出现死叉。 2. Bar柱从红转绿，后市股价有一定幅度的下跌。	正确
2018-12-06	在11月29日的卖点基础上，再次出现死叉，且Bar柱没有向上运动的趋势，所以进一步确认了弱势行情，果然后市股价有一定幅度的下跌。	正确
2019-02-28	1. 虽然出现死叉，但前后几天DIF和DEA均在向上运动。这说明强势并没有结束。 2. Bar柱虽然变绿，但变绿的幅度非常小。 3. 后市价格不是下跌，而是上涨了。	不正确
2019-04-04	1. DIF和DEA在x轴上方出现死叉，说明强势行情有可能即将结束。 2. Bar柱由红开始转绿。 3. 后市价位出现一波短暂反弹，这可以理解成强势的结束，之后出现下跌，且下跌幅度不小。	正确