example_2.m

一本关于dsp的原著教材

colordef none; % Set a black background for the figures
load data; % Load the experimental data
%whos % View the loaded variables
figure;
plot(X1(:,1),X1(:,2))
logProb(1,1)=logfwd(X1,hmm1);
logProb(1,2)=logfwd(X1,hmm2);
logProb(1,3)=logfwd(X1,hmm3);
logProb(1,4)=logfwd(X1,hmm4);
logProb(1,5)=logfwd(X1,hmm5);
logProb(1,6)=logfwd(X1,hmm6);

figure;
plot(X2(:,1),X2(:,2))
logProb(2,1)=logfwd(X2,hmm1);
logProb(2,2)=logfwd(X2,hmm2);
logProb(2,3)=logfwd(X2,hmm3);
logProb(2,4)=logfwd(X2,hmm4);
logProb(2,5)=logfwd(X2,hmm5);
logProb(2,6)=logfwd(X2,hmm6);

figure;
plot(X3(:,1),X3(:,2))
logProb(3,1)=logfwd(X3,hmm1);
logProb(3,2)=logfwd(X3,hmm2);
logProb(3,3)=logfwd(X3,hmm3);
logProb(3,4)=logfwd(X3,hmm4);
logProb(3,5)=logfwd(X3,hmm5);
logProb(3,6)=logfwd(X3,hmm6);

figure;
plot(X4(:,1),X4(:,2))
logProb(4,1)=logfwd(X4,hmm1);
logProb(4,2)=logfwd(X4,hmm2);
logProb(4,3)=logfwd(X4,hmm3);
logProb(4,4)=logfwd(X4,hmm4);
logProb(4,5)=logfwd(X4,hmm5);
logProb(4,6)=logfwd(X4,hmm6);

figure;
plot(X5(:,1),X5(:,2))
logProb(5,1)=logfwd(X5,hmm1);
logProb(5,2)=logfwd(X5,hmm2);
logProb(5,3)=logfwd(X5,hmm3);
logProb(5,4)=logfwd(X5,hmm4);
logProb(5,5)=logfwd(X5,hmm5);
logProb(5,6)=logfwd(X5,hmm6);

figure;
plot(X6(:,1),X6(:,2))
logProb(6,1)=logfwd(X6,hmm1);
logProb(6,2)=logfwd(X6,hmm2);
logProb(6,3)=logfwd(X6,hmm3);
logProb(6,4)=logfwd(X6,hmm4);
logProb(6,5)=logfwd(X6,hmm5);
logProb(6,6)=logfwd(X6,hmm6);

%for i=1:6,
 %   for j=1:6,
  %      stri=num2str(i);
   %     strj=num2str(j);
    %    eval([ 'logProb(', stri, ',', strj, ')=logfwd(X', stri, ',hmm', strj, ');']);
    %end;
    %end;

logProb
[Lik,k]=max(logProb');
MaxLikelihood=[Lik]'
HMM_No=[k]'