新版Matlab 7辅助信号处理技术与应用

				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