看看下面的程序,应该能帮上你的忙,已经通过调试:
Fs=256; %采样频率(Hz)
N=256; %采样点数
t=[0:1/Fs:N/Fs]; %采样时刻
S=2+3*cos(2*pi*10*t+pi*30/180)+cos(2*pi*20*t+pi*90/180);
%我的调试信号,你自己是电流电压数据的话,最开始通过load指令载入就是
Y = fft(S,N); %做FFT变换
Ayy = abs(Y); %取模
Ayy=Ayy/(N/2); %换算成实际的幅度
Ayy(1)=Ayy(1)/2;
F=([1:N]-1)*Fs/N; %换算成实际的频率值,Fn=(n-1)*Fs/N
stem(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果
title('幅度-频率曲线图');
fft(y)