%%?BP-DEMATEL算法設(shè)計(jì)
clear;clc;close?all
%%?數(shù)據(jù)處理
A=xlsread（‘?dāng)?shù)據(jù).xlsx‘‘Sheet2‘）;
%行為指標(biāo)?；?列為樣本
T=xlsread（‘?dāng)?shù)據(jù).xlsx‘‘Sheet2‘）;
%對(duì)應(yīng)的元素個(gè)數(shù)為樣本個(gè)數(shù)
P=（A-mean（A））./std（A）;??%?標(biāo)準(zhǔn)化數(shù)據(jù)
[RRCC]=size（A）;
hiddennum=sqrt（RR+1）;
%%?利用minmax函數(shù)求輸入樣本范圍
net=newff（minmax（P）[hiddennum1]{‘tansig‘‘purelin‘}‘traingdx‘）;
%%?BP神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過(guò)程
net.trainParam.show=50;??????%?顯示訓(xùn)練迭代過(guò)程
net.trainParam.lr=0.05;??????%?學(xué)習(xí)率
net.trainParam.epochs=3000;???%?最大訓(xùn)練次數(shù)
net.trainParam.goal=1e-5;????%?訓(xùn)練要求精度
net.trainParam.mc=0.5;???????%?動(dòng)量因子
[nettr]=train（netPT）;?????%?網(wǎng)絡(luò)訓(xùn)練?，P為輸入數(shù)據(jù)（輸入層）；T為目標(biāo)層（輸出層）；訓(xùn)練記錄為tr
%%?計(jì)算權(quán)值與閾值
W1=?net.iw{11};???????%輸入層到中間層的權(quán)值
B1=?net.b{1};??????????%隱含層神經(jīng)元閾值
W2=?net.lw{21};???????%隱含層到輸出層的權(quán)值
B2=?net.b{