此段代码是基于辛普森公式的积分计算方法 1.代码 %%复合辛普森求积公式 %%Y是数值(attribute=0)或具体表达式(attribute=1),interval是求积区间,n是精度(如果是数值,则为数值长度-1) function CSQF = Compound_Simpson_quadra ...
分类:
其他好文 时间:
2020-01-30 17:22:53
阅读次数:
180
1. 用1阶至4阶Newton-Cotes公式计算积分 程序: function I = NewtonCotes(f,a,b,type) % syms t; t=findsym(sym(f)); I=0; switch type case 1, I=((b-a)/2)*(subs(sym(f),t, ...
分类:
其他好文 时间:
2018-11-10 23:58:24
阅读次数:
683
画了很长时间搞懂常见的灰色模型,以及如何求解。数值计算是一门强大的学科,尤其是结合计算机计算软件(Python,Matlab)。灰色模型中的微分方程求解,涉及数值积分,梯形公式、辛普森公式,继续是复化梯形公式、复化辛普森公式,龙贝格公式(外推),控制误差 全排列算法 无重复元素: 递归实现: (1) ...
分类:
其他好文 时间:
2018-01-25 00:10:23
阅读次数:
141
[吐槽] 嗯一开始接触到这个东西其实是因为某道凸包的题目好像可以用这个奇妙的方法来算 但其实了解也不是很深,只是觉得这个东西十分有意思, 所以先稍微写一下自己的想法,了解更多之后慢慢完善 [正题] 首先这个公式是长这样的: 取一个二次函数上(嗯。。貌似是这样的)的一段区间,左端点为$(l, hl)$ ...
分类:
其他好文 时间:
2017-08-27 22:29:41
阅读次数:
245
对于积分:
只要找到被积公式的原函数F(x),利用牛顿莱普利兹公式有:
但是,实际使用这种求积分的方法往往是有困难的,因为大量的被积函数的原函数是不能用初等函数表示的;另外,当f(x)是由测量或数值计算给出的一张数据表时,牛顿莱普利兹公式也无法直接运用,因此有必要研究积分的数值计算问题。
对于一些理论的推导,大家可以看看维基百科,下面我主要给出牛顿-科特斯公式在n...
分类:
其他好文 时间:
2015-02-02 12:38:58
阅读次数:
243
题目大意:给定一棵由圆台和圆锥构成的柠檬树,月光以α的夹角平行射向地面,求阴影部分面积
补充题目大意:看到这题我产生了心理阴影,求阴影部分面积
题目不好分析,但其实就是求一堆圆和一堆梯形的面积交
样例如图(画的有点烂),将顶点看做半径为0的圆,则图中圆的半径即为给定圆的半径,圆心距为h/tan(α),直线为两圆公切线
这题我们采用辛普森自适应公式
首先辛普森公式见度受百科 htt...
分类:
其他好文 时间:
2014-09-13 17:21:05
阅读次数:
285
一 实验目的1. 掌握复合梯形公式与复合辛普森公式的基本思想。2. 编程实现用复合梯形公式与复合辛普森公式求积分。3. 熟悉matlab软件的使用。二 实验内容1、用复合梯形公式计算积分 I=4/(1+x2)dx ,求它0到1的积分。精确度为10-5.(0.00001),精确到 ●1 计算公式 .....
分类:
其他好文 时间:
2014-06-27 18:13:07
阅读次数:
865
1 //这里f()为被积函数,输入a,b为积分上下限, 2
//eps为计算精度[这里要注意假收敛,一般设小一点好] 3 #include 4 #include 5 #define eps 1e-6 6 using
namespace std; 7 double f(double x){ 8...
分类:
其他好文 时间:
2014-06-07 05:22:57
阅读次数:
239