举例：利用 matlab 将线性系统与非线性系统的加速度时域数据进行频域对比分析
代码部分
close all
%此处为时域数据，例子中即为 线性的数据和非线性的数据
sig1 = feixianxing.data;%非线性
sig2 = xianxing.data;%线性
% 采样频率
fs = 500;
% 信号长度
N1=length(sig1);
N2=length(sig2);
% 使用汉宁窗 窗要与信号（加速度）相同长度
window1= hanning(N1);
window2= hanning(N2);
nfft = 9000;
%周期图法功率谱 periodogram 是求功率谱密度的函数,[a,b]中 a 是功率谱，b 是频率
[Pxx1, f1] =periodogram(sig1, window1, nfft, fs);
[Pxx2, f2] =periodogram(sig2, window2, nfft, fs);
figure
semilogx(f2, Pxx2);
hold on
semilogx(f1, Pxx1);
axis tight
xlabel(' ontname{宋体}频率 ontname{Times New Roman}/(Hz)m ','fontsize',13)
%%注意角加速度单位和加速度单位区别
ylabel(' ontname{ 宋体} 功 率 谱 密 度 ontname{Times New Roman}/(rad^2/s^3)m
','fontsize',13)
legend({' ontname{宋体}线性',' ontname{宋体}非线性'},'fontsize',11)
xlim([10^(-1), 10^2]);
ylim([0 2])
结果图
资源评论
