位置: 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手机模板是什么)

  • 支付宝一起攒钱功能在哪里(支付宝一起攒钱的软件)

    支付宝一起攒钱功能在哪里(支付宝一起攒钱的软件)

  • 荣耀30s能自定义语音唤醒词吗(荣耀30s自定义动态壁纸)

    荣耀30s能自定义语音唤醒词吗(荣耀30s自定义动态壁纸)

  • QQ关联对方为什么都是撤回(qq关联为什么对方能看到我的我看不到他的)

    QQ关联对方为什么都是撤回(qq关联为什么对方能看到我的我看不到他的)

  • iphone6s听筒声音越来越小

    iphone6s听筒声音越来越小

  • 标题艺术字体怎么设置(标题艺术字怎么写)

    标题艺术字体怎么设置(标题艺术字怎么写)

  • 微信拉黑和删除的区别(微信拉黑和删除哪个更狠)

    微信拉黑和删除的区别(微信拉黑和删除哪个更狠)

  • 吸吸在qq中是什么意思(q吸是什么意思)

    吸吸在qq中是什么意思(q吸是什么意思)

  • 微信缩略图和大图区别(微信缩略图和原图不一样)

    微信缩略图和大图区别(微信缩略图和原图不一样)

  • 已上市的4k屏幕手机(2021年4k显示屏)

    已上市的4k屏幕手机(2021年4k显示屏)

  • 惠普打印机加墨后还是提示不足(惠普打印机加墨后怎么操作才能打印)

    惠普打印机加墨后还是提示不足(惠普打印机加墨后怎么操作才能打印)

  • 电脑服务器运行失败怎么办(电脑服务器运行中为什么要切换到激活程序)

    电脑服务器运行失败怎么办(电脑服务器运行中为什么要切换到激活程序)

  • qq群举报会不会显示自己(qq群举报成功会怎么样)

    qq群举报会不会显示自己(qq群举报成功会怎么样)

  • 苹果x和苹果11pro一样大吗(苹果x和苹果11pro参数对比)

    苹果x和苹果11pro一样大吗(苹果x和苹果11pro参数对比)

  • 红米k20长度(红米k20p长度)

    红米k20长度(红米k20p长度)

  • 手机自动播放音乐怎么回事(手机自动播放音乐找不到播放源)

    手机自动播放音乐怎么回事(手机自动播放音乐找不到播放源)

  • 快手作品里已戳怎么取消(快手上已戳别人没事吧)

    快手作品里已戳怎么取消(快手上已戳别人没事吧)

  • vivo怎么自定义微信音(vivo怎么自定义锁屏)

    vivo怎么自定义微信音(vivo怎么自定义锁屏)

  • 如何清理ps暂存盘垃圾(清理ps暂存盘暂存文件)

    如何清理ps暂存盘垃圾(清理ps暂存盘暂存文件)

  • 常见的计算机系统(常见的计算机系统结构分类法有哪三种)

    常见的计算机系统(常见的计算机系统结构分类法有哪三种)

  • 面容id坏了能修吗(面容id坏了能卖多少钱)

    面容id坏了能修吗(面容id坏了能卖多少钱)

  • pc端是电脑还是手机(pc端是电脑端的意思吗)

    pc端是电脑还是手机(pc端是电脑端的意思吗)

  • 已经发过的朋友圈还能设置权限吗(已经发过的朋友圈怎么改可见范围)

    已经发过的朋友圈还能设置权限吗(已经发过的朋友圈怎么改可见范围)

  • 米8闪屏什么原因引起(米8se屏幕老是闪)

    米8闪屏什么原因引起(米8se屏幕老是闪)

  • 怎么把抖音喜欢的清零(怎么把抖音喜欢的作品删除掉)

    怎么把抖音喜欢的清零(怎么把抖音喜欢的作品删除掉)

  • internal/modules/cjs/loader.js:985 throw err; ^ Error: Cannot find module ‘node:util‘(internet笔记)

    internal/modules/cjs/loader.js:985 throw err; ^ Error: Cannot find module ‘node:util‘(internet笔记)

  • Pont Rouge (Red Bridge) over the Palmer River in Sainte-Agathe-de-Lotbinière, Quebec, Canada (© Jean Surprenant/Getty Images)

    Pont Rouge (Red Bridge) over the Palmer River in Sainte-Agathe-de-Lotbinière, Quebec, Canada (© Jean Surprenant/Getty Images)

  • python用户如何自定义异常(python怎么自己写函数)

    python用户如何自定义异常(python怎么自己写函数)

  • 增值税专用发票的税率是多少啊
  • 高新企业享受税收优惠政策备案怎么申请
  • 月末结存材料的实际成本例题
  • 金蝶固定资产累计折旧怎么改
  • 招待费汇算清缴扣除比例
  • 购进的库存商品转自用需要开发票吗
  • 房屋出租收入其税收如何计算呢?
  • 累计折旧是费用岗负责的吗
  • 商贸公司账目
  • 代国外佣金代扣税金账务怎么处理?
  • 重庆房产税政策历史
  • 零售通收银
  • "税务 政策"
  • 用友t3建新账
  • 小规模季度划分
  • 赊销的账务处理流程
  • 10万以下免征增值税 文件
  • 股东认缴和实缴不一致
  • 会计科目发出商品
  • mac 怎么操作
  • 企业所得税汇算清缴调增的项目有哪些
  • 付款凭证怎么填写电脑
  • 库存股的会计科目
  • 民办非企业单位是什么企业类型
  • 什么时编制
  • Vue3+element-plus 后台管理系统(含登陆注册功能页面)
  • 港口机场铁路
  • 今夕七夕
  • 小企业会计准则2023电子版
  • 马拉喀什的历史背景
  • 增值税免征与即征的区别
  • go开发web项目
  • php多维数组合并相同key
  • 自动登录微信号
  • 材料发票和工程发票的区别
  • dedecms默认用户名
  • 灵动cfe
  • mongo聚合查询效率问题
  • 罚款收入算是营收收入吗
  • 法人不发工资可以吗
  • 库存商品和固定成本区别
  • sql server 2008怎么使用sql语句
  • sql函数判断一个值是否是数字
  • 公司承担员工的个税怎么做账
  • 增值是什么意思解释
  • 可转债中签后怎么看涨跌
  • 用友软件接口
  • 短期借款财务处理
  • 有限合伙企业收入怎么做账
  • 只有发票没有银行回单怎么做账
  • 补发以前年度工资的账务处理
  • 结转借贷方向如何计算
  • 旅游发票可以报销入什么明细
  • 第二年发票可以入上年账吗
  • 工程款发票怎么做分录
  • 零星采购多少不需要发票
  • 没收到电费账单怎么办
  • 农民工工资专户销户流程
  • 租赁车子
  • 预付调整到其他应付款
  • mysql mac启动
  • 为什么多出一个系统
  • centos chrony
  • win8应用商店废了
  • xp系统怎么关闭更新系统
  • win8.1设置在哪
  • unity outline性能
  • 我是如何从0开始做到年入1000+万的
  • 如何使用css设置元素的层叠效果?
  • node.js快速入门
  • 3种不同的播种方法
  • android DefaultHttpClient设置setCookieStore
  • android面试2021
  • js如何定义类
  • 河北省税务总局电话
  • 每月个人所得税申报截止时间是几号
  • 福州房管局网签查询
  • 企业怎么成立工会组织
  • 预测收入遵循什么方法
  • 个体户 零申报
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

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

    友情链接: 武汉网站建设