24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 401  |  回复: 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 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cars

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hrbeu2219136

金虫 (小有名气)

引用回帖:
: Originally posted by cars at 2012-03-09 15:38:56:
sin(mpix/R)sin(npix/R)=cos/2-cos/2,可以让双循环变单循环,不过此级数估计不收敛,不知你的截断是否有根据。
此积分实际是可积的,只不过积分结果中包含超越函数cosint(x)

因此没必要用数值积分。

我的截断是根据整个程序得来的。是的,我用mathematica计算出来有超越函数,但是带进数值后计算出来有复数,不知道这个怎么解决?
3楼2012-03-09 22:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hrbeu2219136 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 材料学硕333求调剂 +11 北道巷 2026-03-24 11/550 2026-03-29 08:31 by fmesaito
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-03-26 6/300 2026-03-29 01:11 by hanserlol
[硕博家园] 求调剂 有机化学考研356分 +3 Nadiums 2026-03-25 4/200 2026-03-29 00:50 by 544594351
[考研] 316求调剂 +7 江辞666 2026-03-26 7/350 2026-03-28 21:28 by sanrepian
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +7 幸运哩哩 2026-03-22 11/550 2026-03-28 20:27 by 唐沐儿
[考研] 材料277求调剂 +7 min3 2026-03-24 7/350 2026-03-28 11:39 by xuxiang
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 330一志愿中国海洋大学 化学工程 085602 有读博意愿 求调剂 +3 wywy.. 2026-03-27 4/200 2026-03-28 03:32 by fmesaito
[考研] 一志愿南师大0703化学 275求调剂 +4 Ripcord上岸 2026-03-27 4/200 2026-03-27 17:00 by zhyzzh
[考博] 26申博 +3 加油冲啊! 2026-03-26 3/150 2026-03-27 15:38 by cls512
[考研] 307求调剂 +8 超级伊昂大王 2026-03-24 9/450 2026-03-27 15:34 by 超级伊昂大王
[考研] 一志愿郑大085600,310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[硕博家园] 北京林业大学硕导招生广告 +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
[考研] 一志愿陕师大生物学071000,298分,求调剂 +5 SYA! 2026-03-23 5/250 2026-03-27 09:29 by 不吃魚的貓
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +7 niko- 2026-03-22 7/350 2026-03-25 20:14 by qingfeng258
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
信息提示
请填处理意见