函数(年度)= FFT (y, Fs)

1)计算功率谱密度和振幅谱(P(f)， f (f))
对一维信号y(t)，采样率Fs(奈奎斯特率)，其先验值为%。结果被绘制成3个图，对应于简单
分别为对数PSD、对数PSD (dB)和振幅谱。
＿＿＿＿＿＿＿＿＿＿＿＿＿
振幅(f) = \/ PSD(f)

这个函数的有用之处在于对频率轴的调整。

3)利用Matlab内置函数计算快速傅里叶变换
Fft，但对于长度<1000点的信号，可以使用嵌套的
函数y = Fast_Fourier_Transform (X, N)。

演示:

Fs = 800;
Tf = 2;
t = 0:1 / Fs: Tf;
f = 40 [75];
Amp = (4.5 - 9.22);
σ= 1.33;
y = Amp (1) * exp (j * 2 *π* t * f (1))
+ Amp (2) * exp (j * 2 *π* t * f (2));
N =σ/√(2))* (randn(大小(t)) + j * randn(大小(t)));
y = y + N;
人物,情节(t、y),包含(时间(s)), ylabel(“电压(v)”),
strcat('信号损坏与AWGN， \sigma='，num2str(sigma)))，
财政年度= FFT (y, Fs);

在m文件Demo_FFT中:
第一部分:计算频率为Fc的正弦信号Y(t)的频谱
第二部分:FFT[Y²(t)]

样品率Fs>=4*Fc进行调整。