同花顺通达信公式指标改写成python实现代码，纯numpy，无循环！

myTT库非常精简，70行核心代码就实现同花顺通达信大部分技术指标公式

先看一个最简单的例子

#数字货币行情获取和指标计算演示
from  hb_hq_api import *         #数字货币行情库
from  MyTT import *              #myTT麦语言工具函数指标库

#获取btc.usdt交易对120日的数据
df=get_price('btc.usdt',count=120,frequency='1d');     # ‘4h’是4小时

#-----------df结果如下表(股市也基本一样)-------------------------------------------

	open	close	high	low	vol
2021-05-16	48983.62	47738.24	49800.00	46500.0	1.333333e+09
2021-05-17	47738.24	43342.50	48098.66	42118.0	3.353662e+09
2021-05-18	43342.50	44093.24	45781.52	42106.0	1.793267e+09

#-------有数据了，下面开始正题 -------------
CLOSE=df.close.values;  OPEN=df.open.values
HIGH=df.high.values;    LOW=df.low.values             #基础数据定义

MA5=MA(CLOSE,5)                                       #获取5日均线序列
MA10=MA(CLOSE,10)                                     #获取10日均线序列

print('BTC5日均线', MA5[-1] )                          # 只取最后一个数   
print('BTC10日均线',RET(MA10))                         # RET(MA10) == MA10[-1]
print('今天5日线是否上穿10日线',RET(CROSS(MA5,MA10)))
print('最近5天收盘价全都大于10日线吗？',EVERY(CLOSE>MA10,5) )

MyTT库中的部分工具函数

MA(CLOSE, 5)             # 取得收盘价5日平均线

加权移动平均计算：EMA

EMA(CLOSE, 5)            # 为了精度 ，  EMA至少需要120周期   

中国式的SMA计算：SMA

SMA(CLOSE, 5)            # 为了精度 ，  SMA至少需要120周期   

返回序列标准差：STD

STD(CLOSE, 5)             # 返回收盘价5日内标准差

平均绝对偏差：AVEDEV

AVEDEV(CLOSE, 5)    # 序列与其平均值的绝对差的平均值

金叉判断：CROSS

CROSS(MA(CLOSE, 5), MA(CLOSE, 10))    # 5日均线上穿10日均线

n天内满足条件的天数：COUNT

COUNT(CLOSE > OPEN, 10)           # 最近10天收阳线的天数

n天内全部满足条件的天数：EVERY

EVERY(CLOSE >OPEN, 5)           # 最近5天都是收阳线

n天内是否至少满足条件一次：EXIST

n天内最大值：HHV

HHV(MAX(OPEN, CLOSE), 20)        # 最近20天K线实体的最高价

n天内最小值：LLV

LLV(MIN(OPEN, CLOSE), 60)          # 最近60天K线实体的最低价

求和n日数据 SUM

SUM(CLOSE, 10)                            # 求和10天的收盘价

条件 IF

IF(OPEN > CLOSE, OPEN, CLOSE)        #如果 开盘>收盘  返回OPEN ，否则返回CLOSE

def MACD(CLOSE,SHORT=12,LONG=26,M=9):    # EMA的关系，CLOSE取120日，结果能精确到雪球小数点2位
    DIF = EMA(CLOSE,SHORT)-EMA(CLOSE,LONG);  
    DEA = EMA(DIF,M);      MACD=(DIF-DEA)*2
    return RD(DIF),RD(DEA),RD(MACD)

def KDJ(CLOSE,HIGH,LOW, N=9,M1=3,M2=3):   
    RSV = (CLOSE - LLV(LOW, N)) / (HHV(HIGH, N) - LLV(LOW, N)) * 100
    K = EMA(RSV, (M1*2-1));    D = EMA(K,(M2*2-1));        J=K*3-D*2
    return K, D, J

def RSI(CLOSE, N=24):      
    DIF = CLOSE-REF(CLOSE,1) 
    return RD(SMA(MAX(DIF,0), N) / SMA(ABS(DIF), N) * 100)  

def WR(CLOSE, HIGH, LOW, N=10, N1=6):            #W&R 威廉指标
    WR = (HHV(HIGH, N) - CLOSE) / (HHV(HIGH, N) - LLV(LOW, N)) * 100
    WR1 = (HHV(HIGH, N1) - CLOSE) / (HHV(HIGH, N1) - LLV(LOW, N1)) * 100
    return RD(WR), RD(WR1)

def BIAS(CLOSE,L1=6, L2=12, L3=24):  # BIAS乖离率
    BIAS1 = (CLOSE - MA(CLOSE, L1)) / MA(CLOSE, L1) * 100
    BIAS2 = (CLOSE - MA(CLOSE, L2)) / MA(CLOSE, L2) * 100
    BIAS3 = (CLOSE - MA(CLOSE, L3)) / MA(CLOSE, L3) * 100
    return RD(BIAS1), RD(BIAS2), RD(BIAS3)

def BOLL(CLOSE,N=20, P=2):            #BOLL布林带    
    MID = MA(CLOSE, N); 
    UPPER = MID + STD(CLOSE, N) * P
    LOWER = MID - STD(CLOSE, N) * P
    return RD(UPPER), RD(MID), RD(LOWER)    

def PSY(CLOSE,N=12, M=6):  
    PSY=COUNT(CLOSE>REF(CLOSE,1),N)/N*100
    PSYMA=MA(PSY,M)
    return RD(PSY),RD(PSYMA)

def CCI(CLOSE,HIGH,LOW,N=14):  
    TP=(HIGH+LOW+CLOSE)/3
    return (TP-MA(TP,N))/(0.015*AVEDEV(TP,N))

def ATR(CLOSE,HIGH,LOW, N=20):   #真实波动N日平均值
    TR = MAX(MAX((HIGH - LOW), ABS(REF(CLOSE, 1) - HIGH)), ABS(REF(CLOSE, 1) - LOW))
    return MA(TR, N)

因为语法的问题，条件与是& ，条件或是|

# 收盘价在10日均线上 且10日均线在20日均线上
(C > MA(C, 10)) & (MA(C, 10) > MA(C, 20))

# 收阳线 或 收盘价大于昨收
(C > O) | (C > REF(C, 1))

BOLL带指标数据获取和做图演示 (上证综指)

up,mid,lower=BOLL(CLOSE)                                        #获取布林带数据 

plt.figure(figsize=(15,8))  
plt.plot(CLOSE,label='上证');    plt.plot(up,label='up');        #画图显示 
plt.plot(mid,label='mid');      plt.plot(lower,label='lower');

项目地址：https://github.com/mpquant/MyTT