标签:
Methods
No linear relationship was found between NO2 tropospheric column and NO3-, but both well simulated by y=y0+a*sin(b*x+c)+e*x indicating NO2 tropospheric column could reflect the trend of NO3- but was difficult to quantify NO3-.
data,
(1) User defined function: y=y0+a*sin(b*x+c)+e*x in Sigmaplot
First, generate scatter
Second, define the function "y=y0+a*sin(b*x+c)+e*x"
Equation:
pi=3.14159265358979
f = y0+a*sin(2*pi*x/b+c)
fit f to y
‘‘fit f to y with weight reciprocal_y
‘‘fit f to y with weight reciprocal_ysquare
‘‘fit f to y with weight reciprocal_pred
‘‘fit f to y with weight reciprocal_predsqr
Vavriables:
x = col(1) ‘ {{prevmin: 0.000000}} {{prevmax: 5.000000}} {{intervals: 100}}
y = col(2)
reciprocal_y = 1/abs(y)
reciprocal_ysquare = 1/y^2
reciprocal_pred = 1/abs(f)
reciprocal_predsqr = 1/f^2
‘Automatic Initial Parameter Estimate Functions
F(q) = sinp(x,y)
Initial parameters:
a = F(0)[1] ‘‘Auto {{previous: 95}}
b = if(F(0)[2]<>0, F(0)[2], 1) ‘‘Auto {{previous: 1.8}}
c = F(0)[3] ‘‘Auto {{previous: .5}}
y0 = mean(y) ‘‘Auto {{previous: 5}}
Constraits
a>0
b>0
c>-6.28
c<6.28
Third, gernrate the curve
Additional procedure:
Multiple lines:
Linear trend, then, right click, Curve fit, Line/Scatter
(2) Made "y=y0+a*sin(b*x+c)+e*x" in Matlab
cftool -v1
标签:
原文地址:http://www.cnblogs.com/rockman/p/4860092.html