怎么用matlab做傅里叶变换

这样族仔用matlab做傅里叶变换：

1、我们使用matlab开发的傅立叶变换程序代码如下：

syms x

f = exp(-2*x^2) %our function

ezplot(f,[-2,2])  % plot of our function

FT = fourier(f) % Fourier transform

将其写入到我们的matlab程序模块中。

2、我们运行上面的傅立叶变换程序代码，兄镇将得出运行结果兆尘汪：FT = (2^(1/2)*pi^(1/2)*exp(-w^2/8))/2。

3、如果我们需要更高级的显示，我们修改上述代码即可，如使用ezplot(FT)作傅里叶变换折线图。

以上就是如何用matlab做傅里叶变换的解决步骤。

以下是一份MATLAB程序，使用筛法求含裤得小于10^8的雷打数：

```matlab

n = 1e8% 设置筛法范围

isPrime = true(n, 1)% 初始化数组为全部为true

isPrime(1) = false% 定义1不是素数

% 首先从2开始，将所有的素数的倍数标记为非素数

for i = 2:sqrt(n)

if isPrime(i)

isPrime(i * i : i : n) = false

end

end

% 输出小于10^8的雷打数

for i = 1:length(isPrime)

if mod(i, 4) == 1 &&isPrime(i)

fprintf('%d\n', i)

end

end

```

在程序中，首先定义了一个筛法范围 `n = 1e8`，然后利用逻辑数组 `isPrime` 来记录每个数是否为素数。初始化时，所有的数默认为素数，将数漏者组全部初始化为 `true`；同谈搜简时，定义数字1不是素数。之后，从2开始循环，将所有的素数的倍数标记为非素数。具体实现的方法是，如果检测到一个素数i，则将i * i及其之后的所有i的倍数都标记为非素数。

最后，再次循环判断并输出所有小于10^8的雷打数。因为所有的雷打数的形式为4k+1，所以对于满足此条件的素数，都可以视为雷打数。程序使用 `mod(i, 4) == 1` 检查当前数是否满足4k+1的形式，并使用 `isPrime(i)` 检查该数是否为素数。如果两个条件都满足，则输出该数即可。

需要注意的是，在程序中使用了 `sqrt(n)` 函数，以降低时间复杂度。因为n以内的素数最大可能为 `sqrt(n)`，所以只需要进行到 `sqrt(n)` 的筛法，即可得到全部小于n的素数。

用fft就能画时间数据锋陵猛的频谱图。

比如说你下载的地震波是一个传感器在一段时间内得到的数据，那么这组数据应该是个银桥向量，比如说设成向量x,

假设这组汪裂数据的采样频率是Fs，数据x的长度为L，设其频谱为X。Matlab程序为：

X=fft(x)

X_Mag=abs(X)

f=0:Fs/L:Fs-Fs/L

plot(f,

X_Mag)

这里X_Mag为频谱幅度，用plot画出频谱幅值图。f为对应的频率，画出的频谱从0到Fs,

由于得到的数据是离散的采样数据，所以频率是对称的，你只需要分析频率从0到Fs/2的频率就行了。

---