this thread中,alko发布了一个用于计算多元函数的偏导数的非凡答案.

我现在有一个关于增强此函数以接受输入值数组的后续问题.我有一些代码,我循环遍历一大堆n维点,计算每个变量的偏导数,这在计算上非常昂贵.

使用np.vectorize向所讨论的函数进行矢量化很容易,但它会导致partial_derivative包装器出现问题:

from scipy.misc import derivative

import numpy as np

def foo(x, y):

return(x**2 + y**3)

def partial_derivative(func, var=0, point=[]):

args = point[:]

def wraps(x):

args[var] = x

return func(*args)

return derivative(wraps, point[var], dx=1e-6)

vfoo = np.vectorize(foo)

>>>foo(3,1)

>>>10

>>>vfoo([3,3], [1,1])

>>>array([10,10])

>>>partial_derivative(foo,0,[3,1])

>>>6.0

>>>partial_derivative(vfoo,0,[[3,3], [1,1]])

>>>TypeError: can only concatenate list (not "float") to list

理想情况下,最后一行应返回[6.0,6.0].在这种情况下,提供给向量化函数vfoo的两个数组基本

python求偏导函数_Python中多元函数的向量化偏导数相关推荐

  1. python求偏导_python实现点位精度评定

    点位精度评定,主要是评定点位数据得l离散度.评定主要是一系列得数学值.在python的numpy中有一些列的函数可以用来实现这一系列的评估值. 1.期望 期望表示一点点位最可能出现的位置.一般情况下期 ...

  2. python求乘积内建函数_Python中的内建函数(Built_in Funtions)

    前言 在Python官方文档的标准库章节中,第一节是简介,第二节就是Built_in Functions,可见内建函数是Python标准库的重要组成部分,而有很多内建函数我们平时却很少用到或根本就不知 ...

  3. python算法和数据结构_Python中的数据结构和算法

    python算法和数据结构 To 至 Leonardo da Vinci 达芬奇(Leonardo da Vinci) 介绍 (Introduction) The purpose of this ar ...

  4. python 求导实现_python怎么实现函数求导

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

  5. python决策树 多分类_Python中的决策树分类:您需要了解的一切

    python决策树 多分类 什么是决策树? (What is Decision Tree?) A decision tree is a decision support tool that uses ...

  6. python求众数程序_python求众数问题实例

    本文实例讲述了python求众数问题的方法,是一个比较典型的应用.分享给大家供大家参考.具体如下: 问题描述: 多重集中重数最大的元素称为众数...就是一个可以有重复元素的集合,在这个集合中重复的次数 ...

  7. python数据库模糊查询_python中数据库like模糊查询方式

    python中数据库like模糊查询方式 在Python中%是一个格式化字符,所以如果需要使用%则需要写成%%. 将在Python中执行的sql语句改为: sql = "SELECT * F ...

  8. python求回文_python实现求最长回文子串长度

    给定一个字符串,求它最长的回文子串长度,例如输入字符串'35534321',它的最长回文子串是'3553',所以返回4. 最容易想到的办法是枚举出所有的子串,然后一一判断是否为回文串,返回最长的回文子 ...

  9. python求不规则图形面积_python 微积分之---黎曼和

    黎曼求和 这里有一块形状不规则的土地,要测量它的面积,怎么办呢?一个叫黎曼的德国数学家(Bernhard Riemann, 1826-1866),他想了个办法:将这不规则图形切成一条条的小长条儿,然后 ...

最新文章

  1. 2021-03-21 jdk中体统的URL能访问的协议非常有限(当然可以进行扩展,不过很麻烦),常用的有http,file,ftp等等
  2. InstallShield2013 error 6109
  3. 在博文顶部添加文章字数及阅读时间信息:阅读本文需要xx分钟
  4. 【转】typedef函数指针的用法(C++)
  5. Android实践 -- 对apk进行系统签名
  6. 博士生是大学的廉价劳动力吗
  7. UVa 1153 Keep the Customer Satisfied 【贪心 优先队列】
  8. 使用windows 命令行创建一个空的文本文件
  9. 计算机中缺失freeimage.dll,修复freeimage.dll
  10. Java项目:jsp+servlet网上会议室预约系统
  11. wifi情况下使用fiddler_如何对手机http进行抓包?Fiddler工具超好用
  12. 联想小新i2000装win10教程(u盘装系统)
  13. 每日一结(10.31)
  14. php 替换 全角符号,php如何实现全角符号转半角
  15. XMPP Client 流程及入门
  16. 使用DpInst安装驱动程序
  17. ios开发所需图片尺寸汇总
  18. 集成电路ic是什么?
  19. ubuntu windows双系统默认启动项切换
  20. LINUXUNIX 中文著名网站

热门文章

  1. 《程序员》05期上市,下一代互联网竞争格局
  2. 有哪些简单好用的电脑录音软件
  3. db dbm dbi dbc 的含义
  4. 天文相关Python资料收集
  5. 口袋精灵java源码_“数学口袋精灵”App的第一个Sprint计划----开发日记
  6. 实验吧——安全杂项(四)
  7. 如何压缩图片jpg格式?一分钟学会jpg图片压缩
  8. Mybatis的属性(properties)
  9. 企业数字化转型挑战——知识库管理
  10. README文件编写