dct_m.m

离散信号及其Matlab实现,包括单位冲激序列 ,单位阶跃序列,斜坡序列等

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