Прежде всего, разделите ваш код на маленькие независимые блоки (высокая когезия), поэтому создайте функцию с требуемым вычислением:
def Pmotion(x,t,A,alpha,f):
w=2*np.pi*f
k1 = (w/alpha)
theta = k1*x-w*t
Up = k1*A*complex(-np.sin(theta),np.cos(theta))
Up_real = Up.real
Up_imag = Up.imag
return Up_real, Up_imag
Затем вы можете начать думать о методе построения графика. например,
def plot_Pmotion_t():
t_range = np.arange(0,2,0.1)
reals = [Pmotion(0,t,1,6000,2) for t in t_range]
plt.plot(t_range, reals)
plt.show()
Теперь вы можете свободно изменять или добавлять входы в функцию графика без изменения функции Pmotion.
Примечание. Теперь вы строите как реальные, так и мнимые значения, измените его на reals = [Pmotion(0,t,1,6000,2)[0] for t in t_range]
, чтобы отобразить только действительную часть. Надеюсь, это поможет!