给我一份用DCT压缩jpg图像的代码

这是一个用Matlab实现的DCT压缩JPEG图像的代码：

clear all; close all; clc;

% 读入原始图像
I = imread('lena.jpg');
imshow(I);

% 归一化（0-1范围内）
I = im2double(I);
imwrite(I, 'lena_normal.jpg');

% 色彩空间转换（转换为YUV）
YUV = rgb2ycbcr(I);
imwrite(YUV, 'lena_yuv.jpg');

% 获取YUV中的Y通道
Y = YUV(:, :, 1);

% 设置DCT变换的块大小
N = 8;

% 进行DCT变换
dct_fun = @(block_struct) dct2(block_struct.data);
B = blockproc(Y, [N N], dct_fun);

% 量化
q_mtx = [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];
dq_mtx = bsxfun(@rdivide, q_mtx, N * N);
quant_fun = @(block_struct) round(block_struct.data ./ dq_mtx) .* q_mtx;
QB = blockproc(B, [N N], quant_fun);
imwrite(QB, 'lena_quant.jpg');

% 逆量化
idq_fun = @(block_struct) block_struct.data .* dq_mtx;
IDQB = blockproc(QB, [N N], idq_fun);

% 进行IDCT变换
idct_fun = @(block_struct) idct2(block_struct.data);
IDY = blockproc(IDQB, [N N], idct_fun);

% 恢复图像
IY = Y - IDY;
IYUV = YUV;
IYUV(:, :, 1) = IY;
IRGB = ycbcr2rgb(IYUV);

imshow(IRGB);

注意：此代码仅供学习参考，实际使用时需考虑更多细节和性能优化。