:Ҷ

ʽ:
[er,w,a0,a,b,l]=FourierFitGA(x,y,n,w0,wMax,wMin,Nbird,LoopSum,Error)
[er,w,a0,a,b,l]=FourierFitGA(x,y,n,w0,wMax,wMin,Nbird,LoopSum)
[er,w,a0,a,b,l]=FourierFitGA(x,y,n,w0,wMax,wMin,Nbird)
[er,w,a0,a,b,l]=FourierFitGA(x,y,n,w0,wMax,wMin)
[er,w,a0,a,b,l]=FourierFitGA(x,y,n,w0,wMax)
[er,w,a0,a,b,l]=FourierFitGA(x,y,n,w0)
[er,w,a0,a,b,l]=FourierFitGA(x,y,n)
[er,w,a0,a,b,l]=FourierFitGA(x,y)

x:Ҫϵľ
y:xͬάҪϵľ
n:ϱ,ĬΪ1.庬忴˵.
w0:,ĬΪ6.28318530717958
wMax:ݵ,ĬΪ12.5663706143592
wMin:ݵ,ĬΪ0
Nbird:Ⱥ,Ĭ100
LoopSum:ѭ,Ĭ3000
Error:Ʊ,Ĭ1e-12

er:ص
w:ص
a0:õĳϵ
a:õϵ
b:õϵ
l:سֹʱѭ

˵:

1Ҫ±ʽ,֪(x,y)ϳw,a0,a,b
$$
y = a_0 + \sum_{i=1}^n (a_i\cos(iwx)+ b_i\sin(iwx))
$$
2ڲ{nlinFit2<\nlinFit2>}{SolveEquationGauss<\SolveEquationGauss>}

3עⱾ{FourierFit<\FourierFit>}(δ֪)

4Ϊڲ⴦,ʱֱӵصķϺҪ!

:

x=0:0.2:5;y=4-5.6*sin(3.2*x)+3*cos(3.2*x);//سִб,ͼ֪x,y

[er,w,a0,a,b,l]=FourierFitGA(x,y,1,10)//wʼֵΪ10,(עĬΪ12.5663706143592)سõ½
//õ y=3.99999997420387+3.00000035327529*sin(3.20000002548315*x)-5.59999981846264*cos(3.20000002548315*x)
er =
[ 8.4483621288E-13 ]
w =
[ 3.20000002548315 ]
a0 =
[ 3.99999997420387 ]
a =
[ 3.00000035327529 ]
b =
[-5.59999981846264 ]
l =
[ 10.0000000000000 ]

