在Python中使用Numpy对Hermite_e系列的0轴以上进行整合

当我们使用Python中的Numpy库时,可以很方便地对Hermite_e系列进行整合。下面是详细的攻略:

1. 导入Numpy库

首先,需要在代码中导入Numpy库:

import numpy as np

2. Hermite_e函数

定义Hermite_e函数,可以使用以下代码实现:

def Hermite_e(n, x):
    """
    Compute the nth "probabilists" Hermite polynomial evaluated at x.

    Parameters
    ----------
    n : int
        Degree of the Hermite polynomial. Must be non-negative.
    x : array_like
        A scalar or array of values at which to evaluate the resulting
        polynomial.

    Returns
    -------
    H_n : ndarray
        The value(s) of the nth probabilist's Hermite polynomial at `x`.

    """
    n = np.asarray(n)
    x = np.asarray(x)
    if n.ndim != 0 or x.ndim != 0:
        n, x = np.broadcast_arrays(n, x)
    if np.any(n < 0):
        raise ValueError("Degree of Hermite polynomial must be non-negative.")
    scale = np.sqrt(np.pi)*2**n
    He_0 = 1
    if n == 0:
        return He_0
    He_1 = 2*x
    if n == 1:
        return He_1
    for k in range(1, n):
        He_kplus1 = 2*x*He_1 - 2*k*He_0
        He_0 = He_1
        He_1 = He_kplus1
    return He_1/scale

3. 使用Numpy对Hermite_e的0轴以上进行整合

使用Numpy库的sum函数可以非常方便地整合Hermite_e系列中的0轴以上部分。这里使用以下代码对Hermite_e的0轴以上进行整合:

#计算重积分 cos(x*theta)*Hermite_e(n1, np.sqrt(beta)*x)*Hermite_e(n2, np.sqrt(beta)*x)*np.exp(-0.5*beta*x**2)在区间(-inf, inf)上的值

#首先定义函数f
def f(x, n1, n2, beta, theta):
    return np.cos(x*theta) * Hermite_e(n1, np.sqrt(beta)*x) * Hermite_e(n2, np.sqrt(beta)*x) * np.exp(-0.5*beta*x**2)

#定义待计算的阶数
n1 = 2
n2 = 3
#定义beta和theta的值,这里为了方便直接赋值
beta = 1
theta = 0.3
#使用quad函数计算重积分
result, _ = quad(f, -np.inf, np.inf, args=(n1, n2, beta, theta))
print(result)

这里的示例计算了重积分 cos(x*theta)*Hermite_e(n1, np.sqrt(beta)*x)*Hermite_e(n2, np.sqrt(beta)*x)*np.exp(-0.5*beta*x**2) 在区间 (-inf, inf) 上的值。

另外一个示例是计算Hermite_e 0轴以上阶数为2至5的总和,可以使用以下代码实现:

#定义一个数组保存每个阶数的值
n_array = np.array([2, 3, 4, 5])
#定义beta和x的值,这里为了方便直接赋值
beta = 1
x = 0.5

#使用NumPy的sum函数对0轴以上进行整合求和
result = np.sum(Hermite_e(n_array, np.sqrt(beta)*x))
print(result)

这里的示例计算了Hermite_e阶数在2至5之间的总和,其中 betax 的值分别为 1 和 0.5。

本站文章如无特殊说明,均为本站原创,如若转载,请注明出处:在Python中使用Numpy对Hermite_e系列的0轴以上进行整合 - Python技术站

(0)
上一篇 2023年3月25日
下一篇 2023年3月25日

相关文章

  • 详解Python 实用的WSGI应用程序

    下面详细讲解Python实用的WSGI应用程序的完整攻略。 什么是WSGI WSGI是Web服务器网关接口(Web Server Gateway Interface)的缩写。它是Python Web应用程序和Web服务器之间的一种通用接口,通过该接口,可以使得Python Web应用程序可以被任意一种Web服务器调用和运行。 WSGI框架 在Python中,…

    python-answer 2023年3月25日
    00
  • Python正则表达式中group与groups的用法详解

    以下是“Python正则表达式中group与groups的用法详解”的完整攻略: 一、问题描述 在Python中,我们可以使用正则表达式来匹配字符串。在匹配过程中,我们可以使用group()和groups()方法来获取匹配结果。本文将详细讲解Python正则表达式中group与groups的用法,并提供两个示例说明。 二、解决方案 2.1 group()方法…

    python 2023年5月14日
    00
  • 代码解析python标准库logging模块

    1. 简介 logging是Python标准库中提供的一个标准日志工具模块,其主要的作用是用来记录应用程序的运行时信息。通过使用logging来输出日志可以更加方便的进行开发、调试、测试以及发布运行等的工作。 2. logging模块的基础用法 2.1 创建Logger对象 创建Logger对象是logging模块的第一步,Logger对象主要用于向后端输出…

    python 2023年6月3日
    00
  • 使用Python 自动生成 Word 文档的教程

    请您耐心阅读以下的教程,此教程分为以下几个部分: 介绍Python生成word文档的工具库 安装工具库 创建word文档 添加文本与表格 添加图片与图表 示例说明 总结 1. 介绍Python生成word文档的工具库 目前Python生态圈里提供了多种文档生成的工具库,常用的有:python-docx,python-docx-template和docxtpl…

    python 2023年5月19日
    00
  • Python datetime 格式化 明天,昨天实例

    Python的datetime模块提供了对日期和时间的操作,包括创建日期时间对象、格式化日期时间、日期时间计算等操作。在datetime模块中有一个datetime类,它表示一个日期时间对象。格式化日期时间可以使用strftime()方法来实现。 格式化日期时间 使用strftime()方法可以将datetime对象格式化为字符串。strftime(form…

    python 2023年6月2日
    00
  • Python3.7安装PyQt5 运行配置Pycharm的详细教程

    下面是安装PyQt5并配置PyCharm的详细步骤: 1. 安装Python3.7 首先,你需要在你的电脑上安装Python3.7。你可以在Python官网(https://www.python.org/downloads/)下载最新的Python3.7版本的安装程序,并按照提示进行安装。 2. 利用pip安装PyQt5 在安装完Python3.7后,你需要…

    python 2023年5月18日
    00
  • Python排序算法之堆排序算法

    下面是详细讲解“Python排序算法之堆排序算法”的完整攻略,包含两个示例说明。 堆排序算法 堆排序算法是一种基于二叉堆的排序算法。它的基本思想是将待排序的序列构建成一个二叉堆,然后不断将堆顶元素与堆底元素交换,再重新调整,到整个序列有序为止。 堆排序算法的Python实现 下面是一个示例代码,用于实现堆排序算法: def heap_sort(arr): n…

    python 2023年5月14日
    00
  • Python爬虫入门有哪些基础知识点

    Python爬虫入门有哪些基础知识点 背景介绍 爬虫是一种按照一定规则自动抓取网页信息的程序,近年来日益风行,因其便于获取大量网络数据而受到广泛关注。Python语言作为一种简单易用、生态丰富的编程语言,自然成为了开发爬虫的首选工具。 本文将详细介绍Python爬虫入门所需的基础知识点,旨在帮助初学者快速入门,开启自己的爬虫之路。 知识点一:HTML与CSS…

    python 2023年5月14日
    00
合作推广
合作推广
分享本页
返回顶部