| 查看: 400 | 回复: 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 ] |
回复此楼» 猜你喜欢
352分-085602-一志愿985
已经有3人回复
-
291求调剂
已经有4人回复
-
330分求调剂
已经有4人回复
-
309求调剂
已经有4人回复
-
0856求调剂
已经有13人回复
-
0703化学
已经有11人回复
-
329求调剂
已经有7人回复
-
321求调剂
已经有8人回复
-
0856材料化工调剂 总分330
已经有12人回复
-
070300化学354求调剂
已经有4人回复
【答案】应助回帖
感谢参与，应助指数 +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
6