Using these as sample codes

Perform 8x8 DCT transform on them. Then try the following operations:

a. Leave only the DC values. Delete all AC values.

b. Leave the DC value and 5 AC values (in the zig-zag order) in each block.

c. Quantize the DCT coefficients using the standard JPEG quantization table (Table 9.1 and 9.2 in

your reference book).

% Sampling frequency % sampling period % Length of signal % Time vector 1 Fs 1800; 3 L=1000; t (0:L-1 )*T; = cos(2*p1*50*t); % First row wave % Second row wave % Third row wave 6 x1 x2 x3 = cos(2*p1+150*t); cos(2*pi*300*t); = 8 = 12 for i# 1:3 13 14 15 16 end 17 18 n=2^nextpow2 ( L ) ; 19 20 dim-2; 21 22 Yfft(x,n,dim) 23 24 P2abs(Y/L); 25 P1=P2(:, 1:n/2+1); 26 P1(:,2:end-1)2*P1,2:end-1); 27 figure; 28 for i-1:3 29 30 31 32 end subplot (3,1,1) plot ( t (1:100 ) , X(i,1:100)) title(C'Row num2str(i),' in the Time Domain ]) subplot(3,1,i) plot ( 0: ( Fs): (FS/2-Fs), P1 (, 1:n/2)) title(C'Row ,num2str(i),' in the Frequency Domain']) % Sampling frequency % sampling period % Length of signal % Time vector 1 Fs 1800; 3 L=1000; t (0:L-1 )*T; = cos(2*p1*50*t); % First row wave % Second row wave % Third row wave 6 x1 x2 x3 = cos(2*p1+150*t); cos(2*pi*300*t); = 8 = 12 for i# 1:3 13 14 15 16 end 17 18 n=2^nextpow2 ( L ) ; 19 20 dim-2; 21 22 Yfft(x,n,dim) 23 24 P2abs(Y/L); 25 P1=P2(:, 1:n/2+1); 26 P1(:,2:end-1)2*P1,2:end-1); 27 figure; 28 for i-1:3 29 30 31 32 end subplot (3,1,1) plot ( t (1:100 ) , X(i,1:100)) title(C'Row num2str(i),' in the Time Domain ]) subplot(3,1,i) plot ( 0: ( Fs): (FS/2-Fs), P1 (, 1:n/2)) title(C'Row ,num2str(i),' in the Frequency Domain'])