我在100kHz和200MHz的范围内具有非均匀采样的传递函数数据(仅幅度)。使用此传递函数,我想计算特定输入的输出。
问题:
为了从幅度数据中获得传递函数的相位,我试图在MATLAB中使用希尔伯特变换。
(1)阶段= - imag(希尔伯特(log(TF_magnitude)));
上面的代码是对的吗?我找到了以下计算最小相位的公式。
(2)θ(ω)= - H [ln(G(ω))]
我不知道为什么(1)使用图像功能。
传递函数的采样周期不是线性的,而是对数尺度的均匀。为了使这个线性,我试图使用三次样条插值。但是,传递函数没有任何规律性。
因此,我试图提取与传递函数数据相对应的输入数据,计算输出,然后为线性间隔输出点插值输出。如果我以这种方式工作,这有关系吗?另外,我应该使输出(或传递函数)线性化以在MATLAB中执行IFFT吗?
由于我英语很差,我不确定我的想法是否传达了......但我希望如此。 任何帮助都会很感激。