gmmevaltest.m

一个关于数据聚类和模式识别的程序,在生物化学,化学中因该都可以用到.希望对大家有用,谢谢支持

addpath d:/users/jang/matlab/toolbox/dcpr

fprintf('Traing GMM...\n');
DS=prData('iris');
gmmParam=gmmTrain(DS.input, [8, 1]);
dataNum=100000;
dim=length(gmmParam(1).mu);
data=100*rand(dim, dataNum);

%fprintf('Compiling gmmEvalMex.cpp...\n');
%mex gmmEvalMex.cpp d:/users/jang/c/lib/dcpr.cpp d:/users/jang/c/lib/utility.cpp

fprintf('Testing gmmEvalMex.cpp...\n');
tic, [logProb1, eachProb1]=gmmEval(data, gmmParam); time1=toc;
tic, [logProb2, eachProb2]=gmmEvalMexOld(data, [gmmParam.mu], [gmmParam.sigma], [gmmParam.w]); time2=toc;
tic, [logProb3, eachProb3]=gmmEvalMex(data, gmmParam); time2=toc;

fprintf('Difference in logProb1 and logProb2 = %g\n', sum(abs(logProb1-logProb2)));
fprintf('Difference in eachProb1 and eachProb2 = %g\n', sum(sum(abs(eachProb1-eachProb2))));
fprintf('time1/time2=%g/%g=%g\n', time1, time2, time1/time2);
fprintf('Difference in logProb1 and logProb3 = %g\n', sum(abs(logProb1-logProb3)));
fprintf('Difference in eachProb1 and eachProb3 = %g\n', sum(sum(abs(eachProb1-eachProb3))));