level 3
clc
clear all
a=xlsread('fxxstxy3.xlsx',2);%信号时间6s,6000个点,故采样频率2000,采样点数6000
fs=2000;%采样频率
N=6000;%采样点数
dt=1/fs;%时间间隔
t=0:dt:(N-1)*dt;
b=(fft(a,N)); %a为时域信号,作fft变换
%幅值修正
P2 =2*(b/N);
P1 =P2(1:N/2+1);%此时选取前半部分,因为fft之后为对称的双边谱
%画图
f = fs/N*(0:(N/2));
plot(f,abs(P1))
title('频域信号')
xlabel('f (Hz)')
ylabel('amplitude(m)')
2021年05月01日 12点05分

