% fig33 clear;clc;close all; N=512; err1=[]; err2=[]; n=[0:511]; z=((n-256).*exp(-(n-256).^2/512)); subplot(221);plot(z);axis([0 512 -11 11]); 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))