標(biāo)題: 數(shù)字音頻水印 三級(jí)小波算法及相關(guān)函數(shù)代碼 [打印本頁(yè)]
作者: qizhiweiq22    時(shí)間: 2017-7-13 00:35
標(biāo)題: 數(shù)字音頻水印 三級(jí)小波算法及相關(guān)函數(shù)代碼
數(shù)字音頻水印。三級(jí)小波調(diào)試



源程序如下:
  1. clear all
  2. clc
  3. len=32;ls=60; q=0.00025; N=10;c=[1 1 1 -1 -1 1 -1];    sq=0.05;
  4. thr=6*sq;    %判決門限，這個(gè)數(shù)值可在實(shí)際系統(tǒng)中進(jìn)行調(diào)整%巴克碼
  5. a1=wavread('dsjjj.wav');
  6. figure;
  7. subplot(3,1,1);plot(a1); title('嵌入水印的音頻波形圖');axis([0,1141504,-0.2,0.2]);
  8. lz= length(a1);
  9. lam1=floor((1/N)*lz);
  10. lam=zeros(N-1,1);aa=a1';
  11. out=zeros(1,length(aa));
  12. num=[];   %儲(chǔ)存可能的幀頭

  14. for i=7:N*lam1
  15.     out(i)=sum(aa(i-6:i).*c);    %每7個(gè)數(shù)據(jù)進(jìn)行一次相關(guān)求和
  16.     if(out(i)>thr)           %判斷相關(guān)求和值是否大于門限
  17.         num=[num i-7];
  18.       %  a=a1(num(1):1:1024*N*C);
  19.         %將大于門限的位置進(jìn)行儲(chǔ)存
  20.     end
  21. end
  22. a=a1;

  24. lz= length(a);
  25. lam1=floor((1/N)*lz);
  26. lam=zeros(N-1,1);
  27. lam(1)=lam1;
  28. for i=2:N-1;
  29. lam(i)=i.*lam1;
  30. end
  31. am{1,1}=[a(1:1:lam1)];
  32. for i=2:N-1;
  33. am{i,1}=[a((lam(i-1)+1):1:lam(i))];
  34. end
  35. am{N,1}=[a((lam(N-1)+1):1:lz)];




  40. for j=1:N
  41.     an=am{j,1};an=an';
  42. [c,l]=wavedec(an,3,'haar');%3級(jí)小波分解
  43. ca3=appcoef(c,l,3,'haar');%提取3級(jí)小波分解的最低頻部分
  44. for i=1:len
  45.     ave(i)=sum(ca3(i*ls:(i+1)*ls))/ls;
  46.     wm1(i)=(mod(fix(ave(i)/q+1/2),2));
  47.     wm{j,1}=wm1;
  48. end
  49. end
  50. wm1=wm{1,1};
  51. wm2=wm{2,1};
  52. wm3=wm{3,1};
  53. wm4=wm{4,1};
  54. wm5=wm{5,1};
  55. wm6=wm{6,1};
  56. subplot(3,1,2);plot(wm1);

  58. frr = fopen( 'tiquwave1.txt', 'wb') ;
  59. fwrite( frr, wm1, 'ubit1') ;
  60. fclose( frr);

  62. frr = fopen( 'tiquwave2.txt', 'wb') ;
  63. fwrite( frr, wm2, 'ubit1') ;
  64. fclose( frr);

  66. frr = fopen( 'tiquwave3.txt', 'wb') ;
  67. fwrite( frr, wm3, 'ubit1') ;
  68. fclose( frr);

  70. ……………………

  72. …………限于本文篇幅 余下代碼請(qǐng)從51黑下載附件…………
復(fù)制代碼

所有資料51hei提供下載:
三級(jí)小波算法及相關(guān)函數(shù)代碼.rar (32.7 KB, 下載次數(shù): 11)








歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1