24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1069  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 评审有感 +7 popular289 2026-05-18 12/600 2026-05-19 05:02 by redmoonzpc
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +3 光伏爱好者 2026-05-17 4/200 2026-05-19 01:19 by 风起于浮萍
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 7/350 2026-05-17 19:42 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:37 by ue3ir18jc3
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见