%讀取訓(xùn)練數(shù)據(jù)
[f1f2f3f4class]?=?textread（‘trainData.txt‘??‘%f%f%f%f%f‘150）;

%特征值歸一化
[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?;

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

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

%測(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?）