%%?案例20:神經網絡變量篩選—基于BP的神經網絡變量篩選
%
%
%?
%?			該案例作者申明：				1：本人長期駐扎在此板塊里，對該案例提問，做到有問必答。	2：此案例有配套的教學視頻，配套的完整可運行Matlab程序。						3：以下內容為該案例的部分內容（約占該案例完整內容的1/10）。							4：此案例為原創案例，轉載請注明出處（Matlab中文論壇，《Matlab神經網絡30個案例分析》）。							5：若此案例碰巧與您的研究有關聯，我們歡迎您提意見，要求等，我們考慮后可以加在案例里。							6：您看到的以下內容為初稿，書籍的實際內容可能有少許出入，以書籍實際發行內容為準。					7：此書其他常見問題、預定方式等，請點擊這里。	
%?
%
%%?清空環境變量
clc
clear
%%?產生輸入?輸出數據

%?設置步長
interval=0.01;

%?產生x1?x2
x1=-1.5:interval:1.5;
x2=-1.5:interval:1.5;

%?產生x3?x4（噪聲）
x=rand（1301）;
x3=（x-0.5）*1.5*2;
x4=（x-0.5）*1.5*2;

%?按照函數先求得相應的函數值，作為網絡的輸出。
F?=20+x1.^2-10*cos（2*pi*x1）+x2.^2-10*cos（2*pi*x2）;

%設置網絡輸入輸出值
p=[x1;x2;x3;x4];
t=F;


%%?變量篩選?MIV算法的初步實現（增加或者減少自變量）

p=p‘;
[mn]=size（p）;
yy_temp=p;

%?p_increase為增加10%的矩陣?p_decrease為減少10%的矩陣
for?i=1:n
????p=yy_temp;
????pX=p（:i）;
????pa=pX*1.1;
????p（:i）=pa;
????aa=[‘p_increase‘??int2str（i）?‘=p‘];
????eval（aa）;
end


for?i=1:n
????p=yy_temp;
????pX=p（:i）;
????pa=pX*0.9;
????p（:i）=pa;
????aa=[‘p_decrease‘?int2str（i）?‘=p‘];
????eval（aa）;
end


%%?利用原始數據訓練一個正確的神經網絡
nntwarn?off;

p=p‘;
%?bp網絡建立
net=newff（minmax（p）[81]{‘tansig‘‘purelin‘}‘traingdm‘）;
%?初始化bp網絡
net=init（net）;
%?網絡訓練參數設置
net.trainParam.show=50;
net.trainParam.lr=0.05;
net.trainParam.mc=0.9;
net.trainParam.epochs=2000;

%?bp網絡訓練
net=train（netpt）;


%%?變量篩選?MIV算法的后續實現（差值計算）

%?轉置后sim

for?i=1:n
????eval（[‘p_increase‘num2str（i）‘=transpose（p_increase‘num2str（i）‘）‘]）
end

for?i=1:n
????eval（[‘p_decrease‘num2str（i）‘=transpose（p_decrease‘num2str（i）‘）‘]）
end


%?result_in為增加10%后的輸出?result_de為減少10%后的輸出
for?i=1:n
????eval（[‘result_in‘num2str（i）‘=sim（net‘‘p_increase‘num2str（i）‘）‘]）
end

for?i=1:n
????eval（[‘result_de‘num2str（i）‘=sim（net‘‘p_decrease‘num2str（i）‘）‘]）
end

for?i=1:n
????eval（[‘result_in‘num2str（i）‘=transpose（result_in‘num2str（i）‘）‘]）
end

for?i=1:n
????eval（[‘result_de‘num2str（i）‘=transpose（result_de‘num2str（i）‘）‘]

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????4087??2010-08-26?09:24??matlab神經網絡30個案例分析\chapter20.m

?????文件??????52475??2010-08-25?21:38??matlab神經網絡30個案例分析\chapter3.rar

?????文件?????375712??2010-08-17?19:56??matlab神經網絡30個案例分析\案例1.rar

?????文件???????4832??2010-08-26?09:24??matlab神經網絡30個案例分析\案例10.rar

?????文件???????3000??2010-08-26?09:24??matlab神經網絡30個案例分析\案例11.rar

?????文件??????38241??2010-08-26?09:35??matlab神經網絡30個案例分析\案例12.rar

?????文件?????284336??2010-08-25?21:51??matlab神經網絡30個案例分析\案例13.rar

?????文件?????178523??2010-08-25?14:19??matlab神經網絡30個案例分析\案例14.rar

?????文件?????356881??2010-08-25?14:12??matlab神經網絡30個案例分析\案例15.rar

?????文件??????46036??2010-08-25?14:14??matlab神經網絡30個案例分析\案例16.rar

?????文件???????3711??2010-08-25?14:16??matlab神經網絡30個案例分析\案例17.rar

?????文件???????2186??2010-08-26?09:43??matlab神經網絡30個案例分析\案例18.rar

?????文件???????2495??2010-08-26?09:34??matlab神經網絡30個案例分析\案例19.rar

?????文件??????48076??2010-08-25?14:04??matlab神經網絡30個案例分析\案例2.rar

?????文件??????91506??2010-08-25?22:13??matlab神經網絡30個案例分析\案例21.rar

?????文件????3761746??2010-08-25?22:00??matlab神經網絡30個案例分析\案例22.rar

?????文件???????4571??2010-08-26?09:35??matlab神經網絡30個案例分析\案例23.rar

?????文件??????13923??2010-08-25?14:23??matlab神經網絡30個案例分析\案例24.rar

?????文件??????25751??2010-08-25?21:54??matlab神經網絡30個案例分析\案例25.rar

?????文件???????1956??2010-08-25?21:40??matlab神經網絡30個案例分析\案例26.rar

?????文件??????89795??2010-08-25?21:48??matlab神經網絡30個案例分析\案例27.rar

?????文件???????2850??2010-08-25?21:53??matlab神經網絡30個案例分析\案例28.rar

?????文件??????88036??2010-08-25?21:43??matlab神經網絡30個案例分析\案例29.rar

?????文件?????100738??2010-08-25?21:46??matlab神經網絡30個案例分析\案例4.rar

?????文件??????60252??2010-08-25?14:09??matlab神經網絡30個案例分析\案例5.rar

?????文件??????12649??2010-08-25?14:07??matlab神經網絡30個案例分析\案例6.rar

?????文件???????3183??2010-08-25?14:21??matlab神經網絡30個案例分析\案例7.rar

?????文件???????6877??2010-08-26?09:26??matlab神經網絡30個案例分析\案例8.rar

?????文件???????5879??2010-08-25?14:10??matlab神經網絡30個案例分析\案例9.rar

?????文件???????6918??2010-08-26?09:46??matlab神經網絡30個案例分析\目錄.txt

............此處省略4個文件信息