调制信号识别docxWord文档格式.docx
- 文档编号:18113574
- 上传时间:2022-12-13
- 格式:DOCX
- 页数:10
- 大小:21.37KB
调制信号识别docxWord文档格式.docx
《调制信号识别docxWord文档格式.docx》由会员分享,可在线阅读,更多相关《调制信号识别docxWord文档格式.docx(10页珍藏版)》请在冰豆网上搜索。
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算法的主要程序如下:
functioncoefs=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);
ifwtype==5,psi_integ=conj(psi_integ);
endwtype=5说明如果是没有尺度
函数的复小波,将小波积分值取复共轭
xval=xval-xval
(1);
dx=xval
(2);
xmax=xval(end);
ind=1;
fork=1:
nbscales计算各个尺度的信号的连续小波变换值a=scales(k);
j=[1+floor([0:
a*xmax]/(a*dx))];
设置j,对积分值psi_integ进行采样例a=4,(0:
1:
4*xmax)/4*dx
iflength(j)==1,j=[11];
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)是与(tb)的分段积分进行卷积,所以在程序中出现了一个diffa
运算,对相邻的两个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)可以描述为一个带通
滤波器组的脉冲响应,af0/f,f0是带通滤波器的中心频率,
f是要分析信号
的频率。
随着a的变化,这样的一组滤波器,在时间轴上滑动,即
改变,信号
的不同频率成分将有可能进入其通带,对小波变换的模起到主要作用,当信号的
某个频率不但进入其通带而且其频率恰好等于滤波器组的中心频率f0时,将使
得小波变换在此区域附近产生一个极大值,即W(a,)局部最大。
提取所有时间轴上的模极大值,得到一条脊线,即为小波脊线法。
具体方法
是,对任一固定时刻,遍历小波的尺度a,找到W(a,)在所有尺度上的最大
值。
之后找到每个最大值所对应的尺度,根据尺度和频率的对应关系,af0/f,
将尺度转换成频率,根据极大值的产生原理,这个频率就是输入信号的频率。
对每个时刻进行如此循环操作,便得到输入信号的频率曲线。
问题:
1、主要提取信号的频率特征,通过分析频率曲线的阶数P,可识别FSK信号和ASK、PSK信号。
如果P1,则此信号是ASK或PSK信号;
如果P1,则此信号是FSK信号,并且根据频率曲线可知此信号在某个时刻的频率。
对识
别FSK信号比较有效。
2、当信号的频率比较高时,识别效果比较好。
3、由于cwt变换在信号跳变处的混叠发散现象,在最大值搜索中,搜到一些伪最大值,影响了真实频率的提取。
ASK信号识别,SNR=5.7dB
PSK信号,SNR=3.6dB
FSK信号,SNR=3.8dB
高斯噪声
在通信理论中,最重要的概率密度函数是高斯或正态概率密度函数。
统计学
中的中心极限定理指出:
在非常宽的条件下,大量N个统计独立的随机变量xi之
N
和Zxi的分布律,不管每个xi的分布律如何,在N的极限情况下,趋
i1
于高斯正态分布。
因此,高斯噪声是指其统计分布服从正态分布的噪声。
根据中心极限定理,高斯噪声是普遍存在的一种随机信号,这也是在分析设计中常常采用高斯噪声假设的原因。
七、过零点检测
过零点抽样,在现代模式识别中是一个非常具有吸引力的工具,具有广泛的应用。
当输入信号穿越零值点时,过零点抽样记录下这些时刻。
当接收信号的相位变化时,过零点抽样提供了大量的有效信息,可以进行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(i
1)
x(i)
i
1,2,...,N
z(i)
y(i
y(i)
2
2.相关过零变量的概率密度函数
设接收信号(t)由正弦型信号和噪声组成:
(t)
Acos2fctv(t)
0tT0
第i个过零点为
0.5
i1,2,...,N
(i)
2fc
(i)是由噪声和误差引起的随机变量。
在高
CNR下,
(i)的密度函数是高斯
分布的,其均值为零,方差为
2(2
fc)2
是CNR,定义为
A
。
那么过零间隔y(i)为
(0)
因此,y(i)是关于频率量度的一个序列。
z(i)y(i1)y(i)是对y(i)变换的一个
量度序列。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
上式中
(i
(i),(i)的值依赖于载波的相位。
如果
2fcx(i)在/2
(mod2
)附近,则
1[v(x(i))
v(x(i
1))]
fcA
如果2fcx(i)在3
/2(mod2
[v(x(i))
推导:
是使r(t)为零的点,
当tx(i),r(t)
0。
v[x(i)]
Acos2fc[x(i)]
Acos2fc[i
(i)]
Acos[(i
)2fc
(i)]
A{sin(i
)sin2
fc
(i)}
如果2
fcx(i)在
/2(mod2
)附近,则i
也大约在
)附近。
Asin2
fc
当2fc
(i)很小时,v[x(i)]
A2fc
(i)。
(i)
2A
则(i)
(i1)
{v[x(i
1)]v[x(i)]}
2A
fcx(i)在3
/2(mod2
也大约在3
/2(mod2
v[x(i)]
Asin2fc
1)]v[x(i)]}。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在信号发生变化的时刻,例如幅度,频率,相位发生改变时,所采样的过零
点也会产生不规则的突变,称这些点为码间异变点。
码间异变点影响了对信号频
率和码元速率的估计,因此在估计信号频率的时候,需要剔除这些点。
由于y(i)
序列是关于频率量度的一个序列,将y(i)中过大的点除掉,得到ya(i)序列。
(具体方法是:
首先考察z(i)序列,计算密集部分的z(i)的标准方差za,找
出对应于z(i)
3.034za的y(i)点,将其剔除。
3.034
za为一个适合的门限。
)
,做出fc的直方图,可以分
ya(i)序列的长度为Ny,则信号频率fc
2ya(i)
辨出信号是单频还是多频,如果是多频有哪几个频率组成。
由此可以识别出ASK,PSK和FSK信号,还可以对FSK信号进行类内识别。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 调制 信号 识别 docx