离散信号频谱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('