matlab用辛普森公式求积分_标准正态分布概率密度函数的定积分计算方法及Python实现代码...

最新推荐文章于 2024-04-25 11:47:50 发布
最新推荐文章于 2024-04-25 11:47:50 发布
阅读量3.7k 收藏 5
点赞数 1
文章标签： matlab用辛普森公式求积分 python 计算曲线面积 python解zuobiaoxi方程 python计算卡方分布 利用雅可比方法求线性方程组C语言 如何把密度函数化为标准正态二维分布

efb0bd91594f57ae838c7c96948802c3.png

最近利用碎片时间在读Allen B.Downey的《贝叶斯思维：统计建模的Python学习法》，顺便用手机上的Pythonista写实例。因为Pythonista没有scipy科学计算包，遇上需求标准正态累积分布函数的时候就只能抓瞎，为此决定自己写一个。累积分布函数（Cumulative Distribution Function，CDF）就是概率密度函数（Probability Density Function，PDF）的积分，利用原函数计算定积分的方法建立在牛顿-莱布尼兹公式之上

。然而，原函数可以用初等函数表示的函数为数不多，大部分的可积函数的积分无法用初等函数表示，甚至无法有解析表达式，比如在统计学上很重要的正态分布函数
就是这样一个无法用初等函数直接计算的函数。

在知乎上查找一遍这个积分的计算方法，没找到理想的答案。所以决定写下来和大家一起学习一下。分别用数值积分法的复化辛普森公式和无穷级数的泰勒级数法求解这个积分，包括部分理论、算法和Python代码。

定积分

在几何上可以解释为由
,
,
以及
这四条边所围成的曲边梯形面积。如图1所示。而这个面积之所以难于计算是因为它有一条曲边
。

ffa3db7dd34116d096fe4b901373a920.png
图1：定积分原理

一、数值积分是利用黎曼积分等数学定义，用数值逼近的方法近似计算给定的定积分值。

1. 矩形公式：就是常见的黎曼和，矩形的高由函数值来决定，可选择使用左矩形、右矩形或中矩形。

28a9bf5026113d2c21698aadfd08b850.png
图2

如图2所示，用矩形面积来拟合积分，公式为

。通常将积分区间
最低0.47元/天 解锁文章
weixin_39705193
关注
  • 1
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
python梯形法计算定积分_用矩形法（梯形法）求定积分
weixin_39909859的博客
12-06 3890
分析：高中的时候，我们学习过，可以通过矩形法或者矩形法来求定积分。思路就是将积分区间划分成n等份，然后将这n等份近似看成矩形(或梯形)，然后对所有的矩形(或梯形)的面积进行求和。简单的例子：求函数X^2在的定积分矩形法：#include#includeusing namespace std;int main(){float fun(float x);float a,b;cout<cin&gt...
皮尔逊相关系数计算（Matlab和Python实现）
想握住此生辽阔
11-03 5748
Matlab 注意，数组元素之间用分号隔开的 clear all; a=[ 0.6; 0.9;0.78; 0.28;0.44;0.40;0.69;0.94;0.45] b=[ 31.46;26.75;32.08 ;26.38;29.75;27.57;28.56 ;27.67 ;28.18 ] coeff = corr(a , b,'type','pearson'); Python from math import sqrt def multipl(a,b): sumofab =
数值积分——复化(Simpson)辛普森求积公式 | 北太天元 or matlab
最新发布
以算法实现为主
04-25 3688
数值积分中的复化Simpson求积公式，主要以算法实现为主，使用软件为北太天元
MATLAB 与 C 语言的混合编程实战之辛普森积分法、自适应辛普森积分
菱形继承
11-19 4780
题目要求 题目大意是让你用c系语言实现辛普森积分法对定积分的粗略估计，所谓辛普森积分法即为： 定义：辛普森法则（Simpson's rule）是一种数值积分方法，是牛顿-莱布尼茨公式的特殊形式，以二次曲线逼近的方式取代矩形或梯形积分公式，以求得定积分的数值近似解。其近似值如下： ...
辛普森 matlab,利用MATLAB软件编写辛普森求积公式程序
weixin_28896781的博客
03-23 4291
a=input('a=?');b=input('b=?');x0=(a+b)/2;f=[0 0 0];f(1)=a^2*cos(a);f(2)=x0^2*cos(x0);f(3)=b^2*cos(b);s=((b-a)/6)*(f(1)+4*f(2)+f(3))a=input('a=?');b=input('b=?');x=a:(b-a)/2:bf=(x.^2).*cos(x)s=((b-a)/6...
辛普森求积公式 和 复合辛普森求积公式 Matlab 实现
热门推荐
Gy_1096440940的博客
11-04 5万+
辛普森求积公式 和 复合辛普森求积公式 Matlab 实现辛普森求积公式 利用区间[a,b]的端点及中点计算积分 ∫x2x1f(x)dx≈b−a6×(f(a)+4×f(b−a2)−+f(b)) \int_{x_1}^{x_2} f(x)dx \approx \frac{b-a}{6} \times (f(a) + 4 \times f(\frac{b-a}{2}) -+f(b))
龙贝格求定积分.rar
12-12
《龙贝格求定积分》是一个关于数值积分的MATLAB程序，主要采用了Romberg方法来求解定积分问题。在数值分析中，由于解析求解定积分并不总是可能的，因此我们常常需要借助于数值方法。Romberg方法是其中一种高效且精确...
哈工大计算方法.7z哈工大计算方法.7z
12-03
3. **编程作业**：计算方法往往与编程紧密相连，学生可能需要用MATLAB、Python、C++或其他编程语言实现数值算法。压缩包可能包含示例代码和作业要求。 4. **实验报告**：课程可能包含实验部分，要求学生在实际环境...
合肥工业大学_计算方法实验（完成版，全）.zip
01-12
8. **软件工具的使用**：实验可能需要使用MATLAB、Python等编程语言，教授如何编写代码实现上述算法。 9. **误差分析**：学习如何评估数值方法的精度，理解绝对误差、相对误差和截断误差的概念。 10. **实验报告...
数值积分_积分_数值积分_
09-29
在实验中，学生可以使用各种编程语言（如Python、C++、MATLAB等）实现上述积分方法。例如，Python的`scipy.integrate.quad`函数可以方便地进行数值积分，而自己编写算法则能更好地理解其实现原理。 六、实验练习 ...
实验4 数值积分20160405_mod3_9884010221
08-03
实验的目标是让学习者深入理解数值积分的算法和实现过程，而不是简单地依赖特定的编程环境，如Python、Matlab或Mathematica。虽然C/C++被推荐使用，但实验的目的是掌握算法本身，而非语言的便利性。 首先，实验要求...
复化梯形法和辛普森数值积分的matlab实现程序
06-30
这是一段复化梯形法和辛普森数值积分的matlab实现程序。
自适应Simpson积分公式matlab源代码
04-08
数值积分课程需编制自适应的Simpson公式，此代码采用递归函数，函数中采用了fcnchk函数，matlab6.5及以下版本会报错，只需将函数定义语句改成inline函数即可
17 误差分布曲线的建立 - 辛普森的研究
阿甘兄
09-16 1081
十八世纪，天文学的发展积累了大量的天文学数据需要分析计算，应该如何来处理数据中的观测误差成为一个很棘手的问题。我们在数据处理中经常使用平均的常识性法则，千百年来的数据使用经验说明算术平均能够消除误差，提高精度。平均有如此的魅力，道理何在，之前没有人做过理论上的证明。算术平均的合理性问题在天文学的数据分析工作中被提出来讨论：测量中的随机误差应该服从怎样的概率分布？算术平均的优良性和误差的分布有怎样的...
matlab求条件概率密度_你真的会用程序求多重积分吗？
weixin_39603117的博客
11-22 973
学过高数的同学都应该知道怎么手算多重积分，这里的多一般最多是三重积分，万一是四重或者更多重积分呢？手算肯定不现实了，如何编制程序来计算多重积分呢？其实根据微积分的思想可以把积分转化为离散求和，如下图二重积分示意图，将二重积分转化为一个个小立方体的体积和。这种方法可以推广到任意重积分，但是如果利用这种方法来编制程序求解多重积分效率是十分低的，因为随着积分重数的增加，会调用很多的循环结构，不听劝的同学...
first day!
xchl321的博客
11-19 214
first day! 自定义目录标题) 欢迎使用Markdown编辑器 你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。 新的改变 我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博客： ...
第二单元 用python学习微积分（十六）定积分
牛眼观世界
02-08 1959
本文内容来自于学习麻省理工学院公开课：单变量微积分-定积分-网易公开课 一、定积分的概念 几何意义找到曲线下的面积（另 累积和） ​ 和不定积分的区别在于，不定积分并没有给出上限或下限（a,b） 二、计算面积的步骤 1、切割成多个"矩形" 2、把这些"矩形"的面积累加起来 3、修正前面的结果，（通过让"矩形"变窄来取得极限值，到无穷窄） ​ 三、例子 1、 ，a = 0, b 任意 面积= 这里要计算 ， 老师给了一个算法，把数列中的每个数看作一个高..
矩形法求定积分
weixin_62397764的博客
12-20 1855
、、 说明：sin,cox,exp函数已在系统的数学函数库中，程序开头要用#include<math.h> 矩形法：将一个图形划分为多个小块矩形，求取每个小块矩形的面积，进行相加，得到最终结果 定积分的上下限：0，1 等分区间个数：10000 小块矩形的底部边长：dx=(1-0)/10000 小块矩形的高度：具体使用的函数值 小块矩形的面积：Sn=func(x)*dx S=S1+S2+...+Sn void cal(double (*handle) (double)) {
matlab用辛普森公式求积分_1.1高等数学公式大全
weixin_39658474的博客
11-05 725
摘自各种资源，希望帮助到大家，欢迎关注我！目录导数公式基本积分表三角函数的有理式积分初等函数两个重要极限三角函数公式和差角公式和差化积公式倍角公式半角公式正玄余弦公式莱布尼兹公式、中值定理、曲率定积分应用相关公式空间几何平面方程二次曲面多元函数微分多元复合函数及应用隐函数求导微分在几何上的应用方向导数与梯度多元函数极值重积分及应用柱面坐标和球面坐标曲线积分曲面积分高斯公式斯托克斯公式常数项级数审敛...
math.exp(-x**2)辛普森公式求积分，并求估计误差的python代码
06-13
书，然后修改该图书的存馆数量。注意，在借阅时需要检查该图书的存馆数量是否大于0。 ```c void borrow_book() { char isbn[20]; printf("请输入要借阅的图书的以下是使用辛普森公式求解math.exp(-x**2)的积分以及计算估计误差的书号："); scanf("%s", isbn); int i; for (i = 0; i < num_books; iPython代码： ```python import math # 定义被积函数 def f(x): return math.exp(-x**2) ++) { if (strcmp(isbn, books[i].isbn) == 0) { // 找到了要借阅的图# 定义辛普森公式 def simpson(a, b, n): h = (b - a) / n 书 if (books[i].quantity > 0) { books[i].quantity--; printf("借阅成功\n"); sum1 = 0 sum2 = 0 for i in range(1, n, 2): sum } else { printf("该图书已经全部借出\n"); } break; } } if1 += f(a + i * h) for i in range(2, n, 2): sum2 += f(a + (i == num_books) { printf("未找到该书号对应的图书\n"); } } void return_book i * h) return (b - a) / 3 / n * (f(a) + f(b) + 4 *() { char isbn[20]; printf("请输入要归还的图书的书号："); scanf("%s", isbn sum1 + 2 * sum2) # 计算积分和估计误差 a = 0 b = 1); int i; for (i = 0; i < num_books; i++) { if (strcmp(isbn, books n = 100 integral = simpson(a, b, n) error = abs((1 / math.sqrt(math.pi) - integral[i].isbn) == 0) { // 找到了要归还的图书 books[i].quantity++; printf) / 1 / math.sqrt(math.pi)) # 输出结果 print('积分结果：', integral) print('估计误差("归还成功\n"); break; } } if (i == num_books) { printf("未找到该：', error) ``` 在上面的代码中，我们首先定义了被积函数f(x)，然后定义了辛书号对应的图书\n"); } } ``` 以上是一个初步的设计，可以根据需要进行修改和普森公式simpson(a, b, n)。最后，我们使用定义好的函数计算积分和估计误差，并输出结果。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值