本文整理汇总了Python中vcmq.N.sqrt方法的典型用法代码示例。如果您正苦于以下问题:Python N.sqrt方法的具体用法?Python N.sqrt怎么用?Python N.sqrt使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类vcmq.N
的用法示例。
在下文中一共展示了N.sqrt方法的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: round_date
# 需要导入模块: from vcmq import N [as 别名]
# 或者: from vcmq.N import sqrt [as 别名]
f.close()
# Create hourly time axis
taxi = sp.getTime()
taxi.toRelativeTime("hours since 2000")
ctimesi = taxi.asComponentTime()
ct0 = round_date(ctimesi[0], "hour")
ct1 = round_date(ctimesi[-1], "hour")
taxo = create_time(lindates(ct0, ct1, 1, "hour"), taxi.units)
# Lag error
# - estimation
els = []
lags = N.arange(1, 6)
for lag in lags:
els.append(N.sqrt(((sp[lag:] - sp[:-lag]) ** 2).mean()))
els = N.array(els)
a, b, _, _, _ = linregress(lags, els)
# - plot
P.figure(figsize=(6, 6))
P.subplot(211)
P.plot(lags, els, "o")
P.plot([0, lags[-1]], [b, a * lags[-1] + b], "g")
P.axhline(b, color="0.8", ls="--")
P.ylim(ymin=0)
P.xlabel("Lag [hour]")
P.ylabel("Error [m s-1]")
add_key(1)
P.title("Linear lag error model")
# Interpolation
示例2: cellerr1d
# 需要导入模块: from vcmq import N [as 别名]
# 或者: from vcmq.N import sqrt [as 别名]
# - data
yi = N.array([-2., 2., 6. , 12., 15.])
vari = N.ma.array(yi)
vari[-1] = N.ma.masked
errm = N.ma.array([1, 2, 1., 2., 3.])
yo = N.array([5., 15])
mv = 1e20
vari = vari.reshape(1, -1).filled(mv)
errm = errm.reshape(1, -1).filled(mv)
errl = N.ones(1)
# - interp
varo, erro = cellerr1d(vari, yi, yo, mv, errm, errl)
# - truth
errot = N.array([
1/N.sqrt(
1/(errm[0, 1]**2+N.abs(yo[0]-yi[1])*errl[0]) +
1/(errm[0, 2]**2+N.abs(yo[0]-yi[2])*errl[0])),
N.sqrt(errm[0, 3]**2+N.abs(yo[1]-yi[3])*errl[0])
])
varot = N.array([
vari[0, 1]/(errm[0, 1]**2+N.abs(yo[0]-yi[1])*errl[0]) +
vari[0, 2]/(errm[0, 2]**2+N.abs(yo[0]-yi[2])*errl[0]),
vari[0, 3]/(errm[0, 3]**2+N.abs(yo[1]-yi[3])*errl[0])
])*errot**2
# - check
result = [(N.testing.assert_allclose, [erro.ravel(), errot]),
(N.testing.assert_allclose, [varo.ravel(), varot])]
# 1dx
yi = yi.reshape(1, -1)
varox, errox = cellerr1d(vari, yi, yo, mv, errm, errl)