function?[F]=DCT（A）
???A=im2double（A）;
???[mn]=size（A）;
???F=zeros（mn）;
???S=zeros（1m*n）;
???k=2;
???for?u=1:m
???????for?v=1:n
???????????for?a=1:m
???????????????for?b=1:n
???????????????????S（11）=A（11）*cos（（（2*1+1）*u*pi）/（2*n））*cos（（（2*1+1）*v*pi）/（2*m））;
???????????????????if?k==2
???????????????????????S（1k）=A（ab+1）*cos（（（2*a+1