clc
clear
close all

%---------------------------------------------------
% 產(chǎn)生訓(xùn)練樣本與測試樣本，每一列為一個樣本

n1 = [rand(3,5),rand(3,5)+1,rand(3,5)+2];
x1 = [repmat([1;0;0],1,5),repmat([0;1;0],1,5),repmat([0;0;1],1,5)];

n2 = [rand(3,5),rand(3,5)+1,rand(3,5)+2];
x2 = [repmat([1;0;0],1,5),repmat([0;1;0],1,5),repmat([0;0;1],1,5)]

xn_train = n1;          % 訓(xùn)練樣本
dn_train = x1;          % 訓(xùn)練目標(biāo)

xn_test = n2;           % 測試樣本
dn_test = x2;           % 測試目標(biāo)

%---------------------------------------------------
% 訓(xùn)練與測試

net = newgrnn(xn_train,dn_train);       % 訓(xùn)練
X = sim(net,xn_test);                   % 測試 - 輸出為預(yù)測值
X = full(compet(X))                     % 競爭輸出

%---------------------------------------------------
% 結(jié)果統(tǒng)計

Result = ~sum(abs(X-x2))               % 正確分類顯示為1
Percent = sum(Result)/length(Result)   % 正確分類率