def Tmtx_ri(b_ri, K, D, L):
"""
build convolution matrix associated with b_ri
:param b_ri: a real-valued vector
:param K: number of Diracs
:param D1: expansion matrix for the real-part
:param D2: expansion matrix for the imaginary-part
:return:
"""
b_ri = np.dot(D, b_ri)
b_r = b_ri[:L]
b_i = b_ri[L:]
Tb_r = linalg.toeplitz(b_r[K:], b_r[K::-1])
Tb_i = linalg.toeplitz(b_i[K:], b_i[K::-1])
return np.vstack((np.hstack((Tb_r, -Tb_i)), np.hstack((Tb_i, Tb_r))))
评论列表
文章目录