一维FFT的频率幅值关系

本文的参数变量

信号频率:

本文后续采用Fsign 来表示。

采样频率:

本文后续采用Fsample来表示 必须是Fsign的两倍以上

频率的模:

abs(fft(S))  也就是 转频率后的 根号下(实数的平方+ 虚数的平方)的绝对值。
频率的幅值:
A= abs(fft(S))/(N/2)
n点频率= (n-1)*Fsample/2

执行功能:

信号显示:

%step1 clean all window info
clc
clear
close all;
%step2 gener 
Fsample =500;
Fsign   = 50;%信号1
Tsample  =1/Fsample;
L        = 42;
t = (0:L-1)*Tsample;
%t2 = (0:Tsample:L-1);
S = sin(pi*Fsign*t);

figure(1)
plot(t,S);
grid on
title('Signal AMP')
xlabel('time (sec)')

figure(2)
plot((0:L-1),S);
grid on
title('Signal AMP')
xlabel('sample')

信号转为频率fft:

%step1 clean all window info
clc
clear
close all;
%step2 gener 
Fsample =500;
Fsign   = 10;%信号1
Tsample  =1/Fsample;
L        = 256;
t = (0:L-1)*Tsample;
%t2 = (0:Tsample:L-1);
S = sin(pi*Fsign*t);

figure(1)
plot((0:L-1),S);
grid on
title('Signal AMP')
xlabel('sample')

%step2 fft
NFFT = 2^nextpow2(L); %  1020 =》1024 120=》128
Y = fft(S,FFTN);      
f = Fsample/2*linspace(0,1,FFTN/2+1);
A=2*abs(Y(1:FFTN/2+1))/L; 
figure(2)
plot(f,A)
grid on
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('Amplitude')        

时域转换为频率表示幅值。

  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值