LSB的音频水印嵌入提取算法【Matlab语音处理】.zip

%2010-09-24.隐藏信息为一个文本文件 clear all clc close all [y,fs,nbits,opts]=wavread('C:\Users\baishufeng\Documents\MATLAB\LSB\country1.wav');% 语音载体，16kHz，16bits，单声道y就是音频信号； % fs是采样频率，比如说16000就是每秒16000次； % nbit是采样精度，比如说16就是指16位精度的采样； [m,n]=size(y) wavplay(y,fs);%播放音频文件 id=fopen('test.txt','r'); [msg,msg_len]=fread(id,'ubit1') len=msg_len+2;%为什么？ L=floor(m/len)%计算分帧数，每个帧的长度为len sum=zeros(L*len,1); e=zeros(L,1); for i=1:len:(L-1)*len+1%？？？？？？？ for j=0:len-1 sum(i)=sum(i)+abs(y(i+j,1)); end end for i=1:L e(i)=sum((i-1)*len+1);%计算各个帧的能量 end [y1,inx]=sort(e,'descend'); %对各帧的能量进行从大到小排序 %[sA,index] = sort(A) ，排序后，sA是排序好的向量，index 是 向量sA 中对 A 的索引。 索引使排列逆运算成为可能。 b=(inx(1)-1)*len+1;%能量最大的帧对应的起始位置，为什么？？ c=inx(1)*len;%能量最大的帧对应的结束位置 y2=y(b:c); s=3;%窗口大小，不明白？？窗口有什么用 for i=s:len y2_binary=convertbinary_16(y2(i-s+1:i),s);%十进制转16位二进制 y2(i)=y2(i)+msg(i-s+1)+mod((y2_binary(1,16)+y2_binary(2,16)+y2_binary(3,16)),2); end wavwrite(y2,fs,nbits,'C:\Users\baishufeng\Documents\MATLAB\LSB\country2.wav');%保存文件 [y11,fs,nbits,opts]=wavread('C:\Users\baishufeng\Documents\MATLAB\LSB\country2.wav'); wavplay(y11,fs); figure subplot(2,2,1);%原始音频信号图 plot(y); axis([1,m ,min(y),max(y)]); title('原始音频信号图'); subplot(2,2,2); plot(y11); axis([1,m ,min(y11),max(y11)]);%嵌入水印后的音频信号图 title('嵌入水印后的原始音频信号图');