close all
clear all
clc
t=-1:0.01:1;
x=exp(-0.12*t).*cos(0.25*pi*t);
plot(t,x);
y=fft(x);%matlab自带求频谱函数
%hold on
% figure,plot(y);
y1=fftshift(y);
% figure,plot(y1);
y2=log10(abs(y1));%工程上一般取对数
figure,plot(y2);
没有查相应的公式,就随便回答了,应该是数字信号处理上有相应的求和公式,抱歉手头上没有查阅相应书籍,实现不难,用循环就好了,自带的fft也是求频谱的函数