dan@0: function [M,Mdct]=vec2featmat(D,cep_ord)
dan@0:   
dan@0: b0=[0 0 0 0 1  0  0  0 0];
dan@0: b1=[0 0 2 1 0 -1 -2 0 0]/10;
dan@0: b2=conv(b1(3:7),[2 1 0 -1 -2]/10);
dan@0: 
dan@0: Mdct=[cos(pi*(1:cep_ord)'*((1:D)-0.5)/D);ones(1,D)];
dan@0: 
dan@0: M=kron([b0;-b1;b2],Mdct);