% fig34 clear;clc;close all; N=512; err1=[]; err2=[]; z=zeros(1,N); z(1,[33:96])=1; z(1,[133:260])=2; z(1,[316:411])=4; subplot(221);plot(z);axis([0 512 -1 5]); title('Original Signal'); for K=1:100, % Fourier compression z_hat=fft(z); zz=fliplr(sort(abs(z_hat))); thres=zz(1,K); for k=1:N, if (abs(z_hat(k))