www.pudn.com > yycl.rar > twa.m

%%%获取t_lpc 
function t_lpc=twa(num_group,num_wave) %%time wrapping algorithm 
load('mydata400') 
X=zm_lpc; 
 
num_total=400; 
max_hist=[]; 
for i=1:num_total 
    for j=1:10 
        for k=1:13 
            l=(j-1)*13+k; 
            A(k,j)=X(l,i); %%%%130维分成10帧13维 
        end 
    end 
    n=10; %%%%每个样本开始帧数 
    n_goal=7;  %%%每个样本目标帧数 
 
    while n>n_goal %%%将样本维数压缩到n_goal 
        max=100; 
         for j=1:(n-1) 
            k=j+1; 
                if norm(A(:,j)-A(:,k))max_tag 
               A(:,j)=A(:,j+1); 
            end 
         end 
  end %%%end while 
    b_lpc=[]; 
    for kk=1:n 
        b_lpc=cat(1,b_lpc,A(:,kk)); 
    end    %%%将n_goal帧13维合并称13×n_goal维向量 
    t_lpc(:,i)=b_lpc; 
end %%%end for i=1:total 
save ('twa_lpc.mat','t_lpc','goal_rand');