24小时热门版块排行榜

>论坛更新日志 (4540)
>虫友互识 (516)
>文献求助 (439)
>休闲灌水 (396)
>导师招生 (388)
>考博 (135)
>硕博家园 (131)
>博后之家 (102)
>基金申请 (101)
>招聘信息布告栏 (100)
>论文投稿 (94)
>考研 (85)
>教师之家 (65)
>找工作 (48)
>公派出国 (48)
>绿色求助(高悬赏) (46)

返回列表

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

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有11人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐博导已经有4人回复
求推荐英文EI期刊已经有5人回复
疑惑？已经有5人回复

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

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

返回列表