如何计算非均匀采样传递函数数据的输出?

时间:2016-01-22 06:42:19

标签: matlab fft

我在100kHz和200MHz的范围内具有非均匀采样的传递函数数据(仅幅度)。使用此传递函数,我想计算特定输入的输出。

问题:

  1. 为了从幅度数据中获得传递函数的相位,我试图在MATLAB中使用希尔伯特变换。

    (1)阶段= - imag(希尔伯特(log(TF_magnitude)));

    上面的代码是对的吗?我找到了以下计算最小相位的公式。

    (2)θ(ω)= - H [ln(G(ω))]

    我不知道为什么(1)使用图像功能。

  2. 我没有100kHz以下的传递函数数据,尤其是0Hz。在这种情况下,我应该在100kHz以下推断数据吗?或插入0?插入平均值?这会对产量产生很大影响吗?
  3. 传递函数的采样周期不是线性的,而是对数尺度的均匀。为了使这个线性,我试图使用三次样条插值。但是,传递函数没有任何规律性。

    因此,我试图提取与传递函数数据相对应的输入数据,计算输出,然后为线性间隔输出点插值输出。如果我以这种方式工作,这有关系吗?另外,我应该使输出(或传递函数)线性化以在MATLAB中执行IFFT吗?

  4. 由于我英语很差,我不确定我的想法是否传达了......但我希望如此。 任何帮助都会很感激。

0 个答案:

没有答案