关于音频水印攻击的程序

gemer

有哪位兄弟是做音频水印的，请教一下音频水印的攻击的matlab源程序，先谢了

happy

音频没怎么做过，收集了一些程序，看看下面的是否是你需要的

1. %水印嵌入*********************************************************************************************（1）
   
2. %1同步信号，在音频起始点嵌入10个幅值相同的点
   
3. [x,fs,nbits]=wavread('E:\wave1\0000.wav');
   
4. lenx=length(x);
   
5. w=wavread('E:\wave1\0001.wav');
   
6. y=[0.02; 0.02 ;0.02; 0.02 ;0.02; 0.02 ;0.02; 0.02 ;0.02 ;0.02];
   
7. %input m
   
8. m=1000;
   
9. k=fix(lenx/10/m)*10;
   
10. w1=w(1:m);
    
11. %input a
    
12. a=50;
    
13. %2分段处理,
    
14. for i=1:k:m*k
    
15.    %3离散余弦变换
    
16.     x1=dct(x(i:i+k-1));
    
17.    %排序找出中频系数mid 用公式mid'=mid(1+a*w）
    
18.    [x11,index]=sort(x1);
    
19.    j=fix(i/k)+1;   
    
20.    for t=1:3
    
21.       %4修改中频系数
    
22.       x1(index(t+1))=x11(t+1)*(1+a*w(j));
    
23.       %纪录被改变的系数和序号
    
24.       recordx(j,t)=x11(t+1);
    
25.       recordi(j,t)=index(t+1);
    
26.       %5离散余弦逆变换
    
27.       x(i:i+k-1)=idct(x1);
    
28.    end
    
29. end
    
30. %6合并得到嵌入水印后的音频信号
    
31. x=[y;x];
    
32. %将嵌入水印后的音频信号保存
    
33. wavwrite(x,fs,nbits,'E:\wave1\00000.wav');
    
34. %水印提取*********************************************************************************************（2）
    
35. %1同步检测
    
36. [x2,fs,nbits]=wavread('E:\wave1\00000.wav');
    
37. lenxx=length(x2)-m*k;
    
38. i=1;
    
39. %input b
    
40. b=0.001;
    
41. while and(or(or(abs(x2(i)-0.02)>b,abs(x2(i+1)-0.02)>b),abs(x2(i+2)-0.02)>b),i<lenxx)
    
42.     i=i+1;
    
43. end
    
44. while  and(abs(x2(i+3)-0.02)<b, i<=lenxx)
    
45.    i=i+1;
    
46. end
    
47. %判断是否有同步检测信号
    
48. if i>lenxx
    
49.   error('没有同步检测信号 ');
    
50. end
    
51. %2含有水印的音频信号做分段处理
    
52. x3=x2((i+3):(i+2+m*k));
    
53. length(x3)
    
54. %3原始音频信号做分段处理
    
55. for i=1:k:m*k  
    
56. %4含有水印的音频信号做离散余弦变换
    
57.    x4=dct(x3(i:i+k-1));
    
58.    j=fix(i/k)+1;
    
59. %5水印提取
    
60.    w21=(x4(recordi(j,1))/recordx(j,1)-1)/a;
    
61.    w22=(x4(recordi(j,2))/recordx(j,2)-1)/a;
    
62.    w23=(x4(recordi(j,3))/recordx(j,3)-1)/a;
    
63.    w2(j)=(w21+w22+w23)/3;
    
64. end
    
65. %6求相关系数
    
66. cox1=(w1'*w2')/(norm(w1)*norm(w2))
    
67. %水印攻击*********************************************************************************************（3）
    
68. %1加噪音，    转（2）
    
69. %2滤波   ，  转（2）
    
70. %3剪切，     转（2）
    
71. %4重量化，   转（2）
    
72. %5重采样，   转（2）
    
73. %6MPEG压缩，转（2）

复制代码

gemer

非常感谢，我也是正在做论文，挺着急的，我先看看，再次感谢！！！

liujunjun25

学习一下:lol

liujunjun25

%水印攻击*********************************************************************************************（3）
%1加噪音，    转（2）
%2滤波   ，  转（2）
%3剪切，     转（2）
%4重量化，   转（2）
%5重采样，   转（2）
%6MPEG压缩，转（2）
怎么没攻击的程序呢:@(

lijianzhong85

有没有基于快速傅里叶变换的音频水印源代码 谢了