通达信公式02-攻击&力量Python化代码
1.公式
```python
TYP:=(HIGH+LOW+CLOSE)/3;
攻击:(TYP-MA(TYP,10))/(0.015*AVEDEV(TYP,10)),LINETHICK1;
力量:(TYP-MA(TYP,30))/(0.015*AVEDEV(TYP,30)),COLORYELLOW,LINETHICK2;
通达信翻译：
TYP赋值:(最高价+最低价+收盘价)/3
输出攻击:(TYP-TYP的10日简单移动平均)/(0.015*TYP的10日平均绝对偏差),线宽为1
输出力量:(TYP-TYP的30日简单移动平均)/(0.015*TYP的30日平均绝对偏差),画黄色,线宽为2
```
2.代码化处理
```python
#1.该策略数据来源于baostock
#2.相关的计算与通达信的不同，由于对于计算机等常识的不懂，所以无法调用通达信数据进行分析
#但是实际计算的结果和通达信展示的不同，故而先进行计算，看看效果。
#3.重要的一点是，该数据貌似不是用前复权数据，因为从通达信的公式导出来的数据完全没有。
#4.先用回测看一下，所以将数据导入接入进来了。
import baostock as bs
import pandas as pd
import matplotlib.pyplot as plt
pd.set_option('expand_frame_repr',False)
#py文本名字：AttackPower.py
def dfload(code="sh.600789",startday='1990-01-01',endday='2020-01-01'):
#导入数据
lg=bs.login()
rs = bs.query_history_k_data_plus(code,
"date,code,open,high,low,close,preclose,volume,amount,adjustflag,turn,tradestatus,pctChg,isST",
start_date=startday, end_date=endday,
frequency="d", adjustflag="2")
data_list = []
while (rs.error_code == '0') & rs.next():
# 获取一条记录，将记录合并在一起
data_list.append(rs.get_row_data())
result = pd.DataFrame(data_list, columns=rs.fields)
return result
#5.由于baostock的数据导出均为object，所以需要转化
def astypedata(df):
df['high'] = df['high'].astype(float)
df['low'] = df['low'].astype(float)
df['open'] = df['open'].astype(float)
df['close'] = df['close'].astype(float)
df['date'] = pd.to_datetime(df['date'])
return df
#6.计算绝对误差值
def avedev(df):
'''
TYP:=(HIGH+LOW+CLOSE)/3;
攻击:(TYP-MA(TYP,10))/(0.015*AVEDEV(TYP,10)),LINETHICK1;
力量:(TYP-MA(TYP,30))/(0.015*AVEDEV(TYP,30)),COLORYELLOW,LINETHICK2;
金钻突破:力量-REF(力量,1)>90 OR 攻击-REF(攻击,1)>90;
由于金钻突破只有买点，没有买点，所以选择忽视，只以攻击和力量的交叉点进行判断买点卖点。
:param df:
:return:
'''
df['typ'] = (df['high'] + df['low'] + df['close']) / 3
df['avedev'] = abs(df['typ'] - sum(df['typ']) / len(df))
df['avedev10'] = df['avedev'].rolling(10).mean()
df['avedev30'] = df['avedev'].rolling(30).mean()
df['attack'] = (df['typ'] - df['typ'].rolling(10).mean()) / (0.015 * df['avedev10'])
df['power'] = (df['typ'] - df['typ'].rolling(30).mean()) / (0.015 * df['avedev30'])
df=df.fillna(0)
df=round(df,2)
return df
def sign(df):
for i in range(df.shape[0]):
if df['attack'][i]>df['power'][i] and df['attack'][i-1]<=df['power'][i-1]:
df.ix[i,'signal']=1
if df['attack'][i]<df['power'][i] and df['attack'][i-1]>df['power'][i-1]:
df.ix[i,'signal']=-1
return df
def position(df):
df['position'] = df['signal'].shift()
df['position'].fillna(method='ffill', inplace=True)
# 不能买的信号
cond_cannot_buy = df['open'] > df['close'].shift(1) * 1.097
# 将开盘涨停日、并且当天position为1时的'pos'设置为空值
df.loc[cond_cannot_buy & (df['position'] == 1), 'position'] = None
# 不能卖的信号
# 找出开盘跌停的日期
cond_cannot_sell = df['open'] < df['close'].shift(1) * 0.903 # 今天的开盘价相对于昨天的收盘价下得了9.7%
# 将开盘跌停日、并且当天position为0时的'pos'设置为空值
df.loc[cond_cannot_sell & (df['position'] == -1), 'position'] = None
# position为空的日期，不能买卖。position只能和前一个交易日保持一致。
df['position'].fillna(method='ffill', inplace=True)
# 在position为空值的日期，将position补全为0
df['position'].fillna(value=-1, inplace=True)
return df
#后面就是回测阶段，写到这里就可以了。
```