1.基于压缩感知、高速移动的UFMC系统无线信道估计方法，其特征在于：该方法包括以下步骤：S1：建立系统模型，分析插入PN序列的通用滤波多载波(Universal Filtered Multicarrier,UFMC)系统参数，确定在多径信道时的发送与接收信号模型；

S2：进行多径信道的时间相关特性及稀疏特性分析；

S3：进行基于PN序列自相关的信道粗估计；

S4：进行基于压缩感知算法的信道时延位置精确估计；

S5：进行最小二乘(least square,LS)算法的信道幅度值的精确估计；

S6：性能与仿真结果分析。

2.根据权利要求1所述的基于压缩感知、高速移动的UFMC系统无线信道估计方法，其特征在于：所述步骤S1包括以下步骤：

S101：在通用滤波多载波系统中，设B个子带，每个子带子载波数量为 总的子载波数量为N，采用切比雪夫滤波器h，滤波器长度为Lf，则在一个符号时间间隔内的基带离散时间信号为：其中，输入数据Xi(m)是独立随机变量的比特流，i是子带编号，l是滤波器长度的编号，m是子带载波的编号，xn是经过UFMC调制后得到时域的通用滤波多载波数据信号，n为离散时间索引；

S102：将PN序列cn插入数据信号部分xn得到一个完整的时域训练序列通用滤波多载波TS-UFMC帧结构sn，输出的TS-UFMC符号 由长度为N+Lf-1数据块和长度为M的PN序列c＝[c0,c1,…,cM]T组成，因此，输出的通用滤波多载波符号表示为：

为频域信号；第i个子带的 个复数正交幅度调制符号乘以离散傅里叶反变换(IDFT)矩阵Vi得到子带的时域数据，1≤i≤B；Fi为托普利茨矩阵，包含了有限长单位冲激响应(Finite Impulse Response,FIR)滤波器的冲击响应，将子带的时域数据与滤波器进行线性卷积，每个子带时域数据滤波后再进行累加得到时域的数据信号x；

S103：对长度为L的多径信道，信道离散脉冲响应CIR： 第i个UFMC符号，信道稀疏度为Si的离散脉冲响应hi,n建模为：

其中，为第i个UFMC符号的第l条径的幅度值，为与采样周期进行归一化后的第l条径的时延，为信道的n时刻的冲击响应：路径时延的集合为：

其中， 为了避免2个相邻的数据块之间的干扰(IBI)，假设L≤M；

发送的第i个TS-UMFC时域符号 经过信道 后，接收到的数据信号 其频域的数据符号为

其中， 为零均值方差为σ2的高斯白噪声AWGN， 为第i个符号的第k个子载波的信道频域响应CFR；

S104：根据公式(6)观测信号 根据公式(7)精确估计信道 的情况，并求出发送信号 的值；由于TS-UFMC帧结构没有导频，通过时域插入的PN序列进行信道估计，接收端的PN序列时域符号 表示为：di＝Ψihi+ηi  (8)

其中，ηi为AWGN噪声向量；

在实际的多径信道中，接收到的PN序列时域符号di将与前一个UFMC数据块的拖尾混叠；在当前的数据块的前面数据中同时混叠了当前PN的拖尾，因此，实际的接收信号rin表示为：

因此，根据实际接收信号rin、消除前一数据块拖尾对PN序列的干扰、消除当前PN序列的拖尾对自身数据块内容造成的干扰，准确估计信道；PN序列重构后的表达式 为：

3.根据权利要求1所述的基于压缩感知、高速移动的UFMC系统无线信道估计方法，其特征在于：所述步骤S2具体为：

对于宽带系统，设系统带宽为W，路径时延宽带的分辨率为 同时对h(τ)按照系统时钟 对时间区间[0,τmax]进行均匀采样，从而得到 个采样点；假设i固定不变，无线信道具有较高能量的路径集合Di为：Di＝{n:|hi,n|＞ε}；其中，ε为路径能量的门限值，小于ε的路径均设置为0；当K＝|Di|＜＜L时，则信道的稀疏度为K；系统带宽B越大，信道稀疏特征越明显。

4.根据权利要求1所述的基于压缩感知、高速移动的UFMC系统无线信道估计方法，其特征在于：所述步骤S3具体为：

在第i个TDS-UFMC符号的保护间隔中插入时域PN序列ci，其归一化相关函数表示为：通过时域相关，即接收到的PN序列di和本地PN序列ci相关，得到信道时域冲激响应的粗略估计同时，根据无线信道的时间相关特性，得到在相干时间Tc内，信道的具有相同的冲激响应CIR，且Tc的大小由接收终端的运动速度和载波频率决定；因此，信道的时延在连续2Rd-1个符号内基本不变，其中 同时，在相干时间Tc内，信道增益表示为|αi,p|exp(φ0+2πfdt)，其中，αi,p为第i个符号的第p条径的路径增益；φ0为初始相位，fd为多普勒频率；根据信道增益的表达式，在时间间隔t＜1/2fd的时间内，即 个连续的符号时间内，复信道增益的相位变化不会超过π；同时，假设信道在一个符号内的时延和信道的幅度均不变，为准静态信道；

因此，在充分利用PN序列扩频增益的基础上，平均Rg个连续符号的CIR来降低噪声对信号的影响，进一步提高信道时域冲激响应的粗略估计 的估计精度；令连续的Rg个接收的PN序列与长度为M本地PN序列c与做循环相关；

由于信道的时延在连续2Rd-1个符号内基本不变，因此在2Rd-1个符号范围内共有2Rd-Rg组 值；因此，得到信道粗略估计的结果为：由于信道存在高斯白噪声和多径，以及多径信道的稀疏特性，经过粗略估计得到的信道的时域冲激响应 中的小电平值将被设为0；那么，路径时延集合为： ζ为信道幅度的门限值，其定义为：其中，Ke为线性调整系数用于调整信道幅度的门限值；Lpre为信道主要时延的长度预设值；

那么，得到信道初始稀疏度S0＝||D0||0；同时，考虑部分时延幅度比较小的路径在粗略估计过程中误认为是系统的噪声，设系统的稀疏度S大于S0为S＝S0+Sa；通过PN序列自相关进行信道粗估计得到了信道的主要时延的位置和幅度值，作为稀疏信道估计的预知信息，以提高稀疏信道估计的精度并降低复杂度。

5.根据权利要求的基于压缩感知、高速移动的UFMC系统无线信道估计方法，其特征在于：所述步骤S4具体为：

设x0为信道幅度初始值，r0为残差初始值，dIBI-free为观测向量，Φ为传感矩阵，F0为主要时延位置初始值，D0为信道主要时延位置，Llast为信道最终稀疏度初始值，S0为初始稀疏度，F为主要时延位置，Sstep为步长， 为信道主要时延位置；

S401：初始化信道的主要路径时延： 为D0的矩阵 的主要时延的列；

S402：初始化残差：

S403：初始化主要时延位置：F0＝D0；

S404：初始化最终的稀疏度：Llast＝S0；

S405：主要时延位置判决迭代次数k＝1；

S406：步长Sstep判决迭代次数j＝1；

重复S407-S412；

S407：选择Llast个最大的max(|ΦHrk-1|,Llast)；

S408：Ck＝Fk-1∪Sk；Sk为第k次max(|ΦHrk-1|,Llast)中最大的Llast个值ΦH对应的序列号所组成的集合；Fk-1为第k-1次中集合Ck中前Llast的值；

S409：F＝Ck(1:Llast)，F为第k次循环集合Ck中前Llast的值；

S410：信道粗估计中最大的Llast径幅度值： 为对应集合F的矩阵 的主要时延的列；

S411：计算残差： ΦF为传感矩阵Φ对应集合F所对应的列；

S412：若停止条件为真，则退出迭代；否则，||r||2＞||rk-1||2时，j＝j+1Llast＝Llast+j×Sstep；||r||2≤||rk-1||2时，Fk＝F，rk＝r，k＝k+1，rk为第k次的残差；

S413：若停止条件为真，则输出

6.根据权利要求1所述的基于压缩感知、高速移动的UFMC系统无线信道估计方法，其特征在于：所述步骤S5具体为：

根据信道主要时延位置 得到信道估计 的计算公式：其中，ΦD为矩阵Φ所对应的主要时延的列， 为信道的主要Llast个时延的幅度值，Dc范围的信道的CIR幅度值为0，即 根据公式(18)，采用LS算法得：

7.根据权利要求1所述的基于压缩感知、高速移动的UFMC系统无线信道估计方法，其特征在于：所述步骤S6具体为：

S601：信道估计的误差性能下界CRLB分析

根据得到的AWGN噪声向量ηi，其分布函数为 求出 的在 条件概率密度函数PDF为：

根据向量估计理论，无偏估计 公式(19)的CRLB推导得：当采用LS等常规信道估计算法时，假设完全消除干扰的情况下，得到的最好的均方误差MSE为σ2；当Llast＜LIBI-free时，采用基于PPI-SAMP的信道估计算法得到的MSE值将会更小；

S602：复杂度分析

1)采用基于PN序列自相关进行信道粗估计，将接收到的信号ri和本地PN序列ci相关，对每个TS-UFMC符号，其计算复杂度为O(LIBI-free)；

2)采用PPI-SAMP算法进行信道时延位置精确估计，此步骤计算复杂度为：O(NPPI-SAMP×LIBI-free×L)，其中NPPI-SAMP为余下的稀疏度b所需要循环的次数，L为信道长度；

2

3)采用LS算法精确估计信道幅度；此步的计算复杂度为O(LIBI-free(Llast))；

4)根据分析知，PPI-SAMP算法的总复杂度为O(LIBI-free×(NPPI-SAMP×L+(Llast)2+1))。