音频数字水印matlab仿真实现

clc,clear all;
close all;
%%%%读取音频文件
[waveinfo,Fs,NBIT]=wavread('D:\water\doo.wav');
[rowave,colwave]=size(waveinfo);
lenwave=rowave*colwave;
waveinfo=reshape(waveinfo,1,lenwave);
% for k=1:lenwave
%    waveinfo((k-1)*8+1:k*8)=binwave((k-1)*8+1:k*8);
% end 
%%%%%%%%%%%%读取图像文件
imageinfo=imread('D:\water\soga.bmp');
[rowimage,colimage]=size(imageinfo);
lenimage=rowimage*colimage;
imageinfo=reshape(imageinfo,1,lenimage)
imageinfo=double(imageinfo);%%%%%%%%%%%%%%%%%%%进行数字化处理的时候，比如加减乘除，一定要变为double类型
imageinfo1=zeros(1,8*lenimage);
for k=1:lenimage
    for j=8:-1:1
    imageinfo1((k-1)*8+j)=mod(imageinfo(k),2);
    imageinfo(k)=fix(imageinfo(k)/2);%应用此公式时，数据类型不能为unit
    end
end 
%imageinfo=double(imageinfo1);
%%%%%%%%%%%%%%%%%%%%%%%%
imageinfo=zeros(1,8*lenimage);
imageinfo=imageinfo1;
Lstep=3;
[c,l]=wavedec(waveinfo,3,'db4');%对waveinfo进行‘db4’型3级小波分析。
ca3=appcoef(c,l,'db4',3);