本发明属于射频水印技术领域,更具体地,涉及一种基于k-s检测的射频水印检测方法。
背景技术:
经过数十年的快速发展,数字水印已经不再局限于数字图像、音频、视频等数字产品中,应用于射频系统确保无线系统的使用安全,成为数字水印技术新的发展方向。应用于射频系统的数字水印也称作射频虽然数字水印。如图1所示,典型的基于基带信号水印嵌入算法的射频水印系统,其核心部分包括水印嵌入和水印提取。水印嵌入部分是以帧为单位进行的,嵌入水印信息后,合并为完整信号经上变频、功率放大后发射。水印提取部分在下变频后首先进行帧同步和数据同步,然后进行水印信息的提取。目前的研究主要集中在水印的嵌入和提取上。
常规的水印检测都是基于多媒体信号的,国内外还没有针对调制信号的水印本身的检测。目前几乎所有关于射频水印的研究都是在已知编码方式的情况下进行的。但是,如果接收端仅知道从信道中传过来的一组信号值而对其它的情况一无所知,又该如何判断这一组信号中是否存在射频水印,目前还没有人对这方面的问题作出考虑。
技术实现要素:
针对现有技术的缺陷,本发明的目的在于解决接收端仅知道一组信号值情况如何检测射频水印的技术问题。
为实现上述目的,第一方面,本发明实施例提供了一种基于k-s检测的射频水印检测方法,该方法包括以下步骤:
s1.去除接收到的一组射频信号中的载体信号,得到一组新信号值{z};
s2.计算该组新信号值的累积分布函数fz(z)和噪声累积分布函数fn(z);
s3.对fz(z)和fn(z)进行k-s检测,根据检测结果判断射频水印是否存在。
具体地,步骤s1包括以下子步骤:
s101.估计该组射频信号中每个信号的载体信号,载体信号的计算公式如下:
si=e(|y|)sign(yi)
其中,yi表示第i个射频信号,|·|表示绝对值算符,e(·)表示期望值,sign(yi)为符号函数,当yi大于0时,其值为1,当yi小于0时,其值为-1,i=1,2,…,l,l表示该组射频信号总个数;
s102.将该组射频信号中每个信号去除其载体信号,得到一组新信号值{z}={z1,z2,…,zl},其中,zi=yi-si。
具体地,计算所述该组新信号值中的噪声累积分布函数fn(z)包括以下步骤:
(1)假定信号{z}中只包含噪声,利用该组新的信号值,计算z的功率值
其中,var(z)表示z的方差,σz表示噪声信号的标准差;
(2)基于σz得到噪声累积分布函数fn(z),计算公式如下:
具体地,计算所述该组新信号值的累积分布函数fz(z)包括以下步骤:
(1)对该组新的信号值z升序排列{z′}={z′1,z′2,…z′l},其中,z′1≤z′2≤…≤z′l;
(2)计算fz(z)在l个z′i上的取值,计算公式如下:
其中,i表示z′i排序后的位置,i=1,2,…,l,l表示该组射频信号总个数。
具体地,步骤s3具体包括以下子步骤:
s301.构造统计量
s302.固定d的显著性水平为0.05,当
第二方面,本发明实施例提供了一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述第一方面所述的射频水印检测方法。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
1.本发明针对的载体不是多媒体信号(例如图片、声音等),而是数字调制的射频信号,该信号可在接收端恢复出来,基于该特性,将载体信号解调出来之后去除,再进行k-s检测,消除了载体对水印的干扰,从而提高k-s检测的准确率。
2.本发明将k-s检测应用到调制信号水印检测中,由于嵌入水印前后,信号的统计分布规律会发生变化,k-s检测可以检测到统计规律的变化,判断水印是否存在,从而在在接收端仅知道一组信号值时,发现信号中隐藏的水印信号,该方法应用范围更加广泛。
附图说明
图1为现有技术中典型的射频水印系统结构示意图;
图2为本发明实施例提供的一种基于k-s检测的射频水印检测方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图2所示,一种基于k-s检测的射频水印检测方法,该方法包括以下步骤:
s1.去除接收到的一组射频信号中的载体信号,得到一组新信号值;
s2.计算该组新信号值的累积概率函数fz(z)和噪声累积分布函数fn(z);
s3.对fz(z)和fn(z)进行k-s检测,根据检测结果判断射频水印是否存在。
步骤s1.去除接收到的一组射频信号中的载体信号,得到一组新信号值。
在bpsk调制中,接收到的射频信号可表示为
y=x+s+n
其中,x表示功率为
为了精确分析,本发明实施例选择在接收端每次接收的信号个数为l=500,并将接收到的一组射频信号表示为
{y}={y1,y2,…yl}
s101.估计该组射频信号中每个信号的载体信号。
利用载体信号功率远大于水印信号功率这个特性,将接收到的一组bpsk射频信号进行处理,获得估计的载体信号s。
si=e(|y|)sign(yi)
其中,yi表示第i个射频信号,|·|表示绝对值算符,e(·)表示期望值,sign(yi)为符号函数,当yi大于0时,其值为1,当yi小于0时,其值为-1,i=1,2,…,l。
s102.将该组射频信号中每个信号去除其载体信号,得到一组新信号值。
将接收到的一组信号值y除去估计的载体信号值s,最后得到一组新的信号值z。接收到的射频信号y消除载体信号s后,得到的新信号值可以表示为{z}={z1,z2,…,zl},其中,zi=yi-si。
步骤s2.计算该组新信号值的累积分布函数fz(z)和噪声累积分布函数fn(z)。
s201.假定信号{z}中只包含噪声,利用该组新的信号值,计算z的功率值
其中,var(z)表示z的方差,σz表示噪声信号的标准差。
s202.基于σz得到噪声累积分布函数fn(z)。
s203.对该组新的信号值进行升序排序,利用排序后的信号值计算fz(z)。
首先,将新构造的信号z升序排列{z′}={z′1,z′2,…z′l},其中,z′1≤z′2≤…≤z′l。然后计算fz(z)在z′i(1≤i≤l)上的l个取值如下:
步骤s3.对fz(z)和fn(z)进行k-s检测,根据检测结果判断射频水印是否存在。
本发明构造统计量
以上,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应该以权利要求的保护范围为准。