Python作业5

连续函数的导数定义为

df(x)/dx=lim△x→0(f(x+△x) - f(x))/△x

因此可以用下式近似求取函数在点xi  导数

f’(xi) = (f(xi+1) - f(xi))/△x                       

其中，△x= xi+1 - xi。(1) 自己选择一个函数，编写程序画出函数图形；(2) 根据上式绘制函数导数的图形。将程序代码及图形的截图上传。

我的答案 ：

 

import sympy as sy

import matplotlib.pyplot as plt

import numpy as np

plt.rcParams['font.sans-serif']=['SimHei']

plt.rcParams['axes.unicode_minus']=False

x = sy.symbols('x')

y=5*x**4+6*x

dy = sy.diff(y,x)

x_value = []

y_value = []

dy_value = []

for i in np.arange(-10,10,0.01):

    x_value.append(i)

    y_value.append(y.subs('x',i))

    dy_value.append(dy.subs('x',i))

p1=plt.figure(1)

plt.subplot(221)

plt.title('原函数图像')

plt.xlabel('x')

plt.ylabel('f')

plt.plot(x_value,y_value)

p2=plt.figure(2)

plt.subplot(222)

plt.title('一阶导函数图像')

plt.xlabel('x')

plt.ylabel('df')

plt.plot(x_value,dy_value)