版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>虫友互识 (5791)
>论坛更新日志 (3656)
>留学DIY (3103)
>留学生活 (2462)
>海外博后 (1487)
>育儿交流 (1200)
>海归之家 (1000)
>医药类考试 (1000)
>材料综合 (799)
>人文社科 (798)
>信息科学 (760)
>复合材料 (658)
>地学 (640)
>文学芳草园 (601)
>转基因 (600)
>我与小木虫的故事 (600)

返回列表

Simewe

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 34
帖子: 13
在线: 1.6小时
虫号: 2615539
注册: 2013-08-25
专业: 代数学

[求助] Matlab中自定义的函数怎么能支持向量化操作?

CODE:

function r=primeQ(n)

r=1;

for i=2:sqrt(n);

    if mod(n,i)==0;

        r=0;

        return

    end

end;

上面的primeQ是个判断素数的函数，调用的时候只会这样用
sum( arrayfun(@primeQ, 2:1e6) )
听说arraygfun的速度不快，我希望可以支持向量化操作，就像下面这样
sum(primeQ(2:1e5))
应该怎么做？
注意：我知道判断素数有更好的算法，我关心的不是算法问题

回复此楼

» 猜你喜欢

欢迎发来filecode的Mz6后的代码验证其规律已经有56人回复
有时候，自然基金真的不能太认真（我的申报经验）已经有9人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有3人回复
关于Filecode分析方法已经有13人回复
时间戳又变了8-15 已经有9人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有8人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有9人回复
哪位老哥知道今年的国自然具体哪一天放榜？已经有6人回复

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

MATLAB-pdepe函数使用问题已经有5人回复
用matlab里的nlinfit如何拟合分段函数已经有8人回复
Matlab中如何自己定义一个新函数，然后计算时可以调用？已经有12人回复
如何在matlab模糊逻辑工具箱中创建自定义隶属度函数？已经有4人回复
matlab高手帮忙做个自定义复杂函数拟合，含有复数，已经有25人回复
求matlab画一个函数的图已经有6人回复
MATLAB函数速查手册已经有211人回复
matlab的函数调用问题已经有3人回复
问个关于matlab里find函数的问题！已经有4人回复
请教个matlab中的rand函数已经有3人回复
【讨论】matlab量化程序已经有4人回复
matlab做函数图象的问题已经有6人回复
【求助】求助：matlab下求y=x(e^abs(x)+1.3)的逆函数已经有9人回复
【求助】修改MATLAB支持向量机核函数程序遇到的疑问？求助已经有5人回复
Matlab 的自定义函数拟合已经有3人回复

1楼 2013-09-21 16:29:44

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

prs1022

铁虫 (初入文坛)

应助: 1 (幼儿园)
金币: 330.3
帖子: 42
在线: 17小时
虫号: 2704348
注册: 2013-10-07
性别: GG
专业: 计算机软件

【答案】应助回帖

你首先测试是否为偶数，如果不是，后边的偶数因子可以全部跳过，这样可以减少更多的时间
当然可以也对3做同样的操作，判断数不是2或3的倍数后只要判断那些6k+1和6k+5的因子
可以继续下去，但是越往后提升的效果就越差，所以自行斟酌

对于比较大的数可以考虑用Fermat小定理或者Miller-Rabin检验去判断；双精度最大连续整数上限是2^53

另外可以用matlab coder把m文件转换成mex，这样会更快

赞一下

回复此楼

stayhungry,stayfoolish

2楼2013-10-24 17:20:47

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

相关版块跳转我要订阅楼主 Simewe 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3	k0dTPqJtl0jt 2026-08-14	3/150	2026-08-15 12:41 by KxMI1BYBxWX1
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3	k0dTPqJtl0jt 2026-08-14	3/150	2026-08-15 12:33 by KxMI1BYBxWX1
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	k0dTPqJtl0jt 2026-08-14	5/250	2026-08-15 07:09 by 4wMiSEwB6436
[基金申请] filecode +7	cratir 2026-08-14	11/550	2026-08-15 06:19 by 学员8dgXkO
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	HFw0lei2R37i 2026-08-14	4/200	2026-08-15 04:17 by 4wMiSEwB6436
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4	k0dTPqJtl0jt 2026-08-14	4/200	2026-08-15 02:21 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +11	Tide man 2026-08-10	12/600	2026-08-15 02:12 by home3163
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	k0dTPqJtl0jt 2026-08-14	4/200	2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8	布布和一二 2026-08-10	11/550	2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23	Lanmanbaby 2026-08-09	37/1850	2026-08-14 13:40 by Equinoxhua
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10	家与远方 2026-08-10	15/750	2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4	汪汪锅 2026-08-09	4/200	2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源：本周末出结果 +10	瞬息宇宙 2026-08-12	10/500	2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7	且听虎啸 2026-08-12	9/450	2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4	布布和一二 2026-08-11	5/250	2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11	Tide man 2026-08-11	11/550	2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6	瞬息宇宙 2026-08-10	7/350	2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3	Tide man 2026-08-11	3/150	2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4	Vierhys 2026-08-10	8/400	2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5	西山十月 2026-08-09	7/350	2026-08-10 07:32 by 仁砚薪传