@[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.運行結果
運行結果如下：
感謝閱讀，如有不足之處，歡迎來指正。謝謝！
|
