fft.zip_FFT分析_matlab 频谱分析_matlab频谱分析_spectrum analysis_功率谱分析

clear all; close all; %读取数据 data=xlsread('C:\Users\dc\Desktop\12.xls'); figure(1); subplot(231); plot(data,'-r');%绘制信号波形 xlabel('时间/s');ylabel('幅值/mV'); title('信号时域波形') %进行FFT变换并做频谱图 fs=250; %采样频率 N=15000; %数据长度 n=0:N-1; t=n/fs; %进行FFT变换并做频谱图 y=fft(data,N);%进行fft变换 mag=abs(y);%求幅值 f=(0:N-1)*fs/N;%横坐标频率的表达式为f=(0:M-1)*Fs/M; figure(1); subplot(232); plot(f,mag);%做频谱图 axis([0,20000,0,2000]); xlabel('频率(Hz)'); ylabel('幅值'); title('信号幅频谱图'); %求均方根谱 sq=abs(y); figure(1); subplot(233); plot(f,sq); xlabel('频率(Hz)'); ylabel('均方根谱'); title('信号均方根谱'); grid; %求功率谱 power=sq.^2; figure(1); subplot(234); plot(f,power); xlabel('频率(Hz)'); ylabel('功率谱'); title('信号功率谱'); grid; %求对数谱 ln=log(sq); figure(1); subplot(235); plot(f,ln); xlabel('频率(Hz)'); ylabel('对数谱'); title('信号对数谱'); grid; %用IFFT恢复原始信号 xifft=ifft(y); magx=real(xifft); ti=[0:length(xifft)-1]/fs; figure(1); subplot(236); plot(ti,magx); xlabel('t'); ylabel('y'); title('通过IFFT转换的信号波形');