| 查看: 471 | 回复: 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人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有3人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
【答案】应助回帖
感谢参与，应助指数 +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
20