调制信号识别.docx
调制信号的小波分析
一、小波函数简介
1.Haar 小波
最简单的小波函数, Haar 小波是离散的,与阶跃信号相似,同 Daubechies db1 小波是一样的。
2. Daubechies小波
Daubechies 小波是紧支正则小波,便于进行离散小波分析。这类小波没有显
式的表达式,除了 db1( Haar)。然而它的传递函数的模的平方是有简单的表达
式的。
3.Biorthogonal 小波
此类小波具有线性相位,用于信号和图像重建。
4. Coiflet 小波
这个小波族是 I.Daubechies应 R.Coifman 的要求所创建的, coif N 较 dbN 有更好的对称性。
5.Symlets 小波
此小波由 Daubechies提出,作为对 db 小波族的修正,是一种近似对称小波,它和 db 小波族的性质是近似的。
6.Morlet 小波
x2
其尺度函数不存在,小波函数为 ( x) e 2 cos5x ,Morlet 小波不满足容许性条件。
7.Mexican Hat 小波
小波函数为 ( x) ( 2
1x2
4 )(1
x2 )e 2,它是 Gaussian概率密度函数的二
阶
3
导数,由于它不存在尺度函数,因此不具有正交性。
8.Meyer 小波
Meyer 小波的尺度函数和小波函数都在频域中定义,都具有显式的表达式。
二、连续小波变换
从数学上来说,傅里叶变换就是将信号 f (t ) 乘以一个复指数后在所有的时间
域上求和。变换的结果就是傅里叶系数。
相似的,连续小波变换(CWT)定义为,将信号乘以由尺度和位移确定的
小波函数后,再在整个时间轴上相加。CWT 的变换结果是很多小波系数C,C 是尺度和位移的函数。
大尺度对应于时间上伸展大的小波,小波伸展地越大,所比较的信号段就越长,所以小波系数所量度的信号特征也就越粗糙。
在计算机中,任何实数域的信号处理都是对离散信号的操作,那么, CWT 的连续性及它与 DWT 的区别表现在尺度的选取和对位移的操作。与离散小波变换不同的是,只要在计算机的计算能力之内, CWT 可以在每一个尺度上计算;在位移上连续是指小波可以在待分析函数的整个域上进行平滑的移动。
三、离散小波变换
对于大多数信号来说,低频部分往往是最重要的,给出了信号的特征。而高频部分则与噪音及扰动联系在一起。将信号的高频部分去掉,信号的基本特征仍然可以保留。
信号的概貌主要是系统大的、低频的成分,大尺度;而细节往往是信号局部、高频成分,小尺度。
分解算法: 1.产生两组系数:概貌系数cA1 和细节系数 cD1。通过低通滤波器Lo_D 卷积信号 s 得到 cA1,通过高通滤波器 Hi_D 卷积 s 得到 cD1,之后进行二抽取。每个滤波器的长度是 2N。如果 n = length (s),那卷积后概貌信号和细节信号的长度为 n + 2N - 1,进行二抽取之后 cA1 和 cD1 的长度为 floor ((n-1)/2)+N。
关于 matlab 中 cwt 算法的分析
cwt 算法的主要程序如下:
function coefs = cwt(signal,scales,wname,plotmode,xlim)
precis = 10;
signal = signal(:)'; 输入信号
len= length(signal);
coefs = zeros(length(scales),len);设置小波系数数组
nbscales = length(scales);
[psi_integ,xval] = intwave(wname,precis); 根据不同的小波计算其积分值
wtype = wavemngr('type',wname);
if wtype==5 , psi_integ = conj(psi_integ); end wtype=5 说明如果是没有尺度
函数的复小波,将小波积分值取复共轭xval = xval-xval(1);
dx= xval(2);
xmax = xval(end);
ind = 1;
for k = 1:nbscales计算各个尺度的信号的连续小波变换值
a = scales(k);
j = [1+floor([0:a*xmax]/(a*dx))];设置j,对积分值psi_integ进行采
样例 a=4,(0:1:4*xmax ) /4*dx if length(j)==1 , j = [1 1]; end
f = fliplr(psi_integ(j)); 将积分值即小波滤波器系数反转 coefs(ind,:) = -
sqrt(a)*wkeep(diff(conv(signal,f)),len);
将信号与小波系数 f 进行卷积,再差分,截取中间数值ind= ind+1;
end
dummyCoefs = coefs;
dummyCoefs = abs(dummyCoefs);
plotCOEFS(axeAct,dummyCoefs,plotPARAMS); 可见,cwt 画出的是小波变换系数
的绝对值 dummyCoefs,而返回
值是 coefs,不是绝对值。
算法理论分析:
由于 s(k ) 是与
(
t b
) 的分段积分进行卷积,所以在程序中出现了一个 diff a
运算,对相邻的两个 coefs 值进行相减,因此在变换图中,在不同频率变换处,出现混叠发散现象,难以得到准确清晰的频率分辨。
四、调制信号识别
(一)利用模式识别方法分类调制类型,所用的分类特征归纳起来主要有以下几
种:
1.直方图特征
Liedtke 等人利用幅度、频率和相位的直方图分类通信信号。
2.统计矩特征
由于直方图分类特征的维数太大,现在常用的分类特征是信号瞬时幅度、相位和频率函数的各阶统计矩特征。
3.变换域特征
把信号变换到其它特征空间,利用新特征空间中的特征参数来识别调制类型。
(二)模最大值法
对于 3 种基本的调制信号: ASK,FSK 和 PSK 信号,可以将它们进行小波变换,分析变换后的参数特征来识别。
采用提取模最大值的方法来提取三种信号在小波变换域中的特征进行识别。
模极大值的定义:对 x0邻域内的任意点 x ,若在尺度 s 上满
足Wf ( x, s) Wf ( x0 , s) ,则称( s, x0)为一模极大值点, Wf (x0 , s) 称为在(s, x0)点
的小波变换模极大值。
小波变换模极大值携带了信号的大部分信息,信号的所有奇异值点都被极大值点定位。 Mallat 证明了,通过模极大值可以对原始信号进行重建,得到一个近
似的逼近。因此提取模极大值可以分析信号的特征。
小波变换为什么能产生一个极大值?小波函数a, (t) 可以描述为一个带通滤波器组的脉冲响应, a f 0 / f , f 0是带通滤波器的中心频率, f 是要分析信号
的频率。随着 a 的变化,这样的一组滤波器,在时间轴上滑动,即改变,信号的不同频率成分将有可能进入其通带,对小波变换的模起到主要作用,当信号的
某个频率不但进入其通带而且其频率恰好等于滤波器组的中心频率 f 0时,将使
得小波变换在此区域附近产生一个极大值,即W (a, ) 局部最大。
提取所有时间轴上的模极大值,得到一条脊线,即为小波脊线法。具体方法是,对任一固定时刻,遍历小波的尺度 a ,找到W (a, )在所有尺度上的最大
值。之后找到每个最大值所对应的尺度,根据尺度和频率的对应关系,a f 0 / f ,
将尺度转换成频率,根据极大值的产生原理,这个频率就是输入信号的频率。对
每个时刻进行如此循环操作,便得到输入信号的频率曲线。
问题:1、主要提取信号的频率特征,通过分析频率曲线的阶数P ,可识别FSK 信号和 ASK 、PSK 信号。如果 P 1 ,则此信号是 ASK 或 PSK 信号;如果 P 1, 则此信号是 FSK 信号,并且根据频率曲线可知此信号在某个时刻的频率。对识
别 FSK 信号比较有效。
2、当信号的频率比较高时,识别效果比较好。
3、由于 cwt 变换在信号跳变处的混叠发散现象,在最大值搜索中,搜到一些伪最大值,影响了真实频率的提取。
ASK 信号识别, SNR=5.7dB
PSK 信号, SNR=3.6dB
FSK 信号, SNR=3.8dB
高斯噪声
在通信理论中,最重要的概率密度函数是高斯或正态概率密度函数。统计学中的中心极限定理指出:在非常宽的条件下,大量 N 个统计独立的随机变量x i之
N
和 Z x i的分布律,不管每个x i的分布律如何,在N的极限情况下,趋
i 1
于高斯正态分布。因此,高斯噪声是指其统计分布服从正态分布的噪声。根据中
心极限定理,高斯噪声是普遍存在的一种随机信号,这也是在分析设计中常常采用
高斯噪声假设的原因。
七、过零点检测
过零点抽样,在现代模式识别中是一个非常具有吸引力的工具,具有广泛的应用。当输入信号穿越零值点时,过零点抽样记录下这些时刻。当接收信号的相位变
化时,过零点抽样提供了大量的有效信息,可以进行CW,AM,FSK,PSK 等信号的
识别。
1.3 个序列
利用接收到的信号,可以创建 3 个序列x(i), y(i), z(i )。当接收信号进行过零点
抽样后,过零的时刻组成了一个过零序列{ x(i ), i1,2,..., N } 。为了从 x(i ) 中提取相位和频率信息,又创建了y(i ) 和 z(i ) 两个序列。
y(i )x(i1)x(i )i1,2,..., N1
z(i )y(i1)y(i)i1,2,..., N2
2.相关过零变量的概率密度函数
设接收信号 (t ) 由正弦型信号和噪声组成:
(t) Acos2 f c t v(t )
0 t T 0
第 i 个过零点为
i 0.5 i 1,2,..., N x(i )
(i )
2 f c
(i ) 是由噪声和误差引起的随机变量。在高 CNR 下, (i ) 的密度函数是高斯
分布的,其均值为零,方差为
2
1
2(2 f c ) 2
是 CNR ,定义为
A 2
。那么过零间隔 y(i ) 为
2 (0)
1 (i )
y(i)
2 f c
因此, y(i ) 是关于频率量度的一个序列。 z(i ) y(i 1) y(i) 是对 y(i) 变换的一个
量度序列。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
上式中 (i )
(i
1)
(i ) , (i ) 的值依赖于载波的相位。如果
2 f c x(i ) 在 / 2
( mod2 )附近,则
(i )
2 1
[ v( x(i)) v(x(i
1))]
f c A
如果 2 f c x(i ) 在 3 / 2 (mod 2 )附近,则
(i )
1 v(x(i
1))]
2 [ v( x(i))
f c A
推导:
x(i ) 是使 r (t) 为零的点,
当 t x(i ) , r (t )
0 。
v[ x(i )]
Acos 2 f c [ x(i )]
Acos2 f c [ i
0.5
(i)]
2 f c
Acos[(i) 2 f c(i )]
2
A{ sin(i)sin 2 f c(i)}
2
如果 2 f c x(i ) 在/ 2 ( mod2)附近,则 i
2
也大约在/ 2 ( mod2)附近。
v[ x(i )]Asin 2f c(i )
当 2 f c(i ) 很小时, v[ x(i )] A 2 f c(i ) 。
v[ x(i )]
(i)
f c
2 A
则 (i )(i 1)
1
{ v[ x(i1)] v[ x(i )]}。
(i )
2A f c
如果 2 f c x(i ) 在3/ 2(mod 2)附近,则 i
2
也大约在 3/ 2(mod2)附近。
v[ x(i)]Asin 2 f c(i )
v[ x(i)]
(i )
f c
2 A
则 (i )(i 1)(i )1{ v[ x(i1)] v[ x(i )]} 。
2A f c %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%在信号发生变化的时刻,例如幅度,频率,相位发生改变时,所采样的过零
点也会产生不规则的突变,称这些点为码间异变点。码间异变点影响了对信号频
率和码元速率的估计,因此在估计信号频率的时候,需要剔除这些点。由于y(i)序列是关于频率量度的一个序列,将y(i) 中过大的点除掉,得到y a (i) 序列。
(具体方法是:首先考察z(i ) 序列,计算密集部分的z(i ) 的标准方差za ,找出对应于 z(i ) 3.034 za的 y(i ) 点,将其剔除。 3.034za 为一个适合的门限。)
1
,做出 f c的直方图,可以分y a (i) 序列的长度为 N y,则信号频率 f c
2 y a (i )
辨出信号是单频还是多频,如果是多频有哪几个频率组成。由此可以识别出
ASK,PSK 和 FSK 信号,还可以对 FSK 信号进行类内识别。