CN105743830B - 一种mfsk信号的解调方法 - Google Patents
一种mfsk信号的解调方法 Download PDFInfo
- Publication number
- CN105743830B CN105743830B CN201610116988.7A CN201610116988A CN105743830B CN 105743830 B CN105743830 B CN 105743830B CN 201610116988 A CN201610116988 A CN 201610116988A CN 105743830 B CN105743830 B CN 105743830B
- Authority
- CN
- China
- Prior art keywords
- ramanujan
- summation
- fourier
- mfsk
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/10—Frequency-modulated carrier systems, i.e. using frequency-shift keying
- H04L27/14—Demodulator circuits; Receiver circuits
Abstract
本发明公开了一种MFSK信号解调的方法,其特征在于,该方法包括:对接收MFSK信号做过采样并做短时拉马努金傅里叶变换,得到拉马努金傅里叶谱,取模平方得到拉马努金傅里叶谱图;对拉马努金傅里叶谱图做时间采样,得到码元中心的拉马努金傅里叶谱图;对码元中心的拉马努金傅里叶谱图做谱峰搜索得到谱峰位置频率函数;对谱峰位置频率函数使用时间上逐点的谱峰选择算法得到解调出的数据信息。本发明解决了存在多普勒频移情况下MFSK信号准确解调的问题。
Description
技术领域
本发明涉及信号解调方法,特别涉及一种MFSK信号的解调方法。
背景技术
信号解调是各种通信系统中的一个重要环节,发送端对所欲传送的消息进行调制,产生携带这一消息的信号,接收端必须通过解调才能恢复所传送的消息。解调方法的好坏对通信系统总体性能起着非常重要的作用。在信号传输过程中,由于如加性高斯白噪声AWGN、多普勒频移等影响,如何在接收到的包含噪声的信号中准确恢复调制信息是迫切需要解决的问题。
传统的MFSK信号解调方法包括:包络检波、相干解调和基于短时傅里叶变换ST-DFT的解调方法等。其中,包络检波、相干解调,为MFSK信号解调常用方法,具有比较好的检测性能,但对多普勒频移的处理不够理想;基于短时傅里叶变换ST-DFT的解调方法比较适用于存在较大多普勒频移的情况下,但仍需要提高对抗噪声性能。
发明内容
本发明旨在解决存在多普勒频移情况下信号准确解调的问题。
本发明实施例提供了一种MFSK信号解调方法,其特征在于,该方法包括:
步骤1、对接收到的MFSK信号s(t)过采样得到离散接收信号s(n),其中t为连续的时间,n为离散的时间序号,n=1,2,…,N,对所述离散接收信号s(n)做短时拉马努金傅里叶变换ST-RFT,得到拉马努金傅里叶谱S(m,q),时间m=1,2,…,N,频率q=1,2,…,N,设定所述拉马努金傅里叶谱S(m,q)的拉马努金傅里叶谱图Em(q)=|S(m,q)|2;
进一步地,步骤1中MFSK信号表达式为
其中,Tb为码元周期,g(t)为0到Tb之间幅度为1的门函数,
M个载频设定采样频率fs,则采样周期Ts=1/fs,离散MFSK信号表达式为
短时拉马努金傅里叶变换ST-RFT的表达式为:
为滑动的时间窗,为拉马努金代数和,φ(q)为欧拉函数。
步骤2、对所述拉马努金傅里叶谱图Em(q)做时间采样,得到码元中心处拉马努金傅里叶谱图
Tb为所述接收到的MFSK信号码元周期,Ts为采样周期,round[·]为求圆整,l为码元序号,对所述码元中心处拉马努金傅里叶谱图做谱峰搜索,得到第j高的谱峰位置频率函数qpeak(j,l);
步骤3、设定MFSK载频qc,载频序号c=1,2,…,M,M为MFSK信号阶次,对所述第j高的谱峰位置频率函数qpeak(j,l)与所述MFSK载频qc做时间上逐点的谱峰选择算法,得到解调出的数据信息date(l)。
进一步地,步骤3中所述时间上逐点的谱峰选择算法包括:
a、令l=1;
b、令j=1;
c、计算所述第j高的谱峰位置频率函数qpeak(j,l)与所述MFSK信号载频qc的差qdiff(c,l)=qpeak(j,l)-qc;
d、确定所述拉马努金傅里叶谱图El(q)的判决频率qdec(l),设定频率分辨率qRES=1/I,I为频率插值常数,如果存在一个载频序号c使|qdiff(c,l)|<qRES则qdec(l)=qc;
e、如果任意m都满足|qdiff(m,l)|≥qRES则j=j+1并回到步骤c;
f、根据所述判决频率qdec(l)与数据的映射关系图解调出所述数据信息date(l),l=l+1并回到步骤b。
本发明解决了存在多普勒频移情况下信号准确解调的问题。
附图说明
图1为本发明所述的MFSK解调方法流程图。
图2为本发明实施例中拉马努金傅里叶变换RFT与离散傅里叶变换DFT的性能比较图。
图3为本发明实施例中时间上逐点的谱峰选择算法解调原理图。
图4为本发明实施例中时间上逐点的谱峰选择算法流程图。
图5为本发明实施例中MFSK信号仿真图。
图6为本发明实施例中MFSK信号拉马努金傅里叶谱图仿真图。
图7为本发明实施例中MFSK信号码元中心处拉马努金傅里叶谱图仿真图。
图8为本发明实施例中误码率曲线图。
具体实施方式
以下将参照图1-8对本发明的具体实施方式进行说明。
一种MFSK信号解调方法,如图1所示,该方法包括:
步骤1、对接收到的MFSK信号s(t)过采样得到离散接收信号s(n),其中t为连续的时间,n为离散的时间序号,n=1,2,…,N,对所述离散接收信号s(n)做短时拉马努金傅里叶变换ST-RFT,得到拉马努金傅里叶谱S(m,q),时间m=1,2,…,N,频率q=1,2,…,N,设定所述拉马努金傅里叶谱S(m,q)的拉马努金傅里叶谱图Em(q)=|S(m,q)|2;
进一步地,步骤1中MFSK信号表达式为
其中,Tb为码元周期,g(t)为0到Tb之间幅度为1的门函数,
M个载频设定采样频率fs,则采样周期Ts=1/fs,离散MFSK信号表达式为
短时拉马努金傅里叶变换ST-RFT的表达式为:
为滑动的时间窗,为拉马努金代数和,φ(q)为欧拉函数。
需要说明的是:拉马努金傅里叶变换RFT表达式为
定义拉马努金谱为
P(q)=|S(q)|2
余弦信号的拉马努金傅里叶变换RFT可以提取出余弦信号的频率,设定余弦信号周期为T0,相位为幅度为A,其拉马努金傅里叶变换RFT分两种情况:
当q=T0时
当q≠T0时,Sc(q)趋于0。
因此努金傅里叶变换RFT可以用于信号的频域分析,拉马努金傅里叶变换RFT与离散傅里叶变换DFT在加性高斯白噪声AWGN存在情况下的性能比较如图2所示,图中纵坐标为主峰能量比η,定义为
q3dB(i)为信号的拉马努金谱P(q)中3dB带宽内的频率,i=1,2,…,Nm,总共有Nm个频率。图中曲线显示了在提取信号主频率时,拉马努金傅里叶变换RFT受加性高斯白噪声AWGN影响比离散傅里叶变换DFT小。
步骤2、对所述拉马努金傅里叶谱图Em(q)做时间采样,得到码元中心处拉马努金傅里叶谱图
Tb为所述接收到的MFSK信号码元周期,Ts为采样周期,round[·]为求圆整,l为码元序号,对所述码元中心处拉马努金傅里叶谱图做谱峰搜索,得到第j高的谱峰位置频率函数qpeak(j,l);
需要说明的是:为了保证码元中心处拉马努金傅里叶谱图准确,需要对信号做准确的时间同步,因此本方法需要对MFSK接收信号先做时间同步。
步骤3、设定MFSK载频qc,载频序号c=1,2,…,M,M为MFSK信号阶次,对所述第j高的谱峰位置频率函数qpeak(j,l)与所述MFSK信号载频qc做时间上逐点的谱峰选择算法。
需要说明的是:时间上逐点的谱峰选择算法原理如图3所示,通过比较每个时间l上的第j高的谱峰位置频率函数qpeak(j,l)与MFSK载频qc,确定每个时间上的判决频率qdec(l)。
进一步地,时间上逐点的谱峰选择算法流程图如图4所示,步骤包括:
a、令l=1;
b、令j=1;
c、计算所述第j高的谱峰位置频率函数qpeak(j,l)与所述MFSK载频qc的差qdiff(c,l)=qpeak(j,l)-qc;
d、确定所述拉马努金傅里叶谱图El(q)的判决频率qdec(l),设定频率分辨率qRES=1/I,I为频率插值常数,如果存在一个载频序号c使|qdiff(c,l)|<qRES则qdec(l)=qc;
e、如果任意m都满足|qdiff(m,l)|≥qRES则j=j+1并回到步骤c;
f、根据所述判决频率qdec(l)与数据的映射关系图解调出所述数据信息date(l),l=l+1并回到步骤b。
需要说明的是:所述MFSK载频qc与信号载频fc的对应关系为qm=fs/fm。
下面是一个MFSK信号解调的实例:
设定一个2FSK接收信号,信号载频分别为400MHz、500MHz,码元周期Tb=0.1μs,调制的数据为{1,0,0,1},2FSK接收信号仿真图如图5所示;
对2FSK接收信号进行时域采样,得到离散2FSK信号,采样周期Ts=5×10-5μs,对离散2FSK信号做短时拉马努金傅里叶变换ST-RFT,再计算其拉马努金傅里叶谱图Em(q),如图6所示;
对拉马努金傅里叶谱图做时域采样,得到码元中心处拉马努金傅里叶谱图
如图7所示;
根据时间上逐点的谱峰选择算法得到qdec(l)={50,40,40,50},映射关系图解调出的数据为{1,0,0,1};
做蒙特卡洛实验,得到误码率曲线图,如图8所示,并与基于ST-DFT的解调方法进行了对比。
尽管参考附图详地公开了本发明,但应理解的是,这些描述仅仅是示例性的,并非用来限制本发明的应用。本发明的保护范围由附加权利要求限定,并可包括在不脱离本发明保护范围和精神的情况下针对发明所作的各种变型、改型及等效方案。
Claims (1)
1.一种MFSK信号解调方法,该方法包括:
步骤1、对接收到的MFSK信号s(t)过采样得到离散接收信号s(n),其中t为连续的时间,n为离散的时间序号,n=1,2,…,N,对所述离散接收信号s(n)做短时拉马努金傅里叶变换ST-RFT,得到拉马努金傅里叶谱S(m,q),时间m=1,2,…,N,频率q=1,2,…,N,设定所述拉马努金傅里叶谱S(m,q)的拉马努金傅里叶谱图Em(q)=|S(m,q)|2;MFSK信号表达式为:
其中,Tb为码元周期,g(t)为0到Tb之间幅度为1的门函数,M个载频设定采样频率fs,则采样周期Ts=1/fs,离散MFSK信号表达式为
短时拉马努金傅里叶变换ST-RFT的表达式为:
为滑动的时间窗,为拉马努金代数和,φ(q)为欧拉函数;
步骤2、对所述拉马努金傅里叶谱图Em(q)做时间采样,得到码元中心处拉马努金傅里叶谱图
Tb为所述接收到的MFSK信号码元周期,Ts为采样周期,round[·]为求圆整,l为码元序号,对所述码元中心处拉马努金傅里叶谱图做谱峰搜索,对搜索出的谱峰进行排序,得到第j高的谱峰位置频率函数qpeak(j,l);
步骤3、设定MFSK载频qc,载频序号c=1,2,…,M,M为MFSK信号阶次,对所述第j高的谱峰位置频率函数qpeak(j,l)与所述MFSK载频qc做时间上逐点的谱峰差值选择算法,得到解调出的数据信息date(l);
步骤3中,所述时间上逐点的谱峰差值选择算法包括:
a、令l=1;
b、令j=1;
c、计算所述第j高的谱峰位置频率函数qpeak(j,l)与所述MFSK信号载频qc的差qdiff(c,l)=qpeak(j,l)-qc;
d、确定所述拉马努金傅里叶谱图El(q)的判决频率qdec(l),设定频率分辨率qRES=1/I,I为频率插值常数,如果存在一个载频序号c使|qdiff(c,l)|<qRES则qdec(l)=qc;
e、如果任意m都满足|qdiff(m,l)|≥qRES则j=j+1并回到步骤c;
f、根据所述判决频率qdec(l)与数据的映射关系图解调出所述数据信息date(l),l=l+1并回到步骤b。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610116988.7A CN105743830B (zh) | 2016-03-01 | 2016-03-01 | 一种mfsk信号的解调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610116988.7A CN105743830B (zh) | 2016-03-01 | 2016-03-01 | 一种mfsk信号的解调方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105743830A CN105743830A (zh) | 2016-07-06 |
CN105743830B true CN105743830B (zh) | 2018-09-14 |
Family
ID=56248907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610116988.7A Active CN105743830B (zh) | 2016-03-01 | 2016-03-01 | 一种mfsk信号的解调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105743830B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106789783B (zh) * | 2016-10-18 | 2019-11-08 | 天津理工大学 | 一种利用短时拉马努金傅里叶变换谱图的信号识别方法 |
CN115249486B (zh) * | 2022-07-28 | 2024-04-09 | 哈尔滨工业大学 | 一种旋转机械音异常识别预处理方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103647740A (zh) * | 2013-11-25 | 2014-03-19 | 北京航空航天大学 | 基于拉马努金和的正交非均匀多载波间距的多载波调制解调方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9386385B2 (en) * | 2005-09-27 | 2016-07-05 | Ronald Quan | Method and apparatus to evaluate audio equipment via filter banks for dynamic distortions and or differential phase and frequency modulation effects |
-
2016
- 2016-03-01 CN CN201610116988.7A patent/CN105743830B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103647740A (zh) * | 2013-11-25 | 2014-03-19 | 北京航空航天大学 | 基于拉马努金和的正交非均匀多载波间距的多载波调制解调方法 |
Non-Patent Citations (1)
Title |
---|
Ramanujin Suns and Discrete Fourier Transforms;Saed Samadi等;《IEEE SIGNAL PROCESSING LETTERS》;20050430;第12卷(第4期);第293-296段 * |
Also Published As
Publication number | Publication date |
---|---|
CN105743830A (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | A modified waveform design for radar-communication integration based on LFM-CPM | |
Headley et al. | Asynchronous classification of digital amplitude-phase modulated signals in flat-fading channels | |
KR101462364B1 (ko) | 근거리 자기장 통신 구현 방법 및 시스템 | |
CN103837740A (zh) | 一种高精度数字瞬时测频方法及装置 | |
CN104410590A (zh) | 一种基于压缩感知的短波ofdm抑制干扰联合信道估计方法 | |
US20140376540A1 (en) | Frequency-Domain Symbol And Frame Synchronization In Multi-Carrier Systems | |
US11005525B2 (en) | Methods for improving flexibility and data rate of chirp spread spectrum systems in LoRaWAN | |
JP5740434B2 (ja) | 通信速度およびスペクトル効率を増大させ、他の利益を可能にする共通波通信システムおよび方法ならびに側波帯軽減通信システムおよび方法 | |
CN111935046A (zh) | 一种低复杂度的频移键控信号符号率估计方法 | |
CN105871413B (zh) | 低信噪比直接序列扩频信号检测方法 | |
CN105743830B (zh) | 一种mfsk信号的解调方法 | |
CN105429719A (zh) | 基于功率谱和多尺度小波变换分析强干扰信号检测方法 | |
Macnae | Stripping very low frequency communication signals with minimum shift keying encoding from streamed time-domain electromagnetic data | |
CN103873225B (zh) | 突发通信的定时估计方法 | |
CN110113273B (zh) | 一种基于频谱对称的时频重叠信号盲分离方法 | |
CN106569182B (zh) | 基于最小熵的相位编码信号载频估计方法 | |
CN105656511B (zh) | 一种适应于有频偏和低信噪比环境下的差分相关捕获方法 | |
Cai et al. | Further results on multicarrier MFSK based underwater acoustic communications | |
EP3324590A1 (en) | Signal demodulating device and method applying to closed communication system | |
Kubankova et al. | Algorithms of digital modulation classification and their verification | |
CN110290084A (zh) | 一种基于数据频率能量峰值的短波信道盲符号同步方法 | |
CN107707499B (zh) | 一种Alpha稳定分布噪声下OFDM信号调制参数估计方法 | |
Feng et al. | An accurate modulation recognition method of QPSK signal | |
Molchanov et al. | A novel bispectrum-based frequency encoding technique developed for noisy and multipath fading radio links | |
US7734434B2 (en) | High speed digital waveform identification using higher order statistical signal processing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220311 Address after: No.99 Denglong Road, Fuzhou Economic and Technological Development Zone, Fujian Province 350015 Patentee after: YANGO University Address before: 300384 No. 391 Binshui West Road, Xiqing District, Tianjin Patentee before: TIANJIN University OF TECHNOLOGY |
|
TR01 | Transfer of patent right |