x = 1:64;
a = reshape(x,8,8)';
A1 = floor(2*sqrt(2)*dct(a' - 128))';
A = floor(2*sqrt(2)*dct(A1));


Q = [16 11 10 16 24 40 51 61; 
    12 12 14 19 26 58 60 55;     
    14 13 16 24 40 57 69 56;
    14 17 22 29 51 87 80 62;
    18 22 37 56 68 109 103 77; 
    24 35 55 64 81 104 113 92; 
    49 64 78 87 103 121 120 101; 
    72 92 95 98 112 100 103 99];

R = round(2^15 ./ Q);

Y1 = round((A.*R).*2^(-17));

