一种结合密钥和随机单位正交基的音频伪装及恢复方法与流程

本发明属于信息安全和数字音频信号处理交叉领域，涉及一种音频伪装恢复方法，特别涉及一种结合密钥和随机单位正交基的音频伪装及恢复方法。

背景技术

近年来，伴随着计算机和网络技术的发展，越来越多的数字图像音频在公有信道中传播。对秘密图像音频不加限制地访问和使用，不仅涉及个人隐私、媒体公信和政府诚信，同时也会导致社会动荡和诱发军事冲突。

如何有效地保护公有信道中传输的秘密图像音频，成为目前研究的热点。针对此问题，人们已提出了多种图像音频保护方法，例如将明文图像音频转换为密文的数字图像音频加密技术以及数字图像音频伪装技术等。相对于数字图像音频加密技术，数字图像音频伪装是将机密图像音频伪装成有意义的非机密图像音频，从而在公有信道传输中不易引起拦截者的注意，减少潜在攻击的可能性，进一步保证传输中秘密图像音频的安全性。

tangram方法，也称七巧法或中国拼图方法(丁玮.数字图像信息安全的方法研究[d].北京：中国科学院计算技术研究所，2000.)。其基本思路是选取一幅和秘密图像等大的公开图像，将秘密图像和公开图像分割成不重叠的等大子块。将公开图像子块进行8种等距变换构造字典块，利用最小2乘法寻找秘密图像的每个子块在公开图像字典块中匹配残差最小子块所对应的变换参数，然后利用这些变换参数将公开图像转换为秘密图像。由于tangram方法在匹配过程中需进行全局搜索，编码开销大，严重限制了tangram方法的实际应用价值，同时tangram方法变换过程中仅有8个等距变换，也制约了将秘密图像转变为公开图像的匹配精度，而借助全局搜索找到最合适的等距变换参数也进一步提高了计算复杂性。

为降低经典tangram方法的搜索时间和加快编码速度，吴军和吴秋新等对经典tangram方法进行了改进，添加了小块均化操作，在此基础上给出了基于块均化5元组匹配的tangram方法(吴军，吴秋新.一种基于七巧板游戏的数字图像信息伪装方法[j].计算机应用，2004，24(6):125-128.)。但该方法在减小计算代价的同时，也减少了等距变换数量，从而进一步降低了匹配精度。

为进一步减少搜索时间，加快编码速度，余建德等人对tangram方法进行了改进，利用像素灰度值作图像区域非均匀剖分的思想，提出了基于三角形剖分的tangram方法(余建德，宋瑞霞，齐东旭.基于数字图像三角形剖分的信息伪装方法[j].计算机研究与发展，2009，46(9):1432-1437.)。相对于经典tangram方法，基于三角形剖分的tangram方法不进行全局搜索，只找到对应三角区域的差值，降低了运算代价。但只能对秘密图像的三角形剖分区域进行近似重构，降低了秘密图像重构精度。

为减少计算代价和提高秘密图像的重构质量，在所授权的国家发明专利“基于改进tangram算法和2维双尺度矩形映射的数字图像伪装方法”(邵利平，李苑梦.基于tangram算法和2维双尺度矩形映射的图像伪装及重构方法[p].中国专利，zl201410404838.7，2014-11-12)中，将秘密图像划分的小块作为字典，利用2d双尺度矩形映射确定秘密图像子块和公开图像子块的对应关系，将每个秘密图像子块按8个等距变换直接和对应位置的公开图像子块进行最小2乘法匹配。由于避免了全局匹配，实际编码时间远低于tangram方法。但所提策略同传统tangram方法一样，通过有限的等距变换，由此导致信道中的伪装图像视觉质量不高。

tangram方法在后续的研究中，也被进一步拓展到音频伪装。在所申请的国家发明专利“基于分块序列的数字图像伪装及重构方法”(邵利平，李苑梦，谢贤文.基于分块序列的数字图像伪装及重构方法[p].中国专利:cn104835107a，2015-05-12.)和所授权的国家发明专利“基于分段序列的数字音频伪装及重构方法”(邵利平，谢贤文，李苑梦.基于分段序列的数字音频伪装及重构方法[p].中国专利:zl201510239139.6，2015-05-12.)中，我们分别以图像划分的不重叠图像块和数字音频划分的小段序列构造旋转向量，通过提高等距变换的数量来增加匹配精度，使得公开图像和公开音频伪装质量较好，且重构密图和重构音频的视觉和听觉质量也较为清晰。

但单纯地通过增加等距变换数量来提高可用于最小2乘匹配的向量数量并从中找到最优的变换参数所带来的性能提升依然十分有限。等距变换仅能提供有限数量的元素组合而不是全部组合，因此所找到的最优变换参数并不是所有向量元素组合的最优解，同时由于引入了等距变换，也使得通过向量旋转所产生的多个等距变换向量，都要和公开音频进行最小2乘匹配，从而带来了较高的匹配计算代价，为了减少等距变换数量，所发表的论文(谢贤文,邵利平.结合字典学和排序线性拟合的音频信息隐藏方法.小型微型计算机系统[j],2017,38(12):2658-2667)进一步引入了排序线性拟合,所提方法通过引入了排序拟合使得排序后的公开音频元素和秘密音频划分的小段序列变化趋势保持一致，不仅避免了多个等距变换所带来的高昂匹配代价，也提高了匹配拟合精度。

但无论是①标准tangram算法、还是②块均化5元组匹配tangram方法以及之前给出的③改进tangram算法和2维双尺度矩形映射的数字音频伪装方法、④基于分块序列的数字音频伪装及重构方法、⑤基于分段序列的数字音频伪装及重构方法与⑥结合字典学和排序线性拟合的音频信息隐藏方法，这些方法所基于的变换模型都是仿射变换模型，对于图像，仿射变换模型只有均值块和差异块，对应音频，仿射变换模型只有均值向量和差异向量，且均值块和差异块以及均值向量和差异向量不满足基本的正交关系，导致仿射变换模型的变换精度普遍较低，不能有效的保证秘密图像音频对公开图像音频的拟合精度，从而不能进行有效地信道欺骗，同时也不能有效地保证公开图像音频对秘密图像音频的拟合精度，从而无法保证秘密图像音频的准确重构和获取高精度的秘密图像音频。而另一方面，若仿射模型的差异块和均值块或均值向量和差异向量趋于一致，则对应为恒值块或恒值序列处理，此时则必须添加随机扰动来改善秘密图像音频和公开图像音频划分小块和小段的匹配性能，否则将无法恢复，由此进一步限制了放射变换模型的变换精度。对于⑦基于三角剖分的tangram算法则只能对秘密图像的三角剖分区域进行近似重构，其实际应用价值依然较小。

技术实现要素：

本发明的目的在于克服现有技术缺陷，给出一种结合密钥和随机单位正交基的音频伪装恢复方法，同现有方法相比，本发明可实现不同精度的重构并严格依赖于用户密钥，从而具备较高的安全性。

为达到上述目的，本发明采用以下技术方案予以实现：

一种结合密钥和随机单位正交基的音频伪装方法，包括以下步骤：

第1步：记秘密音频和公开音频分别为长度为l·n和l·n的r位音频序列s＝(si)l·n和p＝(pi)l·n，即si,pi∈{-2^r-1,…,0,…,2^r-1-1}，将s和p分别划分为长度为n和n大小的小段序列，记为其中u＝0,1,…,l-1，由密钥key生成随机序列g＝(gu)l，初始化u＝0；

其中，n,n为正整数，si,pi分别对应为秘密音频和公开音频序列中的第i个元素，r为正整数，对应为秘密音频和公开音频的采样位数，和分别是su和pu中的第i个元素；

第2步：将gu作为密钥，生成n×n的随机矩阵x′，将x′的第0行所有元素置为1，然后对x′进行行单位正交化后得到x；

第3步：记x对应的n个行依次为xf,f＝0,1,…,n-1，对于所有的计算su在单位正交基xf上的投影其中，f为正整数，对应为矩阵x的第f行；

第4步：从中找到除的前k-1个最大幅值系数α0,α1,…,αk-2，通过长度为k-1的索引序列r＝(ri)k-1,ri∈{0,1,…,n-1}来记录α0,…,αk-2在上对应的行索引，然后将记为αk-1；

第5步：将r＝(ri)k-1转化为10进制数记为rn，并进一步通过q进制数进行表达，其中lr为对应的q进制数序列长度；

第6步：将αi,i＝0,1,…,k-1进行q进制数表示；

第7步：记公开音频序列然后将rn和q进制数表示的αi,i＝0,1,…,k-1嵌入到中，从而将pu转换为pu′＝(pi′^u)n，置u＝u+1；

第8步：反复执行第2步到第7步，直至u＝l，然后将所有的pu′,u＝0,1,…,l-1重新拼接为公开音频p，将p作为信道传输音频并输出。

本发明进一步的改进在于：

第2步中，对x′进行行单位正交化的具体方法是按式(1)对x′进行行单位施密特正交化，即将x′的每一行转换为单位正交向量：

x＝schimidt(x′)(1)

式(1)中，函数schimidt()为行单位施密特正交化函数。

第3步中，计算su在向量xf上的投影的具体方法是式(2)：

式(2)中，“||||2”为2范数，其中系数能够按式(3)进行计算：

第5步lr为对应的q进制数序列长度如式(4)所示：

lr＝[logqn^k-1](4)。

第6步中，将αi,i＝0,1,…,k-1进行q进制数表示的具体方法是：

第6.1步：将αk-1按式(5)转换为10进制数lm是其对应的q进制数序列长度，由式(6)确定：

式(5)中，“[]”为四舍五入取整函数

lm＝[logq2^r](6)

第6.2步：将αi,i＝0,1,…,k-1借助进行近似表示，如式(7)所示，然后转换为q进制数表示：

式(7)中，对应为αi的符号部分，分别对应为αi的幂次和有效数字。

第6.2步中，具体的确定方法如式(8)、式(9)和式(10)所示：

式(8)、式(9)和式(10)中，用0或1进行表达，其中0对应为正数，1对应为负数，用位q进制数进行存储，用lint位q进制数进行存储；

第6.2步将转换为q进制数的具体方法是将直接视为2进制数，将其对应的10进制数记为sign2，sign2对应的q进制数长度lsign由式(11)确定：

和对应的q进制数长度分别为lpow,lint，其满足的约束如式(12)所示：

第7步中，将r＝(ri)k-1和q进制数表示的αi,i＝0,1,…,k-1嵌入到的具体方法是式(13)，其中r和对应的10进制数分别为rn和sign2；

式(13)中，函数emd-q()为emd-q全方位扩展嵌入函数，记10进制数t∈[0,1,…,q^d-1]和转换前、后的序列分别为h＝(hi)d和h′＝(h′i)d，则h′＝emd-q(t,h)执行的具体功能是：

第7.1步：按式(14)计算出映射值vmap：

第7.2步：按式(15)计算t与vmap模q^d之差vdiff：

vdiff＝(t-vmap)modq^d(15)

第7.3步：将vdiff按式(16)映射到区间

第7.4步：按式(17)计算hi从而将t嵌入到h中，其中i＝0,1,…,d-1：

式(19)中，参数λi按式(18)计算：

一种结合密钥和随机单位正交基的秘密音频恢复方法，包括以下步骤：

第1步：输入接收到的信道传输音频p′＝(pi′)l·n，将p′划分成长度为n的小段序列p′＝(pu)l，记为其中u＝0,1,…,l-1，由密钥key生成随机序列g＝(gu)l，初始化u＝0；

第2步：将gu作为密钥，依次生成n×n的随机矩阵x′，将x′的第0行置为1，对x′进行行单位正交化后记为x；

第3步：记x的n个行依次为xf,f＝0,1,…,n-1，且

第4步：从pu中提取出隐藏变换参数rn和q进制数表示的αi,i＝0,1,…,k-1，并进一步将其转换为n进制数序列r＝(ri)k-1和10进制数αi,i＝0,1,…,k-1，然后利用这些恢复出的变换参数来重构音频小段序列su，置u＝u+1

第5步：重复第2到4步，直至u＝l，然后将所有的小块su,u＝0,1,…,l-1重新拼接作为解密后的秘密音频s输出。

上述方法进一步的改进在于；

第2步中，对x′分别进行行正交化的具体方法是按式(1)分别对x′进行行单位施密特正交化：

x＝schimidt(x′)(1)

式(1)中，函数schimidt()为单位施密特正交化函数。

第4步中，从pu中提取出的隐藏变换参数提取出隐藏变换参数rm和q进制数表示的αi,i＝0,1,…,k-1，并进一步将其转换为10进制数序列r＝(ri)k-1和10进制数αi,i＝0,1,…,k-1的具体方法是：

第4.1步：按式(19)从中提取出的隐藏变换参数αk-1″、rn和sign2：

第4.2步：对于参数αi,i＝0,1,…,k-2则由式(20)中提取出

第4.3步：将sign2对应的2进制数序列和按式(9)计算α0,α1,…,αk-2：

第4.4步：将rn转化为n进制数序列r＝(ri)k-1；

式(19)和式(20)中，extract()为提取函数，t＝extract(h)表示从序列h＝(hi)d中提取出嵌入的t，具体的提取步骤是：

第4.1.1步：按式(14)计算出映射值vmap：

第4.1.2步：按式(21)提取出秘密值t

第4步中，重构秘密音频小段序列su的具体方法是：

第4.5步：利用序列r＝(ri)k-1从矩阵x中选取对应的行得到k-1个1维序列，记为yi,i＝0,1,…,k-2；

第4.6步：将矩阵x中x0记为yk-1，然后按式(22)计算αk-1：

第4.7步：按式(23)将yi,i＝0,1,…,k-1分别与αi,i＝0,1,…,k-1相乘并累计作为恢复出的小段序列su：

与现有技术相比，本发明具有以下有益效果：

①传统tangram方法通常采用的是仿射变换模型，对于图像，仿射变换模型通常包括均值块和差异块，对于音频，仿射变换模型通常包括均值向量和差异向量，其中均值块和均值向量是由单位元素构成的矩阵小块和矩阵向量所构成，而差异块和差异向量则来源于秘密图像音频或公开图像音频划分的矩阵小块和向量小段。但这些均值块和差异块以及均值向量和差异向量不满足基本的正交关系，存在交叠和冗余，导致仿射线性模型不能对被拟合图像或音频进行有效地表达和拟合，从而不能有效的保证秘密图像音频对公开图像音频的拟合精度，同时也不能有效地保证公开图像音频对秘密图像音频的拟合精度。同传统tangram方法采用的仿射变换模型不同，本发明利用密钥来构造随机单位正交基，从而保证所构造的向量严格正交，避免拟合对象在不同正交向量上的投影系数相互交叠，从而可有效地消除冗余，对待拟合对应进行充分有效地表达，提高拟合精度。

②传统tangram方法通常采用仿射变换模型，对于图像只包含均值块和差异块，对于音频只包含均值向量和差异向量，由此导致仿射变换模型对拟合对象的拟合能力十分有限，因此同传统tangram方法采用的仿射变换模型不同，本发明所构造的音频伪装方法则包含多个正交基，始终通过选取幅度值最大系数对应的正交向量来对秘密音频进行重构，并可通过选取不同数量的正交基来实现秘密音频的不同精度的重构。

③传统tangram方法通常是由秘密图像音频划分的小块和小段来对公开音频划分的小块和小段进行仿射线性拟合，以进行秘密图像音频的信道传输欺骗，或通过公开图像音频划分的小块和小段来对秘密图像音频划分的小块和小段进行重构，但无论是秘密图像音频对公开图像音频的仿射线性拟合还是公开图像音频对秘密图像音频的仿射线性拟合，若差异块和均值块区域一致，则使得恒值块和恒值序列无法对图像音频划分的小块和小段进行正确拟合，同时也不能保证图像音频划分的小块和小段对恒值小块和小段拟合过程中的信息丢失，为避免这点，传统方法通常需增加随机扰动来避免趋于一致，由此会在拟合过程中添加噪声，导致恢复质量下降。同以上方法不同，本发明引入了随机密钥，通过随机密钥生成的正交基向量来对秘密音频进行拟合，并且只有正确的随机密钥才可保证随机正交基可以准确重构，从而不会造成任何信息丢失，同时不同的正交张量基也避免了传统仿射线性拟合模型中均值块和差异块可能导致的一致性。并且所述策略严格依赖于密钥，只有正确的密钥才能对密图进行重构，从而使得该策略具备较高的安全性。

④本发明可直接将变换参数嵌入到公开音频中，为进一步减少嵌入变换参数对公开载体的影响，本发明对嵌入的参数进行了优化，使得能在更大范围内对嵌入的参数进行表示，且对秘密音频的影响较小，本发明并不追求秘密音频和公开音频大小的严格相等，从而可通过选取合适大小的信息隐藏载体来对秘密音频进行高精度的重构。

附图说明

图1是嵌入流程图；

图2是提取流程图；

图3是实施例，选取长度为1128768个字节的16位音频采样序列：陌上花早，作为秘密音频序列；

图4是实施例，选取长度为4797264个字节的16位音频采样序列：倾尽天下，作为公开音频序列；

图5是实施例，取k＝2，按图1嵌入流程嵌入参数后的公开音频1，为长度为4797264个字节的16位音频，相对于图4的snr为67.0138db；

图6是实施例，取k＝3，按图1嵌入流程嵌入参数后的公开音频2，是长度为4797264个字节的16位音频，相对于图4的snr为65.0563db；

图7是实施例，取k＝2，按图2提取流程重构的秘密音频1，是长度为1128768个字节的16位音频，相对于图3的snr为19.9059db；

图8是实施例，取k＝3，按图2提取流程重构的秘密音频2，是长度为1128768个字节的16位音频，相对于图3的snr为21.4306db。

具体实施方式

下面结合附图对本发明做进一步详细描述：

参见图1-2，本发明结合密钥和随机单位正交基的音频伪装方法，包括以下步骤：

第1步：记秘密音频和公开音频分别为长度为l·n和l·n的r位音频序列s＝(si)l·n和p＝(pi)l·n，即将s和p分别划分为长度为n和n大小的小段序列，记为和其中u＝0,1,…,l-1，由密钥key生成随机序列g＝(gu)l，初始化u＝0；

其中，n,n为正整数，si,pi分别对应为秘密音频和公开音频序列中的第i个元素，r为正整数，对应为秘密音频和公开音频的采样位数，和分别是su和pu中的第i个元素；

第2步：将gu作为密钥，生成n×n的随机矩阵x′，将x′的第0行所有元素置为1，然后对x′进行行单位正交化后得到x；其具体方法是按式(1)对x′进行行单位施密特正交化，即将x′的每一行转换为单位正交向量：

x＝schimidt(x′)(1)

式(1)中，函数schimidt()为行单位施密特正交化函数。

第3步：记x对应的n个行依次为xf,f＝0,1,…,n-1，对于所有的计算su在单位正交基xf上的投影其中，f为正整数，对应为矩阵x的第f行；具体的：

计算su在向量xf上的投影的具体方法是式(2)：

式(2)中，“||||2”为2范数，其中系数可按式(3)进行计算：

第4步：从中找到除的前k-1个最大幅值系数α0,α1,…,αk-2，通过长度为k-1的索引序列r＝(ri)k-1,ri∈{0,1,…,n-1}来记录α0,…,αk-2在上对应的行索引，然后将记为αk-1；

第5步：将r＝(ri)k-1转化为10进制数记为rn，并进一步通过q进制数进行表达，其中lr为对应的q进制数序列长度如式(4)所示：

第6步：将αi,i＝0,1,…,k-1进行q进制数表示，其具体方法是：

第6.1步：将αk-1按式(5)转换为10进制数αk-1″，lm是其对应的q进制数序列长度，由式(6)确定：

式(5)中，“[]”为四舍五入取整函数

第6.2步：将αi,i＝0,1,…,k-1借助进行近似表示，如式(7)所示，然后转换为q进制数表示：

式(7)中，对应为αi的符号部分，分别对应为αi的幂次和有效数字。

具体的确定方法如式(8)、式(9)和式(10)所示：

式(8)、式(9)和式(10)中，用0或1进行表达，其中0对应为正数，1对应为负数，用位q进制数进行存储，用lint位q进制数进行存储；

将转换为q进制数的具体方法是将直接视为2进制数，将其对应的10进制数记为sign2，sign2对应的q进制数长度lsign由式(11)确定：

和对应的q进制数长度分别为lpow,lint，其满足的约束如式(12)所示：

第7步：记公开音频序列然后将rn和q进制数表示的αi,i＝0,1,…,k-1嵌入到中，从而将pu转换为pu′＝(pi′^u)n，置u＝u+1；

将r＝(ri)k-1和q进制数表示的αi,i＝0,1,…,k-1嵌入到的具体方法是式(13)，其中r和对应的10进制数分别为rn和sign2；

式(13)中，函数emd-q()为emd-q全方位扩展嵌入函数，记10进制数t∈[0,1,…,q^d-1]和转换前、后的序列分别为h＝(hi)d和h′＝(hi′)d，则h′＝emd-q(t,h)执行的具体功能是：

第7.1步：按式(14)计算出映射值vmap：

第7.2步：按式(15)计算t与vmap模q^d之差vdiff：

vdiff＝(t-vmap)modq^d(15)

第7.3步：将vdiff按式(16)映射到区间

第7.4步：按式(17)计算hi从而将t嵌入到h中，其中i＝0,1,…,d-1：

式(19)中，参数λi按式(18)计算：

第8步：反复执行第2步到第7步，直至u＝l，然后将所有的pu′,u＝0,1,…,l-1重新拼接为公开音频p，将p作为信道传输音频并输出。

一种结合密钥和随机单位正交基的秘密音频恢复方法，包括以下步骤：

第1步：输入接收到的信道传输音频p′＝(pi′)l·n，将p′划分成长度为n的小段序列p′＝(pu)l，记为其中u＝0,1,…,l-1，由密钥key生成随机序列g＝(gu)l，初始化u＝0；

第2步：将gu作为密钥，依次生成n×n的随机矩阵x′，将x′的第0行置为1，对x′进行行单位正交化后记为x；对x′分别进行行正交化的具体方法是按式(1)分别对x′进行行单位施密特正交化：

x＝schimidt(x′)(1)

式(1)中，函数schimidt()为单位施密特正交化函数；

第3步：记x的n个行依次为xf,f＝0,1,…,n-1，且

第4步：从pu中提取出隐藏变换参数rn和q进制数表示的αi,i＝0,1,…,k-1，并进一步将其转换为n进制数序列r＝(ri)k-1和10进制数αi,i＝0,1,…,k-1，然后利用这些恢复出的变换参数来重构密图小块su，置u＝u+1；

从pu中提取出的隐藏变换参数提取出隐藏变换参数rm和q进制数表示的αi,i＝0,1,…,k-1，并进一步将其转换为10进制数序列r＝(ri)k-1和10进制数αi,i＝0,1,…,k-1的具体方法是：

第4.1步：按式(19)从中提取出的隐藏变换参数αk-1″、rn和sign2：

第4.2步：对于参数αi,i＝0,1,…,k-2则由式(20)中提取出

第4.3步：将sign2对应的2进制数序列和按式(9)计算α0,α1,…,αk-2：

第4.4步：将rn转化为n进制数序列r＝(ri)k-1；

式(19)和式(20)中，extract()为提取函数，t＝extract(h)表示从序列h＝(hi)d中提取出嵌入的t，具体的提取步骤是：

第4.1.1步：按式(14)计算出映射值vmap：

第4.1.2步：按式(21)提取出秘密值t

重构秘密音频小段序列su的具体方法是：

第4.5步：利用序列r＝(ri)k-1从矩阵x中选取对应的行得到k-1个一维序列，记为yi,i＝0,1,…,k-2；

第4.6步：将矩阵x中x0记为yk-1，然后按式(22)计算αk-1：

第4.7步：按式(23)将yi,i＝0,1,…,k-1分别与αi,i＝0,1,…,k-1相乘并累计作为恢复出的小段序列su：

第5步：重复第2到4步，直至u＝l，然后将所有的小块su,u＝0,1,…,l-1重新拼接作为解密后的秘密音频s输出。

实施例：

以下以javajdk1.7.0_09为案例实施环境，结合附图对本发明实施方式进行详细说明，其中图1是嵌入流程图，图2是提取流程图。

第1步：选取秘密音频如图3所示，为16位音频序列，即r＝16，长度为1128768，选取公开音频，如图4所示，为16位音频，长度为4797264；取n＝8,n＝34，即将秘密音频和公开音频对应的采样音频序列分别划分为长度为8和34的小段序列，由于是16位音频，因此每一个采样值包含2个字节，总共可将秘密音频和公开音频划分为l＝1128768/(2·8)＝70548个小段，将秘密音频和公开音频对应位置的小段序列分别作为秘密音频小段序列su和公开音频小段序列pu，利用密钥key可生成长度为70548的随机序列g＝(gu)l，初始化u＝0；

第2步：假设由gu可生成1个8×8随机矩阵，将其第0行所有元素置1后按式(1)分别进行行单位正交化，可得到：

第3步：若su＝[15331539160116811762190720232132]，以x的第2行x2＝[-0.37390.47310.28950.09290.4636-0.1712-0.3756-0.3983]为例，计算秘密音频小段序列su在x2正交基上的投影，即将su和x2代入式(2)并进一步通过式(3)可计算出对应的拟合参数

第4步：若取k＝3，su＝[15331539160116811762190720232132]，假设从中找到除的前2个最大幅值系数对应为：由su和x4＝[0.7885-0.07740.0179-0.12790.1309-0.0066-0.5568-0.1686]按式(2)计算出的以及由su和x2按式(2)计算得到的则记录索引序列为r＝(4,2)，由su和x0＝[0.35360.35360.35360.35360.35360.35360.35360.3536]按式(2)可得到

第5步：将r＝(4,2)视为n＝8进制数序列，因此可将其转化为10进制数rn＝4·8⁰+2·8¹＝20，若取q＝5,则按式(4)可计算出

第6步：将α2＝5012按式(5)转化为10进制数由式(6)可得到其对应的q＝5进制数序列长度将α0,α1按式(7)可表示为

由此可得将直接转化为2进制数sign2＝3，sign2对应的5进制数长度由式(11)可得取lr＝3,lm＝7,lsign＝1,n＝34，则由式(12)知：lpow+lint≤11，因此可取对应的5进制数长度分别为lpow＝2,lint＝9；

第7步：假设公开音频小段序列pu扫描的1维序列：

pu＝(2629,2589,2649,2831,2954,3060,2904,2830,2739,2844,2841,2862,2766,2511,2373,2128,2115,2222,2441,2705,2689,2561,2361,2162,2208,2170,2391,2478,2684,2470,2296,2112,2353,2552,2834)

按式(13)将α2″＝34540嵌入到pu[0,…,7]中(取lm＝8)，即取pu[0,…,33]的前8个元素构成：h＝(2629,2589,2649,2831,2954,3060,2904,2830)

将h＝(2629,2589,2649,2831,2954,3060,2904,2830)和q＝5,d＝8代入式(14)计算出vmap＝1881174，将t＝α2＝34540,vmap＝1881174代入式(15)和(16)计算：

vdiff＝-153634

将vdiff＝-153634,q＝5和h＝(2629,2589,2649,2831,2954,3060,2904,2830)代入式(17)和(18)得h′＝(2630,2587,2649,2832,2953,3061,2904,2828)，从而得到p′u[0,…,7]＝(2630,2587,2649,2832,2953,3061,2904,2828)；按式(15)将rn＝20嵌入到pu[8,…,10](取lm＝8,lr＝3,lm+lr-1＝10)，先取pu[8,…,10]构成h＝(2739,2844,2841)，将h＝(2739,2844,2841)和q＝5,d＝3代入式(14)计算出vmap＝109，将t＝rn＝20,vmap＝109代入式(15)和(16)计算vdiff＝-89，将vdiff＝-89,q＝5和h＝(2739,2844,2841)代入式(17)和(18)得h′＝(2740,2846,2842)，从而得到p′u[8,…,10]＝(2740,2846,2842)；按式(15)将sign2＝3嵌入到pu[11](取lm＝8,lr＝3,lsign＝1,lm+lr＝11,lm+lr+lsign-1＝11)，先取pu[11]构成h＝(2862)，将h＝(2862)和q＝5,d＝1代入式(14)计算出vmap＝2，将t＝sign2＝3,vmap＝2代入式(15)和(16)计算vdiff＝1，将vdiff＝1,q＝5和h＝(2862)代入式(17)和(18)得h′＝(2863)，从而得到p′u[11]＝(2863)；按式(15)将嵌入到pu[12,…,15](取lm＝8,lr＝3,lsign＝1,lpow＝2,lm+lr+lsign＝12,lm+lr+lsign+2·lpow-1＝15)，先取pu[12,13]构成h＝(2766,2511)，将h＝(2766,2511)和q＝5,d＝2代入式(14)计算出vmap＝21，将vmap＝21代入式(15)和(16)计算vdiff＝-19，将vdiff＝-19,q＝5和h＝(2766)代入式(17)和(18)得h′＝(2767,2512)，从而得到p′u[12,13]＝(2767,2512)，再取pu[14,15]构成h＝(2373,2128)，将h＝(2373,2128)和q＝5,d＝2代入式(14)计算出vmap＝13，将代入式(15)和(16)计算vdiff＝-11，将vdiff＝-11,q＝5和h＝(2373,2128)代入式(17)和(18)得h′＝(2372,2126)，从而得到p′u[14,15]＝(2372,2126)；按式(15)将嵌入到pu[16,…,33](取lm＝8,lr＝3,lsign＝1,lpow＝2,lint＝9,lm+lr+lsign+2·lpow＝16,lm+lr+lsign+2·lpow+2·lint-1＝33)，先取pu[16,…,24]构成：h＝(2115,2222,2441,2705,2689,2561,2361,2162,2208)

将h＝(2115,2222,2441,2705,2689,2561,2361,2162,2208)和q＝5,d＝9代入式(14)计算出vmap＝268000

将vmap＝268000代入式(15)和(16))计算vdiff＝96480，将vdiff＝vdiff＝96480,q＝5和h代入式(17)和(18)得：

h′＝(2115,2223,2440,2707,2688,2562,2362,2163,2208)

从而得到：

p′u[16,…,24]＝(2115,2223,2440,2707,2688,2562,2362,2163,2208)

再取pu[25,…,33]构成：

h＝(2170,2391,2478,2684,2470,2296,2112,2353,2552)

将h＝(2170,2391,2478,2684,2470,2296,2112,2353,2552)和q＝5,d＝9代入式(14)计算出vmap＝177200，将代入式(15)和(16)计算vdiff＝167027，将vdiff＝167027,q＝5和h＝(2170,2391,2478,2684,2470,2296,2112,2353,2552)代入式(17)和(18)得h′＝(2172,2391,2479,2685,2472,2294,2113,2355,2552)，从而得到p′u[25,…,33]＝(2172,2391,2479,2685,2472,2294,2113,2355,2552)；将α″2,α0,α1,rn隐藏到pu后，得：

pu′＝(2630,2587,2649,2832,2953,3061,2904,2828,2740,2846,2842,2863,2767,2512,2372,2126,2115,2223,2440,2707,2688,2562,2362,2163,2208,2172,2391,2479,2685,2472,2294,2113,2355,2552)

置u＝u+1；

第9步：反复执行上述步骤，直至u＝70548，然后将所有的pu′,u＝0,1,…,l-1重新拼接为公开音频，从而得到嵌有秘密音频的信道传输音频。

秘密音频提取方法具体实施步骤：

第1步：选取图6作为含密载体音频p′，对应为16bit音频，长度为4797264，将其划分为l＝70548个长度为n＝34的小段序列pu，利用密钥key生成长度为70548的随机序列g＝(gu)l，初始化u＝0；

第2步：假设：

pu＝(2630,2587,2649,2832,2953,3061,2904,2828,2740,2846,2842,2863,2767,2512,2372,2126,2115,2223,2440,2707,2688,2562,2362,2163,2208,2172,2391,2479,2685,2472,2294,2113,2355,2552)

以gu作为密钥生成1个8×8随机矩阵，按式(1)将其第0行置为1后进行行单位正交化得到：

第3步：从pu中按式(19)依次提取出隐藏变换参数，其中参数α″2，先取pu[0,…,7]构成h＝(2630,2587,2649,2832,2953,3061,2904,2828)，将h＝(2630,2587,2649,2832,2953,3061,2904,2828)和q＝5,d＝8代入式(14)计算出vmap＝34540，按式(21)得t＝34540，从而得α″2＝t＝34540；其中参数rn，先取pu[8,…,10]构成h＝(2740,2846,2842)，将h＝(2740,2846,2842)和q＝5,d＝3代入式(14)计算出vmap＝20，按式(21)得t＝20，从而得rn＝t＝20；其中参数sign2，先取pu[11]构成h＝(2863)，将h＝(2863)和q＝5,d＝1代入式(14)计算出vmap＝3，按式(21)得t＝3，从而得sign2＝t＝3；其中参数先取pu[12,13]＝(2767,2512)构成h＝(2767,2512)，将h＝(2767,2512)和q＝5,d＝2代入式(14)计算出vmap＝2，按式(21)得t＝2，从而得再取pu[14,15]＝(2372,2126)构成h＝(2372,2126)，将h＝(2372,2126)和q＝5,d＝2代入式(14)计算出vmap＝2，按式(21)得t＝2，从而得其中参数先取pu[16,…,24]＝(2115,2223,2440,2707,2688,2562,2362,2163,2208)构成：

h＝(2115,2223,2440,2707,2688,2562,2362,2163,2208)

将h＝(2115,2223,2440,2707,2688,2562,2362,2163,2208)和q＝5,d＝9代入式(14)计算出vmap＝364480，按式(21)得t＝364480，从而得再取pu[25,…,33]＝(2172,2391,2479,2685,2472,2294,2113,2355,2552)构成h＝(2172,2391,2479,2685,2472,2294,2113,2355,2552)，将h＝(2172,2391,2479,2685,2472,2294,2113,2355,2552)和q＝5,d＝9代入式(14)计算出vmap＝344227，按式(21)得t＝344227，从而得将行号索引rn＝20转化为n＝8进制数序列r＝(4,2)，符号位sign2＝3转化为2进制数序列由式(9)代入和得α0＝(-1)¹×364480×10^-2-1＝-364.48和α1＝(-1)¹×344227×10^-2-1＝-344.227；

第4步：根据行号索引r＝(4,2)从矩阵x得到对应的行

，

记y2＝x0＝[0.35360.35360.35360.35360.35360.35360.35360.3536]，按式(22)代入α2″＝34540得到：

将和y2＝[0.35360.35360.35360.35360.35360.35360.35360.3536]分别与α0＝3-6α14＝.代入式(238)得到提取出的秘密音频小段序列su＝[1613.3251637.3571665.85261786.6471564.7001833.3072104.2581970.550]，置u＝u+1；

第5步：反复执行第2步到第4步，直至u＝70548，然后将所有的小块su,u＝0,1,…,l-1拼接作为解密后的秘密音频s输出。

图5是以图3为秘密音频，图4为公开音频，取k＝2，按图1所示流程嵌入参数后的公开音频1，相对于图4的snr为67.0138db；；

图6是以图3为秘密音频，图4为公开音频，取k＝3，按图1所示流程嵌入参数后的公开音频2，相对于图4的snr为65.0563db；

图7是以图5为含密音频按图2所示提取流程重构的秘密音频1，相对于图3的snr为19.9059db；

图8是以图6为含密音频按图2所示提取流程重构的秘密音频2，相对于图3的snr为21.4306db。

以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。