xxx
add DCT and modify it to be the same with MATLAB except the scale