离散信号matlab频谱分析程序

《离散信号matlab频谱分析程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、离散信号MATLAB频谱分析程序【ZZ】%FFT变换,获得采样数据基本信息，时域图，频域图%这里的向量都用行向量，假设被测变量是速度，单位为m/sclear;closeall;loaddata.txt                      %通过仪器测量的原始数据，存储为data.txt中，附件中有一个模版(该信号极不规则)A=data;                          %将测量数据赋给A，此时A为N×2的数组x=A(:,1);                        %将A中的第一列赋值给x，形成时间序列x=x';             

2、               %将列向量变成行向量y=A(:,2);                        %将A中的第二列赋值给y，形成被测量序列y=y';                            %将列向量变成行向量%显示数据基本信息fprintf('数据基本信息：')fprintf('     采样点数=%7.0f',length(x))                %输出采样数据个数fprintf('     采样时间=%7.3fs',max(x)-min(x))             %输出采样耗时fprintf(' 

3、    采样频率=%7.1fHz',length(x)/(max(x)-min(x)))  %输出采样频率fprintf('     最小速度=%7.3fm/s',min(y))                %输出本次采样被测量最小值fprintf('     平均速度=%7.3fm/s',mean(y))              %输出本次采样被测量平均值fprintf('     速度中值=%7.3fm/s',median(y))            %输出本次采样被测量中值fprintf('     最大速度=%7.3fm/s',max(

4、y))                 %输出本次采样被测量最大值fprintf('     标准方差=%7.3f',std(y))                    %输出本次采样数据标准差fprintf('    协方差=%7.3f',cov(y))                     %输出本次采样数据协方差fprintf('   自相关系数=%7.3f',corrcoef(y))               %输出本次采样数据自相关系数  %显示原始数据曲线图（时域）subplot(2,1,1);plot(x,y)            

5、                                         %显示原始数据曲线图axis([min(x)max(x)1.1*floor(min(y))1.1*ceil(max(y))])        %优化坐标，可有可无xlabel('时间(s)');ylabel('被测变量y');title('原始信号(时域)');gridon;%傅立叶变换y=y-mean(y);                               %消去直流分量，使频谱更能体现有效信息Fs=2000;                                 

6、  %得到原始数据data.txt时，仪器的采样频率。其实就是length(x)/(max(x)-min(x));N=10000;                                离散信号MATLAB频谱分析程序【ZZ】%FFT变换,获得采样数据基本信息，时域图，频域图%这里的向量都用行向量，假设被测变量是速度，单位为m/sclear;closeall;loaddata.txt                      %通过仪器测量的原始数据，存储为data.txt中，附件中有一个模版(该信号极不规则)A=data;                   

7、       %将测量数据赋给A，此时A为N×2的数组x=A(:,1);                        %将A中的第一列赋值给x，形成时间序列x=x';                            %将列向量变成行向量y=A(:,2);                        %将A中的第二列赋值给y，形成被测量序列y=y';                            %将列向量变成行向量%显示数据基本信息fprintf('数据基本信息：')fprintf('     采样点数=%7.0