python numpy用法_Python numpy.polyder() 使用实例

The following are code examples for showing how to use . They are extracted from open source Python projects. You can vote up the examples you like or vote down the exmaples you don’t like. You can also save this page to your account.

Example 1

def calc_approx(time, altitude, div = 5, dt = 1.0, deg = 100, lim = 1):

time_list, altitude_list, velocity_y = [],[],[]

for i in range(div):

f = trendline.approx_func(

time[max(int((i - lim) * len(time) / div), 0):min(len(time), int((i + lim) * len(time) / div))],

altitude[max(int((i - lim) * len(time) / div), 0):min(len(time), int((i + lim) * len(time) / div))], deg)

der = np.polyder(f)

#print(time[int(i * len(time) / div)], min(time[-1], time[int((i + 1) * len(time) / div) - 1]))

for t in np.arange(time[int(i * len(time) / div)], min(time[-1], time[int((i + 1) * len(time) / div) - 1]), dt):

#for t in time:

# if t >= time[int(i * len(time) / div)] and t < min(time[-1], time[int((i + 1) * len(time) / div) - 1]):

time_list.append(t)

altitude_list.append(f(t))

velocity_y.append(der(t))

return time_list, velocity_y, altitude_list

Example 2

def calc_adot(cp,order=1):

a=[]

for i in range(len(cp)):

ptmp = []

tmp = 0

for j in range(len(cp)):

if j != i:

row = []

row.insert(0,1/(cp[i]-cp[j]))

row.insert(1,-cp[j]/(cp[i]-cp[j]))

ptmp.insert(tmp,row)

tmp += 1

p=[1]

for j in range(len(cp)-1):

p = conv(p,ptmp[j])

pder = numpy.polyder(p,order)

arow = []

for j in range(len(cp)):

arow.append(numpy.polyval(pder,cp[j]))

a.append(arow)

return a

Example 3

def t