不知道函数形式，只是到某些(X,Y)，如何先拟合再求积分呢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

b1是响应函数，b2是能量值，第一列是波长(nm)。
如何求响应函数*能量值在400nm到407nm的积分呢?
现在是不知道b1和b2的函数形式

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：函数形式 不知道 响应函数 如何

两种思路：
第一钟是拟合b1，b2两条曲线，然后利用monte carlo投点法计算积分值。
第二种直接利用给定的数据，结合数值积分里面的积分转化为级数求和的方法，两种均可行。

hzcmaster 发表于 2015-5-28 11:28
两种思路：
第一钟是拟合b1，b2两条曲线，然后利用monte carlo投点法计算积分值。
第二种直接利用给定的数 ...

我现在查到的结果是先用三次样条插值，再用数值积分，现在正在写程序，不造对不？

hzcmaster 发表于 2015-5-28 11:28
两种思路：
第一钟是拟合b1，b2两条曲线，然后利用monte carlo投点法计算积分值。
第二种直接利用给定的数 ...

我现在查到的结果是先用三次样条插值，再用数值积分，现在正在写程序，不造对不？

瞌睡的大脑 发表于 2015-5-28 14:29
我现在查到的结果是先用三次样条插值，再用数值积分，现在正在写程序，不造对不？

嗯，可行

hzcmaster 发表于 2015-5-28 14:42
嗯，可行

library(RODBC)

z<-odbcConnectExcel2007("F:\\***\\ESUN2.xlsx")

foo<-sqlFetch(z,'band1')

attach(foo)

close(z)

s<-spline(waveLenght,responseFactor,xout=seq(430,530,0.5))

m<-s$x

n=s$y

f<-function(m) return(n)

source('E:\\desktop\\code\\r\\simpson.R')

simpson(f,440,520)

我这个代码写完以后run不起来，老是说内存不足“Error: cannot allocate vector of size 16.0 Gb”
电脑内存还挺大的唉，24G的，这是怎么回事呢

瞌睡的大脑 发表于 2015-5-28 15:27
library(RODBC)

z

这个问题没遇到过，我建议你如下处理。
不用内置函数积分，得到三次样条函数之后，在400~407之间生成1000个服从均匀分布的随机数，代入到三次样条函数中，得到1000个值，将这1000个值相加乘以7再除以1000就是你的积分值。

hzcmaster 发表于 2015-5-28 16:45
这个问题没遇到过，我建议你如下处理。
不用内置函数积分，得到三次样条函数之后，在400~407之间生成100 ...

嗯，你这个方法比较简单，我也用matlab做了验证是对的结果。matlab没有这个内存限制，里面集成了这些函数稍微好用一下。谢谢你的指点