拟合零 a0 的傅立叶模型

Question

以下代码查找傅里叶级数中的系数

load enso;
f = fit(month,pressure,'fourier2');

代码输出：

 General model Fourier1:
 f(x) =  a0 + a1*cos(x*w) + b1*sin(x*w)
 Coefficients (with 95% confidence bounds):
   a0 =       10.63  (10.22, 11.03)
   a1 =       2.876  (2.187, 3.565)
   b1 =       1.177  (0.09191, 2.261)
   w =      0.5263  (0.5225, 0.5301)

我需要找到函数 f(x) = a1*cos(x*w) + b1*sin(x*w) 的系数，所以 a0 应该是 0 或者根本不应该在公式中。如何做到这一点？

Answer 1

你想要'fourier1' 还是'fourier2'？因为你的问题都有。您可以使用拟合选项将 a0 限制为 0。以下是两者的匹配项，-Inf 到 Inf 是 a、b 和 w 的范围，而 0 到 0 是限制 a0。

load enso;

ft1 = fittype('fourier1');

options1 = fitoptions(ft1);
options1.Lower = [0 -Inf -Inf -Inf];
options1.Upper = [0  Inf  Inf  Inf];

f1 = fit(month,pressure,ft1,options1)

ft2 = fittype('fourier2');

options2 = fitoptions(ft2);
options2.Lower = [0 -Inf -Inf -Inf  -Inf  -Inf];
options2.Upper = [0  Inf  Inf  Inf   Inf   Inf];

f2 = fit(month,pressure,ft2,options2)

编辑：

对我来说运行良好。这很奇怪……

f1 = 

     General model Fourier1:
     f1(x) =  a0 + a1*cos(x*w) + b1*sin(x*w)
     Coefficients (with 95% confidence bounds):
       a0 =           0  (fixed at bound)
       a1 =       2.258  (-1.631, 6.148)
       b1 =       2.406  (-1.317, 6.13)
       w =      0.5311  (0.516, 0.5462)