提升小波的語(yǔ)音端點(diǎn)檢測(cè)算法研究畢業(yè)論文

析為理論基,對(duì)語(yǔ)音端點(diǎn)檢測(cè)中小波系數(shù)方差算法和子噪聲的頻域差,對(duì)這兩種算法進(jìn)行了優(yōu),并應(yīng)差算法耗時(shí)實(shí)時(shí)性差的缺,并克服了子帶平性求:1引言2語(yǔ)音端點(diǎn)的一般原理3語(yǔ)音端點(diǎn)本章提出一種以圖案作水印、且提取時(shí)不需原圖的盲語(yǔ)音端點(diǎn)嵌入算法。這個(gè)算法是基于變換的一種算法，經(jīng)過(guò)研究知道：在小波變換后，紋理和邊緣特征一般集中在高頻子帶LH1、HL1和HH1中，如把水印信息加在低頻系數(shù)LL上，則人眼不易察覺(jué)，可以加大水印的嵌入量。由于經(jīng)過(guò)一般的圖像處理以后，高頻上的信息容易丟失，其魯棒性不夠強(qiáng)。但是，一般情況下在嵌入的水印信息的量是較少的，本文為了增強(qiáng)水印的魯棒性，在嵌入的水印信息量較少的情況下，就把水印信息嵌入到低頻系數(shù)上，這樣水印的魯棒性將會(huì)大大增強(qiáng)。所謂水印攻擊分析，就是對(duì)現(xiàn)有的語(yǔ)音端點(diǎn)系統(tǒng)進(jìn)行攻擊，以檢驗(yàn)其魯棒性，通過(guò)分析其弱點(diǎn)所在與其易受攻擊的原因，以便在以后語(yǔ)音端點(diǎn)系統(tǒng)的設(shè)計(jì)中加以改進(jìn)。攻擊的目的在于使相應(yīng)的語(yǔ)音端點(diǎn)系統(tǒng)的檢測(cè)工具無(wú)確地恢復(fù)水印信號(hào)，或不能檢測(cè)到水印信號(hào)的存在。這和傳統(tǒng)密碼學(xué)中的加密算法設(shè)計(jì)和密碼分析是相對(duì)應(yīng)的。下面我們對(duì)一些典型的攻擊方法進(jìn)行分析[1-9,26-27]。所謂串謀攻擊就是利用同一原始多媒體數(shù)據(jù)集合的不同水印信號(hào)版本，來(lái)生成一個(gè)近似的多媒體數(shù)據(jù)集合，以此來(lái)逼近恢和復(fù)原始數(shù)據(jù)，其目的是使檢測(cè)系統(tǒng)無(wú)法在這一近似的數(shù)據(jù)集合中檢測(cè)出水印信號(hào)的存在，其最簡(jiǎn)單的一種實(shí)現(xiàn)就是平均法。3.3.1小波變換的語(yǔ)音端點(diǎn)算法3.3.2實(shí)驗(yàn)結(jié)果與分析3.4.1DCT域語(yǔ)音端點(diǎn)3.4.2本章算法與實(shí)驗(yàn)結(jié)果分析4小波變換在數(shù)字音頻水印中的應(yīng)用5視頻水印概述與算法分析6程序附錄clearall;closeall;clc;[filename1,pathname]=uigetfile('*.*','selecttheimage');image1=imread(num2str(filename1));subplot(2,2,1);imshow(image1);title('originalimage');%orginalimageforwatermarkingimage1=double(image1);[row,col]=size(image1);imagew=imread('dmg2.tif');subplot(2,2,2);imshow(imagew);title('watermark');%??ó???è?oˉêyμ÷ó?[marked]=blockdwt2(image1,imagew);%generatesthewatermarkedimage%markedmax=max(marked(:));%markscale=marked/markedmax*255;subplot(2,2,3);colormap(gray(256));image(marked);title('Watermarkedimage');%showsthewatermarkedimageimwrite(marked,gray(256),'marked_image.bmp');%savesthewatermarkedimageasabmpfile%*******??ó?ìáè?1y3ì********watermark=marked-image1;%imageadaptivewatermark%watermark=watermark*255/max(watermark(:))fori=1:row%thresholdingforj=1:colifwatermark(i,j)>0watermark(i,j)=255;endifwatermark(i,j)<0watermark(i,j)=0;endendend[w2h2]=size(watermark);[w3h3]=size(imagew);p=1;q=1;mark=zeros(size(imagew));fori=1:w3:w2forj=1:h3:h2mark=mark+watermark(i:p*w3,j:q*h3);q=q+1;endp=p+1;q=1;endmark=mark/8;subplot(2,2,4);colormap(gray(256));image(mark);title('watermark');%showstheimageadaptivewatermarkimwrite(mark,gray(256),'watermark.bmp');%savestheimageadaptivewatermarkasabmpfile%%%%%%%%%%%%D??ü?à??%%%%%%%%%%%%%image1=imread(num2str(filename1));x0=image1;imaged=imread('marked_image.bmp');x2=imaged;%%%%%%%%%%%%%%%%%%%%%%??????è?á???ó?μ?D?o?μ?D???±è%%%%fz=sum(x0.*x0);fm=sum((x0-x2).*(x0-x2));SNR=-10*log(fm/fz)%%%%%%%%%%%%%%%%%%%%%?????ó???ê%%%%%%%%%%%%%%%%%%%%%%mark=imread('dmg2.tif');mark1=imread('watermark.bmp');[row,col]=size(mark);m1=mark(:);mm1=mark1(:);err=sum(sum(xor(m1',mm1')))/row*colfunction[D]=blockdwt2(A,W);[row,col]=size(A);k=0.007;[ca,ch,cv,cd]=dwt2(A,'db1');c1=[chcvcd];[h,w]=size(ca');[m,n]=size(c1');W=dmg(W,A);[caachhcvvcdd]=dwt2(W,'db1');W=caa;size(W);%Addingwatermarkimage.fori=1:hforj=1:wifW(i,j)==0a=-0.1;elsea=0.1;endCa(i,j)=ca(i,j)*(1+k*a);%<--------k*abs(double(c1(i,j)))*W(i,j);deolabilirendendD=double(idwt2(Ca,ch,cv,cd,'db1'));%Damgauretme%%imgew:Damgaimgesi%imge1:Orjinalimge%functionW=dmg(imgew,imge1);[wh]=size(imgew');[w1h1]=size(imge1');a=1;b=1;W=zeros(size(imge1'));fori=1:w:w1forj=1:h:h1W(i:a*w,j:b*h)=imgew';b=b+1;enda=a+1;b=1;endW=W';基于小波變換語(yǔ)音端點(diǎn)（2）clearall;closeall;clc;M=256;%?-í???3¤?èN=64;%??ó?3¤?è[filename1,pathname]=uigetfile('*.*','selecttheimage');image1=imread(num2str(filename1));subplot(2,2,1);imshow(image1);title('originalimage');%orginalimageforwatermarkingimage1=double(image1);imagew=imread('dmg2.bmp');subplot(2,2,2);imshow(imagew);title('originalwatermark');%originalwatermark%??è???ó?[ca,ch,cv,cd]=dwt2(image1,'db1');[cas,chs,cvs,cds]=dwt2(ca,'db1');fori=1:Nforj=1:Nifimagew(i,j)==0a=-1;elsea=1;endCa(i,j)=cas(i,j)*(1+a*0.01);endendIM=idwt2(Ca,chs,cvs,cds,'db1');markedimage=double(idwt2(IM,ch,cv,cd,'db1'));%??ê???è?oó??ó?í???subplot(2,2,3);colormap(gray(256));image(markedimage);title('markedimage');imwrite(markedimage,gray(256),'watermarked.bmp','bmp');%ìáè???ó?image1=imread(num2str(filename1));image1=double(image1);imaged=imread('watermarked.bmp');[ca,ch,cv,cd]=dwt2(image1,'db1');[cas,chs,cvs,cds]=dwt2(ca,'db1');[caa,chh,cvv,cdd]=dwt2(imaged,'db1');[caas,chhs,cvvs,cdds]=dwt2(caa,'db1');forp=1:Nforq=1:Na=caas(p,q)-cas(p,q);ifa<0W(p,q)=0;elseW(p,q)=255;endendend%??ê?ìáè?μ???ó?subplot(2,2,4);colormap(gray(256));image(W);title('′óo???ó?í????Dìáè?μ???ó?');imwrite(W,gray(256),'watermark.bmp','bmp');%%%%%%%%%%%%D??ü?à??image1=imread(num2str(filename1));x0=image1;imaged=imread('watermarked.bmp');x2=imaged;%%%%%%%%%%%%%%%%%%%%%%??????è?á???ó?μ?D?o?μ?D???±è%%%%fz=sum(x0.*x0);fm=sum((x0-x2).*(x0-x2));SNR=-10*log(fm/fz)%%%%%%%%%%%%%%%%%%%%%?????ó???ê%%%%%%%%%%%%%%%%%%%%%%mark=imread('dmg2.bmp');mark1=imread('watermark.bmp');[row,col]=size(mark);m1=mark(:);mm1=mark1(:);err=sum(sum(xor(m1',mm1')))/row*col%??è???ó?μ?3ìDò′ú??M=256;%?-í???3¤?èN=32;%??ó?í???3¤?èK=8;I=zeros(M,M);J=zeros(N,N);BLOCK=zeros(K,K);%??ê??-í???subplot(3,2,1);I=imread('mona.bmp','bmp');imshow(I);title('?-ê?1??aí???');%??ê???ó?í???subplot(3,2,2);J=imread('flag.bmp','bmp');imshow(J);title('??ó?í???');%??è???ó?forp=1:Nforq=1:Nx=(p-1)*K+1;y=(q-1)*K+1;BLOCK=I(x:x+K-1,y:y+K-1);BLOCK=dct2(BLOCK);ifJ(p,q)==0a=-1;elsea=1;endBLOCK(1,1)=BLOCK(1,1)*(1+a*0.03);BLOCK=idct2(BLOCK);I(x:x+K-1,y:y+K-1)=BLOCK;endend%??ê???è???ó?oóμ?í???subplot(3,2,3);imshow(I);title('??è???ó?oóμ?í???');imwrite(I,'watermarked.bmp','bmp');%′ó??è???ó?μ?í????Dìáè???ó?I=imread('mona','bmp');J=imread('watermarked.bmp','bmp');forp=1:Nforq=1:Nx=(p-1)*K+1;y=(q-1)*K+1;BLOCK1=I(x:x+K-1,y:y+K-1);BLOCK2=J(x:x+K-1,y:y+K-1);BLOCK1=idct2(BLOCK1);BLOCK2=idct2(BLOCK2);a=BLOCK2(1,1)/BLOCK1(1,1)-1;ifa<0W(p,q)=0;elseW(p,q)=1;endendend%??ê?ìáè?μ???ó?subplot(3,2,4);imshow(W);title('′óo???ó?í????Dìáè?μ???ó?');%??ó?1￥?÷2aê?M=256;N=32;K=8;I=zeros(M,M);J=zeros(M,M);w=zeros(N,N);BLOCK1=zeros(K,K);BLOCK2=zeros(K,K);%????è???ó?oóμ?í?????DDJPEG?1??L=imread('watermarked.bmp','bmp');imwrite(L,'attack.jpg','jpeg','Quality',45);J=imread('attack.jpg','jpeg');subplot(3,2,5);imshow(J);title('?1??oóμ?í???');I=imread('mona','bmp');%′ó?1??μ?í????Dìáè???ó?forp=1:Nforq=1:Nx=(p-1)*K+1;y=(q-1)*K+1;BLOCK1=I(x:x+K-1,y:y+K-1);BLOCK2=J(x:x+K-1,y:y+K-1);BLOCK1=idct2(BLOCK1);BLOCK2=idct2(BLOCK2);a=BLOCK2(1,1)/BLOCK1(1,1)-1;ifa<0W(p,q)=0;elseW(p,q)=1;endendend%??ê?ìáè?μ???ó?subplot(3,2,6);imshow(W);title('′ó?-1y?1??μ?í????Dìáè?μ???ó?');clearall;closeall;clc;[x,fs,nbits]=wavread('F:\?Dò?×÷òμ?à1?\ò??μ\2blues1.wav');[c,l]=wavedec(x,3,'db4');%èy??D?2¨·??aca3=appcoef(c,l,'db4',3);cd3=detcoef(c,l,3);cd2=detcoef(c,l,2);cd1=detcoef(c,l,1);x=ca3;%ìáè?μí?μ?μêylenx=length(x);%×÷?aí?2?D?o?%y=[0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02];%?áè???ó?D??￠mark=imread('F:\?Dò?×÷òμ?à1?\ò??μ\name.bmp');[row,col]=size(mark);w=mark(:);m=length(w');k=fix(lenx/10/m)*10;subplot(2,1,1);imshow(mark);title('??ó?í???','Fontsize',16,'color','blue');w1=w(1:m);a=0.01;fori=1:k:m*kx1=x(i:i+k-1);j=fix(i/k)+1;fort=1:3recordx(j,t)=x1(t);x1(t)=x1(t)*(1+a*w(j));x(i:i+k-1)=x1;endendc1=[x',cd3',cd2',cd1'];s1=waverec(c1,l,'db4');%6o?2￠μ?μ???è???ó?oóμ?ò??μD?o?%marked=[y;s1];wavwrite(s1,fs,nbits,'F:\?Dò?×÷òμ?à1?\ò??μ\em2blues1.wav');%??ó?ìáè?[x2,fs,nbits]=wavread('F:\?Dò?×÷òμ?à1?\ò??μ\em2blues1.wav');[cr,lr]=wavedec(x2,3,'db4');car3=appcoef(cr,lr,'db4',3);cdr3=detcoef(cr,lr,3);cdr2=detcoef(cr,