﻿功能:傅里叶级数拟合

格式:
[f,c,a,b] = FourierFit(x,y,n,s)
[f,c,a,b] = FourierFit(x,y,n)

x: 矩阵变量存储的自变量
y: 与x维数对应的因变量
n: 拟合的最高次数,此参数大于0
s: 与周期相关的系数,此参数默认为1

f: 返回拟合后的三角函数表达式
c: 拟合的常系数
a: 拟合的系数
b: 拟合的系数

说明:
1、本函数主要是针对如下关系式进行拟合得到a,b,c
$$
\begin{cases}
s &=\dfrac{2\pi}{t}
\\
y &= c + \sum_{k=1}^n(a_k\cos(ksx)+b_k\sin(ksx))
\end{cases}
$$
上面的t为循环周期

2. 如果不确定周期相关系数s,建议使用{FourierWFit<矩阵运算/FourierWFit>}函数进行拟合

例子:

x=0:001:10;y=4+5*sin(x)+6*cos(2*x);//先产生测试数据

[f,c,a,b] = FourierFit(x,y,3)//这里使用3次拟合系数,执行得到如下结果
f =
{4+5*Sin(x)+6*Cos(2*x) }
c =
[ 4.00000000000000 ]
a =
[ 0.00000000000000
  6.00000000000000
  0.00000000000000 ]
b =
[ 5.00000000000000
  0.00000000000000
  0.00000000000000 ]
