import numpy as np import pandas as pd from sklearn.model_selection import train_test_split from keras.models import Sequential from keras.layers import Dense from pyswarm import pso import matplotlib.pyplot as plt from sklearn.preprocessing import StandardScaler from sklearn.metrics import mean_absolute_error from sklearn.metrics import mean_squared_error from sklearn.metrics import r2_score file = "zhong.xlsx" data = pd.read_excel(file) #reading file X=np.array(data.loc[:,'种植密度':'有效积温']) y=np.array(data.loc[:,'产量']) y.shape=(185,1) # 将数据集分为训练集和测试集 X_train, X_test, y_train, y_test = train_test_split(X,y, test_size=0.25, random_state=10) SC=StandardScaler() X_train=SC.fit_transform(X_train) X_test=SC.fit_transform(X_test) y_train=SC.fit_transform(y_train) y_test=SC.fit_transform(y_test) print("X_train.shape:", X_train.shape) print("X_test.shape:", X_test.shape) print("y_train.shape:", y_train.shape) print("y_test.shape:", y_test.shape) # 定义BP神经网络模型 def nn_model(X): model = Sequential() model.add(Dense(8, input_dim=X_train.shape[1], activation='relu')) model.add(Dense(12, activation='relu')) model.add(Dense(1)) model.compile(loss='mean_squared_error', optimizer='adam') return model # 定义适应度函数 def fitness_func(X): model = nn_model(X) model.fit(X_train, y_train, epochs=60, verbose=2) score = model.evaluate(X_test, y_test, verbose=2) print(score) # 定义变量的下限和上限 lb = [5, 5] ub = [30, 30] # 利用PySwarm库实现改进的粒子群算法来优化BP神经网络预测模型 result = pso(fitness_func, lb, ub) # 输出最优解和函数值 print('最优解:', result[0]) print('最小函数值:', result[1]) mpl.rcParams["font.family"] = "SimHei" mpl.rcParams["axes.unicode_minus"] = False # 绘制预测值和真实值对比图 model = nn_model(X) model.fit(X_train, y_train, epochs=60, verbose=2) y_pred = model.predict(X_test) y_true = SC.inverse_transform(y_test) y_pred=SC.inverse_transform(y_pred) plt.figure() plt.plot(y_true,"bo-",label = '真实值') plt.plot(y_pred,"ro-", label = '预测值') plt.title('神经网络预测展示') plt.xlabel('序号') plt.ylabel('产量') plt.legend(loc='upper right') plt.show() print("R2 = ",r2_score(y_test, y_pred)) # R2 # 绘制损失函数曲线图 model = nn_model(X) history = model.fit(X_train, y_train, epochs=60, validation_data=(X_test, y_test), verbose=2) plt.plot(history.history['loss'], label='train') plt.plot(history.history['val_loss'], label='test') plt.legend() plt.show() mae = mean_absolute_error(y_test, y_pred) print('MAE: %.3f' % mae) mse = mean_squared_error(y_test, y_pred) print('mse: %.3f' % mse)

时间: 2024-03-21 12:40:29 浏览: 179
RAR

first_NN_keras_keras_神经网络_多层感知机_糖尿病数据集_糖尿病_

你的代码中存在一些问题。首先,在导入模块时,你应该将每个模块分开导入,而不是在同一行导入。例如,应该将`import numpy as np import pandas as pd`分开成两行导入。其次,在进行数据标准化时,你应该使用同一个`StandardScaler`对象对训练集和测试集进行转换,而不是分别创建两个不同的对象进行转换。应该如下修改: ``` SC=StandardScaler() X_train=SC.fit_transform(X_train) X_test=SC.transform(X_test) y_train=SC.fit_transform(y_train) y_test=SC.transform(y_test) ``` 另外,你在训练模型时,对于`y_train`和`y_test`也进行了标准化,这是不正确的。你应该只对输入变量进行标准化,而不是对输出变量进行标准化。应该修改为: ``` SC=StandardScaler() X_train=SC.fit_transform(X_train) X_test=SC.transform(X_test) ``` 最后,你在定义适应度函数时,应该将`X`作为参数传递给`nn_model`函数,而不是直接使用全局变量`X_train`和`y_train`。应该修改为: ``` def fitness_func(X): model = nn_model(X) model.fit(X_train, y_train, epochs=60, verbose=2) score = model.evaluate(X_test, y_test, verbose=2) print(score) ``` 除了以上问题,你的代码看起来应该能够正常运行。
阅读全文

相关推荐

import numpy as np import tensorflow as tf from keras.models import Sequential from keras.layers import Dense, Activation, Dropout, Flatten from keras.layers.convolutional import Conv2D, MaxPooling2D from keras.utils import np_utils from keras.datasets import mnist from keras import backend as K from keras.optimizers import Adam import skfuzzy as fuzz import pandas as pd from sklearn.model_selection import train_test_split # 绘制损失曲线 import matplotlib.pyplot as plt from sklearn.metrics import accuracy_score data = pd.read_excel(r"D:\pythonProject60\filtered_data1.xlsx") # 读取数据文件 # Split data into input and output variables X = data.iloc[:, :-1].values y = data.iloc[:, -1].values X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=42) # 导入MNIST数据集 # 数据预处理 y_train = np_utils.to_categorical(y_train, 3) y_test = np_utils.to_categorical(y_test, 3) # 创建DNFN模型 model = Sequential() model.add(Dense(64, input_shape=(11,), activation='relu')) model.add(Dense(128, activation='relu')) model.add(Dropout(0.5)) model.add(Dense(3, activation='softmax')) # 编译模型 model.compile(loss='categorical_crossentropy', optimizer=Adam(), metrics=['accuracy']) # 训练模型 history = model.fit(X_train, y_train, validation_data=(X_test, y_test), epochs=10, batch_size=128) # 使用DNFN模型进行预测 y_pred = model.predict(X_test) y_pred= np.argmax(y_pred, axis=1) print(y_pred) # 计算模糊分类 fuzzy_pred = [] for i in range(len(y_pred)): fuzzy_class = np.zeros((3,)) fuzzy_class[y_pred[i]] = 1.0 fuzzy_pred.append(fuzzy_class) fuzzy_pred = np.array(fuzzy_pred) print(fuzzy_pred)获得其运行时间

mport numpy as np import tensorflow as tf from keras.models import Sequential from keras.layers import Dense, Activation, Dropout, Flatten from keras.layers.convolutional import Conv2D, MaxPooling2D from keras.utils import np_utils from keras.datasets import mnist from keras import backend as K from keras.optimizers import Adam import skfuzzy as fuzz import pandas as pd from sklearn.model_selection import train_test_split # 绘制损失曲线 import matplotlib.pyplot as plt import time from sklearn.metrics import accuracy_score data = pd.read_excel(r"D:\pythonProject60\filtered_data1.xlsx") # 读取数据文件 # Split data into input and output variables X = data.iloc[:, :-1].values y = data.iloc[:, -1].values X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=42) # 导入MNIST数据集 # 数据预处理 y_train = np_utils.to_categorical(y_train, 3) y_test = np_utils.to_categorical(y_test, 3) # 创建DNFN模型 start_time=time.time() model = Sequential() model.add(Dense(64, input_shape=(11,), activation='relu')) model.add(Dense(128, activation='relu')) model.add(Dropout(0.5)) model.add(Dense(3, activation='softmax')) # 编译模型 model.compile(loss='categorical_crossentropy', optimizer=Adam(), metrics=['accuracy']) # 训练模型 history = model.fit(X_train, y_train, validation_data=(X_test, y_test), epochs=10, batch_size=128) # 使用DNFN模型进行预测 y_pred = model.predict(X_test) y_pred= np.argmax(y_pred, axis=1) print(y_pred) # 计算模糊分类 fuzzy_pred = [] for i in range(len(y_pred)): fuzzy_class = np.zeros((3,)) fuzzy_class[y_pred[i]] = 1.0 fuzzy_pred.append(fuzzy_class) fuzzy_pred = np.array(fuzzy_pred) end_time = time.time() print("Total time taken: ", end_time - start_time, "seconds")获得运行结果并分析

import pandas as pd data = pd.read_csv(C:\Users\Administrator\Desktop\pythonsjwj\weibo_senti_100k.csv') data = data.dropna(); data.shape data.head() import jieba data['data_cut'] = data['review'].apply(lambda x: list(jieba.cut(x))) data.head() with open('stopword.txt','r',encoding = 'utf-8') as f: stop = f.readlines() import re stop = [re.sub(' |\n|\ufeff','',r) for r in stop] data['data_after'] = [[i for i in s if i not in stop] for s in data['data_cut']] data.head() w = [] for i in data['data_after']: w.extend(i) num_data = pd.DataFrame(pd.Series(w).value_counts()) num_data['id'] = list(range(1,len(num_data)+1)) a = lambda x:list(num_data['id'][x]) data['vec'] = data['data_after'].apply(a) data.head() from wordcloud import WordCloud import matplotlib.pyplot as plt num_words = [''.join(i) for i in data['data_after']] num_words = ''.join(num_words) num_words= re.sub(' ','',num_words) num = pd.Series(jieba.lcut(num_words)).value_counts() wc_pic = WordCloud(background_color='white',font_path=r'C:\Windows\Fonts\simhei.ttf').fit_words(num) plt.figure(figsize=(10,10)) plt.imshow(wc_pic) plt.axis('off') plt.show() from sklearn.model_selection import train_test_split from keras.preprocessing import sequence maxlen = 128 vec_data = list(sequence.pad_sequences(data['vec'],maxlen=maxlen)) x,xt,y,yt = train_test_split(vec_data,data['label'],test_size = 0.2,random_state = 123) import numpy as np x = np.array(list(x)) y = np.array(list(y)) xt = np.array(list(xt)) yt = np.array(list(yt)) x=x[:2000,:] y=y[:2000] xt=xt[:500,:] yt=yt[:500] from sklearn.svm import SVC clf = SVC(C=1, kernel = 'linear') clf.fit(x,y) from sklearn.metrics import classification_report test_pre = clf.predict(xt) report = classification_report(yt,test_pre) print(report) from keras.optimizers import SGD, RMSprop, Adagrad from keras.utils import np_utils from keras.models import Sequential from keras.layers.core import Dense, Dropout, Activation from keras.layers.embeddings import Embedding from keras.layers.recurrent import LSTM, GRU model = Sequential() model.add(Embedding(len(num_data['id'])+1,256)) model.add(Dense(32, activation='sigmoid', input_dim=100)) model.add(LSTM(128)) model.add(Dense(1)) model.add(Activation('sigmoid')) model.summary() import matplotlib.pyplot as plt import matplotlib.image as mpimg from keras.utils import plot_model plot_model(model,to_file='Lstm2.png',show_shapes=True) ls = mpimg.imread('Lstm2.png') plt.imshow(ls) plt.axis('off') plt.show() model.compile(loss='binary_crossentropy',optimizer='Adam',metrics=["accuracy"]) model.fit(x,y,validation_data=(x,y),epochs=15)

#导入所需库 import numpy as np import pandas as pd from sklearn.preprocessing import LabelEncoder from keras.utils import to_categorical from keras.models import Sequential from keras.layers import Dense from sklearn.model_selection import KFold #读入数据 train_data = pd.read_csv('ProSeqs_Train.txt', delimiter=' ', header=None) test_data = pd.read_csv('ProSeqs_Test.txt', delimiter=' ', header=None) #预处理训练集数据 X = train_data.iloc[:, 2:].values y = train_data.iloc[:, 1].values le = LabelEncoder() y = le.fit_transform(y) y = to_categorical(y) #定义模型 model = Sequential() model.add(Dense(64, input_dim=X.shape[1], activation='relu')) model.add(Dense(32, activation='relu')) model.add(Dense(2, activation='softmax')) model.compile(loss='categorical_crossentropy', optimizer='adam', metrics=['accuracy']) #K折交叉验证训练模型 kf = KFold(n_splits=5, shuffle=True, random_state=42) fold_scores = [] for train_index, valid_index in kf.split(X): train_X, train_y = X[train_index], y[train_index] valid_X, valid_y = X[valid_index], y[valid_index] model.fit(train_X, train_y, validation_data=(valid_X, valid_y), epochs=50, batch_size=32, verbose=2) fold_scores.append(model.evaluate(valid_X, valid_y, verbose=0)[1]) print('KFold cross-validation accuracy: {:.2f}%'.format(np.mean(fold_scores) * 100)) #预处理测试集数据 test_X = test_data.iloc[:, 1:].values #预测测试集结果 preds = model.predict(test_X) preds = np.argmax(preds, axis=1) #保存预测结果至文件中 np.savetxt('preds.txt', preds, fmt='%d') #输出预测结果 print('Predictions:') print(preds)该蛋白质功能预测实验涉及分类模型的理论基础

最新推荐

recommend-type

Font Awesome图标字体库提供可缩放矢量图标,它可以被定制大小、颜色、阴影以及任何可以用CSS的样式

Font Awesome图标字体库提供可缩放矢量图标,它可以被定制大小、颜色、阴影以及任何可以用CSS的样式
recommend-type

EDAfloorplanning

介绍了physical design的floorplanning问题
recommend-type

数学建模培训资料 数学建模实战题目真题答案解析解题过程&论文报告 最低生活保障问题的探索 共20页.pdf

数学建模培训资料 数学建模实战题目真题答案解析解题过程&论文报告 最低生活保障问题的探索 共20页.pdf
recommend-type

俄罗斯RTSD数据集实现交通标志实时检测

资源摘要信息:"实时交通标志检测" 在当今社会,随着道路网络的不断扩展和汽车数量的急剧增加,交通标志的正确识别对于驾驶安全具有极其重要的意义。为了提升自动驾驶汽车或辅助驾驶系统的性能,研究者们开发了各种算法来实现实时交通标志检测。本文将详细介绍一项关于实时交通标志检测的研究工作及其相关技术和应用。 ### 俄罗斯交通标志数据集(RTSD) 俄罗斯交通标志数据集(RTSD)是专门为训练和测试交通标志识别算法而设计的数据集。数据集内容丰富,包含了大量的带标记帧、交通符号类别、实际的物理交通标志以及符号图像。具体来看,数据集提供了以下重要信息: - 179138个带标记的帧:这些帧来源于实际的道路视频,每个帧中可能包含一个或多个交通标志,每个标志都经过了精确的标注和分类。 - 156个符号类别:涵盖了俄罗斯境内常用的各种交通标志,每个类别都有对应的图像样本。 - 15630个物理符号:这些是实际存在的交通标志实物,用于训练和验证算法的准确性。 - 104358个符号图像:这是一系列经过人工标记的交通标志图片,可以用于机器学习模型的训练。 ### 实时交通标志检测模型 在该领域中,深度学习模型尤其是卷积神经网络(CNN)已经成为实现交通标志检测的关键技术。在描述中提到了使用了yolo4-tiny模型。YOLO(You Only Look Once)是一种流行的实时目标检测系统,YOLO4-tiny是YOLO系列的一个轻量级版本,它在保持较高准确率的同时大幅度减少计算资源的需求,适合在嵌入式设备或具有计算能力限制的环境中使用。 ### YOLO4-tiny模型的特性和优势 - **实时性**:YOLO模型能够实时检测图像中的对象,处理速度远超传统的目标检测算法。 - **准确性**:尽管是轻量级模型,YOLO4-tiny在多数情况下仍能保持较高的检测准确性。 - **易集成**:适用于各种应用,包括移动设备和嵌入式系统,易于集成到不同的项目中。 - **可扩展性**:模型可以针对特定的应用场景进行微调,提高特定类别目标的检测精度。 ### 应用场景 实时交通标志检测技术的应用范围非常广泛,包括但不限于: - 自动驾驶汽车:在自动驾驶系统中,能够实时准确地识别交通标志是保证行车安全的基础。 - 智能交通系统:交通标志的实时检测可以用于交通流量监控、违规检测等。 - 辅助驾驶系统:在辅助驾驶系统中,交通标志的自动检测可以帮助驾驶员更好地遵守交通规则,提升行驶安全。 - 车辆导航系统:通过实时识别交通标志,导航系统可以提供更加精确的路线规划和预警服务。 ### 关键技术点 - **图像处理技术**:包括图像采集、预处理、增强等步骤,为后续的识别模型提供高质量的输入。 - **深度学习技术**:利用深度学习尤其是卷积神经网络(CNN)进行特征提取和模式识别。 - **数据集构建**:构建大规模、多样化的高质量数据集对于训练准确的模型至关重要。 ### 结论 本文介绍的俄罗斯交通标志数据集以及使用YOLO4-tiny模型进行实时交通标志检测的研究工作,显示了在该领域应用最新技术的可能性。随着计算机视觉技术的不断进步,实时交通标志检测算法将变得更加准确和高效,进一步推动自动驾驶和智能交通的发展。
recommend-type

管理建模和仿真的文件

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

预测区间与置信区间:机器学习中的差异与联系

# 1. 机器学习中的统计基础 在当今数据驱动的时代,机器学习已经成为了理解大数据的关键途径。在这一章节中,我们将探索机器学习与统计学之间密不可分的关系,重点介绍统计学在机器学习中的核心地位及其应用。我们将从最基本的统计概念入手,为读者建立起机器学习中的统计基础。 ## 1.1 统计学的核心概念 统计学为我们提供了一套强大的工具,用以描述、分析以及从数据中得出结论。核心概念包括均值、方差、标准差等描述性统计指标,它们是理解数据集基本特征的关键。 ## 1.2 统计推断基础 统计推断是建立在概率论基础上的,允许我们在有限的数据样本上做出关于整体的结论。我们将解释置信区间和假设检验等基本概念
recommend-type

基于KNN通过摄像头实现0-9的识别python代码

基于KNN(K-Nearest Neighbors,最近邻算法)实现摄像头实时抓取图像并识别0-9数字的Python代码需要几个步骤,包括数据预处理、训练模型和实际应用。这里是一个简化版本的示例: ```python # 导入必要的库 import cv2 from sklearn.neighbors import KNeighborsClassifier import numpy as np # 数据预处理:假设你已经有一个包含手写数字的训练集 # 这里只是一个简化的例子,实际情况下你需要一个完整的图像数据集 # X_train (特征矩阵) 和 y_train (标签) X_train
recommend-type

易语言开发的文件批量改名工具使用Ex_Dui美化界面

资源摘要信息:"文件批量改名工具-易语言"是一个专门用于批量修改文件名的软件工具,它采用的编程语言是“易语言”,该语言是为中文用户设计的,其特点是使用中文作为编程关键字,使得中文用户能够更加容易地编写程序代码。该工具在用户界面上使用了Ex_Dui库进行美化,Ex_Dui是一个基于易语言开发的UI界面库,能够让开发的应用程序界面更美观、更具有现代感,增加了用户体验的舒适度。 【易语言知识点】: 易语言是一种简单易学的编程语言,特别适合没有编程基础的初学者。它采用了全中文的关键字和语法结构,支持面向对象的编程方式。易语言支持Windows平台的应用开发,并且可以轻松调用Windows API,实现复杂的功能。易语言的开发环境提供了丰富的组件和模块,使得开发各种应用程序变得更加高效。 【Ex_Dui知识点】: Ex_Dui是一个专为易语言设计的UI(用户界面)库,它为易语言开发的应用程序提供了大量的预制控件和风格,允许开发者快速地制作出外观漂亮、操作流畅的界面。使用Ex_Dui库可以避免编写繁琐的界面绘制代码,提高开发效率,同时使得最终的软件产品能够更加吸引用户。 【开源大赛知识点】: 2019开源大赛(第四届)是指在2019年举行的第四届开源软件开发竞赛活动。这类活动通常由开源社区或相关组织举办,旨在鼓励开发者贡献开源项目,推广开源文化和技术交流,提高软件开发的透明度和协作性。参与开源大赛的作品往往需要遵循开放源代码的许可协议,允许其他开发者自由使用、修改和分发代码。 【压缩包子文件的文件名称列表知识点】: 文件名称列表中包含了几个关键文件: - libexdui.dll:这显然是一个动态链接库文件,即DLL文件,它是由Ex_Dui库提供的,用于提供程序运行时所需的库函数和资源。DLL文件可以让程序调用相应的函数,实现特定的功能。 - 文件批量改名工具.e:这可能是易语言编写的主程序文件,带有.e扩展名,表明它是一个易语言源代码文件。 - Default.ext:这个文件名没有给出具体扩展名,可能是一个配置文件或默认设置文件,用户可以通过修改它来自定义软件的行为。 - Source:这可能是一个包含易语言源代码的目录,里面应该包含了文件批量改名工具的源代码,供开发者阅读和学习。 - Res:这个目录通常用于存放资源文件,如图形、声音等。在易语言项目中,Res目录下可能存放了程序运行所需的各种资源文件。 通过对标题、描述、标签以及文件名列表的分析,我们可以了解到这款文件批量改名工具采用了易语言编程,并且界面通过Ex_Dui库进行美化。它可能被提交到了2019年第四届开源大赛中,是开发者为用户提供的一个实用工具,用于提高处理文件时的效率。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

【机器学习预测区间入门】:从概念到实现

![【机器学习预测区间入门】:从概念到实现](https://img-blog.csdnimg.cn/img_convert/6158c68b161eeaac6798855e68661dc2.png) # 1. 机器学习预测区间的概念理解 在数据科学和机器学习中,预测区间是衡量模型预测不确定性和精确性的重要工具。**预测区间**是一个围绕预测值的范围,它提供了一个概率区间,旨在包含未来观测值的概率,表明模型预测的可信度。 预测区间的概念易于理解,比如在天气预报中,预报员会给出一个温度预测范围,而不是单一的数字,这个范围就是一种预测区间。它表明了在一定置信水平下,未来观测值可能落在的区间内。