%%?第26章?有導(dǎo)師學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)的分類(lèi)——鳶尾花種類(lèi)識(shí)別
%?
%?			該案例作者申明：		1：本人長(zhǎng)期駐扎在此板塊里，對(duì)該案例提問(wèn)，做到有問(wèn)必答。	2：此案例有配套的教學(xué)視頻，視頻下載請(qǐng)點(diǎn)擊http://www.matlabsky.com/forum-91-1.html。?								3：此案例為原創(chuàng)案例，轉(zhuǎn)載請(qǐng)注明出處（《MATLAB智能算法30個(gè)案例分析》）。							4：若此案例碰巧與您的研究有關(guān)聯(lián)，我們歡迎您提意見(jiàn)，要求等，我們考慮后可以加在案例里。						5：以下內(nèi)容為初稿，與實(shí)際發(fā)行的書(shū)籍內(nèi)容略有出入，請(qǐng)以書(shū)籍中的內(nèi)容為準(zhǔn)。		
%?

%%?清空環(huán)境變量
clear?all
clc

%%?訓(xùn)練集/測(cè)試集產(chǎn)生

%?導(dǎo)入數(shù)據(jù)
load?iris_data.mat
%?隨機(jī)產(chǎn)生訓(xùn)練集和測(cè)試集
P_train?=?[];
T_train?=?[];
P_test?=?[];
T_test?=?[];
for?i?=?1:3
????temp_input?=?features（（i-1）*50+1:i*50:）;
????temp_output?=?classes（（i-1）*50+1:i*50:）;
????n?=?randperm（50）;
????%?訓(xùn)練集——120個(gè)樣本
????P_train?=?[P_train?temp_input（n（1:40）:）‘];
????T_train?=?[T_train?temp_output（n（1:40）:）‘];
????%?測(cè)試集——30個(gè)樣本
????P_test?=?[P_test?temp_input（n（41:50）:）‘];
????T_test?=?[T_test?temp_output（n（41:50）:）‘];
end

%%?模型建立?
result_grnn?=?[];
result_pnn?=?[];
time_grnn?=?[];
time_pnn?=?[];
for?i?=?1:4
????for?j?=?i:4
????????p_train?=?P_train（i:j:）;
????????p_test?=?P_test（i:j:）;
???????%%?GRNN創(chuàng)建及仿真測(cè)試
????????t?=?cputime;
????????%?創(chuàng)建網(wǎng)絡(luò)
????????net_grnn?=?newgrnn（p_trainT_train）;
????????%?仿真測(cè)試
????????t_sim_grnn?=?sim（net_grnnp_test）;
????????T_sim_grnn?=?round（t_sim_grnn）;
????????t?=?cputime?-?t;
????????time_grnn?=?[time_grnn?t];
????????result_grnn?=?[result_grnn?T_sim_grnn‘];
???????%%?PNN創(chuàng)建及仿真測(cè)試
????????t?=?cputime;
????????Tc_train?=?ind2vec（T_train）;
????????%?創(chuàng)建網(wǎng)絡(luò)
????????net_pnn?=?newpnn（p_trainTc_train）;
????????%?仿真測(cè)試
????????Tc_test?=?ind2vec（T_test）;
????????t_sim_pnn?=?sim（net_pnnp_test）;
????????T_sim_pnn?=?vec2ind（t_sim_pnn）;
????????t?=?cputime?-?t;
????????time_pnn?=?[time_pnn?t];
????????result_pnn?=?[result_pnn?T_sim_pnn‘];
????end
end

%%?性能評(píng)價(jià)

%?正確率accuracy
accuracy_grnn?=?[];
accuracy_pnn?=?[];
time?=?[];
for?i?=?1:10
????accuracy_1?=?length（find（result_grnn（:i）?==?T_test‘））/length（T_test）;
????accuracy_2?=?length（find（result_pnn（:i）?==?T_test‘））/length（T_test）;
????accuracy_grnn?=?[accuracy_grnn?accuracy_1];
????accuracy_pnn?=?[accuracy_pnn?accuracy_2];
end
%?結(jié)果對(duì)比
result?=?[T_test‘?result_grnn?result_pnn]
accuracy?=?[accuracy_grnn;accuracy_pnn]
time?=?[time_grnn;time_pnn]
%%?繪圖
figure（1）
plot（1:30T_test‘bo‘1:30result_grnn（:4）‘r-*‘1:30result_pnn（:4）‘k:^‘）
grid?on
xlabel（‘測(cè)試集樣本編號(hào)‘）
ylabel（‘測(cè)試集樣本類(lèi)別‘）
string?=?{‘測(cè)試集預(yù)測(cè)結(jié)果對(duì)比（GRNN?v

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????1105??2010-10-17?14:51??chapter26\iris_data.mat
?????文件????????4233??2010-11-30?20:45??chapter26\main.m
?????目錄???????????0??2018-01-02?21:38??chapter26\