python如何使用三角函数_python中三角函数的积分

这可能是您想要的：

{1美元^import numpy as np

import matplotlib.pyplot as plt

from scipy import integrate

# step Function

def jerkstep(xprime,k):

return (1 - 2*( int(xprime)%2 == 1 ) )*k

# Triangle function ( step integration)

def jerkint(x,k,ac):

return integrate.quad( jerkstep ,0, x, args =(k) )[0] + ac*x

# Integrate of triangle

def jerkintint(x,k,ac,vc):

return integrate.quad( jerkint, 0 , x , args = (1,0) )[0] + vc*x

# linearly spaced time samples (~ x axis)

t = np.linspace(2.0, 10.0, num = 100)

#optional arguments

k = 1

ac = 0

vc = 0

#output functions

jerk = np.array(t);

triangle = np.array(t);

curve = np.array(t);

# Warning : it can be long ( integration not optimized )

for indice,time in zip(range(len(t)),t): # Double indexation on t indices and t values

print "Time, in seconds : ", time

jerk[indice] = jerkstep(time,k)

print "Step function : " , jerk[indice]

triangle[indice] =jerkint(time,k,ac)

print "Step function integrate ( = triangle ) : " , triangle[indice]

curve[indice] = jerkintint(time,k,ac,vc)

print "Step function double integrate ( = triangle integrate ) : " , curve[indice]

print "\n"

# plot figure in mathplotlib

fig, axs = plt.subplots(3, 1, sharex=True, sharey=False)

axs[0].plot(t, jerk, 'o-', label=str("Step Function"))

axs[1].plot(t, triangle, 'o-' , label=str("Triangle Function"))

axs[2].plot(t, curve, 'o-' , label=str("Curve Function"))

fig.show()

plt.show()