24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

求助用matlab软件编一个计算器已经有6人回复
求教，matlab并行运算加速已经有12人回复
求助使用MATLAB计算绿色叶片的覆盖度的问题已经有24人回复
matlab 计算过程中工作空间变量保存及读取求解决内存不足的方法已经有5人回复
如何用MATLAB 实现化学反应方程式（写程序代码）? 已经有20人回复
【素材】常用数值计算Matlab程序（M文件）已经有238人回复
看到讨论的多是Matlab,计算程序为什么不远c++编写呀？已经有4人回复
matlab工具箱中curve fitting中的图如何添加横纵坐标轴的名称已经有5人回复
关于MATLAB和OpenCV对相同算法的运算速度对比问题已经有20人回复
【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？已经有3人回复
【求助】MATLAB编写最大熵法计算概率密度的程序已经有8人回复

1楼2011-06-15 09:37:10

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

superstat54

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 3809.1
散金: 136
红花: 10
帖子: 256
在线: 146.8小时
虫号: 1327326
注册: 2011-06-20
专业: 催化化学

【答案】应助回帖

ll20100996(金币+2): 我用fortran编了，快一倍。呵呵 2011-06-22 08:56:39

高效的Matlab程序，尽量避免使用循环语句，更多的采用向量化的编程方式，即使用矩阵。修改如下：
clear;
tic
N=10;tn=4E4;a=3.905E-10;b=3.905E-10;c=4.2E-10;K0=0; n3=4-1;
         for n1=1:tn
            n2=1:n1-1;
            K01=2*a*b*c*(3*n3^2*c^2./(n1^2*a^2+n2.^2*b^2+n3^2*c^2)-1)./((n1^2*a^2+n2.^2*b^2+n3^2*c^2).^1.5*2*pi*8.85E-12);
            K0 = K0 - sum(K01);
            K0=K0-a*b*c*(3*n3^2*c^2/(n1^2*a^2+n1^2*b^2+n3^2*c^2)-1)/((n1^2*a^2+n1^2*b^2+n3^2*c^2).^1.5*2*pi*8.85E-12);
         end
         n1=1:tn;
         K02=a*b*c*(3*n3^2*c^2./(n1.^2*a^2+n3^2*c^2)-1)./((n1.^2*a^2+n3^2*c^2).^1.5*2*pi*8.85E-12);%n2=0 Or b=0
         K0 = K0 - sum(K02);
         K0=K0-a*b*c*(3*n3^2*c^2/(n3^2*c^2)-1)/((n3^2*c^2)^1.5*8*pi*8.85E-12)
toc;
对比：
tn=4E3，
原程序：
时间：2.550174 seconds，K0=-6.8388e+006；
修改后，时间：2.190174 seconds，K0=-6.8388e+006；

tn=4E4，
原程序，运行N久，没有等到结果就放弃了
修改后，时间155 s, K0 = -6.85e+005
tn=4E5时，你自己尝试一下吧，运行时间估计随tn成平方关系增长。