谱分析,色谱分析,matlab

%*******************************清除环境变量******************************** clc; clear; close all %********************************读取数据******************************** load 105.mat %更改数据名； 因为一个mat文件种包含三个数据，每次选取一个进行分析 y=X105_BA_time(1:12000); %读取其中12000个点分析 %********************************相关参数******************************** N=length(y); %采样点数 Fs=12000; %采样频率 t=1/Fs:1/Fs:N/Fs; %采样时间 f=Fs*(1:N)/N; %频域点展开 %%********************************频谱分析******************************** y_quzhi=y-mean(y); %去除直流分量 Mag=abs(fft(y_quzhi)); %能量谱 Mag_real=2*abs(Mag)/N; %能量的真实幅值（频谱图） %%********************************细化谱分析******************************** f1=1200; %细化频率段起点 f2=1600; %细化频率段终点 M=4000; %细化频段的频点数，决定细化精度 w=exp(-1i*2*pi*(f2-f1)/(Fs*M)); %细化频段的跨度（步长） a=exp(1i*2*pi*f1/Fs); %细化频段的起始点 h=0:1:M-1; %细化频点序列 fx=(f2-f1)/M*h+f1; %频域点展开 Mag_x=czt(y,M,w,a); %细化谱 Mag_xreal=2*Mag_x/N; %细化谱的真实幅值 %********************************画图******************************** figure; plot(t,y) xlabel('时间(s)'); ylabel('幅值'); title('振动信号') figure; subplot(2,1,1); plot(f(1:N/2),Mag_real(1:N/2),'r') xlabel('频率(Hz)') ylabel('幅值') title('振动信号频谱图(显示范围：1200-1600Hz)') xlim([f1,f2]);%要显示的频段 subplot(2,1,2) plot(fx,abs(Mag_xreal),'r'); xlabel('频率(Hz)') ylabel('幅值') title('振动信号细化谱(细化范围：1200-1600Hz)');