有一列数据，现在想用R求曲线下某段区间的定积分该怎么办呢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如题，附件为文件，第一列Time，第二列Fz，图片为Fz~Time，想要计算曲线下某一段的面积(数据的第10行到180行)，就是求定积分，该怎么做呢谢谢各位。
谢谢各位

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：怎么办

你因该是用R

1. f <- approxfun(x = Time, y = Fz)
   
2. integrate(f, Time[10], Time[180])

复制代码

deem 发表于 2017-8-3 12:49
你因该是用R

好像不能用呢，提示 non-finite function value

deem 发表于 2017-8-3 12:49
你因该是用R

如果计算的是整个曲线下的面积呢，谢啦!!☆⌒(*＾-゜)v

deem 发表于 2017-8-3 12:49
你因该是用R

整明白了，试了试，发现上下界限值不能为端点，否则会提示maximum number of subdivisions reached但是如果是计算整个曲线下面积，上下界限到底怎么设定呢？

苏谷壳 发表于 2017-8-3 14:47
整明白了，试了试，发现上下界限值不能为端点，否则会提示maximum number of subdivisions reached但是如 ...

可能是端点处有点问题，你试试

1. integrate(f, Time[0]+1e-8, Time[180]-1e-8)

复制代码

deem 发表于 2017-8-3 15:27
可能是端点处有点问题，你试试

谢谢呀，请问1e-8这主要起什么作用呢

苏谷壳 发表于 2017-8-3 15:39
谢谢呀，请问1e-8这主要起什么作用呢

积分区间缩小一个非常小的数字

deem 发表于 2017-8-3 15:45
积分区间缩小一个非常小的数字

谢谢啦!!☆⌒(*＾-゜)v用上边附件数据，仍然出现那个问题，我再学习学习

大家还哟什么意见吗