写一段python代码,读取一个csv文件,该文件的第一列为时间编号,编号以1,2,3……这样自然数的形式出现,第二列为充电量数值,将该数据集的前70%划分为训练集,其余部分设置为测试集,将该数据集所构成的时间序列进行EMD分解,分解出的子序列中低频的部分用svr模型进行训练,高频的部分用lstm模型进行训练,然后用stacking方法的集成学习对总体进行训练,并对总序列进行预测,并打印rmse作为评价指标,并分别绘制训练集的真实值和预测值图片以及测试集的真实值和预测值图片,两张图片的横坐标单位为time/h,纵坐标单位为kwh

时间: 2023-10-30 16:08:08 浏览: 285
# 导入需要的库 import pandas as pd import numpy as np from PyEMD import EMD from sklearn.svm import SVR from sklearn.preprocessing import StandardScaler from keras.models import Sequential from keras.layers import Dense, LSTM from sklearn.metrics import mean_squared_error import matplotlib.pyplot as plt # 读取csv文件 df = pd.read_csv('data.csv') # 划分训练集和测试集 train_size = int(len(df) * 0.7) train, test = df.values[:train_size, :], df.values[train_size:, :] # 对时间序列进行EMD分解 emd = EMD() imfs = emd(df.iloc[:, 1].values) # 分离出低频部分和高频部分 low_freq = imfs[0] high_freq = imfs[1:] # 划分训练集和测试集 train_low, test_low = low_freq[:train_size], low_freq[train_size:] train_high, test_high = high_freq[:, :train_size], high_freq[:, train_size:] # 对低频部分进行SVR模型训练 scaler = StandardScaler() train_low_scaled = scaler.fit_transform(train_low.reshape(-1, 1)) test_low_scaled = scaler.transform(test_low.reshape(-1, 1)) svr = SVR(kernel='rbf', C=100, gamma=0.1, epsilon=.1) svr.fit(train_low_scaled, train[:, 1]) # 对高频部分进行LSTM模型训练 train_high = train_high.T test_high = test_high.T def create_lstm_model(train_X, train_y, n_steps, n_features): model = Sequential() model.add(LSTM(50, activation='relu', input_shape=(n_steps, n_features))) model.add(Dense(1)) model.compile(optimizer='adam', loss='mse') model.fit(train_X, train_y, epochs=200, verbose=0) return model n_steps = 3 n_features = len(high_freq) train_X, train_y = [], [] test_X, test_y = [], [] for i in range(n_steps, len(train_high)): train_X.append(train_high[i-n_steps:i, :]) train_y.append(train[i, 1]) for i in range(n_steps, len(test_high)): test_X.append(test_high[i-n_steps:i, :]) test_y.append(test[i, 1]) train_X, train_y = np.array(train_X), np.array(train_y) test_X, test_y = np.array(test_X), np.array(test_y) lstm = create_lstm_model(train_X, train_y, n_steps, n_features) # 进行集成学习 def stacked_dataset(train_X, train_y, test_X): stackX = None for model in [svr, lstm]: model.fit(train_X, train_y) yhat = model.predict(test_X) if stackX is None: stackX = yhat else: stackX = np.column_stack((stackX, yhat)) return stackX stacked_train_X = stacked_dataset(train_high.T, train[:, 1], train_X.reshape((train_X.shape[0], n_steps*n_features))) stacked_test_X = stacked_dataset(test_high.T, test[:, 1], test_X.reshape((test_X.shape[0], n_steps*n_features))) model = SVR(kernel='linear') model.fit(stacked_train_X, train[:, 1]) yhat = model.predict(stacked_test_X) print('RMSE:', np.sqrt(mean_squared_error(test[:, 1], yhat))) # 绘制训练集和测试集的真实值和预测值图片 plt.plot(train[:, 0], train[:, 1], label='True Train') plt.plot(train[:, 0], svr.predict(train_low_scaled), label='Predicted Train - SVR') plt.plot(train[:, 0][n_steps:], lstm.predict(train_X).flatten(), label='Predicted Train - LSTM') plt.plot(test[:, 0], test[:, 1], label='True Test') plt.plot(test[:, 0], yhat, label='Predicted Test - Stacking') plt.xlabel('Time/h') plt.ylabel('kwh') plt.legend() plt.show()
阅读全文

相关推荐

最新推荐

recommend-type

aiohttp-3.7.3-cp36-cp36m-win_amd64.whl.rar

python whl离线安装包 pip安装失败可以尝试使用whl离线安装包安装 第一步 下载whl文件,注意需要与python版本配套 python版本号、32位64位、arm或amd64均有区别 第二步 使用pip install XXXXX.whl 命令安装,如果whl路径不在cmd窗口当前目录下,需要带上路径 WHL文件是以Wheel格式保存的Python安装包, Wheel是Python发行版的标准内置包格式。 在本质上是一个压缩包,WHL文件中包含了Python安装的py文件和元数据,以及经过编译的pyd文件, 这样就使得它可以在不具备编译环境的条件下,安装适合自己python版本的库文件。 如果要查看WHL文件的内容,可以把.whl后缀名改成.zip,使用解压软件(如WinRAR、WinZIP)解压打开即可查看。 为什么会用到whl文件来安装python库文件呢? 在python的使用过程中,我们免不了要经常通过pip来安装自己所需要的包, 大部分的包基本都能正常安装,但是总会遇到有那么一些包因为各种各样的问题导致安装不了的。 这时我们就可以通过尝试去Python安装包大全中(whl包下载)下载whl包来安装解决问题。
recommend-type

基于Java中的swing类的图形化飞机游戏的开发练习.zip

基于Java中的Swing类开发的图形化飞机游戏练习包,为初学者和进阶学习者提供了实践Java GUI编程的绝佳机会。通过本资源,开发者可以利用Java语言和Swing库构建一个用户交互式的2D游戏,深入理解图形用户界面(GUI)编程和事件处理机制。该游戏的核心包括玩家飞机的控制、敌机的生成与移动、子弹发射与碰撞检测以及游戏胜负判定等逻辑。玩家通过鼠标移动控制己方飞机,实现平滑的移动和连续的子弹发射;而敌方飞机则按照一定算法无规律出现,随着游戏进程难度逐渐增加。游戏中还引入了特殊NPC,增加了额外的挑战和乐趣。为了提高游戏体验,游戏还包含了开始背景、结束背景以及背景音乐等元素。当玩家击毁敌机时,会有相应的得分计算和展示;若被敌机击中,则游戏结束并显示最终得分。此外,游戏还提供了查看历史前十记录、帮助和退出等选项,方便玩家进行游戏设置和了解游戏玩法。本资源适用于计算机科学与技术、软件工程、信息管理及相关专业的课程设计、毕业设计等环节,为学生提供实践操作的机会,帮助他们巩固Java编程知识,提高动手能力和发散思维。同时,也为希望学习不同技术领域的学习者提供了一个优秀的入门项目。
recommend-type

SQLite:SQLite数据库创建与管理.docx

SQLite:SQLite数据库创建与管理
recommend-type

【完整源码+数据库】SpringBoot 集成 Spring Security短信验证码登录

Spring Security 默认是账号和密码登录,现在是对 Spring Security 进行扩展,来实现短信验证码方式登录。 SpringBoot 集成 Spring Security短信验证码登录【完整源码+数据库】
recommend-type

去年和朋友一起做的java小游戏.游戏具体界面在readme中,游戏设计的uml图在design.pdf中.zip

本资源是一个Java小游戏项目,由我和我的朋友在去年共同开发。这个项目不仅包含了完整的游戏代码,还有详细的设计文档和UML图,适合作为学习和参考的素材。游戏的界面设计简洁明了,玩法有趣且富有挑战性,能够让玩家在游戏中体验到乐趣。在readme文件中,你可以找到游戏的具体界面展示,让你对游戏的外观有一个直观的了解。而design.pdf中则包含了游戏的UML图,详细展示了游戏的设计结构和各个模块之间的关系,对于理解游戏的整体架构非常有帮助。这个Java小游戏项目是一个非常好的学习资源,无论是对于初学者还是有一定经验的开发者来说,都可以通过这个项目来提升自己的编程技能和游戏设计能力。通过阅读代码和设计文档,你可以了解到如何构建一个功能完整的游戏,并且可以根据自己的需要进行修改和扩展。总之,这个Java小游戏项目是一个值得学习和探索的资源,希望对你有所帮助!
recommend-type

探索AVL树算法:以Faculdade Senac Porto Alegre实践为例

资源摘要信息:"ALG3-TrabalhoArvore:研究 Faculdade Senac Porto Alegre 的算法 3" 在计算机科学中,树形数据结构是经常被使用的一种复杂结构,其中AVL树是一种特殊的自平衡二叉搜索树,它是由苏联数学家和工程师Georgy Adelson-Velsky和Evgenii Landis于1962年首次提出。AVL树的名称就是以这两位科学家的姓氏首字母命名的。这种树结构在插入和删除操作时会维持其平衡,以确保树的高度最小化,从而在最坏的情况下保持对数的时间复杂度进行查找、插入和删除操作。 AVL树的特点: - AVL树是一棵二叉搜索树(BST)。 - 在AVL树中,任何节点的两个子树的高度差不能超过1,这被称为平衡因子(Balance Factor)。 - 平衡因子可以是-1、0或1,分别对应于左子树比右子树高、两者相等或右子树比左子树高。 - 如果任何节点的平衡因子不是-1、0或1,那么该树通过旋转操作进行调整以恢复平衡。 在实现AVL树时,开发者通常需要执行以下操作: - 插入节点:在树中添加一个新节点。 - 删除节点:从树中移除一个节点。 - 旋转操作:用于在插入或删除节点后调整树的平衡,包括单旋转(左旋和右旋)和双旋转(左右旋和右左旋)。 - 查找操作:在树中查找一个节点。 对于算法和数据结构的研究,理解AVL树是基础中的基础。它不仅适用于算法理论的学习,还广泛应用于数据库系统、文件系统以及任何需要快速查找和更新元素的系统中。掌握AVL树的实现对于提升软件效率、优化资源使用和降低算法的时间复杂度至关重要。 在本资源中,我们还需要关注"Java"这一标签。Java是一种广泛使用的面向对象的编程语言,它对数据结构的实现提供了良好的支持。利用Java语言实现AVL树,可以采用面向对象的方式来设计节点类和树类,实现节点插入、删除、旋转及树平衡等操作。Java代码具有很好的可读性和可维护性,因此是实现复杂数据结构的合适工具。 在实际应用中,Java程序员通常会使用Java集合框架中的TreeMap和TreeSet类,这两个类内部实现了红黑树(一种自平衡二叉搜索树),而不是AVL树。尽管如此,了解AVL树的原理对于理解这些高级数据结构的实现原理和使用场景是非常有帮助的。 最后,提及的"ALG3-TrabalhoArvore-master"是一个压缩包子文件的名称列表,暗示了该资源是一个关于AVL树的完整项目或教程。在这个项目中,用户可能可以找到完整的源代码、文档说明以及可能的测试用例。这些资源对于学习AVL树的实现细节和实践应用是宝贵的,可以帮助开发者深入理解并掌握AVL树的算法及其在实际编程中的运用。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

【ggplot2绘图技巧】:R语言中的数据可视化艺术

![【ggplot2绘图技巧】:R语言中的数据可视化艺术](https://www.lecepe.fr/upload/fiches-formations/visuel-formation-246.jpg) # 1. ggplot2绘图基础 在本章节中,我们将开始探索ggplot2,这是一个在R语言中广泛使用的绘图系统,它基于“图形语法”这一理念。ggplot2的设计旨在让绘图过程既灵活又富有表现力,使得用户能够快速创建复杂而美观的图形。 ## 1.1 ggplot2的安装和加载 首先,确保ggplot2包已经被安装。如果尚未安装,可以使用以下命令进行安装: ```R install.p
recommend-type

HAL库怎样将ADC两个通道的电压结果输出到OLED上?

HAL库通常是指硬件抽象层(Hardware Abstraction Layer),它是一个软件组件,用于管理和控制嵌入式系统中的硬件资源,如ADC(模拟数字转换器)和OLED(有机发光二极管显示屏)。要将ADC读取的两个通道电压值显示到OLED上,你可以按照以下步骤操作: 1. **初始化硬件**: 首先,你需要通过HAL库的功能对ADC和OLED进行初始化。这包括配置ADC的通道、采样速率以及OLED的分辨率、颜色模式等。 2. **采集数据**: 使用HAL提供的ADC读取函数,读取指定通道的数据。例如,在STM32系列微控制器中,可能会有`HAL_ADC_ReadChannel()
recommend-type

小学语文教学新工具:创新黑板设计解析

资源摘要信息: 本资源为行业文档,主题是设计装置,具体关注于一种小学语文教学黑板的设计。该文档通过详细的设计说明,旨在为小学语文教学场景提供一种创新的教学辅助工具。由于资源的标题、描述和标签中未提供具体的设计细节,我们仅能从文件名称推测文档可能包含了关于小学语文教学黑板的设计理念、设计要求、设计流程、材料选择、尺寸规格、功能性特点、以及可能的互动功能等方面的信息。此外,虽然没有标签信息,但可以推断该文档可能针对教育技术、教学工具设计、小学教育环境优化等专业领域。 1. 教学黑板设计的重要性 在小学语文教学中,黑板作为传统而重要的教学工具,承载着教师传授知识和学生学习互动的重要角色。一个优秀的设计可以提高教学效率,激发学生的学习兴趣。设计装置时,考虑黑板的适用性、耐用性和互动性是非常必要的。 2. 教学黑板的设计要求 设计小学语文教学黑板时,需要考虑以下几点: - 安全性:黑板材质应无毒、耐磨损,边角处理要圆滑,避免在使用中造成伤害。 - 可视性:黑板的大小和高度应适合小学生使用,保证最远端的学生也能清晰看到上面的内容。 - 多功能性:黑板除了可用于书写字词句之外,还可以考虑增加多媒体展示功能,如集成投影幕布或电子白板等。 - 环保性:使用可持续材料,比如可回收的木材或环保漆料,减少对环境的影响。 3. 教学黑板的设计流程 一个典型的黑板设计流程可能包括以下步骤: - 需求分析:明确小学语文教学的需求,包括空间大小、教学方法、学生人数等。 - 概念设计:提出初步的设计方案,并对方案的可行性进行分析。 - 制图和建模:绘制详细的黑板平面图和三维模型,为生产制造提供精确的图纸。 - 材料选择:根据设计要求和成本预算选择合适的材料。 - 制造加工:按照设计图纸和材料标准进行生产。 - 测试与评估:在实际教学环境中测试黑板的使用效果,并根据反馈进行必要的调整。 4. 教学黑板的材料选择 - 传统黑板:传统的黑板多由优质木材和专用黑板漆制成,耐用且书写流畅。 - 绿色环保材料:考虑到环保和学生健康,可以选择无毒或低VOC(挥发性有机化合物)排放的材料。 - 智能材料:如可擦洗的特殊漆料,使黑板表面更加光滑,便于擦拭。 5. 教学黑板的尺寸规格 黑板的尺寸规格应根据实际教室空间和学生的平均身高来设计。一般来说,小学教室的黑板高度应设置在120cm至150cm之间,长度则根据教室墙壁的长度而定,但至少应保证可以容纳整页A4纸的书写空间。 6. 教学黑板的功能性特点 - 书写性能:黑板表面应具备良好的书写性能,使粉笔或马克笔的书写和擦拭都十分顺畅。 - 可视化辅助:集成的可视化工具,如辅助灯、放大镜等,可以帮助教师更有效地展示教学内容。 - 互动性设计:考虑增加互动性元素,例如磁性或可擦写的表面,可以提高学生参与度。 7. 教学黑板的互动功能 随着信息技术的发展,教学黑板可以集成多媒体技术,如触摸屏功能、电子白板功能、互联网接入等,实现与电子设备的互动,从而丰富教学手段,提高教学的趣味性和效率。 综上所述,本资源提供的设计装置文档,聚焦于一种小学语文教学黑板的设计,涵盖了从设计理念到功能实现的全方位内容,旨在通过创新的设计提升小学语文教学的品质和效率。