信号fft变换

信号频谱图

fun = @(x)0.*(x<0 |x>5)+x.*(x>=0 & x<1)+1.*(x>=1 & x<4) +(5-x).*(x>=4 & x<5);%分段函数
x = -5:0.2:10;  % x自变量
y = fun(x);  % 得到输入的平稳信号y
Dfy = fft(y);  % 离散Fourier变换
Dfy_shift = fftshift(Dfy);  % 对称变换得到对称的Fourier频谱
figure
plot(y); axis([0,80,-0.05,1.05]),xlabel('n'),ylabel('y');  % 原始平稳信号
figure
plot(abs(Dfy)),xlabel('n'),ylabel('|Dfy|');  % |Dft|为离散频谱幅值信息
figure
plot(abs(Dfy_shift)),xlabel('n'),ylabel('|Dfy|');  % 对称变换后的频谱