【信号隐藏】基于奇异值分解SVD实现音频水印嵌入提取含加噪滤波检验鲁棒性含误码率附Matlab代码

天天Matlab代码科研顾问

于 2024-04-05 23:44:49 发布

阅读量852

点赞数 19

分类专栏：信号处理文章标签：音视频 matlab 开发语言

本文链接：https://blog.csdn.net/m0_57702748/article/details/137412649

版权

信号处理专栏收录该内容

370 篇文章 57 订阅

订阅专栏

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，代码获取、论文复现及科研仿真合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab完整代码及仿真定制内容点击👇

智能优化算法神经网络预测雷达通信无线传感器电力系统

信号处理图像处理路径规划元胞自动机无人机

物理应用机器学习

🔥 内容介绍

本文提出了一种基于奇异值分解（SVD）的音频水印嵌入提取方法，该方法具有较好的鲁棒性，能够抵抗加噪和滤波处理。该方法首先对音频信号进行SVD分解，然后将水印信息嵌入到奇异值中。在提取阶段，通过计算水印信息与提取信号的奇异值之间的相关性，可以提取出嵌入的水印信息。实验结果表明，该方法能够有效地抵抗高斯白噪声、粉红噪声和低通滤波处理，误码率较低。

引言

音频水印技术是一种将版权信息或其他信息嵌入到音频信号中，而不影响音频信号的感知质量的技术。音频水印技术在版权保护、内容认证和信息隐藏等领域有着广泛的应用。

传统的音频水印技术大多基于时域或频域处理，这些技术在抵抗加噪和滤波处理方面表现不佳。奇异值分解（SVD）是一种信号处理技术，它可以将信号分解为奇异值、左奇异向量和右奇异向量的乘积。奇异值表示信号的能量分布，奇异向量表示信号的结构信息。

本文提出了一种基于SVD的音频水印嵌入提取方法，该方法将水印信息嵌入到奇异值中。该方法具有较好的鲁棒性，能够抵抗加噪和滤波处理。

嵌入算法

该方法的嵌入算法如下：

对音频信号进行SVD分解，得到奇异值矩阵S、左奇异向量矩阵U和右奇异向量矩阵V。
将水印信息转换为二进制序列。
将二进制序列嵌入到奇异值矩阵S中。
重构水印嵌入后的音频信号。

提取算法

该方法的提取算法如下：

对水印嵌入后的音频信号进行SVD分解，得到奇异值矩阵S'、左奇异向量矩阵U'和右奇异向量矩阵V'。
计算S'和S之间的相关性矩阵R。
将R中的最大相关性值对应的二进制序列提取出来。
将二进制序列转换为水印信息。

实验结果

为了验证该方法的鲁棒性，我们进行了以下实验：

**加噪实验：**向水印嵌入后的音频信号添加高斯白噪声和粉红噪声，信噪比分别为20dB、15dB和10dB。
**滤波实验：**对水印嵌入后的音频信号进行低通滤波处理，截止频率分别为1kHz、2kHz和3kHz。

实验结果表明，该方法能够有效地抵抗加噪和滤波处理，误码率较低。

结论

本文提出了一种基于SVD的音频水印嵌入提取方法，该方法具有较好的鲁棒性，能够抵抗加噪和滤波处理。该方法将水印信息嵌入到奇异值中，通过计算水印信息与提取信号的奇异值之间的相关性，可以提取出嵌入的水印信息。实验结果表明，该方法能够有效地抵抗高斯白噪声、粉红噪声和低通滤波处理，误码率较低。

📣 部分代码

function out=zigzag(in)% Zig-zag scanning% IN specifies the input matrix of any size% OUT is the resulting zig-zag scanned (1-D) vector% having length equal to the total number of elements in the 2-D input matrix%% As an example,% IN = [1  2  6  7%    3  5  8  11%    4  9  10  12];% OUT = ZIGZAG(IN)% OUT=%  1     2     3     4     5     6     7     8     9    10    11    12[num_rows num_cols]=size(in);% Initialise the output vectorout=zeros(1,num_rows*num_cols);cur_row=1;  cur_col=1;  cur_index=1;% First element%out(1)=in(1,1);while cur_row<=num_rows & cur_col<=num_cols  if cur_row==1 & mod(cur_row+cur_col,2)==0 & cur_col~=num_cols    out(cur_index)=in(cur_row,cur_col);    cur_col=cur_col+1;              %move right at the top    cur_index=cur_index+1;      elseif cur_row==num_rows & mod(cur_row+cur_col,2)~=0 & cur_col~=num_cols    out(cur_index)=in(cur_row,cur_col);    cur_col=cur_col+1;              %move right at the bottom    cur_index=cur_index+1;      elseif cur_col==1 & mod(cur_row+cur_col,2)~=0 & cur_row~=num_rows    out(cur_index)=in(cur_row,cur_col);    cur_row=cur_row+1;              %move down at the left    cur_index=cur_index+1;      elseif cur_col==num_cols & mod(cur_row+cur_col,2)==0 & cur_row~=num_rows    out(cur_index)=in(cur_row,cur_col);    cur_row=cur_row+1;              %move down at the right    cur_index=cur_index+1;      elseif cur_col~=1 & cur_row~=num_rows & mod(cur_row+cur_col,2)~=0    out(cur_index)=in(cur_row,cur_col);    cur_row=cur_row+1;    cur_col=cur_col-1;  %move diagonally left down    cur_index=cur_index+1;      elseif cur_row~=1 & cur_col~=num_cols & mod(cur_row+cur_col,2)==0    out(cur_index)=in(cur_row,cur_col);    cur_row=cur_row-1;    cur_col=cur_col+1;  %move diagonally right up    cur_index=cur_index+1;      elseif cur_row==num_rows & cur_col==num_cols  %obtain the bottom right element        out(end)=in(end);              %end of the operation    break                    %terminate the operation    endend