强化学习是机器学习的方式之一,它与监督学习、无监督学习并列,是三种机器学习训练方法之一。
在围棋上击败世界第一李世石的 alphago、在《星际争霸2》中以 10:1 击败了人类顶级职业玩家的alphastar,他们都是强化学习模型。诸如此类的模型还有 alphago zero 等。
强化学习的原理非常简单,它非常像心理学中新行为主义派的斯金纳发现的操作性条件反射。
操作性条件反射是什么?当年斯金纳做了一个箱子,进行了两次实验。
第一次实验,箱子里放了一只饥饿的老鼠,在箱子的一边有一个可供按压的杠杆,在杠杆旁边有一个放置食物的小盒子。动物在箱内按下杠杆,食盒就会释放食物进入箱内,动物可以取食。结果:小鼠自发学会了按按钮。这是积极强化。
另一次实验是,每次小白鼠不按下按钮,则给箱子通电,小白鼠因此学会了按按钮以防自己遭受电击。这是消极强化(负向强化)。
这就是斯金纳发现的操作性条件反射,当行为得到奖励或惩罚时出现刺激,反过来控制这种行为。
强化学习与操作性条件反射有异曲同工之妙,以人类玩游戏为例,如果在游戏中采取某种策略购买某种类型的装备可以取得较高的得分,那么就会进一步“强化”这种策略,以期继续取得较好的结果。
网上有不少强化学习的例子,鉴于读者中对股票感兴趣的同学比较多,我们以股票预测为例,实验一下 wangshub 的 rl-stock 项目。
下面就试一下这个强化学习项目,前往github下载 rl-stock:
https://github.com/wangshub/rl-stock/
如果你无法使用github,也可以在python实用宝典公众号后台回复:股票强化学习1 下载全文完整代码,这份代码包括第三部分的多进程优化逻辑。
1.准备
开始之前,你要确保python和pip已经成功安装在电脑上,如果没有,可以访问这篇文章:超详细python安装指南 进行安装。
**(可选1) **如果你用python的目的是数据分析,可以直接安装anaconda:python数据分析与挖掘好帮手—anaconda,它内置了python和pip.
**(可选2) **此外,推荐大家用vscode编辑器,它有许多的优点:python 编程的最好搭档—vscode 详细指南。
请注意,由于tensorflow版本限制,这个强化学习项目只支持 python3 以上,python3.6 及以下的版本,因此我建议使用anaconda创建一个新的虚拟环境运行这个项目:
conda create -n rlstock python=3.6另外,实测依赖需要改动 requirements.txt 的tensorflow-gpu版本至1.14:
windows环境下打开cmd(开始—运行—cmd),苹果系统环境下请打开terminal(command+空格输入terminal),进入 rl-stock 项目文件夹输入命令安装依赖:
pip install -r requirements.txt2.小试强化学习
运行rl-stock项目前,需要下载数据。进入刚创建的虚拟环境,运行get_stock_data.py代码会自动下载数据到 stockdata 目录中:
python get_stock_data.py如果你使用的是在github上下载的代码而不是python实用宝典后台改好的代码,请注意 get_stock_data.py 的第46行,必须对 row[code_name] 去除 * 号,否则windows系统下可能会存在问题:
df_code.to_csv(f'{self.output_dir}/{row[code]}.{row[code_name].strip(*)}.csv', index=false)数据下载完成后就可以运行 main.py 执行强化学习训练和测试,不过在训练之前,我们先简单了解下整个项目的输入状态、动作、和奖励函数。
输入状态(观测 observation)
策略网络观测的就是一只股票的各类数据,比如开盘价、收盘价、成交量等,它可以由许多因子组成。为了训练时网络收敛,观测状态数据输入时必须要进行归一化,变换到 ** [-1, 1] ** 的区间内。rl-stock输入的观测数据字段如下:
动作 action
共有 买入 、卖出和**持有 **3 种动作,定义动作( ** action ** )为长度为 2 的数组
**action[0] **为操作类型;action[1] 表示买入或卖出百分比;动作类型action[0]说明
1 买入action[1]
2 卖出action[1]
3 持有
注意,当动作类型 action[0] = 3 时,表示不买也不抛售股票,此时 action[1] 的值无实际意义,网络在训练过程中,agent 会慢慢学习到这一信息。agent,实称代理,在我们的上下文中,你可以视其为策略。
奖励 reward
奖励函数的设计,对强化学习的目标至关重要。在股票交易的环境下,最应该关心的就是当前的盈利情况,故用当前的利润作为奖励函数。
# profitsreward = self.net_worth - initial_account_balancereward = 1 if reward > 0 else -100为了使网络更快学习到盈利的策略,当利润为负值时,给予网络一个较大的惩罚 (-100)。
梯度策略
作者采用了基于策略梯度的ppo 算法,openai 和许多文献已把 ppo 作为强化学习研究中首选的算法。ppo 优化算法 python 实现参考 stable-baselines。
数据集及自定义
在数据集上,作者使用了1990年至2019年11月作为训练集,2019年12月作为测试集。
1990-01-01~2019-11-292019-12-01~2019-12-31
训练集 测试集
如果你要调整这个训练集和测试集的时间,可以更改 get_stock_data.py 的以下部分:
if __name__ == '__main__': # 获取全部股票的日k线数据 # 训练集 mkdir('./stockdata/train') downloader = downloader('./stockdata/train', date_start='1990-01-01', date_end='2019-11-29') downloader.run() # 测试集 mkdir('./stockdata/test') downloader = downloader('./stockdata/test', date_start='2019-12-01', date_end='2019-12-31') downloader.run()训练并测试
首先,我们尝试一下单一代码的训练和测试,修改main.py里的股票代码,比如我这里修改为601919中远海控:
if __name__ == '__main__': # multi_stock_trade() test_a_stock_trade('sh.601919') # ret = find_file('./stockdata/train', '601919') # print(ret)运行下面的命令,执行此深度学习模型的训练和测试。
python main.py训练完成后,会自动进行模拟操作测试集这20个交易日,然后会输出这20个交易日的测试结果:
------------------------------step: 20balance: 0.713083354256014shares held: 2060 (total sold: 2392)avg cost for held shares: 5.072161917927474 (total sales value: 12195.091008936648)net worth: 10930.56492977963 (max net worth: 10930.56492977963)profit: 930.5649297796299------------------------------step: 21balance: 0.713083354256014shares held: 2060 (total sold: 2392)avg cost for held shares: 5.072161917927474 (total sales value: 12195.091008936648)net worth: 10815.713083354256 (max net worth: 10930.56492977963)profit: 815.713083354256利润图如下:
然后我们看一下中远海控2019年12月的走势:
可以看到这个月的中远海控是一个上升趋势,一共上涨了12%,而这个模型捕捉到其中8%左右的利润,还是相当不错的。
当然,凡事不能只看个体,下面我们修改下作者的源代码,试一下其在市场里的整体表现。
3.强化学习模型整体表现
由于作者原有的模型是单进程的计算,为了测试全市场的表现,我进行了多进程改造。
我将作者的训练及测试任务集成到一个函数中,并使用celery做并行:
@app.taskdef multi_stock_trade(code): stock_file = find_file('./stockdata/train', str(code)) if stock_file: try: profits = stock_trade(stock_file) with open(f'result/code-{code}.pkl', 'wb') as f: pickle.dump(profits, f) except exception as err: print(err)将训练集的时间范围设置为1990年到2021年,测试集的测试周期改为最近一个月:
1990-01-01~2021-11-252021-11-26~2021-12-25
训练集 测试集
开启redis-server 及 celery worker:
# redis-server 独占一个进程,所以需要另开一个窗口celery -a tasks worker -l info遍历所有的股票代码做并发测试:
files = os.listdir(stockdata/train)files_test = os.listdir(stockdata/test)all_files_list = list(set(files) & set(files_test))for i in all_files_list: # 使用celery做并发 code = ..join(i.split(.)[:2]) # multi_stock_trade.apply_async(args=(code,)) multi_stock_trade(code)再对生成的结果进行统计,测试结果如下:
对这个模型在2021-11-26到2021-12-25的测试结果表明,有40.8%的股票进行了交易并且获利,有49.9%的股票没有进行操作,有9.4%的股票进行了交易并亏损。平均每次交易利润为445元,作为一个测试策略,这个结果已经很不错了。
由于只是一个测试策略,这里就不做详细的风险分析了,实际上我们还需要观察这个策略的最大回撤率、夏普率等指标才能更好地评判此策略的好坏。
我认为这个项目还有很大的改造空间,原逻辑中只观察了ohlc等基本数据,我们还可以增加很多指标,比如基于ta-lib,算出macd、rsi等技术指标,再将其加入observation中,让模型观察学习这些数据的特征,可能会有不错的表现。
热阻相关的JEDEC标准介绍
思特威推出全新0.7μm像素尺寸5200万CMOS图像传感器
采用AT89C52单片机实现彩色像片洗印机温控系统的设计方案
【号外】科陆电子中标深圳充电桩建设项目大单
嵌入式环形队列和消息队列的实现
什么是强化学习
介绍基于不确定的语法条件生成类似Java的强类型程序
2019年华为5G手机总售出690万台,智能手机发货量超过2.4亿台
微型电池需求爆发 小电池企业利好
关于外转子无刷电机电磁设计案例的详细说明
AMD桌面处理器份额3年来首次环比下滑
使用TDA2050构建一个2x32瓦音频放大器电路
六氟化硫气体检测——SF6气体红外泄漏报警监测系统
浅谈建筑防火监控系统在民用住宅建筑设计中的应用
怎样设计出接近要求规格的DC/DC转换器电路
什么是LED光色质量
关于人工智能最佳交互界面的分析介绍
为什么我抓不到baidu的数据包
Amber Solutions 专利支持电网的数字化管理
Renesas推出SH-Mobile高端产品