源代码-MIDI信息隐藏技术研究.docVIP

水印嵌入 若（s[i]0xf0）=0x90，则在第i+2个字节s[i+2]中嵌入1比特水印信息，表示修改声音开启的最低位。若（s[i]0xf0）=0xc0，则在第i+1个字节s[i+1]中嵌入1比特水印信息，表示修改乐器编号的最低位比特。若在第i+1个字节s[i+1]中嵌入4比特水印信息，表示修改通道触动压力的低4比特位。 源代码midihiding.m如下： clc; clear; wm=randsrc(1,3000,[0 1]); %产生长度为3000的随机水印信息 save(‘watermark.mat’,’wm’); wm=uint8(wm); fid=fopen(‘love.mid’,’r’); %读入音频文件 wa=fread(fid,inf); [length,col]=size(wa); fclose(fid); modcount=1; watermarklen=3000; j=1; for i=1:length if(bitand(uint8(wa(i,1)),uint8(240))= =144)(j=watermarklen-1) wa((i+2),1)=bitand(wa(i+2),254); wa((i+2),1)=bitset(wa((i+2),1),1,wm(1,j)); modcount=modcount+1; i=i+3; j=j+1; elseif (bitand(uint8(wa(i,1)),uint8(240))= =176)(j=watermarklen-1) wa((i+1),1)=bitand(wa(i+1),254); wa((i+1),1)=bitset(wa((i+1),1),1,wm(1,j)); modcount=modcount+1; i=i+2; j=j+1; elseif (bitand(uint8(wa(i,1)),uint8(240))= =208)(j=watermarklen-4) wa((i+1),1)=bitand(wa(i+1),254); wa((i+1),1)=bitset(wa((i+1),1),1,wm(1,j)); wa((i+1),1)=bitset(wa((i+1),1),2,wm(1,j+1)); wa((i+1),1)=bitset(wa((i+1),1),3,wm(1,j+2)); wa((i+1),1)=bitset(wa((i+1),1),4,wm(1,j+3)); i=i+2; j=j+4; modcount=modcount+1; else i=i+1; end end fid=fopen(‘watermarked.mid’,’wb’); fwrite(fid,wa); fclose(fid); 2、水印提取 源代码midiextract.m如下： clc; clear; fid=fopen(‘watermarked.mid’,’r’); %读入携密音频文件 wa=fread(fid,inf); [length,col]=size(wa); fclose(fid); watermarklen=3000; for j=1: watermarklen wm1(j)=0; end j=1; for i=1:length if(bitand(uint8(wa(i,1)),uint8(240))= =144)(j=watermarklen) wm1(1,j)=bitget(wa(i+2),1); i=i+2; j=j+1; elseif (bitand(uint8(wa(i,1)),uint8(240))= =176)(j=watermarklen-1) wm1(1,j)=bitget(wa(i+1),1); i=i+1; j=j+1; elseif (bitand(uint8(wa(i,1)),uint8(240))= =208)(j=watermarklen-4) wm1(1,j)=bitget(wa(i+2),1); wm1(1,j+1)=bitget(wa(i+2),2); wm1(j+2)=bitget(wa(i+2),3); wm1(j+3)=bitget(wa(i+2),4); i=i+2; j=j+4; end end load(‘watermark.mat’); diff=wm-wm1; 通信0801班 ***