用matlab代码实现JPEG图像压缩

好的，以下是一个简单的 MATLAB 代码实现 JPEG 图像压缩的过程：

% 读取图像
img = imread('image.jpg');

% 将图像从RGB格式转为YCbCr格式
ycbcr = rgb2ycbcr(img);

% 将YCbCr格式的图像分成3个通道
y = ycbcr(:,:,1);
cb = ycbcr(:,:,2);
cr = ycbcr(:,:,3);

% 对Y通道进行离散余弦变换（DCT）
dct_y = dct2(y);

% 对DCT系数进行量化
quant_y = round(dct_y ./ [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]);

% 对Cb和Cr通道进行4:2:0采样
sub_cb = cb(1:2:end,1:2:end);
sub_cr = cr(1:2:end,1:2:end);

% 对Cb和Cr通道进行离散余弦变换（DCT）
dct_cb = dct2(sub_cb);
dct_cr = dct2(sub_cr);

% 对DCT系数进行量化
quant_cb = round(dct_cb ./ [17 18 24 47 99 99 99 99; 18 21 26 66 99 99 99 99; 24 26 56 99 99 99 99 99; 47 66 99 99 99 99 99 99; 99 99 99 99 99 99 99 99; 99 99 99 99 99 99 99 99; 99 99 99 99 99 99 99 99; 99 99 99 99 99 99 99 99]);
quant_cr = round(dct_cr ./ [17 18 24 47 99 99 99 99; 18 21 26 66 99 99 99 99; 24 26 56 99 99 99 99 99; 47 66 99 99 99 99 99 99; 99 99 99 99 99 99 99 99; 99 99 99 99 99 99 99 99; 99 99 99 99 99 99 99 99; 99 99 99 99 99 99 99 99]);

% 将量化后的DCT系数转为 Zigzag 顺序
zigzag_y = Zigzag(quant_y);
zigzag_cb = Zigzag(quant_cb);
zigzag_cr = Zigzag(quant_cr);

% 将Zigzag顺序的系数进行哈夫曼编码
huff_y = huffman(zigzag_y);
huff_cb = huffman(zigzag_cb);
huff_cr = huffman(zigzag_cr);

% 将哈夫曼编码后的系数写入文件
fid = fopen('compressed.jpg', 'w');
fwrite(fid, huff_y);
fwrite(fid, huff_cb);
fwrite(fid, huff_cr);
fclose(fid);

其中，`Zigzag` 和 `huffman` 都是自定义的函数，用于将系数转为 Zigzag 顺序和进行哈夫曼编码。这里就不一一展示代码了。