clc
clear?all
%****************************
%BP神經網絡，6-20-1網絡結構
%預測未來10年的各水質河段比例
%****************************
tic?
first10_data=load（‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\前十年排污量.txt‘）;
first10_rate=load（‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\前十年6類水質河長比例.txt‘）;
load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\forecast_data.mat‘;

%***************歸一化
PW_DATA=zeros（201）;%20年的排污量
PW_DATA（1:101）=first10_data;
PW_DATA（11:201）=forecast_data;
[first10_scal_ratefs]?=?mapminmax（first10_rate01）;?%前10年6種水質河段比例歸一化?
[PW_scal_DATAds]?=?mapminmax（?PW_DATA‘02）;%后20年排污量歸一化?
PW_scal_DATA=PW_scal_DATA‘;

%*****訓練初始化**********
yeta=0.2;%學習率
yeta1=0.055;%防止激活函數進入飽和區
train_time=20000;%最大訓練次數
Input_rate=zeros（206）;
Input_rate（1:10:）=first10_scal_rate;?%存貯各種水質20年的比例

for?class=1:6????%6種水質
????if?class==1
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\W1.mat‘;
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\V1.mat‘;??
????elseif?class==2
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\W2.mat‘;
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\V2.mat‘;
????elseif?class==3
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\W3.mat‘;
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\V3.mat‘;
?????elseif?class==4
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\W4.mat‘;
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\V4.mat‘;
????elseif?class==5
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\W5.mat‘;
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\V5.mat‘;
????elseif?class==6
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\W6.mat‘;
??????load?‘E:\研究生階段\研一時期\下學期\學習與實驗\課程\模式識別\作業二\mat\V6.mat‘;?
????end
????
??????%***********預測第?forecast_index?年****************
??for?forecast_index=11:20???%第11年到第20?年
????%得到隱層輸出
????y0=[PW_scal_DATA（forecast_index1）;Input_rate（forecast_index-5:forecast_index-1class）]‘*V;%
????%激活
????y=1./（1+exp（-y0*yeta1））;
????
????%得到輸出層輸出
????output0=y*W;%
????output=1./（1+exp（-output0*yeta1））;%1*1
????Input_rate（forecast_indexclass）=output;%得到預測值
????
?????????%*************訓練開始*********************
???%**以上述預測值為樣本再次進行訓練，更新權值
???%**為預測下一年做準備
????V=double（rands（620））;%初始化權重（輸入→隱含層）
????W=double（rands（201））;%初始化權重（隱含層→輸入）
???for?time=?1:train_time
??????for?years=1:forecast_index-5???%1?-?15
?????????
????????%得到隱層輸入
?????????y0=[PW_scal_DATA（years+51）;Input_rate（years:years+4class）]‘*V;
????????%激活
?????????y=1./（1+exp（-y0*yeta1））;%1*20
????
????????%得到輸出層輸入
????????output0=y*W;
????????%激活
????????output=1./（1+exp（-output0*yeta1））;%1*1
????
????????%誤差反傳
???????delet=（Input_rate（years+5class）-output）*output*（1-output）;%1*1誤差??
???????delet_W=yeta*y‘.*repmat（delet201）;%20*1?
???????tmp=（W.*repmat（delet201））‘;%1*20??
???????tmp=tmp.*y.*（1-y）;?%1*20?
???????delet_V=yeta*repmat（（[PW_scal_DATA（years+51）;Input_rate（years:years+4class）]‘）‘120）.*repmat（tmp61）;%?
???????
????%更新權值?
????W=W+delet_W;??
????V=V+delet_V;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-14?16:22??長江水質預測（2005年數學建模）+神經網絡\
?????文件?????1275904??2018-04-12?18:27??長江水質預測（2005年數學建模）+神經網絡\1703080202002_王齊勝；1708085202002馬衛衛.doc
?????文件??????455321??2018-04-09?10:05??長江水質預測（2005年數學建模）+神經網絡\1704081002005_王程偉—1704085208013_嚴微.doc
?????文件??????418575??2012-07-09?17:24??長江水質預測（2005年數學建模）+神經網絡\2005年中國大學生數學建模競賽論文（長江水質的評價和趨勢分析模型）II.pdf.pdf
?????文件????????3604??2018-04-14?11:20??長江水質預測（2005年數學建模）+神經網絡\BP_forecas_train_1.m
?????文件????????3916??2018-04-06?22:21??長江水質預測（2005年數學建模）+神經網絡\BP_forecast_test.m
?????文件????????3876??2018-04-08?19:26??長江水質預測（2005年數學建模）+神經網絡\BP_forecast_test_1.m
?????文件????????3522??2018-04-08?19:28??長江水質預測（2005年數學建模）+神經網絡\BP_forecast_train.m
?????文件????????1997??2018-04-08?21:36??長江水質預測（2005年數學建模）+神經網絡\BP_pw_test.m
?????文件????????2400??2018-04-08?21:34??長江水質預測（2005年數學建模）+神經網絡\BP_pw_train.m
?????文件????????1597??2018-04-08?21:25??長江水質預測（2005年數學建模）+神經網絡\GM.m
?????目錄???????????0??2018-04-06?22:10??長江水質預測（2005年數學建模）+神經網絡\mat\
?????文件?????????963??2018-04-06?18:12??長江水質預測（2005年數學建模）+神經網絡\mat\Input_rate.mat
?????文件?????????445??2018-04-06?14:44??長江水質預測（2005年數學建模）+神經網絡\mat\OUTPUT.mat
?????文件?????????298??2018-04-06?16:14??長江水質預測（2005年數學建模）+神經網絡\mat\OUTPUT1.mat
?????文件?????????682??2018-04-06?18:10??長江水質預測（2005年數學建模）+神經網絡\mat\Reversed.mat
?????文件?????????467??2018-04-06?16:16??長江水質預測（2005年數學建模）+神經網絡\mat\Reversed1.mat
?????文件????????1160??2018-04-06?14:43??長江水質預測（2005年數學建模）+神經網絡\mat\V1.mat
?????文件????????1158??2018-04-06?14:43??長江水質預測（2005年數學建模）+神經網絡\mat\V2.mat
?????文件????????1163??2018-04-06?14:43??長江水質預測（2005年數學建模）+神經網絡\mat\V3.mat
?????文件????????1163??2018-04-06?14:43??長江水質預測（2005年數學建模）+神經網絡\mat\V4.mat
?????文件????????1163??2018-04-06?14:44??長江水質預測（2005年數學建模）+神經網絡\mat\V5.mat
?????文件????????1163??2018-04-06?14:44??長江水質預測（2005年數學建模）+神經網絡\mat\V6.mat
?????文件????????1146??2018-04-06?16:14??長江水質預測（2005年數學建模）+神經網絡\mat\VV1.mat
?????文件????????1163??2018-04-06?16:14??長江水質預測（2005年數學建模）+神經網絡\mat\VV2.mat
?????文件????????1163??2018-04-06?16:14??長江水質預測（2005年數學建模）+神經網絡\mat\VV3.mat
?????文件????????1163??2018-04-06?16:14??長江水質預測（2005年數學建模）+神經網絡\mat\VV4.mat
?????文件????????1163??2018-04-06?16:14??長江水質預測（2005年數學建模）+神經網絡\mat\VV5.mat
?????文件????????1157??2018-04-06?16:14??長江水質預測（2005年數學建模）+神經網絡\mat\VV6.mat
?????文件????????1163??2018-04-06?19:07??長江水質預測（2005年數學建模）+神經網絡\mat\V_1.mat
?????文件????????1153??2018-04-06?19:07??長江水質預測（2005年數學建模）+神經網絡\mat\V_10.mat
............此處省略53個文件信息