基于DCT算法的音频信息隐藏研究.docx

基于DCT算法的音频信息隐藏研究基于DCT算法的音频信息隐藏研究引言随着通信技术的发展和全球信息化进程的推进，信息技术向各个领域不断延伸，使得信息的安全和合法运用成为当前研究的热点和难点，数字版权保护的相关技术日益重要并提上日程。国内外研究机构借助图像实现信息的隐藏研究比较成熟。由于人的听觉要比视觉敏感得多，如何借助音频信息实现信息隐藏的研究是国内外近年来探讨的主题，音频信息隐藏技术根据将隐藏信息嵌入的位置可分为时域与变换域。到目前为止，公认比较成熟的时域音频信息隐藏技术有4种：最不重要位法、回声隐藏法、相位编码法、扩频法。与时域相比，变换域因其较好的鲁棒性而被广泛关注，应用较多的离散傅立叶变换法、离散余弦变换法、离散小波变换法和倒谱法等。本文在介绍基于DCT算法的音频信息隐藏原理的基础上，将欲隐藏信息的音频部分载体的低频部分提取出来并分段处理，然后用隐藏的信息替换掉载体音频低频部分的奇数段中的DCT系数，进而实现音频信息的隐藏，最后通过信息的还原算法实现隐藏信息的还原。1基于DCT算法的音频信息隐藏原理离散余弦变换(DCT)运算简单，有很好的能量压缩性能，音频信号经过DCT变换后只有实部，没有虚部，便于保密信息的嵌入和嵌入强度的控制。在DCT变换域上，信号的能量主要集中在几个变化系数上，仅用少数几个变换系数就可表示信号的总体，这一特点是用DCT变换进行数据压缩的基本依据。长度为N的音频信号f（x）的一维DCT变换和IDCT变换公式如下：正变换：F(0)=1N∑N-1x=0f(x),u=0（1）F(u)=2N∑N-1u=1F(x)cos［π2N(2x+1)u］,u=1,2,…,N-1（2）反变换：f(x)=1NF(0)+2N∑N-1u=1F(u)cos［π2N(2x+1)u］,x=0,1,…,N-1（3）F(0)为称为DC系数，其余N-1个为AC系数，修改DCT系数实现信息的隐藏。2基于DCT算法的音频信息隐藏的MATLAB仿真2.1隐藏信息的预处理设隐藏信息的文件名“hs.bmp”，是256*256的二值图像，若图像的行和列用P和Q表示，则隐藏信息共有P*Q=65536个像素（如图1所示），再把该二值图像进行置乱压缩处理。图像置乱压缩方法有很多,如基于Arnold变换、幻方变换、魔方变换和基于面包师的图像置乱技术等，这里定义一个修改图片函数functionb=modifyImage(hs)，在二值图像信息不变的情况下产生相同大小的图片从而实现置乱。为匹配一维音频信息载体，还需要将置乱压缩后的二值图像的像素点的位置和数据进行改变，二维二值图像降维操作，两次置乱后的图像如图2所示。C=reshape(HSZ,1,P*Q);%两次置乱后的图片HSZ降维成一维序列存入矩阵C中。图1置乱前的隐藏信息图2置乱后的隐藏信息2.2载体音频信息的预处理FDATool工具是MATLAB信号处理工具箱里专用的滤波器设计分析工具，可以通过调整滤波器的类型和参数获得所需的滤波器的幅频相频特性设计。载体音频信息的低通滤波器的设计步骤如下：①在MATLAB的命令窗口中输入：“fdatool”，启动FilterDesign&AnalyzeTool（FDATool）；②FilterDesign&AnalyzeTool（FDATool）的主界面总共分两大部分，一部分是特性区，在界面的上半部分，用来显示滤波器的各种特性；另一部分是DesignFilter，在界面的下半部，用来设置滤波器的设计参数；③在主界面下半部的参数设计界面中根据需要选择合适的选项，本文设计一个低通滤波器，ResponseType选择Lowpass，DesignMethod选择IIR中的Butterworth这个传递函数，设定filterorder项中的specifyorder，specifyorder是在条件的允许下使设计的滤波器尽可能的接近理想滤波器，一般specifyorder越高越好，相对来说消耗的资源多些，指定specifyorder＝5。设置完以后点击DesignFilter即可得到所设计的载体音频信息IIR低通滤波器；④把滤波器导入到空间变量里面，滤波器参数名为Hd1。使用滤波器时，调用空间变量Hd1。载体音频信息通过低通滤波器进行滤波得到滤波后载体音频信息。A=wavread('E:\j.wav');%读取存在E盘的载体音频j.wav到数组A。AL=filter(Hd1,A);%A通过低通滤波器后的参数存到矩阵AL。2.3隐藏信息嵌入载体音频的DCT算法实现（1）载体音频信息分段。首先对滤波后载体音频信息AL分段处理，滤波后载体音频信息AL的长度与原始音频的长度一样，滤波前后不会改变载体音频信号的长度，下面一条语句得到