我有等式 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),所以我从上面创建了两个值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
?