python计算导数并绘图的实例
作者:福斯基 发布时间:2023-01-19 21:16:51
标签:python,导数,绘图
我就废话不多说了,直接上代码吧!
import math
import numpy as np
import matplotlib.pyplot as plt
from sympy import * #用于求导积分等科学计算
def dif(left,right,step):#求导 左右区间以及间隔
x,y = symbols('x y')#引入x y变量
expr = pow(x,5)#计算表达式
x_value = [] #save x value
y_value = [] #save x f(x) value
y_value_dif = [] #save x f(x)_dot value
y_value_dif2 = [] #save x f(x)_dot2 value
y_value_dif3 = [] #save x f(x)_dot3 value
y_value_dif4 = [] #save x f(x)_dot4 value
#print(expand(exp(I*x), complex=True))#将复指数展开成实部虚部形式
expr_dif = diff(expr,x,1)
expr_dif2 = diff(expr,x,2)
expr_dif3 = diff(expr,x,3)
expr_dif4 = diff(expr,x,4)
for i in np.arange(left,right,step):
x_value.append(i)
y_value.append(expr.subs('x',i))#将i值代入表达式
y_value_dif.append(expr_dif.subs('x',i))#将i值代入求导表达式
y_value_dif2.append(expr_dif2.subs('x',i))#将i值代入2阶求导表达式
y_value_dif3.append(expr_dif3.subs('x',i))#将i值代入3阶求导表达式
y_value_dif4.append(expr_dif4.subs('x',i))#将i值代入4阶求导表达式
draw_plot_set()#设置画图格式
plt.plot(x_value,y_value,"b-",linewidth=1,label='f(x)='+str(expr)) #画图
plt.plot(x_value,y_value_dif,"r-",linewidth=1,label='f(x)_prim') #画图
plt.plot(x_value,y_value_dif2,"y-",linewidth=1,label='f(x)_prim2') #画图
plt.plot(x_value,y_value_dif3,"g-",linewidth=1,label='f(x)_prim3') #画图
plt.plot(x_value,y_value_dif4,"b-",linewidth=1,label='f(x)_prim4') #画图
plt.legend()#显示图例
plt.show()#显示图像
def draw_plot_set():#设置画图格式
plt.figure()
ax = plt.gca()
#改变坐标轴位置
ax.spines['right'].set_color('none')#删除原来轴
ax.spines['top'].set_color('none')#删除原来轴
ax.xaxis.set_ticks_position('bottom')#在0点处增加轴
ax.spines['bottom'].set_position(('data',0))
ax.yaxis.set_ticks_position('left')#在0点处增加轴
ax.spines['left'].set_position(('data',0))
#设置坐标名
plt.ylabel('f(x)')
plt.xlabel('x')
plt.grid(True)#打开网格
if __name__ == '__main__':
dif(-5,5,0.01)
补充拓展:python利用sympy库对某个函数求导,numpy库使用该求导结果计算的程序
在python数据处理过程中,我们经常会遇见这样一种情况。需要对一个函数表达式求偏导,并将具体数值代入导数式。
而python中通常可用于函数求导的函数是sympy库中的diff()函数。
但他通常所求得的导数只是一个符号表达式。不能直接带入数据使用。
如下例:
import sympy as sp
import numpy as np
x,y = sp.symbols('x y')
z = sp.sin(2*sp.pi*x+2*y/5)
zx = sp.diff(z,x)
zy = sp.diff(z,y)
print(zx)
print(zy)
其输出为:
2*pi*cos(2*pi*x + 2*y/5)
2*cos(2*pi*x + 2*y/5)/5
那么该如何解决这个问题呢?
对x,y使用evalf()函数分别赋值后,用float进行类型转换后,才能利用numpy进行数值计算。
如下例:
import sympy as sp
import numpy as np
x,y = sp.symbols('x y')
z = sp.sin(2*sp.pi*x+2*y/5)
zx = sp.diff(z,x)
zy = sp.diff(z,y)
x1 = 10
y1 = 5
z_x1 = float(zx.evalf(subs={x:x1,y:y1}))
z_y1 = float(zy.evalf(subs={x:x1,y:y1}))
print(z_x1)
print(z_y1)
其输出结果:
-2.61472768902227
-0.16645873461885696
那如果我的x或y不是单一的值呢?而是一个数组。
我们可以利用一个循环来完成。
如下例:
import sympy as sp
import numpy as np
x,y = sp.symbols('x y')
z = sp.sin(2*sp.pi*x+2*y/5)
zx = sp.diff(z,x)
zy = sp.diff(z,y)
x_array = np.linspace(-5, 5, 10)
y_array = np.linspace(-5, 5, 10)
temp_x = []#先定义一个用于存储x偏导的空列表
temp_y = []#先定义一个用于存储y偏导的空列表
for i in range(10):
z_x = float(zx.evalf(subs={x:x_array[i],y:y_array[i]}))
temp_x.append(z_x)#将计算得到的偏导值一一添加到列表中
z_y = float(zy.evalf(subs={x:x_array[i],y:y_array[i]}))
temp_y.append(z_y)
zx_array = np.array(temp_x)#将列表转换为数组
zy_array = np.array(temp_y)
print(zx_array)
print(zy_array)
输出结果为:
[-2.61472769 4.11163864 6.02946289 0.89585862 -5.2854481 -5.2854481
0.89585862 6.02946289 4.11163864 -2.61472769]
[-0.16645873 0.26175505 0.38384753 0.05703213 -0.33648208 -0.33648208
0.05703213 0.38384753 0.26175505 -0.16645873]
由此便实现了由sympy得到求导结果,到numpy库进行数值计算。
来源:https://blog.csdn.net/cool_bot/article/details/90768063
0
投稿
猜你喜欢
- 一、前言确保安装scikit-imagenumpy二、Dataset一个例子:# 导入需要的包import torchimport torc
- base64模块是用来作base64编码解码的。这种编码方式在电子邮件中是很常见的。它可以把不能作为文本显示的二进制数据编码为可显示的文本信
- 本文实例讲述了python用来获得图片exif信息的库用法。分享给大家供大家参考。具体分析如下:exif-py是一个纯python实现的获取
- 为了防止再次被攻击,做个验证码过滤程序是必要的。我在网上找了一些资料,觉得用别人做好的代码总是很不爽,自己做麻又不会写复杂的代码,特别是生成
- matplotlib官方除了提供了鼠标十字光标的示例,还提供了同一图像内多子图共享光标的示例,其功能主要由widgets模块中的MultiC
- 我就废话不多说了,大家还是直接看代码吧!#coding=utf-8import threadingimport timeimport cx_
- Python编写从ZabbixAPI获取信息此脚本用Python3.6执行是OK的。# -*- coding: utf-8 -*-impor
- 自己搭建后台网站,需求:实现类似django Admin站点对每一张表的增删改查权限控制。实现步骤:1.权限控制Django框架已自带,共6
- 本文实例为大家分享了python爬取哈尔滨天气信息的具体代码,供大家参考,具体内容如下环境:windows7python3.4(pip in
- 本文实例讲述了PHP变量作用域(全局变量&局部变量)&global&static关键字用法。分享给大家供大家参考,具
- 比如我们要读取一个桌面路径下的文件 设计 一个函数 怎么写才能正确?注意以下两点就可以了1、对于一个路径 例如C:\Users\xiaomi
- 一、图的定义图是一种比树更复杂的一种数据结构,在图结构中,结点之间的关系是任意的,任意两个元素之间都可能相关,因此,它的应用极广。图中的数据
- python中字典和列表的使用,在数据处理中应该是最常用的,这两个熟练后基本可以应付大部分场景了。不过网上的基础教程只告诉你列表、字典是什么
- 概述从今天开始我们将开启一段自然语言处理 (NLP) 的旅程. 自然语言处理可以让来处理, 理解, 以及运用人类的语言, 实现机器语言和人类
- 前言我们在日常开发中,经常会遇到一些关于文件的操作,例如,实现查看目录内容的功能。类似Linux下的tree命令。统计目录下指定后缀文件的行
- 伙伴遇到一个关于excel导入数据到python中,日期变成数字而不是日期格式的问题。第一反应这个数字应该是excel里面的时间戳类似的,所
- pandas每次多Sheet写入文件,只能一次性存入,不然每次会重写文件,最后只保留最后一次的写入。# !usr/bin env pytho
- using System; using System.Data; using System.Configuration; using Sys
- 脚本内容代码如下:from mitmproxy import http, ctxfrom multiprocessing import Lo
- 概要 “SQL Server 桌面引擎”(也叫 MSDE 2000)没有自己的用户界面,因为它主要设计为在后台运行。用户通过 MSDE 20