y1dct.m

来自「自己手写的图像的FFT变换和DCT变换的MATLAB代码」· M 代码 · 共 16 行

%******Y1DCT() 一维的快速正余弦变换
%***********************************
function f=Y1DCT(x,N)
f=zeros(1,N);
f1=zeros(1,2*N);
 for i=1:N
     f1(i)=x(i);
 end;
c=(-1)^(0.5);
f1=Y1FFT(f1,N*2);
k=(2/N)^0.5;
f(1)=(1/N)^(0.5)*f1(1);
for i=2:N
f(i)=k*(real(f1(i))*cos((i-1)*pi/(N*2))+imag(f1(i))*sin((i-1)*pi/(N*2)));
end
end