24小时热门版块排行榜

返回列表

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))
应该怎么做？
注意：我知道判断素数有更好的算法，我关心的不是算法问题

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复
浙江师范大学是怎么坑我的已经有4人回复
闲聊已经有6人回复
江苏省自然基金什么时候出结果已经有7人回复
浙江师范大学是怎么坑我的已经有10人回复
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？已经有6人回复
面上有专家说实验设备不是我们单位的已经有8人回复
科研求助已经有4人回复

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

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 的主题更新

返回列表