用matlab做谱分析程序,求一个能谱分析的matlab程序

MATLAB频谱分析程序

%FFT变换,获得采样数据基本信息，时域图，频域图

%这里的向量都用行向量，假设被测变量是速度，单位为m/s

clear;

close all;

load data.txt              %通过仪器测量的原始数据，存储为data.txt中，附件中有一个模版(该信号极不规则)

A=data;                  %将测量数据赋给A，此时A为N×2的数组

x=A(:,1);                 %将A中的第一列赋值给x，形成时间序列

x=x';                     %将列向量变成行向量

y=A(:,2);                  %将A中的第二列赋值给y，形成被测量序列

y=y';                     %将列向量变成行向量

%显示数据基本信息

fprintf('\n数据基本信息：\n')

fprintf('        采样点数 = %7.0f \n',length(x))        %输出采样数据个数

fprintf('        采样时间 = %7.3f s\n',max(x)-min(x))   %输出采样耗时

fprintf('        采样频率 = %7.1f Hz\n',length(x)/(max(x)-min(x)))   %输出采样频率

fprintf('        最小速度 = %7.3f m/s\n',min(y))      %输出本次采样被测量最小值