matlab下信息隐藏算法

clc; clear; %读入音频文件 fid=fopen('1.wav','rb'); oa=fread(fid,inf,'uchar'); n=length(oa)-44;%wav文件从文件头到数据前面共有44字节的文件格式说明，读取的时候跳过这44个字节的文件格式说明 fclose(fid); %读取二值图像水印文件 io=imread('bupt.bmp'); [row col]=size(io); if row*col>n error('载体太小，请更换载体'); end wi=io(:); [Y0,I0]=sort(wi); %对原始音频进行DCT变换 dctv1=dct(oa); [m,n]=size(dctv1); A=dctv1(:); [Y1,I1]=sort(A); x=m*n; LENGTH=row*col; k=LENGTH; M=zeros(x,1); alfa=.1; for i=1:x if k>=1 M(x)=Y1(x)*(1+alfa*Y0(k));%嵌入水印后的音频 k=k-1; else M(x)=Y1(x); end x=x-1; end N=zeros(x,1); x=m*n; for i=1:x N(I1(i))=M(i); end a=1; for j=1:n, for i=1:m, dctv2(i,j)=N(a); a=a+1; end end %DCT反变换，得到嵌入水印的音频 idctv1=idct(dctv2); figure; subplot(3,1,1);plot('oa');title('original audio'); subplot(3,1,2);plot(io);title('original image'); subplot(3,1,3);plot(idctv1,[]);title('qianruhou audio ');