matlab时频域分析代码

matlab时频域分析代码

clear;

clc;

load('f:\a.txt'); x=a(1:length(a),1); y=a(1:length(a),2); fs=length(x)/(max(x)-min(x));

n=length(x)-1;

t=n/fs;

N = 2^nextpow2(n);

z=fft(y,N);

mag=2*abs(z)/N;

f=(0:length(z)-1)'*fs/length(z);

figure(1);

plot(f(1:N/2),mag(1:N/2));

grid on

axis ([0 6000 0 0.15]) xlabel('频率/Hz')

ylabel('幅值/V')

title('幅频谱')

figure(2);

z2=rceps(y);

plot(x,z2)

grid on

axis([0 0.1 -0.2 0.2]) xlabel('时间/s')

ylabel('幅值/V')

title('倒频谱') figure(3);

plot(x,y);

grid on;

xlabel('时间/s') ylabel('幅值/V') title('时间曲线')