%--------------------------------------------------------------------------
% exa080901_dct1.m,  for example 8.9.1
% to test dct.m;
% 調用該程序的主程序是exa080901_dct1_test.m,
%--------------------------------------------------------------------------
function x1=exa080801_dct1(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

%subplot(211)
%plot(x)
%subplot(212)
%plot(x1)