[DSP学习笔记]基于TMS320F28335的FFT及加窗函数实现

一、新建工程

首先我们先将C2000 wave中我们所需要的函数库给导入到我们的新建工程中(可见于我实现FIR滤波文章中导入函数库的操作)。

工程中可见有以下文件。(仅FFT可删除滤波有关文件)

接着新建我们所需的文件INPUT.c于Window.c。

二、输入信号 

借助MATLAB我们来生成所需的输入信号于窗函数

生成输入信号的MATLAB程序如下:

%生成固定点数固定频率的正弦信号
%fs为采样率
%f为多频信号频率
%N为采样点数
%amp为信号幅值
%phase为信号初相位
function [input,t] = mult_freq_signal_test(fs,f,N,amp,phase)
    dt = 1/fs;        %采样周期
    t = dt*N;         %采样时间
    temp = (0:dt:(N-1)*dt)';
    y = amp.*sin(2*pi*f.*temp+phase/180*pi);
    Num = size(f,2);    %获取f频率个数
    input = y*ones(Num,1);
    subplot(1,2,1);
    plot(temp/dt,input);
    title('信号时域');
    xlabel('点数 N');ylabel('幅值 A');
    subplot(1,2,2);
    Y = fftshift(fft(input));%做傅里叶变换并把零频点移到频谱中心
    F = (-N/2:N/2-1)*(fs/N);%设置横坐标
    plot(F,abs(Y)/max(abs(Y)));
    title('信号频域');
    xlabel('频率 f');ylabel('幅值 A')
end
%与mult_freq_signal_test.m函数式文件配合使用,将生成的输出信号保存到本地
fs=8e4;             %设采样率为80K
f = [1e4,1.5e4,3e4];  %频率矩阵(1*N)测试信号频率为10K,15K,30K
phase = [20,0,40];    %初相位分别为为20,0,40
amp = [0.5/3*4095,0.3/3*4095,0.4/3*4095];%对应幅值为0.5V,0.3V,0.4V信号的采集值(ADC为0-3V)
N=1024;              %采样点数

y=mult_freq_signal_test(fs,f,N,amp,phase);

y=reshape(y,[4,N/4]);
y=y+1.2/3*4095;       %用于增加直流偏置,无需直流偏置时可注释,直流偏置为1V
y=y';

%writematrix(y,'record.txt','delimiter',',');%将生成的初始信号保存到.m文件路径下的record.txt

fileID = fopen('record.txt','w');
for i=1:1:N/4
    fprintf(fileID,'%.6f,%.6f,%.6f,%.6f,\n',y(i,1),y(i,2),y(i,3),y(i,4));
end
fclose(fileID);%将生成的初始信号保存到.m文件路径下的record.txt

将record.txt的文件内容复制到INPUT.c中,如下:

//#############################################################################
//! \file   input.c
//! \brief  Input Vector (2048)
//! \author Vishal Coelho
//! \date   19-Jan-2023
//!
//
//  Group:          C2000
//  Target Family:  $DEVICE$
//
//#############################################################################
//
//
// $Copyright: Copyright (C) 2020 Texas Instruments Incorporated -
//             http://www.ti.com/ ALL RIGHTS RESERVED $
//#############################################################################

#include <stdint.h>


#pragma DATA_SECTION(FFT_input, "FFT_buffer3")


// Matrix is 2 times its required sizeas it will be used to store
// 1024 floats at a later point
uint16_t FFT_input[2048] = {
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
                            2222.390783,2682.470769,2150.640175,2313.650403,
                            644.109217,815.149034,1704.480279,1183.969400,
                            2222.390783,1925.813432,1571.519721,2627.068134,
                            1463.109217,1128.566765,1125.359825,427.312063,
};




// End of File

三、窗函数

接下我们进行窗函数的生成。

1.加窗之前,首先我们要明白我们为什么要加窗。

在实际应用中,我们采集到的信号通常是较长的非周期时间序列,在进行FFT时,需要对信号进行截断。截断后得到的信号通常是非周期信号,所以FFT结果就变成了连续的谱线,周期信号的频率-幅值不再是一条垂直的直线,而是一条曲线,同时幅值下降,泄露到曲线其他频率点,即我们常说的频谱泄露。而窗函数一个周期函数,窗函数与截断信号相乘后,会变成比较接近周期函数的波形(实际上就是利用窗函数进行调制)。加窗后的信号进行FFT后,得到的频谱曲线会更窄,趋向于垂直直线。因此对频谱泄露有一定改善。

2.知道为什么要加窗后,我们再来根据实际情况来选择我们所需要的窗函数。

通常情况下,我们不对输入信号进行处理,而直接截断情况下相当于给信号加了一个矩形窗。

如果我们对于FFT后的结果要求精确的幅值,平顶窗是我们的首选。

如果我们对应FFT后的结果要求精确的频率点,则可以选择汉宁窗。

如果检测两个频率相近、幅度不同的信号,可以选用布莱克曼窗。

如果对应频率于幅值均有一定要求,则可以选用凯赛窗。

3.选择完窗函数,我们还要对FFT后的结果乘以他们对应的恢复系数(窗函数对信号有衰减),如下表:

 了解完以上知识,我们便可以继续我们的实验。

在MATLAB命令行窗口敲入window,进入下列界面并完成下列配置。(本次测试选用汉宁窗)

点击文件,导出

导出到工作区,我使用默认的变量名。我们便在工作区得到了从我们所需要的窗函数。

 但是,我们工程实际应用,还需对得到的窗函数进行处理,将其转化为Q16型。

官方例程中提出,如下

通过编写MATLAB文件实现,并保存到本地的window.txt文件中。

%将窗设计器生成的窗系数打印到window.txt文件中

fileID = fopen('window.txt','w');
N=1024;%生成窗点数

Win=reshape(Win',N,1);
Windows=uint16(floor(Win.*(2^16))); % convert to Q16

Windows=reshape(Windows,[4,N/4]);
Windows=Windows';
for i=1:1:N/8
    fprintf(fileID,'%d,%d,%d,%d,\n',Windows(i,1),Windows(i,2),Windows(i,3),Windows(i,4));
end
fclose(fileID);%将生成的初始信号保存到.m文件路径下的window.txt

 将window.txt的内容复制到window.c中,如下:

//#############################################################################
//! \file   window.c
//! \brief  window Vector (512)
//! \author yang
//! \date   19-Jan-2023
//! 此窗函数为汉宁窗,由MATLAB生成,加窗后FFT结果需乘以恢复系数,汉宁窗恢复系数为2
//
//  Group:          C2000
//  Target Family:  $DEVICE$
//
//#############################################################################
#include "dsp.h"

const uint16_t Window[512]={
                            0,0,2,5,
                            9,15,22,30,
                            39,50,61,74,
                            88,104,121,138,
                            158,178,200,222,
                            246,272,298,326,
                            355,385,416,449,
                            483,518,554,592,
                            630,670,711,754,
                            797,842,888,935,
                            983,1033,1084,1136,
                            1189,1243,1299,1355,
                            1413,1472,1533,1594,
                            1657,1720,1785,1851,
                            1919,1987,2057,2128,
                            2199,2273,2347,2422,
                            2499,2576,2655,2735,
                            2816,2898,2982,3066,
                            3152,3238,3326,3415,
                            3505,3596,3688,3782,
                            3876,3972,4068,4166,
                            4265,4364,4465,4567,
                            4670,4774,4880,4986,
                            5093,5201,5311,5421,
                            5532,5645,5758,5873,
                            5988,6105,6222,6341,
                            6460,6581,6702,6825,
                            6948,7072,7198,7324,
                            7451,7580,7709,7839,
                            7970,8102,8235,8369,
                            8504,8640,8776,8914,
                            9052,9192,9332,9473,
                            9615,9758,9901,10046,
                            10191,10338,10485,10633,
                            10782,10931,11082,11233,
                            11385,11538,11692,11846,
                            12002,12158,12315,12472,
                            12631,12790,12950,13110,
                            13272,13434,13597,13760,
                            13924,14090,14255,14422,
                            14589,14757,14925,15094,
                            15264,15434,15606,15777,
                            15950,16123,16296,16471,
                            16646,16821,16997,17174,
                            17351,17529,17708,17887,
                            18066,18246,18427,18608,
                            18790,18972,19155,19339,
                            19522,19707,19891,20077,
                            20263,20449,20636,20823,
                            21010,21198,21387,21576,
                            21765,21955,22145,22336,
                            22527,22718,22910,23102,
                            23295,23487,23681,23874,
                            24068,24262,24457,24652,
                            24847,25042,25238,25434,
                            25630,25827,26024,26221,
                            26418,26616,26813,27011,
                            27210,27408,27607,27806,
                            28005,28204,28403,28603,
                            28802,29002,29202,29402,
                            29603,29803,30003,30204,
                            30405,30606,30806,31007,
                            31208,31409,31611,31812,
                            32013,32214,32415,32617,
                            32818,33019,33220,33422,
                            33623,33824,34025,34226,
                            34427,34628,34829,35030,
                            35231,35431,35632,35832,
                            36033,36233,36433,36633,
                            36832,37032,37232,37431,
                            37630,37829,38028,38226,
                            38425,38623,38821,39018,
                            39216,39413,39610,39807,
                            40003,40199,40395,40591,
                            40786,40981,41176,41370,
                            41564,41758,41951,42144,
                            42337,42529,42721,42913,
                            43104,43294,43485,43675,
                            43864,44054,44242,44431,
                            44618,44806,44993,45179,
                            45365,45551,45736,45921,
                            46105,46288,46471,46654,
                            46836,47017,47198,47379,
                            47559,47738,47917,48095,
                            48272,48449,48626,48802,
                            48977,49151,49325,49499,
                            49672,49844,50015,50186,
                            50356,50526,50694,50862,
                            51030,51197,51363,51528,
                            51693,51857,52020,52182,
                            52344,52505,52665,52825,
                            52984,53142,53299,53455,
                            53611,53766,53920,54073,
                            54226,54378,54529,54679,
                            54828,54976,55124,55271,
                            55416,55561,55706,55849,
                            55991,56133,56273,56413,
                            56552,56690,56827,56963,
                            57099,57233,57366,57499,
                            57630,57761,57891,58020,
                            58147,58274,58400,58525,
                            58649,58772,58894,59015,
                            59135,59254,59372,59489,
                            59605,59720,59834,59947,
                            60058,60169,60279,60388,
                            60496,60602,60708,60813,
                            60916,61019,61120,61221,
                            61320,61418,61515,61611,
                            61706,61800,61893,61985,
                            62075,62165,62253,62340,
                            62426,62511,62595,62678,
                            62760,62840,62919,62998,
                            63075,63151,63226,63299,
                            63372,63443,63513,63582,
                            63650,63717,63782,63847,
                            63910,63972,64033,64092,
                            64151,64208,64264,64319,
                            64373,64425,64477,64527,
                            64576,64624,64670,64716,
                            64760,64803,64844,64885,
                            64924,64962,64999,65035,
                            65069,65102,65134,65165,
                            65195,65223,65250,65276,
                            65301,65324,65346,65367,
                            65387,65406,65423,65439,
                            65454,65467,65480,65491,
                            65501,65509,65517,65523,
                            65528,65532,65534,65535,
};

四、FFT以及计算幅值

接下来进行编写对FFT进行初始化的函数。

void RFFT_Init()
{
    RFFT_f32_setInputPtr(hnd_rfft, (float *)FFT_input);
    RFFT_f32_setOutputPtr(hnd_rfft, FFT_output);
    RFFT_f32_setStages(hnd_rfft,  FFT_STAGES);
    RFFT_f32_setFFTSize(hnd_rfft, FFT_SIZE);
    RFFT_f32_setMagnitudePtr(hnd_rfft, FFT_magnitude_phase);
    RFFT_f32_setPhasePtr(hnd_rfft, FFT_magnitude_phase);

    RFFT_f32_setTwiddlesPtr(hnd_rfft, FFT_twiddleFactors);
    RFFT_f32_sincostable(hnd_rfft);

    RFFT_ADC_f32_setTailPtr(hnd_rfft_adc, &(hnd_rfft->OutBuf));
    RFFT_ADC_f32_setInBufPtr(hnd_rfft_adc, FFT_input);

    RFFT_adc_f32_win(FFT_input, Window, FFT_SIZE );
}

进行FFT以及对FFT结果进行计算的函数。

void RFFT_Caculate()
{
    float* mag;

    RFFT_adc_f32(hnd_rfft_adc);
    RFFT_f32_mag(hnd_rfft);     //计算幅度
    mag=RFFT_f32_getMagnitudePtr(hnd_rfft);

    RFFT_Analyse(mag);
}

需注意的是,FFT的返回值不是幅值,我们需对其进行转换。

首先对FFT的结果进行取模,使用RFFT_f32_mag()函数,得到幅度。

接着对幅度再进行相应的处理,如下表:

 即,直流分量的幅值等于幅度除以采样点数,其余的的幅值等于幅值除以采样点数的一半。

对幅度进行分析的函数

void RFFT_Analyse(float *src)
{
    Uint16 i;
    float FFT_Value;        //FFT_Value为该频率点的FFT后的幅度
    Uint16 peak_index[RFFT_Analyse_MaxNum];//存放的峰值个数为RFFT_Analyse_MaxNum个

    memset(peak_index, 0U, RFFT_Analyse_MaxNum*sizeof(Uint16));     //初始化将指定区块为0值
    memset(rfft_analyse.freq, 0U, RFFT_Analyse_MaxNum*sizeof(float));
    memset(rfft_analyse.Amp, 0U, RFFT_Analyse_MaxNum*sizeof(float));

    rfft_analyse.Wave_Num=RFFT_Peak_Find(src,FFT_SIZE/2,peak_index);

    for(i=0;i<rfft_analyse.Wave_Num;i++)
    {
        rfft_analyse.freq[i]=(float) Fs / FFT_SIZE * peak_index[i];//计算频率,已调试过,为准确频率
        FFT_Value=*(src+peak_index[i]);
        rfft_analyse.Amp[i] = 2 * FFT_Value / FFT_SIZE * 3 * 2;   //将幅度转换为幅值公式:幅值=(幅度/(采样点数/2))该点不为FFT后结果的第一个点和最后一个点,乘3是ADC的转化范围为0-3V,乘2是汉宁窗的幅度恢复系数
    }
    rfft_analyse.DC_Amp=src[0]/FFT_SIZE * 3 * 2;                  //直流分量的幅值=幅度/采样点数,乘2是汉宁窗的幅度恢复系数
}

在主函数中,对RFFT_Init()和RFFT_Caculate()完成调用,即可。

五、核验结果

准备阶段结束,进入仿真。

仿真结果如下:

频率点:

 对应幅值:

 直流分量:

 可以看到加窗后FFT的结果于输入信号相吻合,测试结束。

需要工程文件的可以私信我,本文如有错误,也请多多指出。

本文编写时,参考了以下相关文章:

【1】窗函数总结_wxyczhyza的博客-CSDN博客_窗函数衰减​​​​​​【2】TMS320F28335调用官方库进行FFT频谱分析_PeepFuture橙子的博客-CSDN博客

  • 5
    点赞
  • 44
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
No.4 简单的MFC多对话框演示程序 “MultiDialog” 演示如何在工程中拥有多个对话框,及如何在一个对话框中调用另一个对话框。 重点:1、多个对话框类的建立;2、对话框的模式(Modal)调用方法;3、Spin控件的使用。 新建一个基于对话框的MFC工程, 通过菜单"Insert>>Resource"打开添资源对话框, 在其中选择Dialog后点击New创建一个新的对话框 按Ctrl+W打开ClassWizard,系统会提示刚才创建了一个新的对话框资源,是否建立对应的类, 选择建立,然后在New Class窗口中Class Name栏输入它的名称:CSubClass1, 确定后系统会自动生成SubClass1.h和SubClass1.cpp并入工程中,其中有已经创建好的CSubClass1的类的基本代码。 把这个对话框的Caption属性改为“难度选择”,在它上面画三个Radio“简单”、“标准”、“困难”, 并建立相关联的变量m_Option1。(要注意的是在ClassWizard中注意Class Name中应该选CSubClass1而不是之前的主对话框类) 按照相同方法建立第二个新对话框,类名“CSubClass2”,Caption为“关卡选择”。 在上面画一个Edit和一个Spin,注意先画Edit后画Spin,将Spin的Auto buddy和Set buddy integer勾上。 按Ctrl+W打开ClassWizard,为Edit建立关联变量,不过注意是int型而不是CString型, 也为Spin建立关联变量m_Spin1,注意这次是Control型变量CSpinButtonCtrl。 下面为这两个对话框添代码。 双击“难度选择”对话框的OK按钮,建立对话框的OnOK映射。 在其中入(在CDialog::OnOK();之前): UpdateData(TRUE); if ((m_Option1>2) ||(m_Option1<0)) { MessageBox("错误的选择!","提示",MB_OK); return; } 下面对“关卡选择”对话框添初始化代码, 由于关卡的有效值只有1到6,因此需要在初始化时设置Spin控件的有效值范围。 按Ctrl+W打开ClassWizard,在左侧列表选择这个对话框类CSubDialog2, 在右侧列表中选择WM_INITDIALOG,点击右边的“Add Function...”按钮, 接着点击右边的“Edit Code”按钮,在其中中入(在CDialog::OnInitDialog()那句之后,在return TRUE那句之前): m_Spin1.SetRange(1,6); m_Text1=1; m_Spin1.SetPos(1); 其中CSpinButtonCtrl::SetRange()函数的作用是设置和他关联的Spin控件的范围,两个参数分别是下界和上界。 而CSpinButtonCtrl::SetPos()是设定Spin的当前位置。 两个新的对话框都已建立完毕,下来是如何在主对话框中使用的问题。 首先,两个新对话框都有各自的类,分别在SubDialog1.h和SubDialog2.h中有定义。(类的细节则在对应的cpp中定义) 因此,主对话框想要调用这两个新对话框,需要先包含这两个头文件, 在你要使用的地方(本例是MultiDialogDlg.cpp中)文件前面上 #include "SubDialog1.h" #include "SubDialog2.h" 然后在想要调用的地方就可以使用了。 本例中,首先为主窗口的两个Edit建立CString型关联变量m_Text1和m_Text2, 然后在两个按钮的消息映射函数中分别入: CSubDialog1 dialog1; //定义CSubDialog1型对话框的一个新对象 dialog1.DoModal(); //使用“模式”调用,显示对话框 m_Text1.Format("%d",dialog1.m_Option1); //此句在上面对话框没有关闭前不会执行到 UpdateData(FALSE); 和 CSubDialog2 dialog1; dialog1.DoModal(); m_Text2.Format("%d",dialog1.m_Text1); UpdateData(FALSE); 其中第一句均为定义对话框新实例的语句,定义一个你想要的类型的对话框。 第二句是通过调用CDialog::DoModal()方法,来显示这个对话框,并进入“模式”(Modal)状态 在“模式”状态,当子对话框没有关闭之前,调用它的父对话框不能被响应, 并且其语句执行会停留在刚才的DoModal语句上等待,直到子对话框关闭才接着执行下一个语句。 第三第四句将子对话框得到的数据(即类的成员变量)显示在父对话框的Edit上。 四句执行完后退出该函数,这时刚才定义的CSubDialog1等对话框类变量被销毁,因此创建的话框也被销毁。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值