soberprogress 发表于 2008-9-21 10:40

求助:一个matlab广义积分的编程?

求助:一个matlab广义积分的编程?我遇到比较复杂的函数,为了清晰,我把被积函数f 分为几个部分,假定:
f1=f1(k,s1,s2);f2=f2(k,s1,s2);f3=f3(k,s1,s2);
被积函数 f=f1*f2*f3,
要求完成f对s1,s2在0~无穷大区间上的2重积分。其中k是参数,作后要求做出积分结果和k的曲线。

希望matlab高手给予指点,已经考虑很久了!谢谢!如果有相似的实例当然更好!

sigma665 发表于 2008-9-21 10:51

dblquad 应该可以解决

soberprogress 发表于 2008-9-21 10:57

回复 沙发 sigma665 的帖子

关键是积分区间是无穷大啊?:funk:

sigma665 发表于 2008-9-21 14:09

无穷大用inf或者取个足够大的值

soberprogress 发表于 2008-9-21 14:24

回复 地板 sigma665 的帖子

好像这不够严谨!

sigma665 发表于 2008-9-21 17:49

回复 5楼 soberprogress 的帖子

数值解,只要有足够的精度就行了

想严谨,用迭代
对积分区间中的b进行迭代
每次的b得到一个积分值
而后b+1,得到新积分值
比较2个差别,满足精度即认为求解

soberprogress 发表于 2008-9-21 18:59

回复 6楼 sigma665 的帖子

谢谢你,思想我倒是比较清楚的,我希望有个恰当的例子存在!谢谢!

sigma665 发表于 2008-9-22 07:39

回复 7楼 soberprogress 的帖子

已经说到这了,其他的自己试试吧

soberprogress 发表于 2008-9-22 09:31

回复 8楼 sigma665 的帖子

谢谢你的答复!不好意思占用了你很多时间!谢谢!

solitarytear 发表于 2008-12-25 14:45

可否利用FFT算法,先将k进行离散,再利用2维FFT计算得到f,最后画图。
不知道这样是不是可行,希望讨论!
页: [1]
查看完整版本: 求助:一个matlab广义积分的编程?