:Ŵ㷨Ȼָ

ʽ:
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0,bMax,bMin,cMax,cMin,Nbird,LoopSum,Error)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0,bMax,bMin,cMax,cMin,Nbird,LoopSum)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0,bMax,bMin,cMax,cMin,Nbird)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0,bMax,bMin,cMax,cMin)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0,bMax,bMin,cMax)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0,bMax,bMin)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0,bMax)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0,c0)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n,b0)
[y0,a,b,c,er,m]=ExpFitGA(x,y,n)
[y0,a,b,c,er,m]=ExpFitGA(x,y)

x:ҪϵԱ()
y:Ҫϵ()
n:Ҫϵĺ,ĬΪ1,庬忴˵
b0:ݻ߾,ʾʽbĳʼֵ,Ϊ,ԪظΪn
c0:ݻ߾,ʾʽcĳʼֵ,Ϊ,ԪظΪn
bMax:ݻ߾,ʾʽbĵֵ,ĬΪ100,Ϊ,ԪظΪn
bMin:ݻ߾,ʾʽbĵֵ,ĬΪ-100,Ϊ,ԪظΪn
cMax:ݻ߾,ʾʽcĵֵ,ĬΪ100,Ϊ,ԪظΪn
cMin:ݻ߾,ʾʽcĵֵ,ĬΪ-100,Ϊ,ԪظΪn
Nbird:Ⱥ,Ĭ200
LoopSum:,Ĭ3000
Error:Ʊ,Ĭ1E-16

y0:صϵy0ֵ
a:صϵaֵ
b:صϵbֵ
c:صϵcֵ
er:صľ
m:سֹĵ

˵:

1Ҫһຯ,ʽеnӦʽĲn
$$
y = y_0 +\sum_{i=1}^na_i\exp(b_ix+c_ix^2)
$$
2Ϊڲ㷨⴦,ʱֱӵصķϺҪ!

3עڱʽ exp(b*x+c*x^2+d) = exp(b*x+c*x^2)*exp(d) = exp(b*x+c*x^2)*a , a=exp(d)

4ֻ y = y0 + a*exp(b*x) ֱʽ,ǿԽcMax,cMinΪ0

:

//Ϊ˱ʾ,Ȳһ

x=0:0.1:5; y=5+6*exp(4*x-8*x^2);//سһ(x,y),y0=5,a=6,b=4,c=-8

[y0,a,b,c,er,m]=ExpFitGA(x,y)//Ϊֻһ,˿ֱӲĬϷʽسõ½,ȻҲӦķΧ
y0 =
[ 5.00000000017300 ]
a =
[ 5.99999999508852 ]
b =
[ 4.00000000660456 ]
c =
[-8.00000000997583 ]
er =
[ 4.9607443875E-17 ]
m =
[ 14.0000000000000 ]

//By 2015/2/2