iydct.m

自己手写的图像的FFT变换和DCT变换的MATLAB代码

%************************************
%*****    快速余弦反变换
%************************************
function f=iyDCT(x,N)
f=zeros(1,N);
f1=zeros(1,2*N);
c=(-1)^(0.5);
 for i=1:N
     f1(i)=x(i)*exp(c*(i-1)*pi/2/N);
 end;
f1=Y1iFFT(f1,N*2);
for i=2:N+1
f(i-1)=(1-sqrt(2))*sqrt(1/N)*x(1)+sqrt(2/N)*real(f1(i))*N*2;
end
end