由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Computation版 - 定积分的求解
相关主题
有没有什么积分的快速算法做计算的大牛应该有好的机器吧?
一个奇怪的MATLAB数值积分问题Hamiltonian System有什么好的方法吗?
问一个关于使用Matlab的积分函数quad求扇形面积的问题integration question
如何高速计算这个积分MATLAB积分求教
求教一个积分int( (1/t)*exp(t), t )请教一个Gaussian quadrature的问题
how to remove a singularity in an integr感觉这儿有很多大牛
$20K for a computer这里有搞积分方程(Integral equation)的吗?
Sci computation core 2 Duo 8400 Vs. Core 2 Quad 9400求解两个联立的二阶常微分方程 (我的方程)
相关话题的讨论汇总
话题: end话题: 积分话题: spacing话题: integral话题: 简单
进入Computation版参与讨论
1 (共1页)
q***z
发帖数: 543
1
x 和 y 都已经是数值了,想算一个定积分的值。试了一下quad,这个用于函数的,用
不起来。不知道还有什么可以算,谢谢!
r****y
发帖数: 1437
2
最简单的就是composite trapezoid rule啊,用matlab一行搞定
integral = 0.5 * sum((y(1:end-1) + y(2:end)) .* diff(x))
上底加下底乘高除二,直观明了,管你x是不是equal spacing都适用.虽然
讲运算速度不是最优的,写起来最简单,也最便于检查。
I assume your values in x-vector is monotonically increasing

要是x equal spacing, 更简单了
dx = x(2)-x(1);
integral = (sum(y(2:end-1)) + 0.5 * y(1) + 0.5 * y(end) ) * h

【在 q***z 的大作中提到】
: x 和 y 都已经是数值了,想算一个定积分的值。试了一下quad,这个用于函数的,用
: 不起来。不知道还有什么可以算,谢谢!

1 (共1页)
进入Computation版参与讨论
相关主题
求解两个联立的二阶常微分方程 (我的方程)求教一个积分int( (1/t)*exp(t), t )
[转载] 会数值积分的请帮忙看一下how to remove a singularity in an integr
有人紧急问一个问题$20K for a computer
[转载] 多重积分一问Sci computation core 2 Duo 8400 Vs. Core 2 Quad 9400
有没有什么积分的快速算法做计算的大牛应该有好的机器吧?
一个奇怪的MATLAB数值积分问题Hamiltonian System有什么好的方法吗?
问一个关于使用Matlab的积分函数quad求扇形面积的问题integration question
如何高速计算这个积分MATLAB积分求教
相关话题的讨论汇总
话题: end话题: 积分话题: spacing话题: integral话题: 简单