Znn3bq.jpeg
查看: 1040  |  回复: 3

jollage

银虫 (正式写手)

[求助] 如何让mathematica把这个程序运行快一点

我想试图求文件里的alpha_2。但是好像在那一个积分的步骤,程序运行很慢。我没有等到结果,哪位高手能不能看一下什么方法可以加快运行的速度。题目应该是对的。谢谢。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : xiaomuchong.nb
  • 2013-06-23 16:10:06, 2.83 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

walk1997

金虫 (著名写手)

★ ★ ★ ★ ★ ★
jollage(xzhdty代发): 金币+6, 代楼主发放金币 2013-06-25 07:15:09
数值结果 不知道对否
Clear["Global`*"];
a1 = 1/9 (7 - 12 \[Alpha]2);
a2 = 1/36 + (19 \[Alpha]2)/6;
\[Alpha]1 = 1/3 + 4 \[Alpha]2;
fint = (k - (2 a2*Sin[2 k] + 2 a1*Sin[k])/(1 + 2 \[Alpha]1*Cos[k] +
        2 \[Alpha]2*Cos[2 k]))^2;
temp = D[fint, \[Alpha]2];
test[\[Alpha]\[Alpha]2_?NumericQ] :=
  temp /. {\[Alpha]2 -> \[Alpha]\[Alpha]2};
test2[\[Alpha]\[Alpha]2_?NumericQ] :=
NIntegrate[test[\[Alpha]\[Alpha]2], {k, 0, Pi/2}]
FindRoot[test2[\[Alpha]2], {\[Alpha]2, 0.1}]
FindRoot[test2[\[Alpha]2], {\[Alpha]2, -0.1}]
2楼2013-06-23 21:11:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jollage

银虫 (正式写手)

引用回帖:
2楼: Originally posted by walk1997 at 2013-06-23 21:11:13
数值结果 不知道对否
Clear;
a1 = 1/9 (7 - 12 \2);
a2 = 1/36 + (19 \2)/6;
\1 = 1/3 + 4 \2;
fint = (k - (2 a2*Sin + 2 a1*Sin)/(1 + 2 \1*Cos +
        2 \2*Cos))^2;
temp = D;
test :=
  temp /. ...

好的,谢了,好像是对的,这个没有正确答案,但是和别的结果是吻合的。

是因为我其实是在用symbolic的方法去求所以求的很慢?我没有意识到这一点,用数值就快很多了。而且你把积分和微分的顺序换了,这样也令问题好解了一点。
谢了

你要这些金币吗?前面一个帖子里你也没有应助。
3楼2013-06-23 22:07:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
想速度快,除了尽量避免用符号运算,你还可以看看Compile这个函数。
4楼2013-06-24 17:02:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jollage 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 请问地理、遥感方面,可以做哪些横向项目啊,纵向完不成考核啊 +3 锦衣卫寒战 2026-04-07 5/250 2026-04-11 20:51 by 豫椒
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 295求调剂 +3 桂秋二十八 2026-04-05 5/250 2026-04-11 11:36 by zhq0425
[考研] 化学308分求调剂 +22 你好明天你好 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[硕博家园] 新一代电子信息294求调剂 不挑学校 +6 Ytyt11 2026-04-09 7/350 2026-04-11 10:52 by AA小小木虫
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 求调剂 +15 张zic 2026-04-05 16/800 2026-04-10 08:12 by kangsm
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 278求调剂 +27 范婷娜 2026-04-07 31/1550 2026-04-09 20:49 by zhouxiaoyu
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 本科生物信息学,总分362 求07 08调剂 +6 q小倩1210 2026-04-06 6/300 2026-04-07 19:40 by macy2011
[考研] 一志愿南科大生物学297分,求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 319分085702安全工程求调剂 +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
信息提示
请填处理意见