@[TOC](Matlab——快速查找素数(0-50))
# 1.分析
**<font color='red' size=4>素数**就是质数,素数有无限个。一个大于1的自然数,且除了1和它本身外,不能被其他自然数整除的数叫素数。换句话说就是,除了1和该数本身以外不再有其他的因数的数被称为素数。
# 2.源程序代码
代码如下:
````matlab
%% 自定义一段算法代码, 求 50 以内的素数
clear;clc;close all
ss=[]; % ss 用来标记是否为素数
pr=[2 3]; % pr 用来存放素数, 2 和 3 均为素数, 先存放在 pr 中
for k=4:50 % 查找 4~50 里的素数
for m=2:fix(sqrt(k))
if mod(k,m)==0 % 判断 k 是否可以被 m 整除, 如果能被整除, 说明 k 不是素数, 用 ss=0 来标记
ss=0;
break % 跳出内循环
else
ss=1;
end
end
if ss==1 % 如果 k 是素数, 则保存至 pr
pr=[pr,k];
end
end
%% 使用 MATLAB 自带的 primes 函数, 直接求 50 以内的素数
p=primes(50);
%% 比较两种方法的结果是否一致
pr
p
````
# 3.运行结果
运行结果如下:
感谢阅读,如有不足之处,欢迎来指正。谢谢!
|