在利用 fft 变换分析频谱时,常要考虑将采样时间和信号的周期设为互质的关系。为了方便的找到所选的质数,我们可以利用 matlab 自带的质数表函数 primes() 来加以查找。
下面所示为利用 primes 函数求得100以内的质数:
>> primes(100)
ans =
Columns 1 through 19
2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67
Columns 20 through 25
71 73 79 83 89 97
除了质数表函数 primes() 之外,另一种更为灵活的方式是利用检测质数的函数 isprime(),下面的方法就算出了 150 到 200 的质数:
>> A=150:200;
>> B=A(isprime(A))
B =
151 157 163 167 173 179 181 191 193 197 199