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

  • 苹果手机如何进入补电模式(苹果手机如何进入快捷指令)

    苹果手机如何进入补电模式(苹果手机如何进入快捷指令)

  • 微信小程序定位不了位置怎么办(微信小程序定位权限怎么打开)

    微信小程序定位不了位置怎么办(微信小程序定位权限怎么打开)

  • mssqlserver服务无法启动(mssqlserver服务无法启动17051 sql2008)

    mssqlserver服务无法启动(mssqlserver服务无法启动17051 sql2008)

  • 塞班系统还能用吗(塞班系统还能用qq和微信)

    塞班系统还能用吗(塞班系统还能用qq和微信)

  • 微信删除账单怎么关闭人脸识别(微信删除账单怎么恢复人脸识别)

    微信删除账单怎么关闭人脸识别(微信删除账单怎么恢复人脸识别)

  • 微信 完成了群公告是什么意思(微信完成了群公告什么意思,怎么做)

    微信 完成了群公告是什么意思(微信完成了群公告什么意思,怎么做)

  • 微信视频号怎么关注别人(微信视频号怎么直播)

    微信视频号怎么关注别人(微信视频号怎么直播)

  • iphone11高德地图用不了(苹果11高德地图信号弱怎么解决)

    iphone11高德地图用不了(苹果11高德地图信号弱怎么解决)

  • 韩剧tv屏幕自动变暗(为啥韩剧tv会自动播放)

    韩剧tv屏幕自动变暗(为啥韩剧tv会自动播放)

  • 计算机主要有哪几部分组成(计算机主要有哪些应用领域)

    计算机主要有哪几部分组成(计算机主要有哪些应用领域)

  • 抖音直播结束后如何看重播(抖音直播结束后能看回放吗)

    抖音直播结束后如何看重播(抖音直播结束后能看回放吗)

  • qq后面有个帆船标志(qq后面有个帆船是什么意思)

    qq后面有个帆船标志(qq后面有个帆船是什么意思)

  • 安卓抖音怎么投屏到电视(安卓抖音怎么投影电视上)

    安卓抖音怎么投屏到电视(安卓抖音怎么投影电视上)

  • 华为p40pro如何截屏(mate50怎么长截图)

    华为p40pro如何截屏(mate50怎么长截图)

  • k歌红包怎么没了(k歌红包没领完可以退回来嘛)

    k歌红包怎么没了(k歌红包没领完可以退回来嘛)

  • ip协议位于osi什么层(ip协议对应于osi7层模型中的第几层)

    ip协议位于osi什么层(ip协议对应于osi7层模型中的第几层)

  • 手机勿扰模式怎么关闭(手机勿扰模式怎么打进去电话)

    手机勿扰模式怎么关闭(手机勿扰模式怎么打进去电话)

  • airpods亮红灯是什么意思(airpods指示灯红色)

    airpods亮红灯是什么意思(airpods指示灯红色)

  • iphone11什么时候开始卖(iPhone11什么时候生产的)

    iphone11什么时候开始卖(iPhone11什么时候生产的)

  • id不在激活状态怎么办(id不在激活状态能不能退出)

    id不在激活状态怎么办(id不在激活状态能不能退出)

  • 小米钱包安装别的手机(钱包原小米钱包)

    小米钱包安装别的手机(钱包原小米钱包)

  • 对方手机关机微信语音会显示什么(对方手机关机微信语音能打通吗)

    对方手机关机微信语音会显示什么(对方手机关机微信语音能打通吗)

  • Nvidia GPU 最新计算能力表(CUDA Compute Capability)(nvidia gpu系列)

    Nvidia GPU 最新计算能力表(CUDA Compute Capability)(nvidia gpu系列)

  • YOLOv5-v6.0学习笔记(yolov5m)

    YOLOv5-v6.0学习笔记(yolov5m)

  • 劳务公司在异地做项目需要提供当地完税证明
  • 投标报名费开什么类别发票
  • 政府专项补助资金需要交所得税吗
  • 不得免征和抵扣税额抵减额
  • 应付账款和长期负债
  • 手工账写错字怎么改呢
  • 计提分红款体现在利润表中哪一个
  • 合同履约成本合同结算如何处理抵消
  • 支付利息的诉讼请求
  • 金蝶软件如何设置单价小数点
  • 小规模企业怎么报税
  • 分公司开票总公司收款行吗
  • 净利润亏损怎么结转
  • 无形资产租赁税率多少
  • 税务局退的教育费附加怎么做账
  • 销售坚果税率多少
  • 公司投资理财需要纳税吗
  • 本月无票收入到下月开具如何申报
  • 低值易耗品怎么结转
  • 暂估在建工程会计科目
  • 应交税费已交税金借方有余额
  • 外包公司代缴的税怎么算
  • 研发产品对外销售研发材料不得加计扣除分录
  • 购买办公用品发放给各部门使用
  • 金税卡服务费抵扣政策
  • rftray.exe - rftray是什么进程 有什么用
  • .inc是什么文件
  • 普通发票主营业务收入销项负数发票怎么做账
  • 电脑文件删除如何找到
  • PHP:mcrypt_enc_is_block_algorithm()的用法_Mcrypt函数
  • 摊销租赁费的会计分录
  • 通行费发票怎么做分录
  • LNMP部署laravel以及xhprof安装使用教程
  • php有面向对象吗
  • yii2中dropDownList实现二级和三级联动写法
  • php扫二维码
  • 固定资产为什么提折旧,有何实际意义
  • vue区别
  • 深入理解linux内核
  • js调用函数的几种方法
  • 餐饮小微企业优惠政策2021
  • 企业应收票据贴现,应按实际收到的金额计入
  • viite
  • php去掉最后一个逗号
  • RabbitMQ个人实践
  • mongodb replica set 配置高性能多服务器详解
  • sql server 2008 2014
  • 固定资产盘点盘什么
  • 本年利润的会计分录
  • 增值税进项税销项税月末需要结转吗
  • 股权收购协议书实际案例
  • 金银首饰以旧换新所得税处理
  • 暂估入库成本的会计分录
  • 跨年发票应该如何处理
  • 生产车间包括什么
  • 预付账款摊销会计分录
  • 申报个税劳务报酬按20%扣
  • 借款和工程款可以分开吗
  • 项目完工后员工工作总结
  • 公司增资怎么办理手续
  • sqlserver获取数据库名
  • ubuntu命令行中文乱码怎么解决
  • ubuntu更换版本
  • RedHat 9.0下Apache+PHP+MySQL服务器安装配置
  • windows7的使用方法
  • win7待机时间在哪调
  • excel嵌入图片变成代码怎么办
  • flash和javascript
  • android系统介绍
  • shell脚本 su
  • js实现双击屏幕放大
  • js中截取字符串的方法
  • android 自定义style
  • python调大小
  • javascript 类
  • 天津武清房管局电话
  • 税务开票系统对账流程
  • 残疾人交房产契税有优惠吗
  • 泰兴市地方税务局
  • 12月份医保没有到账
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

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

    友情链接: 武汉网站建设