|
|
@[TOC](Matlab——快速查找素?cái)?shù)（0-50）)
# 1.分析
**<font color='red' size=4>素?cái)?shù)**就是質(zhì)數(shù)，素?cái)?shù)有無(wú)限個(gè)。一個(gè)大于1的自然數(shù)，且除了1和它本身外，不能被其他自然數(shù)整除的數(shù)叫素?cái)?shù)。換句話(huà)說(shuō)就是，除了1和該數(shù)本身以外不再有其他的因數(shù)的數(shù)被稱(chēng)為素?cái)?shù)。
# 2.源程序代碼
代碼如下：
````matlab
%% 自定義一段算法代碼, 求 50 以?xún)?nèi)的素?cái)?shù)
clear;clc;close all
ss=[]; % ss 用來(lái)標(biāo)記是否為素?cái)?shù)
pr=[2 3]; % pr 用來(lái)存放素?cái)?shù), 2 和 3 均為素?cái)?shù), 先存放在 pr 中
for k=4:50 % 查找 4~50 里的素?cái)?shù)
for m=2:fix(sqrt(k))
if mod(k,m)==0 % 判斷 k 是否可以被 m 整除, 如果能被整除, 說(shuō)明 k 不是素?cái)?shù), 用 ss=0 來(lái)標(biāo)記
ss=0;
break % 跳出內(nèi)循環(huán)
else
ss=1;
end
end
if ss==1 % 如果 k 是素?cái)?shù), 則保存至 pr
pr=[pr,k];
end
end
%% 使用 MATLAB 自帶的 primes 函數(shù), 直接求 50 以?xún)?nèi)的素?cái)?shù)
p=primes(50);
%% 比較兩種方法的結(jié)果是否一致
pr
p
````
# 3.運(yùn)行結(jié)果
運(yùn)行結(jié)果如下：
感謝閱讀，如有不足之處，歡迎來(lái)指正。謝謝！
|
|
