人脸识别HMM

clc;
A=[1 3 2; 
    4 8 3;
    10 22 3;
    21 34 33;
    9 3 9;
    5 3 9;
    10 3 8;
    ]
H = 7;   %%  图片高度
W = 3;   %%  图片宽度
L = 5;   %%  窗口宽度
P = 3;   %%  连续窗口重合宽度
T = floor((H-L)/(L-P));      %%  观察序列个数 (这里和论文不一致，因为下降几个，有几个观察)
%% 上面的T和下面的top，bottom和论文不太一致，都是因为MATLAB数组下标从1开始
%%%%%%%
%

D=L-P;     %% 每次下滑距离


%% 向量化编程，令T为对应序列号，则top和bottom对应每次观察窗口的上下边缘
data=zeros(t,P)
for t=1:T
    %% 获取窗口上下边缘
    top = 1 + t * D;
    bottom = top + P - 1;
    
    %% 提取窗口，特征值
    [U,S,V]= svd( A(top:bottom,:) );
    
    data(t,:)=diag(S)';
    
end