使用方法:
1、首先打开MATLAB ,点击 导入
在弹出的界面选择文件后弹出如下界面
然后点击 右上的 绿√ 即可导入所选内容。这里注意数值必须不能出现文字,否则计算会出错。
最后在命令界面输入如下代码即可。
%% FFT
%clear;clc;close all
% import ACC_X
Fs=100; % 采集频率
T=1/Fs; % 采集时间间隔
N=length(ACC_X); % 采集信号的长度--采样点数
f1=33; % 第一个余弦信号的频率
f2=200; % 第二个余弦信号的频率
t=(0:1:N-1)*T; % 定义整个采集时间点
t=t'; % 转置成列向量
%y=1.2+2.7*cos(2*pi*f1*t+pi/4)+5*cos(2*pi*f2*t+pi/6); % 时域信号
y=ACC_X-min(ACC_X);
%% 绘制时域信号
figure
plot(t,y)
xlabel('时间')
ylabel('信号值')
title('时域信号')
%% fft变换
Y=fft(y); % Y为fft变换的结果,为复数向量
A=abs(Y); % 复数的幅值(模)
A(find(A==max(A)))=0;
RE=real(Y); % 复数的实部
IM=imag(Y); % 复数的虚部
%% 绘制fft变换结果(幅值,实部,虚部)
figure
subplot(3,1,1)
plot(0:1:N-1,A)
xlabel('序号 0 ~ N-1')
ylabel('幅值')
grid on