位置: IT常识 - 正文

python怎么实现函数求导(python如何编写函数)

编辑:rootadmin

推荐整理分享python怎么实现函数求导(python如何编写函数),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:python中编写函数,用python编写函数,用python写函数,python编写函数,python怎么运用函数,python怎么运用函数,python怎么运用函数,python怎么运用函数,内容如对您有帮助,希望把文章链接给更多的朋友!

python实现函数求导的方法是:1、利用sympy库中的symbols方法传入x和y变量;2、利用sympy库中的diff函数传入需要求导的函数即可返回求导之后的结果。

python利用sympy库对某个函数求导,numpy库使用该求导结果计算的程序

在python数据处理过程中,我们经常会遇见这样一种情况。需要对一个函数表达式求偏导,并将具体数值代入导数式。

而python中通常可用于函数求导的函数是sympy库中的diff()函数。但他通常所求得的导数只是一个符号表达式。不能直接带入数据使用。

如下例:

importsympyasspimportnumpyasnpx,y=sp.symbols('xy')z=sp.sin(2*sp.pi*x+2*y/5)zx=sp.diff(z,x)zy=sp.diff(z,y)print(zx)print(zy)python怎么实现函数求导(python如何编写函数)

其输出为:

2*pi*cos(2*pi*x+2*y/5)2*cos(2*pi*x+2*y/5)/5

那么该如何解决这个问题呢?

对x,y使用evalf()函数分别赋值后,用float进行类型转换后,才能利用numpy进行数值计算。

如下例:

importsympyasspimportnumpyasnpx,y=sp.symbols('xy')z=sp.sin(2*sp.pi*x+2*y/5)zx=sp.diff(z,x)zy=sp.diff(z,y)x1=10y1=5z_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不是单一的值呢?而是一个数组。

我们可以利用一个循环来完成。

如下例:

importsympyasspimportnumpyasnpx,y=sp.symbols('xy')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偏导的空列表foriinrange(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.614727694.111638646.029462890.89585862-5.2854481-5.28544810.895858626.029462894.11163864-2.61472769][-0.166458730.261755050.383847530.05703213-0.33648208-0.336482080.057032130.383847530.26175505-0.16645873]

由此便实现了由sympy得到求导结果,到numpy库进行数值计算。

本人还是python初学者,有什么错误恳请各位大佬及时指正~学习路上共同进步~

本文链接地址:https://www.jiuchutong.com/zhishi/304565.html 转载请保留说明!

上一篇:phpcms v9无法更新url(phpcms v9官网)

下一篇:帝国cms手机模板怎么用(帝国cms手机模板是什么)

  • 微信公众号怎样推广出去(微信公众号怎样一天发多次)

    微信公众号怎样推广出去(微信公众号怎样一天发多次)

  • 苹果13怎么设置蜂窝号码标签(苹果13怎么设置壁纸和锁屏)

    苹果13怎么设置蜂窝号码标签(苹果13怎么设置壁纸和锁屏)

  • 华为微信聊天记录怎么迁移到另一个手机(华为微信聊天记录删除了怎么恢复)

    华为微信聊天记录怎么迁移到另一个手机(华为微信聊天记录删除了怎么恢复)

  • 电话号码被标记了怎么取消(电话号码被标记房产中介怎么解除)

    电话号码被标记了怎么取消(电话号码被标记房产中介怎么解除)

  • 美团兑换的店铺红包能换回来吗(美团兑换的店铺红包怎么用)

    美团兑换的店铺红包能换回来吗(美团兑换的店铺红包怎么用)

  • 世界上第几台计算机名为ENIAC(世界上第几台计算机叫eniac)

    世界上第几台计算机名为ENIAC(世界上第几台计算机叫eniac)

  • qq没有消息有却有红点(qq没有消息但是显示有消息)

    qq没有消息有却有红点(qq没有消息但是显示有消息)

  • 接受cookie是什么意思(接收cookie)

    接受cookie是什么意思(接收cookie)

  • 显卡掉了一个很小的电阻(显卡掉了一个很小的电容)

    显卡掉了一个很小的电阻(显卡掉了一个很小的电容)

  • 苹果se2系统是ios多少(苹果se2是ios几)

    苹果se2系统是ios多少(苹果se2是ios几)

  • 苹果手机系统触感反馈是什么意思(苹果手机系统触感反馈怎么不能用了)

    苹果手机系统触感反馈是什么意思(苹果手机系统触感反馈怎么不能用了)

  • 苹果8日版和国行有什么区别(苹果8日版和国行哪个好)

    苹果8日版和国行有什么区别(苹果8日版和国行哪个好)

  • rx460相当于gtx什么卡

    rx460相当于gtx什么卡

  • 微信有自动回复功能吗(微信有自动回复功能在哪里)

    微信有自动回复功能吗(微信有自动回复功能在哪里)

  • windows的桌面是指什么(windows的桌面是一个文件夹)

    windows的桌面是指什么(windows的桌面是一个文件夹)

  • 华为agsw09是什么型号(agsw09是华为哪一款)

    华为agsw09是什么型号(agsw09是华为哪一款)

  • 坚果pro3是什么处理器(坚果pro3是什么指纹传感器)

    坚果pro3是什么处理器(坚果pro3是什么指纹传感器)

  • csdn怎么发布博客(csdn如何发表博文)

    csdn怎么发布博客(csdn如何发表博文)

  • 微信怎么加密码锁屏(微信怎么加密码不让别人看)

    微信怎么加密码锁屏(微信怎么加密码不让别人看)

  • 爱奇艺扫一扫登陆在哪(爱奇艺扫一扫登录会员)

    爱奇艺扫一扫登陆在哪(爱奇艺扫一扫登录会员)

  • 怎样找回id密码(怎样找回id密码 苹果手机7p la s)

    怎样找回id密码(怎样找回id密码 苹果手机7p la s)

  • 苹果自定义运营商名字(苹果自定义运营商名称个性)

    苹果自定义运营商名字(苹果自定义运营商名称个性)

  • 13306怎么添加新乘客(怎么添加新的)

    13306怎么添加新乘客(怎么添加新的)

  • Ghyakar村,尼泊尔上木斯塘 (© Frank Bienewald/Alamy)(尼泊尔乡村)

    Ghyakar村,尼泊尔上木斯塘 (© Frank Bienewald/Alamy)(尼泊尔乡村)

  • React入门笔记(react 入门)

    React入门笔记(react 入门)

  • 10-Ajax&Vue

    10-Ajax&Vue

  • 吸收合并的税务问题
  • 普通增值税税率多少
  • 维修费税率是13的话归口什么类别
  • 增值税延期滞纳金是多少
  • 本年利润借方余额意味什么
  • 买手机手续费怎么算
  • 运输费全额抵扣什么意思
  • 税控盘服务商
  • 个税抵扣租房能抵扣多少
  • 现代服务业包括哪些行业税率
  • 旧房转让土地增值税计算案例
  • 购买资产佣金应该怎么算
  • 费用怎样确认是可抵扣暂时性差异?
  • 2018年收入等级划分图
  • 顺丰的收派服务费可以抵扣吗
  • 租赁类专票不含增值税吗
  • 广告传媒公司的名字
  • 如何处理走逃失的人
  • 营业外支出调减应该填在纳税调整项目明细表哪一项
  • 住房公积金证书更新
  • 上月计提的工资,下月做调整
  • 党员经费的使用范围
  • 调整以前年度应付工资怎么调
  • 华硕笔记本电脑售后维修服务网点
  • 其他权益工具投资公允价值变动怎么计算
  • pb估值是什么意思
  • 5000以下的固定资产当月一次性扣除怎么申报
  • php serialize()与unserialize() 不完全研究
  • 招商银行直播在哪个平台
  • 其他综合收益和其他收益的区别
  • typecho插件开发教程
  • 外资企业采购退税流程
  • 如何分清福利性劳动
  • wordpress自定义函数
  • 无形资产的摊销应计入什么科目
  • 哈士奇新手礼包
  • 一文看懂华为新品发布会
  • 现代c++教程
  • 企业所得税季报资产总额怎么填
  • 财政拨入的科技费是什么
  • spring 门面模式
  • 数据库平移
  • 企业单方面调整员工的工作岗位
  • 一般纳税人内账考虑税额吗
  • 节能控制模块
  • 耕地占用税如何计算
  • 装修店面装修费的会计分录
  • 百旺如何作废发票
  • sql server定时作业
  • 应补退税额是什么意思 法律
  • 债券到期的账面价值计算
  • 预提费用汇算清缴前取得发票
  • 收到银行承兑如何做账
  • 以公允价值计量的金融资产
  • 制造费用如何控制
  • 发票抬头注意事项
  • 如何进行固定资产审计
  • 接手前往年的银行余额不平怎么办?
  • 私企需要计提盈余公积吗
  • mac80211源码分析
  • linux查看du
  • ubuntu20.04命令
  • windows7开机
  • win10系统如何设置不休眠
  • windows8怎么装
  • win8设置无线网络
  • win8怎么设置时间自动更新
  • 程序员用linux可以干什么
  • android 数据库app
  • linux安装unixodbc
  • js开发android应用
  • NGUI之UIRoot
  • webpack循环引用
  • 哪个机场离东莞樟木头最近
  • 党建服务中心属于哪个区
  • 自然人办税服务平台怎么登录不进?
  • 2021河南税务政策
  • 税务备案表付汇日期是否可以提前
  • 耕地占用税免税项目需要计入计税依据吗
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设