Matlab 用失真的时基拟合噪声正弦波

我最近问了一个关于用扭曲的时基拟合正弦波的问题，在下面的链接中得到了一些很好的答案：

但是，我想执行类似的拟合，但现在信号有噪声。我假设我必须在-1和1之间进行一些归一化，并消除缓慢漂移

原问题如下:

我想知道在Matlab中拟合具有扭曲时基的正弦波的最佳方法

时间上的失真由一个n阶多项式（n~10）给出，其形式为

t\u扭曲=p（t）

例如，考虑失真<代码> Ty-扭曲＝8 +12t+6t^ 2 +t^ 3 < /代码>（这只是<代码>（t2）^ 3 <代码>的幂级数展开）。 这将使正弦波失真，如下所示：

我希望能够找到这个扭曲的正弦波的失真。（也就是说，我想找到函数

t=G（t_扭曲）

，但是

t_扭曲=P（t）

是未知的。）

现在，当失真的正弦波在-1和1之间完全正常化时（除了时基之外没有其他失真），我不想找到失真的正弦波的失真，它是有噪声的，也有漂移。（注意，我根本不关心这个漂移，最好是将其移除）。这样的信号如下：

---

你想要失真还是正弦波？最终是正弦波，但我很确定我可以从失真中得到。失真发生的频率高于正弦波。只需应用高通滤波器，并降低截止频率，直到更高的频率消失。