Question

我正在使用doatools.py库（https://github.com/morriswmz/doatools.py）现在，我的代码如下：

import numpy as np
from scipy import constants as const
import math
import doatools.model as model
import doatools.estimation as estimation



def calculate_wavelength(frequency):
    return const.speed_of_light / frequency


# Uniform circular array
#        X
#        |
#   X---------X
#        |
#        X
NUMBER_OF_ELEMENTS = 4  # elements are shown as "X"
RADIUS = 0.47 / 2
FREQ_MHZ = 315
freq = FREQ_MHZ * const.mega
wavelength = calculate_wavelength(freq)

antenna_array = model.UniformCircularArray(NUMBER_OF_ELEMENTS, RADIUS)

# Create a MUSIC-based estimator.
grid = estimation.FarField1DSearchGrid()
estimator = estimation.MUSIC(antenna_array, wavelength, grid)


R = np.array([[1.5, 2, 3, 4], [4, 5, 6, 5], [45, 5, 5, 6], [5, 1, 0, 5]])

_, estimates = estimator.estimate(R, 1, return_spectrum=False, refine_estimates=True)
print('Estimates: {0}'.format(estimates.locations))

我可以使用此库生成信号，但是如何使用我自己的库？例如，来自ADC的信号（例如：

-> Switching to antenna 0 : [0, 4, 7, 10]
-> Switching to antenna 1 : [5, 6, 11, 83] 
-> Switching to antenna 2 : [0, 23, 2, 34]
-> Switching to antenna 3 : [23, 105, 98, 200]

）

Answer 1

我认为您的问题是如何从天线馈送真实数据，对吗？假设您的数据应该随时间变化。我的意思是在“天线0：[0、4、7、10]”的情况下，0是第一个数据，顺序是4、7，而10是时间中的最后一个。如果是，则可以将它们保留为简单的矩阵，就像上面键入的内容一样：

r = matrix 4x4 of
  0, 4, 7, 10
  5, 6, 11, 83 
  0, 23, 2, 34
  23, 105, 98, 200

  //===============

  r(0,0) = 0,   r(0,1) = 4,  r(0,2) = 7,   r(0,3) = 10

  r(1,0) = 5,   r(1,1) = 6,  ... etc.

  r(2,0) = 0,   ...etc.

  //==============

R = r及其厄米矩阵的乘积（在python中为r.h）。

R = r @ r.h

这是您需要在函数中作为第一个参数填写的协方差矩阵。

DOATools.py-使用我自己的信号源（未生成）

1 个答案: