faceprob.m

It is for Face Recognition

%function [p]=faceProb(T1, T2, T3, mT1Face, mT2Face, mT3Face, x, styl)
function [p]=faceProb(T1, T2, T3, mT1Face, mT2Face, mT3Face, x, styl)
 % const
 stdT1Face = [258.6135   99.7751   92.8625   79.4257   60.8324   54.7351   53.1098   46.3724   45.6027   41.6553   40.5613   38.7417]';
 stdT2Face = [191.6423   92.2207   75.4775   58.2254   49.4999   41.9451   39.8454   37.1355]';
 stdT3Face = [187.1485   67.5272   54.5260   49.3976   46.3355   40.8735   34.2275   30.2965]';

 y1 = T1*x          - mT1Face;
 y2 = T2*x(1:96)    - mT2Face;
 y3 = T3*x(97:192)  - mT3Face;

 % face: pca: 5-9
 if y1(5) < 121 & y1(9) < 95
    p1 = 1;
 else
    p1 = 0;
 end
  
 % upper face: pca 4-7: 
 if y2(4) > -100 & y2(7) < 85
    p2 = 1;
 else
    p2 = 0;
 end
 
 % lower face: pca 1:6
 if y3(1)^2 + y3(6)^2 - 185^2 < 0 
    p3 = 1;
 else
    p3 = 0;
 end
 
 p = p1*p2*p3;

 d1 = y1./stdT1Face;
 d2 = y2./stdT2Face;
 d3 = y3./stdT3Face;

 d = [d1' d2' d3'];
 styl = 'xx';
 if styl ~= 'xx'
     hold on;
     stem(d, styl);
 end