clc;
clear all;
gray=imread('lena.jpg');
temp=double(gray);
[m,n]=size(temp);
T=dctmtx(8);%产生二维DCT变换矩阵
Q_table=[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];
B=blkproc(temp,[8,8],'P1*x*P2',T,T');%blkproc分块处理函数
B=blkproc(B,[8,8],'round(x./P1)',Q_table);%round四舍五入函数，得到DCT系数
C=zeros(1,m*n-m*n/64);%创建一个1行m*n-m*n/64列的矩阵,存储量化后的DCT系数
K=1;
for i=1:m
for j=1:n
if (mod(i,8)~=1)||(mod(j,8)~=1)%去掉DC(直流)值，每一块左上角的第一个数是DC系数。
C(1,K)=B(i,j);%提取AC值到C矩阵中
K=K+1;
end
end
end
%计算AC系数不等于0且不等于1的个数
v=0;
for i=1:m*n-m*n/64%m*n-m*n/64是指用所有像素个数减去DC值，就剩下AC值
if(C(1,i)~=0)&&(C(1,i)~=1)%系数不等于0且不等于1
v=v+1;
end
end
rate=0.8;
A=floor(v*rate);%向下取整，取一部分
M=rand(1,A)>0.5;%随机产生矩阵
%从C中取一不大于A的矩阵
j=0;
for i=1:m*n-m*n/64
if(C(1,i)~=0)&&(C(1,i)~=1)
j=j+1;
end
if j==A
break;
end
end%在长度为A的隐写率下，C应该取到第i位
S=zeros(1,i);%构造C的一个子集，i=j+AC位0或1的个数
%隐写过程
l=1;
for r=1:i
S(1,r)=C(1,r);
end
for r=1:i
if (S(1,r)~=0)&&(S(1,r)~=1)
if (S(1,r)>0)%正数隐写
S(1,r)=S(1,r)-mod(S(1,r),2)+uint8(M(1,l));%隐写过程，M是随机产生的秘密信息
elseif (S(1,r)<0)%负数隐写
if (mod(S(1,r),2)==1)%负奇数项
if (M(1,l)==0)
S(1,r)=S(1,r)-1;
else
S(1,r)=S(1,r);
end
else%负偶数项
if (M(1,l)==1)
S(1,r)=S(1,r)+1;%-4变-3，-2变-1
else
S(1,r)=S(1,r);
end
end
end
l=l+1;
if l>A
break;%隐写结束
end
end
end
for r=1:i
C(1,r)=S(1,r);%隐写之后的C矩阵
end
p = JSteg(C,m,n,rate);
