[发明专利]信号重建方法及基于该方法的土壤湿度干涉测量方法

权利要求书

1.一种信号重建方法，其特征在于：所述的方法包括以下步骤：

步骤1：原始数据采集，

使用接收机以固定采样间隔采集原始数据；

步骤2：卫星仰角、方位角计算，

利用步骤1中卫星星历信息结合采样时间与接收天线位置，计算卫星的仰角θ(n)与方位角，其中n为采样点；

步骤3：卫星数据筛选，

根据步骤2中卫星的仰角和方位角，结合天线架设高度H，计算镜面反射点位置，筛选低仰角以及镜面反射点落在需测量的土壤区域内的卫星的信噪比数据；

步骤4：刻度转换，

对步骤3筛选出的卫星的信噪比数据进行刻度转换，从对数刻度转化为线性刻度，此时信噪比数据用下式描述：

其中S_d(θ(n))为直射信号信噪比，S_r(θ(n))为反射信号信噪比，二者构成了信噪比数据的趋势项；f(n)为信噪比的震荡频率，与天线相对于有效反射面的高度H有关；为信噪比数据的初始相位；

步骤5：中心化，

利用低阶多项式对信噪比数据进行拟合，得到信噪比数据的趋势项P₀(θ(n))，并将其从信噪比数据中提剔除得到以0为中心变化的信噪比数据SNR_c(θ(n))，用下式表示：

其中P₁(θ(n))为由低阶多项式的建模误差导致的残余项；P₂(θ(n))为信噪比的时变震荡幅度；f(n)为信噪比的时变振荡频率；为信噪比的初始相位；

步骤6：时频分析，

对中心化后的信噪比数据进行时频分析得到信噪比的时变频率f(n)；

步骤7：利用自适应滤波器估计信噪比的残余项与时变震荡幅度，

采用自适应滤波器进行估计，对式(2)进行改写：

其中[1 cos(2πf(n)sin(θ(n))) -sin(2πf(n)sin(θ(n)))]^T为自适应滤波器的输入信号向量也即待滤波的信号，利用步骤2计算出的卫星仰角θ(n)与步骤6得到的时变频率f(n)生成；c₀(n),c₁(n),c₂(n)为滤波器系数；SNR_f(θ(n))则为滤波器输出；滤波器的自适应调整算法不断调整滤波器系数，使得滤波器输出SNR_f(θ(n))与SNR_c(θ(n))的均方根误差最小，此时滤波器系数c₀(n),c₁(n),c₂(n)分别为对P₁(θ(n)),的估计；

步骤8：直射信号与反射信号重建，

利用趋势项P₀(θ(n))以及步骤7得到的滤波系数构建关于直射信号与反射信号的方程组，如下：

其中S_d(θ(n)),S_r(θ(n))分别为直射信号与反射信号的信噪比。求解该方程组得到：

。

2.根据权利要求1所述的信号重建方法，其特征在于：所述的步骤1中的接收机为GNSS-IR接收机。

3.根据权利要求1所述的信号重建方法，其特征在于：所述的步骤1中的原始数据包括各个卫星的信噪比观测值与星历信息。