| 查看: 404 | 回复: 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 ] |
回复此楼» 猜你喜欢
085601材料工程找调剂
已经有11人回复
-
总分293求调剂
已经有9人回复
-
334分 一志愿武理 材料求调剂
已经有5人回复
-
327求调剂
已经有6人回复
-
329求调剂
已经有10人回复
-
275求调剂
已经有15人回复
-
291求调剂
已经有6人回复
-
309求调剂
已经有7人回复
-
297求调剂
已经有11人回复
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历
已经有6人回复
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
