dct_m.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 22 行

function x1=dct_m(N,k)

% 调出信号 noissin; 
load noissin; 
x=noissin(1:N); 
y=zeros(size(x));

% 对该信号作DCT;
for i=1:N
   if rem(i,8)==0
      y(i-8+1:i)=dct(x(i-8+1:i)); 
      y(i-k+1:i)=0;
   end; 
end

% 重建原信号；
for i=1:N
   if rem(i,8)==0
      x1(i-8+1:i)=idct(y(i-8+1:i)); 
   end; 
end