1.一种基于共生矩阵分析的MP3音频隐写检测方法，其特征在于包括以下步骤：

①选取N个变化风格各不相同且未压缩的WAV音频，其中，N≥100；

②利用MP3音频8HZ编码器对每个未压缩的WAV音频进行压缩编码，得到每个未压缩的WAV音频对应的未隐写MP3压缩音频，由得到的N个未隐写MP3压缩音频构成音频负样本库；

并利用MP3Stego音频隐写算法对每个未压缩的WAV音频隐写不同长度、不同内容的隐秘信息，得到每个未压缩的WAV音频对应的隐写MP3压缩音频，由得到的N个隐写MP3压缩音频构成第一类音频正样本库；利用基于窗口类型选择的隐写算法对每个未压缩的WAV音频隐写不同长度、不同内容的隐秘信息，得到每个未压缩的WAV音频对应的隐写MP3压缩音频，由得到的N个隐写MP3压缩音频构成第二类音频正样本库；利用基于Hufffman码表索引选择的隐写算法对每个未压缩的WAV音频隐写不同长度、不同内容的隐秘信息，得到每个未压缩的WAV音频对应的隐写MP3压缩音频，由得到的N个隐写MP3压缩音频构成第三类音频正样本库；

再由音频负样本库、第一类音频正样本库、第二类音频正样本库和第三类音频正样本库构成的一个样本库；

③利用MP3音频lame解码器对样本库中的每个样本进行解压缩，得到样本库中的每个样本对应的WAV音频；然后利用MP3音频lame编码器对样本库中的每个样本对应的WAV音频进行压缩编码，得到样本库中的每个样本重压缩后对应的载体估计；

④利用MP3音频lame解码器对样本库中的每个样本进行解压缩，提取出样本库中的每个样本中的每帧的576个量化后的MDCT系数，以每个样本中的每帧的576个量化后的MDCT系数作为一行，将样本库中的每个样本对应的所有量化后的MDCT系数构成一个系数矩阵，将样本库中的第i个样本对应的所有量化后的MDCT系数构成的系数矩阵记为Xi，其中，1≤i≤4N，Xi的维数为 表示样本库中的第i个样本中包含的帧的总数，x1,1、x1,2、x1,575、x1,576对应表示样本库中的第i个样本中的第1帧的第1个、第2个、第575个、第576个量化后的MDCT系数，x2,1、x2,2、x2,575、x2,576对应表示样本库中的第i个样本中的第2帧的第1个、第2个、第575个、第576个量化后的MDCT系数， 对应表示样本库中的第i个样本中的第帧的第1个、第2个、第575个、第576个量化后的MDCT系数；

⑤对样本库中的每个样本对应的所有量化后的MDCT系数构成的系数矩阵中值大于

300的系数进行修正，得到样本库中的每个样本对应的新系数矩阵，将样本库中的第i个样本对应的新系数矩阵记为Xi'；

然后根据样本库中的每个样本对应的新系数矩阵，构造样本库中的每个样本对应的新系数矩阵的水平方向共生矩阵、垂直方向共生矩阵、45度角方向共生矩阵和135度角方向共生矩阵，将Xi'的水平方向共生矩阵、垂直方向共生矩阵、45度角方向共生矩阵和135度角方向共生矩阵对应记为Pi,0、Pi,90、Pi,45和Pi,135，将Pi,0中下标为(p,q)处的元素记为Pi,0(p,q)， 将Pi,90中下标为(p,q)处的元素记为Pi,90(p,q)， 将Pi,45中下标

为(p,q)处的元素记为Pi,45(p,q)，

将Pi,135中下标为 (p,q) 处的元素记为 P i,135(p,q)，

其中，1≤p≤300,1≤q≤300，

d表示共生矩阵的步长，x'u,v表示Xi'中下标为(u,v)处的系数，x'u,v+d表示Xi'中下标为(u,v+d)处的系数，x'u+d,v表示Xi'中下标为(u+d,v)处的系数，x'u+d,v+d表示Xi'中下标为(u+d,v+d)处的系数；

⑥根据样本库中的每个样本对应的新系数矩阵的水平方向共生矩阵，构造样本库中的每个样本对应的新系数矩阵的水平方向共生矩阵的包含有12个高阶统计特征的特征向量，将Pi,0的包含有12个高阶统计特征的特征向量记为Fi,0；

根据样本库中的每个样本对应的新系数矩阵的垂直方向共生矩阵，构造样本库中的每个样本对应的新系数矩阵的垂直方向共生矩阵的包含有12个高阶统计特征的特征向量，将Pi,90的包含有12个高阶统计特征的特征向量记为Fi,90；

根据样本库中的每个样本对应的新系数矩阵的45度角方向共生矩阵，构造样本库中的每个样本对应的新系数矩阵的45度角方向共生矩阵的包含有12个高阶统计特征的特征向量，将Pi,45的包含有12个高阶统计特征的特征向量记为Fi,45；

根据样本库中的每个样本对应的新系数矩阵的135度角方向共生矩阵，构造样本库中的每个样本对应的新系数矩阵的135度角方向共生矩阵的包含有12个高阶统计特征的特征向量，将Pi,135的包含有12个高阶统计特征的特征向量记为Fi,135；

⑦根据样本库中的每个样本对应的新系数矩阵的水平方向共生矩阵、垂直方向共生矩阵、45度角方向共生矩阵、135度角方向共生矩阵各自的特征向量，获取样本库中的每个样本的包含有48个特征的特征向量，将样本库中的第i个样本的包含有48个特征的特征向量记为Fi，Fi为Fi,0、Fi,90、Fi,45、Fi,135四个特征向量共包含的48个高阶统计特征按序构成；

⑧按照步骤④至步骤⑦的操作，以相同的方式获取样本库中的每个样本重压缩后对应的载体估计的包含有48个特征的特征向量，将样本库中的第i个样本重压缩后对应的载体*估计的包含有48个特征的特征向量记为Fi；

⑨计算样本库中的每个样本的特征向量与样本库中的每个样本重压缩后对应的载体估计的特征向量中对应的两个特征的差值的绝对值，得到样本库中的每个样本对应的包含有48个元素的隐写分析特征向量，将样本库中的第i个样本对应的包含有48个元素的隐写分析特征向量记为Yi，将Yi中的第t个元素记为Yi(t)， 其中，* *

1≤t≤48，Fi(t)表示Fi中的第t个特征，Fi(t)表示Fi中的第t个特征，符号“| |”为取绝对值符号；

⑩对样本库中的每个样本对应的隐写分析特征向量进行归一化处理，得到样本库中的每个样本对应的归一化后的隐写分析特征向量；

将样本库中属于音频负样本库的所有样本各自对应的归一化后的隐写分析特征向

量标记为-1，将样本库中属于第一类音频正样本库的所有样本各自对应的归一化后的隐写分析特征向量标记为+1，然后将标记后的属于音频负样本库的所有样本各自对应的归一化后的隐写分析特征向量和属于第一类音频正样本库的所有样本各自对应的归一化后的隐写分析特征向量输入SVM分类器进行训练，得到MP3Stego音频隐写算法的检测模板，记为M1；

同样，将样本库中属于音频负样本库的所有样本各自对应的归一化后的隐写分析特征向量标记为-1，将样本库中属于第二类音频正样本库的所有样本各自对应的归一化后的隐写分析特征向量标记为+1，然后将标记后的属于音频负样本库的所有样本各自对应的归一化后的隐写分析特征向量和属于第二类音频正样本库的所有样本各自对应的归一化后的隐写分析特征向量输入SVM分类器进行训练，得到基于窗口类型选择的隐写算法的检测模板，记为M2；

将样本库中属于音频负样本库的所有样本各自对应的归一化后的隐写分析特征向量标记为-1，将样本库中属于第三类音频正样本库的所有样本各自对应的归一化后的隐写分析特征向量标记为+1，然后将标记后的属于音频负样本库的所有样本各自对应的归一化后的隐写分析特征向量和属于第三类音频正样本库的所有样本各自对应的归一化后的隐写分析特征向量输入SVM分类器进行训练，得到基于Hufffman码表索引选择的隐写算法的检测模板，记为M3；

对于任意一个待检测的MP3压缩音频，按照步骤③至步骤⑩的过程，以相同的方式获取该MP3压缩音频对应的归一化后的隐写分析特征向量；然后分别利用MP3Stego音频隐写算法的检测模板M1、基于窗口类型选择的隐写算法的检测模板M2和基于Hufffman码表索引选择的隐写算法的检测模板M3对该MP3压缩音频对应的归一化后的隐写分析特征向量进行检测，以确定该MP3压缩音频是否经过隐写及隐写所采用的音频隐写算法。

2.根据权利要求1所述的一种基于共生矩阵分析的MP3音频隐写检测方法，其特

征在于所述的步骤⑤中样本库中的第i个样本对应的新系数矩阵Xi'的获取过程为：将Xi'中下标为(u,v)处的系数记为x'u,v， 其中，

1≤v≤576，xu,v表示Xi中下标为(u,v)处的系数。

3.根据权利要求1或2所述的一种基于共生矩阵分析的MP3音频隐写检测方法，其

特征在于所述的步骤⑥中Fi,0的获取过程为：Fi,0＝[Fi,0(1),Fi,0(2),Fi,0(3),Fi,0(4),Fi,0(5),Fi,0(6),Fi,0(7),Fi,0(8),Fi,0(9),Fi,0(10),Fi,0(11),Fi,0(12)]，其中，Fi,0(1),Fi,0(2),…,Fi,0(12)对应表示Fi,0中的第1个至第12个高阶统计特征，Fi,0(2)＝max(Pi,0)，

max()表示取最大值函数，

Pi,0(p,p)表示Pi,0中下标为(p,p)处的元素，Pi,0(p,300-p)表示Pi,0中下标为(p,300-p)处的元素；

按照Fi,0的获取过程，以相同的方式分别获取Fi,90、Fi,45和Fi,135。