%%?該代碼為基于模糊神經(jīng)網(wǎng)絡(luò)的水質(zhì)評(píng)價(jià)代碼
%
%?
%?			該案例作者申明：		1：本人長(zhǎng)期駐扎在此板塊里，對(duì)該案例提問，做到有問必答。本套書籍官方網(wǎng)站為：video.ourmatlab.com		2：點(diǎn)此從當(dāng)當(dāng)預(yù)定本書：《Matlab神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析》。	3：此案例有配套的教學(xué)視頻，視頻下載方式video.ourmatlab.com/vbuy.html。?
								4：此案例為原創(chuàng)案例，轉(zhuǎn)載請(qǐng)注明出處（《Matlab神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析》）。							5：若此案例碰巧與您的研究有關(guān)聯(lián)，我們歡迎您提意見，要求等，我們考慮后可以加在案例里。			
%?
%%?清空環(huán)境變量
clc
clear

%%?參數(shù)初始化
xite=0.001;
alfa=0.05;

%網(wǎng)絡(luò)節(jié)點(diǎn)
I=6;???%輸入節(jié)點(diǎn)數(shù)
M=12;??%隱含節(jié)點(diǎn)數(shù)
O=1;???%輸出節(jié)點(diǎn)數(shù)

%系數(shù)初始化
p0=0.3*ones（M1）;p0_1=p0;p0_2=p0_1;
p1=0.3*ones（M1）;p1_1=p1;p1_2=p1_1;
p2=0.3*ones（M1）;p2_1=p2;p2_2=p2_1;
p3=0.3*ones（M1）;p3_1=p3;p3_2=p3_1;
p4=0.3*ones（M1）;p4_1=p4;p4_2=p4_1;
p5=0.3*ones（M1）;p5_1=p5;p5_2=p5_1;
p6=0.3*ones（M1）;p6_1=p6;p6_2=p6_1;

%參數(shù)初始化
c=1+rands（MI）;c_1=c;c_2=c_1;
b=1+rands（MI）;b_1=b;b_2=b_1;

maxgen=100;?%進(jìn)化次數(shù)

%網(wǎng)絡(luò)測(cè)試數(shù)據(jù)，并對(duì)數(shù)據(jù)歸一化
load?data1?input_train?output_train?input_test?output_test

%選連樣本輸入輸出數(shù)據(jù)歸一化
[inputninputps]=mapminmax（input_train）;
[outputnoutputps]=mapminmax（output_train）;
[nm]=size（input_train）;

%%?網(wǎng)絡(luò)訓(xùn)練
%循環(huán)開始，進(jìn)化網(wǎng)絡(luò)
for?iii=1:maxgen
????iii;
????for?k=1:m????????
????????x=inputn（:k）;
????????
????????%輸出層結(jié)算
????????for?i=1:I
????????????for?j=1:M
????????????????u（ij）=exp（-（x（i）-c（ji））^2/b（ji））;
????????????end
????????end
????????
????????%模糊規(guī)則計(jì)算
????????for?i=1:M
????????????w（i）=u（1i）*u（2i）*u（3i）*u（4i）*u（5i）*u（6i）;
????????end????
????????addw=sum（w）;
????????
????????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=yi*w‘;
????????%網(wǎng)絡(luò)預(yù)測(cè)計(jì)算
????????yn（k）=addyw/addw;
????????e（k）=outputn（k）-yn（k）;
????????
????????%計(jì)算p的變化值
????????d_p=zeros（M1）;
????????d_p=xite*e（k）*w./addw;
????????d_p=d_p‘;
????????
????????%計(jì)算b變化值
????????d_b=0*b_1;
????????for?i=1:M
????????????for?j=1:I
????????????????d_b（ij）=xite*e（k）*（yi（i）*addw-addyw）*（x（j）-c（ij））^2*w（i）/（b（ij）^2*addw^2）;
????????????end
????????end??
????????
????????%更新c變化值
????????for?i=1:M
????????????for?j=1:I
????????????????d_c（ij）=xite*e（k）*（yi（i）*addw-addyw）*2*（x（j）-c（ij））*w（i）/（b（ij）*addw^2）;
????