方程式 ln(c)=-1/2k^2 * z^2 があります。y = ln(c)、x = z^2、a = -1/2k^2 です。
を推定したいaので、次のようにします。
a = polyfit(z.^2, log(abs(c)), 1)
(初期) 方程式 c = exp(-z^2/2k^2) があるため、上から の 2 つの値をa見つけて、k (k1) を推定したいので、次のようにします。
k1 = sqrt(-1/2*a(1))
ここで、k1 と z の値を使用して c とエラーを予測したいと思います。私もです：
c_predict = polyval(a,z)
c1 = exp((-z.^2)/2*k1^2)
error = c_predict - c1
あるいは単に：
c1 = exp((-z.^2)./2*s1^2)
error1 = c - c1
何が正しいですか？
error = c_predict - c1
また
error = c - c1
?