复合五点高斯公式计算（Python实现）

刺骨的言语ヽ痛彻心扉 2022-05-23 08:46 165阅读 0赞

## 被积函数 ##

y=11\+x2  y = 1 1 + x 2

x∈\[−1,1\]  x ∈ \[ − 1 , 1 \]

## 算法分析 ##

这里使用分段的方式，这里的节点做了变换之后记得对于结果记得要做类似的变换，来确保这个数值没有问题。  
**注意到，这里本来想用跟之前的节点数目更少的方法来做分析的时候，发现解方程，算不出对应的系数Ai值**。  
会陷入开顿当中，或者是报错。  
所以，这里就直接找了对应的数据做伸缩变换。直接求解。

## 代码 ##

import numpy as np
    from sympy import *
    
    
    def t(begin, end):
        xs = np.array([-0.9061798, -0.5384693, 0, 0.5384693, 0.9061798]) * (end - begin) / 2 + (begin + end) / 2
        if point == 5:
            temp = [0.2369269, 0.4786287, 0.5688889, 0.4786287, 0.2369269]
            T = 0
            for i in range(point):
                T += (Y.subs(x, xs[i]) * temp[i])
            return T * (end - begin) / 2
        return 0
    
    
    def loss(begin=-1, end=1):
        T = sum([t(xl[i], xl[i + 1]) for i in range(n)])
        I = integrate(Y, (x, begin, end))
        print('%.18f' % (I - T).evalf())
    
    
    if __name__ == '__main__':
        x = symbols('x')
        point = 5
        n = 4
        xl = np.linspace(-1, 1, n + 1)
        Y = 1 / (1 + x ** 2)
        loss()

发表评论取消回复

表情：

评论列表（有 0 条评论，165人围观）

还没有评论，来说两句吧...

相关阅读

相关高等数学积分关系定理（格林公式、高斯公式、斯托克斯公式）的理解

1 格林公式、高斯公式、斯托克斯公式 1.1 格林公式（Green formula） ![在这里插入图片描述][77cd6a3a76044890a03d6f3c82

女爷i/ 2024年02月17日 12:32/ 0 赞/ 72 阅读

相关高斯核函数 python实现

高斯过程 Gaussian Processes 是概率论和数理统计中随机过程的一种，是多元高斯分布的扩展，被应用于机器学习、信号处理等领域。本文对高斯过程进行公式推导、原理阐述

矫情吗；*/ 2023年09月26日 15:41/ 0 赞/ 140 阅读

相关贝叶斯公式

讲解的不错尤其是`现象`和`规律`的说法 ![1515886-20190928121823736-2138558902.png][] [贝叶斯公式由浅入深大讲解—A

清疚/ 2023年06月01日 05:48/ 0 赞/ 44 阅读

相关高数高斯公式、泰勒公式、麦克劳伦公式、欧拉公式、傅里叶公式

![gif.latex][]![gif.latex][]一、高斯公式参考文章：[图像处理高斯滤波（带权值的均值滤波）][Link 1] 1、一维高斯滤波。

妖狐艹你老母/ 2023年05月30日 03:46/ 0 赞/ 50 阅读

相关计算公式

股票今日涨幅如何计算涨幅=(现价-上一个交易日收盘价）/上一个交易日收盘价\100% 涨幅=（今天收盘价 / 昨天收盘价 - 1）\ 100 换手率怎么算

桃扇骨/ 2022年08月02日 04:54/ 0 赞/ 317 阅读

相关复合高斯积分（节点数小于等于3的版本Python实现）

被积函数 y=11\+x2 y = 1 1 + x 2 x∈\[−1,1\] x ∈ \[ − 1 , 1 \] 算法描述缺点：在下面代码

- 日理万妓/ 2022年05月23日 08:47/ 0 赞/ 105 阅读

相关复合五点高斯公式计算（Python实现）

被积函数 y=11\+x2 y = 1 1 + x 2 x∈\[−1,1\] x ∈ \[ − 1 , 1 \] 算法分析这里使用分段的方式

刺骨的言语ヽ痛彻心扉/ 2022年05月23日 08:46/ 0 赞/ 166 阅读

相关贝叶斯公式

例已知概率： P(X = A) –原因为A的概率 P(Y=B | X=A) –在原因为A的前提下，结果为B的条件概率需要求的

一时失言乱红尘/ 2022年05月19日 05:08/ 0 赞/ 245 阅读

相关使用递归高斯滤波器实现快速高斯模糊

高斯窗常用于对图像进行模糊或低通滤噪，但是随着高斯半径的增加，时间消耗会逐级增加如高斯半径为N时，计算每个输出采样点需要计算的乘法次数为（2N+1）\模糊方向数，加法次数为

电玩女神/ 2021年11月23日 21:50/ 0 赞/ 402 阅读

相关贝叶斯公式

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ub

﹏ヽ暗。殇╰゛Y/ 2021年11月19日 13:06/ 0 赞/ 381 阅读