关于分块kalman 滤波

[ydaj]

请问博主，在频域分块kalman算法里，更新滤波器系数时，对于constrain matrix GL您在代码里好像没有体现出来啊，另外，为啥计算X(k)P(k)XH(k)时候先计算X2之后要按行求和啊，X2 = np.sum(np.abs(self.X) ** 2, axis=0)
Re = 0.5 * self.Rm * X2 + np.abs(E) ** 2 / self.N

[nicheng0019]

请问博主，在频域分块kalman算法里，更新滤波器系数时，对于constrain matrix GL您在代码里好像没有体现出来啊，另外，为啥计算X(k)P(k)XH(k)时候先计算X2之后要按行求和啊，X2 = np.sum(np.abs(self.X) ** 2, axis=0) Re = 0.5 * self.Rm * X2 + np.abs(E) ** 2 / self.N

请问频域分块kalman算法在计算mu = self.P / (Pe + 1e-10)的时候为什么没有乘以0.5？这个代码实现参照的论文是什么？谢谢

[Fei Jiang]

@ydaj constrain matrix的实现在这呢

pyaec/frequency_domain_adaptive_filters/pfdkf.py

Lines 61 to 70 in 5b9c02c

 if self.partial_constrain: 
   h = ifft(self.H[self.p]) 
   h[self.M:] = 0 
   self.H[self.p] = fft(h) 
   self.p = (self.p + 1) % self.N 
 else: 
   for p in range(self.N): 
     h = ifft(self.H[p]) 
     h[self.M:] = 0 
     self.H[p] = fft(h)