xes1);plot(y);gridon;axistight;title('原始音频信号的时域波形');xlabel('time(s)');ylabel('幅度');axes(handles.axes2);plot(y_fft_f);gridon;axistight;title('原始音频信号的频域波形');xlabel('f(Hz)');ylabel('幅度');sound(y,Fs);%播放原始语音3.2读取原始音频clareset;Fs=44100;FILE1='Rihanna-TakeABow.wav';[y,Fs,bits]=wavread(FILE1);%音频信号y,采样率fs,采样精度bits%用db4小波对原始音频信号进行3级小波分解[c,l]=wavedec(y,3,'db4');%3级小波分解,低频部分为相似,高频部分为细节%提取3级小波分解的低频系数和高频系数ca3=appcoef(c,l,'db4',3);%提取三级小波分解的最低频分cd3=detcoef(c,l,3);%提取三级小波分解的次低频分cd2=detcoef(c,l,2);cd1=detcoef(c,l,1);x=ca3;%提取三级小波分解的最低频部分%找到插入位置,检测特征点s=max(abs(x))*0.2;i=find(abs(x)>s);lx=length(x(i));%读取水印音频FILE2='test_new.wav';[mark,Fs,bits]=wavread(FILE2);mark=mark(1:lx);mark_fft=fft(mark,Fs);mark_fft_f=2*sqrt(mark_fft.*conj(mark_fft));axes(handles.axes1);plot(mark);gridon;axistight;title('水印音频信号的时域波形');
全文预览
