福利到!用 Python 实现海龟交易系统

2016-09-05 14:43:21 +08:00
 mushroomqiu
前言

海龟交易系统本质上是一个趋势跟随的系统,但是最值得我们学习的,是资金管理尤其是分批建仓及动态止损的部分

一、趋势捕捉

唐奇安通道

该指标是有 Richard Donchian 发明的,是有 3 条不同颜色的曲线组成的,该指标用周期(一般都是 20 )内的最高价和最低价来显示市场价格的波动性,当其通道窄时表示市场波动较小,反之通道宽则表示市场波动比较大。 如图所示:
该具体分析为:
当价格冲冲破上轨是就是可能的买的信号;反之,冲破下轨时就是可能的卖的信号。
该指标的计算方法为:

上线=Max (最高低, n )
下线=Min (最低价, n )
中线=(上线+下线)/2
高清源代码请移步: https://uqer.io/community/share/57bd5864228e5b79a575a9b2
海龟交易就是利用唐奇安通道的价格突破来捕捉趋势。
不过我们在向下突破 10 日唐奇安下沿卖出。
二、资金管理

2.1 、 N 值计算

N 值是仓位管理的核心,涉及加仓及止损。另外, N 值与技术指标平均真实波幅 ATR 很相似
首先介绍真实波幅: 真实波幅是以下三个值中的最大值
1 、当前交易日最高价和最低价的波幅
2 、前一交易日的收盘价与当前交易日最高价的波幅
3 、前一交易日的收盘价与当前交易日最低价的波幅
用公式写就是:
TrueRange=Max(High−Low,High−PreClose,PreClose−Low)
接下来, N 值计算公式为:
N=PreN[−19 :]+TrueRange20
其中 preN 为前面 N 值, TrueRange 为当前的真实波幅,此公式的真是含义为计算之前 20 天(包括今天在内)的 N 的平均值
另外,有些海龟交易系统用的是 ATR 来代替 N 值, ATR 为真实波幅的 20 日平均。
2.2 买卖单位及首次建仓
先给出公式:
Unit=1%∗Account/N
首次建仓的时候,当捕捉到趋势,即价格突破唐奇安上轨时,买入 1 个 unit 。
其意义就是,让一个 N 值的波动与你总资金 1%的波动对应,如果买入 1unit 单位的资产,当天震幅使得总资产的变化不超过 1%。例如:
现在你有 10 万元资金, 1%波动就是 1000 元。假如标 X 的 N 值为 0.2 元, 1000 元÷0.2 元=5000 股。也就是说,你的第一笔仓位应该是在其突破上轨(假设为 5 元)时立刻买入 5000 股,耗资 25000 元。
2.3 加仓

若股价在上一次买入(或加仓)的基础上上涨了 0.5N ,则加仓一个 Unit 。
接上面的例子:假如 N 值仍为 0.2 。
价格来到 5 + 0.2*0.5 = 5.1 时,加仓 1 个 Unit ,买入 5000 股,耗资 25500 元,剩余资金 49500 元
价格来到 5.1 + 0.2*0.5 = 5.2 时再加仓 1 个 unit 。买入 5000 股,耗资 26000 元,剩余资金 23500 元
2.4 动态止损

当价格比最后一次买入价格下跌 2N 时,则卖出全部头寸止损。
接上面的例子,最后一次加仓价格为 5.2 。假如此时 N 值 0.2 元。 当价格下跌到 5.2 - 2*0.2 = 4.8 元时,清仓。
持仓成本为 ( 5+5.1+5.2 )*5000/15000 = 5.1 元。 此时亏损 ( 5.1-4.8 )*15000 = 4500 元 对于 10 万来说 这波亏损 4.5%
2.5 止盈

当股价跌破 10 日唐奇安通道下沿,清空头寸结束本次交易
三、代码实现

本代码用 ATR 代替 N 值进行计算,其他逻辑不变:
ATR=MA(TrueRange,20)
我们以单只股票为标,建立海龟交易系统,当然,可以将总资产均分为 n 份,同时交易 n 个标。

计算 ATR 值用日线数据,监控价格突破采用分钟线
0 初始化参数,在 initialize(account)写入
def initialize(account):
account.last_buy_prcie = 0 #上一次买入价
account.hold_flag = False # 是否持有头寸标志
account.limit_unit = 4 # 限制最多买入的单元数
account.unit = 0 # 现在买入 1 单元的股数
1 唐奇安通道计算及判断入场离场:

我们设计个函数,传入值为回测中 account.get_history()取得的某单个股票的历史数据、股票现价、 T 为计算唐奇安通道的数据长度,转化为 dataframe 格式
def IN_OR_OUT(data,price,T):
up = max(data['highPrice'].iloc[-T:])
down = min(data['lowPrice'].iloc[-int(T/2):]) # 这里是 10 日唐奇安下沿
if price>up:
return 1
elif price<down:
return -1
else:
return 0
2. ATR 值计算:


1
def CalcATR(data):
2
TR_List = []
3
for i in range(1,21):
4
TR = max(data['highPrice'].iloc[i]-data['lowPrice'].iloc[i],data['highPrice'].iloc[i]-data['closePrice'].iloc[i-1],data['closePrice'].iloc[i-1]-data['lowPrice'].iloc[i])
5
TR_List.append(TR)
6
ATR = np.array(TR_List).mean()
7
return ATR
3. 计算 unit,注意股数为 100 的整数倍


1
def CalcUnit(perValue,ATR):
2
return int((perValue/ATR)/100)*100
查看全部
4.判断是否加仓或止损:

当价格相对上个买入价上涨 0.5ATR 时,再买入一个 unit
当价格相对上个买入价下跌 2ATR 时,清仓

1
def Add_OR_Stop(price,lastprice,ATR):
2
if price >= lastprice + 0.5*ATR:
3
return 1
4
elif price <= lastprice - 2*ATR:
5
return -1
6
else:
7
return 0
查看全部
5 判断上次卖出操作是否成功(可能出现当日买进,之后却判断需要卖出)


1
def SellComplete(hold_flag,security_position):
2
if len(security_position)>0 and hold_flag==False:
3
return True
4
else:
5
return False
查看全部
构建策略

分钟线回测时间略长啊~
先把上面写的函数集中下,方便微核充启后运行函数

1
################################################### 计算、判断函数 #####################################################################
2
def IN_OR_OUT(data,price,T):
3
up = max(data['highPrice'].iloc[-T:])
4
down = min(data['lowPrice'].iloc[-int(T/2):]) # 这里是 10 日唐奇安下沿
5
if price>up:
6
return 1
7
elif price<down:
8
return -1
9
else:
10
return 0
11

12
def CalcATR(data):
13
TR_List = []
14
for i in range(1,21):
15
TR = max(data['highPrice'].iloc[i]-data['lowPrice'].iloc[i],data['highPrice'].iloc[i]-data['closePrice'].iloc[i-1],data['closePrice'].iloc[i-1]-data['lowPrice'].iloc[i])
16
TR_List.append(TR)
17
ATR = np.array(TR_List).mean()
18
return ATR
19

20
def CalcUnit(perValue,ATR):
21
return int((perValue/ATR)/100)*100
22

23
def Add_OR_Stop(price,lastprice,ATR):
24
if price >= lastprice + 0.5*ATR:
25
return 1
26
elif price <= lastprice - 2*ATR:
27
return -1
28
else:
29
return 0
30

31
def SellComplete(hold_flag,security_position):
32
if len(security_position)>0 and hold_flag==False:
33
return True
34
else:
35
return False
查看全部

1
import numpy as np
2
import pandas as pd
3
from __future__ import division
4
from CAL.PyCAL import *
5
import matplotlib.pyplot as plt
6

7
start = '2012-01-01' # 回测起始时间
8
end = '2016-01-01' # 回测结束时间
9
benchmark = '000001.XSHE'
10
universe = ['000001.XSHE']
11
capital_base = 100000 # 起始资金
12
freq = 'm' # 策略类型,'d'表示日间策略使用日线回测,'m'表示日内策略使用分钟线回测
13
refresh_rate = 1 # 调仓频率,表示执行 handle_data 的时间间隔,若 freq = 'd'时间间隔的单位为交易日,若 freq = 'm'时间间隔为分钟
14

15

16
#----------------------------------- 记录部分数据 -----------------------------
17
global record
18
record = {'break_up':{},'break_down':{},'stop_loss':{},'position':{},'ATR':{}} # 记录入场、离常、止损点、持仓比、 ATR
19
#---------------------------------------------------------------------------------------
20

21
高清源代码请移步: https://uqer.io/community/share/57bd5864228e5b79a575a9b2
2855 次点击
所在节点    Python
0 条回复

这是一个专为移动设备优化的页面(即为了让你能够在 Google 搜索结果里秒开这个页面),如果你希望参与 V2EX 社区的讨论,你可以继续到 V2EX 上打开本讨论主题的完整版本。

https://www.v2ex.com/t/304039

V2EX 是创意工作者们的社区,是一个分享自己正在做的有趣事物、交流想法,可以遇见新朋友甚至新机会的地方。

V2EX is a community of developers, designers and creative people.

© 2021 V2EX