版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上该代码为基于模糊神经网络的水质评价代码清空环境变量参数初始化网络训练网络预测嘉陵江实际水质预测清空环境变量clcclear参数初始化xite=0.001;alfa=0.05;%网络节点I=6; %输入节点数M=12; %隐含节点数O=1; %输出节点数%系数初始化p0=0.3*ones(M,1);p0_1=p0;p0_2=p0_1;p1=0.3*ones(M,1);p1_1=p1;p1_2=p1_1;p2=0.3*ones(M,1);p2_1=p2;p2_2=p2_1;p3=0.3*ones(M,1);p3_1=p3;p3_2=p3_1;p4=0.3*ones(M,1
2、);p4_1=p4;p4_2=p4_1;p5=0.3*ones(M,1);p5_1=p5;p5_2=p5_1;p6=0.3*ones(M,1);p6_1=p6;p6_2=p6_1;%参数初始化c=1+rands(M,I);c_1=c;c_2=c_1;b=1+rands(M,I);b_1=b;b_2=b_1;maxgen=100; %进化次数%网络测试数据,并对数据归一化load data1 input_train output_train input_test output_test%选连样本输入输出数据归一化inputn,inputps=mapminmax(input_train);outp
3、utn,outputps=mapminmax(output_train);n,m=size(input_train);网络训练%循环开始,进化网络for iii=1:maxgen iii; for k=1:m x=inputn(:,k); %输出层结算 for i=1:I for j=1:M u(i,j)=exp(-(x(i)-c(j,i)2/b(j,i); end end %模糊规则计算 for i=1:M w(i)=u(1,i)*u(2,i)*u(3,i)*u(4,i)*u(5,i)*u(6,i); end addw=sum(w); for i=1:M yi(i)=p0_1(i)+p1_1
4、(i)*x(1)+p2_1(i)*x(2)+p3_1(i)*x(3)+p4_1(i)*x(4)+p5_1(i)*x(5)+p6_1(i)*x(6); end addyw=yi*w' %网络预测计算 yn(k)=addyw/addw; e(k)=outputn(k)-yn(k); %计算p的变化值 d_p=zeros(M,1); d_p=xite*e(k)*w./addw; d_p=d_p' %计算b变化值 d_b=0*b_1; for i=1:M for j=1:I d_b(i,j)=xite*e(k)*(yi(i)*addw-addyw)*(x(j)-c(i,j)2*w(i)
5、/(b(i,j)2*addw2); end end %更新c变化值 for i=1:M for j=1:I d_c(i,j)=xite*e(k)*(yi(i)*addw-addyw)*2*(x(j)-c(i,j)*w(i)/(b(i,j)*addw2); end end p0=p0_1+ d_p+alfa*(p0_1-p0_2); p1=p1_1+ d_p*x(1)+alfa*(p1_1-p1_2); p2=p2_1+ d_p*x(2)+alfa*(p2_1-p2_2); p3=p3_1+ d_p*x(3)+alfa*(p3_1-p3_2); p4=p4_1+ d_p*x(4)+alfa*(p
6、4_1-p4_2); p5=p5_1+ d_p*x(5)+alfa*(p5_1-p5_2); p6=p6_1+ d_p*x(6)+alfa*(p6_1-p6_2); b=b_1+d_b+alfa*(b_1-b_2); c=c_1+d_c+alfa*(c_1-c_2); p0_2=p0_1;p0_1=p0; p1_2=p1_1;p1_1=p1; p2_2=p2_1;p2_1=p2; p3_2=p3_1;p3_1=p3; p4_2=p4_1;p4_1=p4; p5_2=p5_1;p5_1=p5; p6_2=p6_1;p6_1=p6; c_2=c_1;c_1=c; b_2=b_1;b_1=b; en
7、d E(iii)=sum(abs(e);endfigure(1);plot(outputn,'r')hold onplot(yn,'b')hold onplot(outputn-yn,'g');legend('实际输出','预测输出','误差','fontsize',12)title('训练数据预测','fontsize',12)xlabel('样本序号','fontsize',12)ylabel('水质等级
8、39;,'fontsize',12)Warning: Ignoring extra legend entries. 网络预测%数据归一化inputn_test=mapminmax('apply',input_test,inputps);n,m=size(inputn_test)for k=1:m x=inputn_test(:,k); %计算输出中间层 for i=1:I for j=1:M u(i,j)=exp(-(x(i)-c(j,i)2/b(j,i); end end for i=1:M w(i)=u(1,i)*u(2,i)*u(3,i)*u(4,i)*u
9、(5,i)*u(6,i); end addw=0; for i=1:M addw=addw+w(i); end for i=1:M yi(i)=p0_1(i)+p1_1(i)*x(1)+p2_1(i)*x(2)+p3_1(i)*x(3)+p4_1(i)*x(4)+p5_1(i)*x(5)+p6_1(i)*x(6); end addyw=0; for i=1:M addyw=addyw+yi(i)*w(i); end %计算输出 yc(k)=addyw/addw;end%预测结果反归一化test_simu=mapminmax('reverse',yc,outputps);%作图f
10、igure(2)plot(output_test,'r')hold onplot(test_simu,'b')hold onplot(test_simu-output_test,'g')legend('实际输出','预测输出','误差','fontsize',12)title('测试数据预测','fontsize',12)xlabel('样本序号','fontsize',12)ylabel('水质等级',
11、'fontsize',12)n = 6m = 50Warning: Ignoring extra legend entries. 嘉陵江实际水质预测load data2 hgsc gjhy dxg%-红工水厂-zssz=hgsc;%数据归一化inputn_test =mapminmax('apply',zssz,inputps);n,m=size(zssz);for k=1:1:m x=inputn_test(:,k); %计算输出中间层 for i=1:I for j=1:M u(i,j)=exp(-(x(i)-c(j,i)2/b(j,i); end end
12、for i=1:M w(i)=u(1,i)*u(2,i)*u(3,i)*u(4,i)*u(5,i)*u(6,i); end addw=0; for i=1:M addw=addw+w(i); end for i=1:M yi(i)=p0_1(i)+p1_1(i)*x(1)+p2_1(i)*x(2)+p3_1(i)*x(3)+p4_1(i)*x(4)+p5_1(i)*x(5)+p6_1(i)*x(6); end addyw=0; for i=1:M addyw=addyw+yi(i)*w(i); end %计算输出 szzb(k)=addyw/addw;endszzbz1=mapminmax(&
13、#39;reverse',szzb,outputps);for i=1:m if szzbz1(i)<=1.5 szpj1(i)=1; elseif szzbz1(i)>1.5&&szzbz1(i)<=2.5 szpj1(i)=2; elseif szzbz1(i)>2.5&&szzbz1(i)<=3.5 szpj1(i)=3; elseif szzbz1(i)>3.5&&szzbz1(i)<=4.5 szpj1(i)=4; else szpj1(i)=5; endend% %-高家花园-zssz
14、=gjhy;inputn_test =mapminmax('apply',zssz,inputps);n,m=size(zssz);for k=1:1:m x=inputn_test(:,k); %计算输出中间层 for i=1:I for j=1:M u(i,j)=exp(-(x(i)-c(j,i)2/b(j,i); end end for i=1:M w(i)=u(1,i)*u(2,i)*u(3,i)*u(4,i)*u(5,i)*u(6,i); end addw=0; for i=1:M addw=addw+w(i); end for i=1:M yi(i)=p0_1(i)
15、+p1_1(i)*x(1)+p2_1(i)*x(2)+p3_1(i)*x(3)+p4_1(i)*x(4)+p5_1(i)*x(5)+p6_1(i)*x(6); end addyw=0; for i=1:M addyw=addyw+yi(i)*w(i); end %计算输出 szzb(k)=addyw/addw;endszzbz2=mapminmax('reverse',szzb,outputps);for i=1:m if szzbz2(i)<=1.5 szpj2(i)=1; elseif szzbz2(i)>1.5&&szzbz2(i)<=2
16、.5 szpj2(i)=2; elseif szzbz2(i)>2.5&&szzbz2(i)<=3.5 szpj2(i)=3; elseif szzbz2(i)>3.5&&szzbz2(i)<=4.5 szpj2(i)=4; else szpj2(i)=5; endend% %-大溪沟水厂-zssz=dxg;inputn_test =mapminmax('apply',zssz,inputps);n,m=size(zssz);for k=1:1:m x=inputn_test(:,k); %计算输出中间层 for i=1:
17、I for j=1:M u(i,j)=exp(-(x(i)-c(j,i)2/b(j,i); end end for i=1:M w(i)=u(1,i)*u(2,i)*u(3,i)*u(4,i)*u(5,i)*u(6,i); end addw=0; for i=1:M addw=addw+w(i); end for i=1:M yi(i)=p0_1(i)+p1_1(i)*x(1)+p2_1(i)*x(2)+p3_1(i)*x(3)+p4_1(i)*x(4)+p5_1(i)*x(5)+p6_1(i)*x(6); end addyw=0; for i=1:M addyw=addyw+yi(i)*w(
18、i); end %计算输出 szzb(k)=addyw/addw;endszzbz3=mapminmax('reverse',szzb,outputps);for i=1:m if szzbz3(i)<=1.5 szpj3(i)=1; elseif szzbz3(i)>1.5&&szzbz3(i)<=2.5 szpj3(i)=2; elseif szzbz3(i)>2.5&&szzbz3(i)<=3.5 szpj3(i)=3; elseif szzbz3(i)>3.5&&szzbz3(i)<=4.5 szpj3(i)=4; else szpj3(i)=5; endendfigure(3)plot(szzbz1,'o-r')h
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 一次性剃须刀行业市场发展分析及发展前景与投资机会研究报告
- 2024-2030年棕榈酸甲酯搬迁改造项目可行性研究报告
- 2024-2030年住宅用水表行业市场现状供需分析及重点企业投资评估规划分析研究报告
- 2024-2030年中小学教材行业市场发展现状及竞争形势与投资前景研究报告
- 2024-2030年中国音乐制作软件行业市场现状供需分析及市场深度研究发展前景及规划战略投资分析研究报告
- 2024-2030年中国阻燃复合材料行业供需态势及投资效益预测报告
- 2024-2030年中国针炙针行业市场现状分析及竞争格局与投资发展研究报告
- 2024-2030年中国重组人表皮生长因子溶液行业市场现状分析及竞争格局与投资发展研究报告
- 2024-2030年中国通信系统及软件行业市场前瞻与未来投资战略分析报告
- 2024-2030年中国远程数字视频采集软件行业发展策略与需求趋势预测研究报告
- 减速器毕业设计答辩(课堂PPT)
- 离心泵主要零部件的强度计算
- 制浆洗漂详细过程工艺
- iatf16949纠正与预防措施管理程序-
- 国家开放大学《管理英语1》边学边练参考答案
- 钳工用扁锉热处理工艺(T12)碳素钢
- 小型河流开挖穿越施工方案_0
- 不锈钢管道焊接工艺标准规范标准
- 2021《反假货币知识竞赛》题库(试题82道含答案)
- 异氰酸酯的化学反应.PPT
- 2022年2022年北师大版二年级下册数学竞赛题
评论
0/150
提交评论