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；

步骤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)。

2.如权利要求1所述的方法，其特征在于，步骤1中：

短时拉马努金傅里叶变换ST-RFT的表达式为：

为滑动的时间窗， 为拉马努金代数和，φ(q)为欧

拉函数。

3.如权利要求1所述的方法，其特征在于，步骤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)|

e、如果任意m都满足|qdiff(m,l)|≥qRES则j＝j+1并回到步骤c；

f、根据所述判决频率qdec(l)与数据的映射关系图解调出所述数据信息date(l)，l＝l+

1并回到步骤b。