| 查看: 361 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
hrbeu2219136金虫 (小有名气)
|
[求助]
积分求助
|
|
|
如附件所示的积分，在一个计算程序中需要多次计算，m和n的上限为截断值，但截断值不一样的时候，调用计算积分点次数不一样，从而导致计算时间的长短，现在但m和n取为20的时候计算时间就很长了，请问大家有什么好的办法进行计算吗？ 现在采用这样的方式计算的 function y = fx1(m,m1,R,r0,t) FF = @(x)cos(lam(R,m)*x).*cos(lam(R,m1)*x)./(x+r0).^t; y = quad(FF,0,R); 请大家多多指点，谢谢  [ Last edited by hrbeu2219136 on 2012-3-9 at 11:30 ] |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有4人回复
-
基金委咋了？2026年的指南还没有出来？
已经有9人回复
-
基金申报
已经有5人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
hrbeu2219136
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 434.1
- 散金: 1888
- 帖子: 133
- 在线: 65.2小时
- 虫号: 636691
- 注册: 2008-10-25
- 专业: 固体力学
3楼2012-03-09 22:04:51
【答案】应助回帖
感谢参与，应助指数 +1
|
sin(mpix/R)sin(npix/R)=cos[(m+n)*pi*x/R]/2-cos[(m-n)*pi*x/R]/2,可以让双循环变单循环，不过此级数估计不收敛，不知你的截断是否有根据。 此积分实际是可积的，只不过积分结果中包含超越函数cosint(x) [见http://www.mathworks.cn/help/toolbox/symbolic/cosint.html]。 因此没必要用数值积分。 |
2楼2012-03-09 15:38:56