[求助]数据处理的问题

最近用mathematica做一些引力波的简单计算，遇到了一个头疼的问题。
根据某个双星的运动计算其引力波时，我们可以通过mathematica把引力波的两个模式hxx和hxy以插值函数的形式解出来，现在问题是怎么分析插值函数的性质。
代码如下：
Clear["Global`*"]
a = 1.9943`*10^9;
e = 0.617; M = m = 2.07*10^3; \[Mu] = (M m)/(M + m); rd =
2.4*10^20; P = 8.372*10^12; \[Omega] = (2 \[Pi])/P;
r[t_] := (a (1 - e^2))/(1 + e Cos[\[Theta][t]])
sol = NDSolve[{\[Theta]'[t] == (a (1 - e^2) (M + m))^(1/2)/
r[t]^2, \[Theta][0] == 0}, \[Theta][t], {t, 0, 4 P}]
\[Theta]1[t_] := \[Theta][t] /. sol[[1]];
\[Theta]1[t];
r1[t_] := (a (1 - e^2))/(1 + e Cos[\[Theta]1[t]]);
Ixx[t_] := \[Mu] r1[t]^2 (Cos[\[Theta]1[t]]^2 - 1/3);
Iyy[t_] := \[Mu] r1[t]^2 (Sin[\[Theta]1[t]]^2 - 1/3);
Ixy[t_] := \[Mu] r1[t]^2 Sin[\[Theta]1[t]] Cos[\[Theta]1[t]];
hxx[t_] := -(1/rd) (D[Ixx[t], {t, 2}] - D[Iyy[t], {t, 2}]);
hxy[t_] := 2/rd D[Ixy[t], {t, 2}];
Plot[{Evaluate[hxx[t]], Evaluate[hxy[t]]}, {t, 0, 2 P},
PlotPoints -> 500, PlotRange -> All]
计算出插值函数并作图后，想进一步分析插值函数的性质，比如其傅里叶系数以及模型拟合时都遇到了很大的问题，比如FindFit计算出来的参数作图后跟原图像差别非常大，请问大大们该怎么处理呢？