1.一种基于压缩感知的信号恢复方法,其特征在于,所述方法包括:
对满足稀疏性的卫星上行连续时间信号进行次奈奎斯特采样,以获得采样观测数据;
基于所述次奈奎斯特采样的采样方式确定采样矩阵;
根据所述采样矩阵利用压缩感知理论从所述采样观测数据中恢复出原始信号。
2.如权利要求1所述的方法,其特征在于,所述次奈奎斯特采样的采样方式为多陪集采样;
相应地,所述对满足稀疏性的卫星上行连续时间信号进行次奈奎斯特采样,以获得采样观测数据,进一步包括:
对满足稀疏性的卫星上行连续时间信号按照预设时间周期T进行奈奎斯特采样,以获得离散数据;
对所述离散数据进行多陪集采样,以获得所述采样观察数据。
3.如权利要求2所述的方法,其特征在于,所述采样矩阵A中的第i行第k列的元素Aik通过下式确定,
其中,L为多陪集采样分组时的数组长度,ci为常数集合C中的第i个参数,所述参数用于反映其在所述常数集合C中的位置。
4.如权利要求1~3中任一项所述的方法,其特征在于,所述根据所述采样矩阵利用压缩感知理论从所述采样观测数据中恢复出原始信号,进一步包括:
A101:利用正交性,通过下式寻找最优的指标λv,
其中,*为共轭转置符,Aj为采样矩阵A中的第j列元素组成的向量,rv-1为上一次计算的数据残差,d为采样矩阵A中元素的列数;
A102:根据所述指标λv通过下式对数据矩阵进行处理,
其中,为本次计算的数据矩阵,为上一次计算的数据矩阵,为采样矩阵A中的第λt列元素组成的向量,为将添加至中的计算符;
A103:利用最小二乘法,通过下式求解最优化问题,以获得本次计算的恢复信号
其中,y为采样观测数据;
A104:通过下式计算本次的数据残差,
rv=y-yv,
其中,yv为观测近似值;
A105:对v进行自增运算,判断v是否达到预设值,若否,则返回步骤A101,否则将本次计算的恢复信号作为恢复的原始信号。
5.如权利要求2所述的方法,其特征在于,所述卫星上行连续时间信号满足以下3个条件时满足稀疏性:
(1)所述卫星上行连续时间信号x(t)为带限信号,其频率区间为[0,1/T];
(2)所述卫星上行连续时间信号x(t)在频域上具有有限个频段,且各频段之间不相交;
(3)所述卫星上行连续时间信号x(t)在频域上所有频段的带宽均为有限宽度,且小于1/T。
6.一种基于压缩感知的信号恢复装置,其特征在于,所述装置包括:
信号采样单元,用于对满足稀疏性的卫星上行连续时间信号进行次奈奎斯特采样,以获得采样观测数据;
矩阵确定单元,用于基于所述次奈奎斯特采样的采样方式确定采样矩阵;
信号恢复单元,用于根据所述采样矩阵利用压缩感知理论从所述采样观测数据中恢复出原始信号。
7.如权利要求6所述的装置,其特征在于,所述次奈奎斯特采样的采样方式为多陪集采样;
相应地,所述信号采样单元,进一步用于对满足稀疏性的卫星上行连续时间信号按照预设时间周期T进行奈奎斯特采样,以获得离散数据;对所述离散数据进行多陪集采样,以获得所述采样观察数据。
8.如权利要求7所述的装置,其特征在于,所述采样矩阵A中的第i行第k列的元素Aik通过下式确定,
其中,L为多陪集采样分组时的数组长度,ci为常数集合C中的第i个参数,所述参数用于反映其在所述常数集合C中的位置。
9.如权利要求6~8中任一项所述的装置,其特征在于,所述信号恢复单元,进一步包括:
指标查找子单元,用于利用正交性,通过下式寻找最优的指标λv,
其中,*为共轭转置符,Aj为采样矩阵A中的第j列元素组成的向量,rv-1为上一次计算的数据残差,d为采样矩阵A中元素的列数;
矩阵处理子单元,用于根据所述指标λv通过下式对数据矩阵进行处理,
其中,为本次计算的数据矩阵,为上一次计算的数据矩阵,为采样矩阵A中的第λt列元素组成的向量,为将添加至中的计算符;
求解恢复子单元,用于利用最小二乘法,通过下式求解最优化问题,以获得本次计算的恢复信号
其中,y为采样观测数据;
残差计算子单元,用于通过下式计算本次的数据残差,
rv=y-yv,
其中,yv为观测近似值;
自增判断子单元,用于对v进行自增运算,判断v是否达到预设值,若否,则调用所述指标查找子单元,否则将本次计算的恢复信号作为恢复的原始信号。
10.如权利要求7所述的装置,其特征在于,所述卫星上行连续时间信号满足以下3个条件时满足稀疏性:
(1)所述卫星上行连续时间信号x(t)为带限信号,其频率区间为[0,1/T];
(2)所述卫星上行连续时间信号x(t)在频域上具有有限个频段,且各频段之间不相交;
(3)所述卫星上行连续时间信号x(t)在频域上所有频段的带宽均为有限宽度,且小于1/T。