如何理解FFT

Posted 2023-04-11

老师让我用MATLAB来理解FFT，请问哪位大侠讲一下，谢谢

fft在matlab中可以按原理采样来运算，也可以用快速fft算法。
clear
fs=1000
t=0:1/fs:0.6;
f1=100;
f2=300;
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);
subplot(711)
plot(x);
title('f1（100Hz）\f2(300Hz)的正弦信号，初相0')
xlabel('序列（n）')
grid on

number=512

y=fft(x,number);
n=0:length(y)-1;
f=fs*n/length(y);
subplot(713)
plot(f,abs(y));
title('f1\f2的正弦信号的FFT（512点）')
xlabel('频率Hz')
grid on

x=x+randn(1,length(x));
subplot(715)
plot(x);
title('原f1\f2的正弦信号（含随机噪声）')
xlabel('序列（n）')
grid on

y=fft(x,number);
n=0:length(y)-1;
f=fs*n/length(y);
subplot(717)
plot(f,abs(y));
title('原f1\f2的正弦信号（含随机噪声）的FFT（512点）')
xlabel('频率Hz')
grid on 参考技术A 快速傅氏变换?
function X=myfft(x)
%myfft函数 用递归实现
N=length(x);
t=log2(N);
t1=floor(t);
t2=ceil(t);
if t1~=t2; %若x的长度N不为2的整数次幂，则补0至最接近的2的整数次幂
x=[x zeros(1,2^t2-N)];
N=2^t2;
end
w0=exp(-j*2*pi/N);
X=zeros(1,N);
if N==2
X(1)=x(1)+x(2);
X(2)=x(1)-x(2);
else
n=1:N/2;
xe(n)=x(2*n-1);
xo(n)=x(2*n);
XE=myfft(xe); %递归调用
XO=myfft(xo);
for n=1:N/2
X(n)=XE(n)+XO(n)*(w0^(n-1));
X(n+N/2)=XE(n)-XO(n)*(w0^(n-1));
end
end

如何使用 NAudio 对波形文件执行 FFT

【中文标题】如何使用 NAudio 对波形文件执行 FFT

【英文标题】：How to perform the FFT to a wave-file using NAudio

【发布时间】：2013-01-23 12:10:26

【问题描述】：

我正在使用 NAudio 库，并希望将快速傅立叶转换为 WaveStream。我看到 NAudio 已经内置了 FFT 但我该如何使用呢？

我听说我必须使用 SampleAggregator 类。

【参考方案1】：

您需要阅读整篇 blog article 以更好地理解我提取的以下代码示例，以确保即使文章没有保留该示例：

using (WaveFileReader reader = new WaveFileReader(fileToProcess))

    IWaveProvider stream32 = new Wave16toFloatProvider(reader);
    IWaveProvider streamEffect = new AutoTuneWaveProvider(stream32, autotuneSettings);
    IWaveProvider stream16 = new WaveFloatTo16Provider(streamEffect);
    using (WaveFileWriter converted = new WaveFileWriter(tempFile, stream16.WaveFormat))
    
        // buffer length needs to be a power of 2 for FFT to work nicely
        // however, make the buffer too long and pitches aren't detected fast enough
        // successful buffer sizes: 8192, 4096, 2048, 1024
        // (some pitch detection algorithms need at least 2048)
        byte[] buffer = new byte[8192]; 
        int bytesRead;
        do
        
            bytesRead = stream16.Read(buffer, 0, buffer.Length);
            converted.WriteData(buffer, 0, bytesRead);
         while (bytesRead != 0 && converted.Length < reader.Length);
    

但简而言之，如果您创建了 WAV 文件，则可以使用该示例将其转换为 FFT。