r*y
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，这个用于函数的，用
: 不起来。不知道还有什么可以算，谢谢！
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，这个用于函数的，用
: 不起来。不知道还有什么可以算，谢谢！
相关阅读
懂fortran的帮忙看一下恳求一文 http://iopscience.iop.org/0256-307X/27/11/110301/Chase信用卡$500刀大放送 双倍Bonus 限时抢！未名有申请研究生的版面吗？请教如何将.m文件和matlab关联急问电脑硬盘坏了如何恢复数据是当初出茅庐的30多岁导师的大弟子还是当泰山北斗级的60多岁大牛的关门弟子?给大家拜个早年！问个数据库的入门问题Computational Chemistry PhD找软件工程的工作有戏吗？问一个有关Fortran的问题问个很2的问题，我装的Parallel Studio XE 2011怎么找不到？ (招聘SQA Instructor。 Great Pay！！！user agent questionc++, sort() 为啥显示 结果不对求助：一个C语言code在自己机器上跑正常，在cluster上出现segmentation fault，可能是什么原因。请大家推荐一下湾区biotech的工作机会（计算或编程方向)刚捐了100到版面，发几个包子Amex Blue Cash preferred 6% Cashback+$150bonus超市专用信用卡★★★有个问题请教大家 (5个包子感谢帮忙)★★★