24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3193  |  回复: 7

fanfanjia

新虫 (初入文坛)

[求助] 请高手指教如何提高这个BP神经网络的预测精度

请高人指点,以下为我为毕业论文设计用MATLAB做的BP神经网络,目前存在两个问题:
1.模型预测精度特别低。
我已经尝试改变隐含层节点数(从20到100都一一试过了),并将传递函数和训练函数都做了替换,但均无明显效果。不知该从何处着手提高模型的预测和训练精度。
2.每次程序运行完,都要提示错误,貌似是运算符合不对,但是改过之后依然有误。不知该怎么改?


clc;
clear;
close all;
warning off;
tic;
%用2012年数据创建模型,2011年数据验证模型精确度。 自动读取文本文档,第一行为农学参数,第二行至最后一行都为光谱值。
data1=importdata('E:\2012SR.txt');
data2=importdata('E:\2011SR.txt');

M1 = size(data1,2);%行      
Y1 = data1(1,1:M1);%第一行所有列,农学参数,2012

M2 = size(data2,2);      
Y2 = data2(1,1:M2);%第一行所有列,农学参数,2011

P =data1(2:size(data1,1),: );      %2012年光谱数据,文本里第二行开始,所有列
T = Y1 ;    %2012年农学参数,第一行所有列
N2011_validation=data2(2:size(data2,1),: ) ;  %2011年光谱数据,文本里第二行开始,所有列
N2012_data=P;       %2012年光谱数据

net=newff(P,T,[61 1],{'tansig' 'purelin'},'traingd')
t1                      = clock;                              %计时开始
net                     = fitnet(70);
net.trainParam.epochs   = 5000;                               %设置训练次数
net.trainParam.goal     = 0.01;                            %设置性能函数(训练要求精度)
net.trainParam.show     = 1;                                  %每10显示
net.trainParam.Ir       = 0.005;                              %设置学习速率
net                     = train(net,P,T);                  %训练BP网络(P为输入,T为输出)
datat                   = etime(clock,t1);
Nets                    = net;
view(Nets);

%2011年数据进行验证
y = sim(net,N2011_validation);
figure;
plot(y,'k-x');hold on
plot(N2011_validation(1:126),'.k');hold on
ymin=-1;
ymax=4;
axis([0 length(N2011_validation) floor(ymin) ceil(ymax)]);
legend('Pridicted values','Measured values');
title('模型预测结果, 烟叶N含量');
%2012年数据进行测试
y1=sim(net,N2012_data);
figure;
plot(y1,'k-x');
hold on;
plot(N2012_data(1:102),'.k');
ymin=0;
ymax=4;
axis([0 length(N2012_data) floor(ymin) ceil(ymax)]);
legend('Pridicted values','Measured values');
title('模型训练结果, 烟叶N含量');
%导入到文本文件
fid=fopen('E:\2012SR1.txt','wt');
if fid == -1
    error('文件打开失败');
end
fprintf(fid,'实际值 预测值\n');
A=[N2012_data;y1];
fprintf(fid,'%f %f\n',A);
fclose(fid);
fid=fopen('E:\2011SR1.txt','wt');
if fid == -1
    error('文件打开失败');
end
fprintf(fid,'实际值 预测值\n');
A=[N2011_validation;y];
fprintf(fid,'%f %f\n',A);
fclose(fid);
%预测效果分析
for i = 1:length(N2011_validation);
    K(i)  = abs(y(i) - N2011_validation(i))/(N2011_validation(i));
end
RMSE = sqrt(sum((y-N2011_validation)*(y-N2011_validation)')/length(N2011_validation))
E= 1-sum((y-N2011_validation)*(y-N2011_validation)')/sum((N2011_validation-mean(N2011_validation))*(N2011_validation-mean(N2011_validation))')
k= average(sum(abs((y-N2011_validation)/(N2011_validation))))
%训练效果分析
for i = 1:length(N2012_data);
    K_1(i)  = average(sum(abs((y1(i)-N2012_data(i))/(N2012_data(i)))));
end
RMSE_1=sqrt(sum((y1-N2012_data)*(y1-N2012_data)')/length(N2012_data))
E_1= 1-sum((y1-N2012_data)*(y1-N2012_data)')/sum((N2012_data-mean(N2012_data))*(N2012_data-mean(N2012_data))')
k_1 = average(sum(abs((y1-N2012_data)/(N2012_data))))


预测出来的结果总是整体偏高,不知道跟我没有进行数据归一化有没有关系?如果进行数据归一化和反归一化,该在哪个地方怎么添加程序呢?请高手指教[ Last edited by jjdg on 2013-1-23 at 13:03 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 2012SR.txt
  • 2012-12-31 17:40:22, 23.52 K
  • 附件 2 : 2011SR.txt
  • 2012-12-31 17:40:33, 30.54 K

» 收录本帖的淘帖专辑推荐

matlab

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanfanjia

新虫 (初入文坛)

我自己看书已经解决了
2楼2013-01-10 21:32:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

脱皮的洋葱

木虫 (小有名气)

引用回帖:
2楼: Originally posted by fanfanjia at 2013-01-10 21:32:02
我自己看书已经解决了

楼主 能分享一下经验吗 最近也在接触 bp 神经网络。。。也有预测精度这方面的问题。。。
3楼2013-01-10 23:22:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wgq840720

新虫 (初入文坛)

【答案】应助回帖

班主你好 我想写关于带式焙烧机热工方面的建模  得用到bp神经网络和MATLAB,请求助,如果能帮忙写2张 我会给予相关报酬
4楼2013-01-23 12:15:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zoro1101

禁虫 (小有名气)

本帖内容被屏蔽

5楼2013-09-27 16:02:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lanxinldw

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by fanfanjia at 2013-01-10 21:32:02
我自己看书已经解决了

楼主,我想问下你是怎么解决的啊?我现在的预测精度也不高,你是用改变归一化函数解决的么?
6楼2014-10-30 10:52:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ning_fen

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by fanfanjia at 2013-01-10 21:32:02
我自己看书已经解决了

楼主,请指点下您的解决办法吧,我是菜鸟一只,刚刚接触这方面的知识,请指教。
7楼2015-05-19 11:16:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小红叶

新虫 (初入文坛)

楼主,求相助  如何做到了
8楼2015-07-24 15:11:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fanfanjia 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 266分,求材料冶金能源化工等调剂 +7 哇呼哼呼哼 2026-03-27 9/450 2026-03-28 12:22 by zllcz
[考研] 289求调剂 +10 新时代材料 2026-03-27 10/500 2026-03-28 09:53 by 无际的草原
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +5 dongfang59 2026-03-27 5/250 2026-03-28 04:53 by wxiongid
[考研] 330一志愿中国海洋大学 化学工程 085602 有读博意愿 求调剂 +3 wywy.. 2026-03-27 4/200 2026-03-28 03:32 by fmesaito
[考研] 求调剂 +4 零八# 2026-03-27 4/200 2026-03-27 18:07 by yu221
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 322求调剂 +4 我真的很想学习 2026-03-23 4/200 2026-03-27 13:51 by 杨杨杨紫
[考研] 324求调剂 +8 hanamiko 2026-03-26 10/500 2026-03-27 08:06 by hypershenger
[考研] 求调剂 一志愿 本科 北科大 化学 343 +6 13831862839 2026-03-24 7/350 2026-03-26 22:57 by 不吃魚的貓
[考研] 294分080500材料科学与工程求调剂 +4 柳溪边 2026-03-26 4/200 2026-03-26 21:14 by XPU李庆
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 085600 材料与化工 329分求调剂 +9 Mr. Z 2026-03-25 9/450 2026-03-26 10:36 by baoball
[考研] 材料与化工328分调剂 +6 。,。,。,。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +4 我爱学电池 2026-03-23 4/200 2026-03-25 00:59 by 1027_324
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
信息提示
请填处理意见