s089.m

来自「简单的数字图像处理程序」· M 代码 · 共 28 行

 % 矩阵作二维离散余弦变换，对应64个基本函数的频率示意  s089 矩阵作二维离散余弦变换，对应64个基本函数的频率示意  s089
a=zeros(8);
b=zeros(71);
for u=1:8
   for v=1:8
      for x=1:8
         for y=1:8
             a(u,v)=cos((2*(x-1)+1)*(u-1)*pi/16)*cos((2*(y-1)+1)*(v-1)*pi/16);
             b(9*(u-1)+x,9*(v-1)+y)=a(u,v);
         end
      end
  end
end
for u=1:8
  for v=1:8
   for x=1:8
     for y=1:8
     if u==1&&v>1 b(9*(u-1)+x,9*(v-1)+y)=sqrt(1/8)*sqrt(2/8)*b(9*(u-1)+x,9*(v-1)+y);
      elseif u>1&&v==1 b(9*(u-1)+x,9*(v-1)+y)=sqrt(1/8)*sqrt(2/8)*b(9*(u-1)+x,9*(v-1)+y);
      elseif u==1&&v==1 b(9*(u-1)+x,9*(v-1)+y)=sqrt(2)*(1/8)*b(9*(u-1)+x,9*(v-1)+y);
      elseif u>1&&v>1 b(9*(u-1)+x,9*(v-1)+y)=sqrt(2/8)*sqrt(2/8)*b(9*(u-1)+x,9*(v-1)+y);
     end
    end
   end
  end
end
figure,imshow(b,[],'notruesize')