应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » BP神经网络进行水质评价遇到问题
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1231  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kathy2836

铁杆木虫 (小有名气)

[求助] BP神经网络进行水质评价遇到问题 已有1人参与

用BP神经网络进行水质评价，选取6个水质因子进行评价，通过rand函数在水质标准各等级标准值区间分别随机生成400个数值。因为标准中添加了劣五类水质，所以共随机生成2400组数据，其中每个等级区间内选取300个用作训练样本，其余100个为测试样本。输出目标格式为：
T=[1 0 0 0 0 0；0 1 0 0 0 0；0 0 1 0 0 0；0 0 0 1 0 0；0 0 0 0 1 0；0 0 0 0 0 1]’。
输入神经元为6个，输出神经元为6个。通过试错法得到隐含层神经元数量为10的时候误差最小。训练网络的时候也能达到设定的目标误差，但是再用水质标准值归一化后通过训练好的网络得到的输出却和预计输出差很远。代码如下，请高手指点一下问题出在哪里，十分感谢！

****************************************************************************************************
% 6个因子的水质标准，包含劣五类水质
P=[7.5 15 2 0.015 0.2 0.01;6 15 4 0.5 0.5 0.025;5 20 6 1 1 0.05;
    3 30 10 1.5 1.5 0.1;2 40 15 2 2 0.2;0 200 25 8 40 1]';
%随机生成训练样本
P1=[(15-P(1,1))*rand(1,400)+P(1,1);P(2,1)*rand(1,400);P(3,1)*rand(1,400);
    P(4,1)*rand(1,400);P(5,1)*rand(1,400);P(6,1)*rand(1,400)];
P2=[(P(1,1)-P(1,2))*rand(1,400)+P(1,2);P(2,2)*rand(1,400);(P(3,2)-P(3,1))*rand(1,400)+P(3,1);
    (P(4,2)-P(4,1))*rand(1,400)+P(4,1);(P(5,2)-P(5,1))*rand(1,400)+P(5,1);(P(6,2)-P(6,1))*rand(1,400)+P(6,1)];
P3=[(P(1,2)-P(1,3))*rand(1,400)+P(1,3);(P(2,3)-P(2,2))*rand(1,400)+P(2,2);(P(3,3)-P(3,2))*rand(1,400)+P(3,2);
    (P(4,3)-P(4,2))*rand(1,400)+P(4,2);(P(5,3)-P(5,2))*rand(1,400)+P(5,2);(P(6,3)-P(6,2))*rand(1,400)+P(6,2)];
P4=[(P(1,3)-P(1,4))*rand(1,400)+P(1,4);(P(2,4)-P(2,3))*rand(1,400)+P(2,3);(P(3,4)-P(3,3))*rand(1,400)+P(3,3);
    (P(4,4)-P(4,3))*rand(1,400)+P(4,3);(P(5,4)-P(5,3))*rand(1,400)+P(5,3);(P(6,4)-P(6,3))*rand(1,400)+P(6,3)];
P5=[(P(1,4)-P(1,5))*rand(1,400)+P(1,5);(P(2,5)-P(2,4))*rand(1,400)+P(2,4);(P(3,5)-P(3,4))*rand(1,400)+P(3,4);
    (P(4,5)-P(4,4))*rand(1,400)+P(4,4);(P(5,5)-P(5,4))*rand(1,400)+P(5,4);(P(6,5)-P(6,4))*rand(1,400)+P(6,4)];
P6=[(P(1,5)-P(1,6))*rand(1,400)+P(1,6);(P(2,6)-P(2,5))*rand(1,400)+P(2,5);(P(3,6)-P(3,5))*rand(1,400)+P(3,5);
    (P(4,6)-P(4,5))*rand(1,400)+P(4,5);(P(5,6)-P(5,5))*rand(1,400)+P(5,5);(P(6,6)-P(6,5))*rand(1,400)+P(6,5)];
%随机生成的样本矩阵（6*2400）
PP=[P1 P2 P3 P4 P5 P6];
%各选取300个做为训练样本
PI=[P1(:,1:300) P2(:,1:300) P3(:,1:300) P4(:,1:300) P5(:,1:300) P6(:,1:300)];
%其余100个做为测试样本
PT=[P1(:,301:400) P2(:,301:400) P3(:,301:400) P4(:,301:400) P5(:,301:400) P6(:,301:400)];

%读取目标输出矩阵文件（6*2400）
T=xlsread('C:\Users\Kathy\Desktop\ANN\T.xlsx','Sheet1');
%对应训练样本和测试样本的输出矩阵
TI=[T(:,1:300) T(:,401:700) T(:,801:1100) T(:,1201:1500) T(:,1601:1900) T(:,2001:2300)];
TTEST=[T(:,301:400) T(:,701:800) T(:,1101:1200) T(:,1501:1600) T(:,1901:2000) T(:,2301:2400)];

%将训练样本归一化至[-1,1]区间
[PN,ps]=mapminmax(PI);

%建立训练网络。这里我开始用的是tansig激活函数，但是三条线的收敛效果很差，所以改成了logsig。不知道有没有问题？是不是归一化区间为[-1,1]要对应使用tansig函数？
net6=newff(minmax(PN),TI,[10,6],{'logsig','purelin'},'trainlm');
%设定输入层权重和阈值
inputWeight6=net6.IW{1,1};
inputbias6=net6.b{1};
%设定当前层权重和阈值
layerWeight6=net6.LW{2,1};
layerbias6=net6.b{2};
% 设定训练参数
net6.trainParam.show=5;
net6.trainParam.lr=0.05;
net6.trainParam.mc=0.95;
net6.trainParam.epochs=3000;
net6.trainParam.goal=0.0001;
%对网络进行训练
net6=train(net6,PN,TI);

A6=sim(net6,PN);

E6=TI-A6;
MSE6=mse(E6);
******************************************************************************
通过以上步骤之后，得到的A6矩阵很满意，和预计输出差不多，但是当我用这个网络对水质标准值（就是上面代码中的P矩阵）归一化之后进行操作时，得到的结果却和预计输出差很多，

对30个站位的水质监测值进行同样操作，结果也是有一些站位的值异常。

请高手给予指点，问题出在哪里呢？困扰我很久了，十分感谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-02-04 11:28:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xl2420

银虫 (小有名气)
楼主犯了个小错误，P=[7.5 15 2 0.015 0.2 0.01;6 15 4 0.5 0.5 0.025;5 20 6 1 1 0.05;
    3 30 10 1.5 1.5 0.1;2 40 15 2 2 0.2;0 200 25 8 40 1]';中第一行第四列应该为0.15，不是0.015
赞一下 回复此楼
xiaoling
3楼2015-07-15 17:12:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

hakuna

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
kathy2836: 金币+5, ★有帮助, 谢谢你！ 2014-04-26 22:53:04
没看太明白，但 大致知道你做的是个分类问题。
考虑一下两个方面：
1.你的描述水质的指标是否合理？难道你的训练集不是测定的数据样本？训练集是否有代表性？
2.判断水质类别与定量预测不同，不许要过度追求寻两误差有多小，只要训练分类合理即可，过度要求小的训练误差会导致训练过度，降低网络的泛化能力。
赞一下 回复此楼
2楼2014-02-06 20:29:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伟长

新虫 (初入文坛)
楼主，你解决了吗？求助求助

发自小木虫Android客户端
赞一下 回复此楼
4楼2018-01-28 07:25:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 材料专硕调剂 +9 哈哈哈吼吼吼哈 2026-04-07 10/500 2026-04-07 22:48 by 来看流星雨10
[考研] 一志愿085404，总分291，四级已过，求调剂 +7 阿俊阿俊阿俊 2026-04-04 9/450 2026-04-07 22:25 by dick_runner
[考研] 0703调剂，一志愿天津大学319分 +19 haaaabcd 2026-04-05 21/1050 2026-04-07 22:18 by dick_runner
[考研] 304求调剂（085602，过四级，一志愿985） +25 化工人999 2026-04-04 26/1300 2026-04-07 22:06 by hemengdong
[基金申请] 请问共同通讯和共同一作的认可度问题 10+5 psa1234 2026-04-01 11/550 2026-04-07 16:57 by caiqing
[考研] 一志愿华东理工085601材料工程303分求调剂 +8 a1708 2026-04-06 8/400 2026-04-07 11:20 by 诗与自由
[考研] 085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干 +3 1021075758 2026-04-06 4/200 2026-04-07 09:23 by 蓝云思雨
[硕博家园] 0856材料化工求调剂，一志愿211，初试成绩349 +4 江淮北月 2026-04-05 4/200 2026-04-06 22:44 by chenzhimin
[考研] 材料工程调剂 +3 小刘同学吖吖 2026-04-06 4/200 2026-04-06 20:56 by lbsjt
[考研] 第一志愿东南大学物理313，有科研竞赛获奖经历，希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +15 路痴小琪 2026-04-05 15/750 2026-04-06 09:27 by cql1109
[考研] 一志愿哈尔滨工业大学085600英一数二337分求调剂 +12 lyz0427 2026-04-03 12/600 2026-04-06 06:37 by houyaoxu
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 一志愿北京科技大学材料工程085601，求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
信息提示
关闭
请填处理意见
关闭