| 查看: 474 | 回复: 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 ] |
回复此楼» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有13人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
有没有大神帮我看看基金代码
已经有9人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
【答案】应助回帖
感谢参与，应助指数 +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
hrbeu2219136
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 434.1
- 散金: 1888
- 帖子: 133
- 在线: 65.2小时
- 虫号: 636691
- 注册: 2008-10-25
- 专业: 固体力学
3楼2012-03-09 22:04:51
28