clc;
clear;
load（‘irisdata.mat‘）;
%讀取訓(xùn)練數(shù)據(jù)
%[f1f2f3f4class]?=?textread（‘trainData.txt‘??‘%f%f%f%f%f‘150）;
f1=[setosa（:1）;versicolour（:1）;virginica（:1）];
f2=[setosa（:2）;versicolour（:2）;virginica（:2）];
f3=[setosa（:3）;versicolour（:3）;virginica（:3）];
f4=[setosa（:4）;versicolour（:4）;virginica（:4）];
class=[ones（501）;ones（501）*2;ones（501）*3];

%特征值歸一化
[inputminImaxI]?=?premnmx（?[f1??f2??f3??f4?]‘）??;

%構(gòu)造輸出矩陣
s?=?length（?class）?;
output?=?zeros（?s??3??）?;
for?i?=?1?:?s?
???output（?i??class（?i?）??）?=?1?;
end

%創(chuàng)建神經(jīng)網(wǎng)絡(luò)
net?=?newff（?minmax（input）??[10?3]?）;?%?{?‘logsig‘?‘purelin‘?}??‘traingdx‘?）?;?

%設(shè)置訓(xùn)練參數(shù)
net.trainparam.show?=?50?;
net.trainparam.epochs?=?500?;
net.trainparam.goal?=?0.01?;
net.trainParam.lr?=?0.01?;

%開始訓(xùn)練
net?=?train（?net?input??output‘?）?;

%讀取測(cè)試數(shù)據(jù)
%?[t1?t2?t3?t4?c]?=?textread（‘testData.txt‘??‘%f%f%f%f%f‘150）;
t1=f1;
t2=f2;
t3=f3;
t4=f4;
c=class;

%測(cè)試數(shù)據(jù)歸一化
testInput?=?tramnmx?（?[t1t2t3t4]‘??minI?maxI?）?;

%仿真
Y?=?sim（?net??testInput?）?;

%統(tǒng)計(jì)識(shí)別正確率
[s1??s2]?=?size（?Y?）?;
hitNum?=?0?;
for?i?=?1?:?s2
????[m??Index]?=?max（?Y（?:???i?）?）?;
????if（?Index??==?c（i）???）?
????????hitNum?=?hitNum?+?1?;?
????end
end
sprintf（‘識(shí)別率是?%3.3f%%‘100?*?hitNum?/?s2?）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1332??2015-01-13?17:14??BP_demo.m
?????文件????????2210??2015-01-13?16:54??irisdata.mat