车辆信号试验数据分析

%原始信号
load data3.dat; %读入数据
data=data3; %将数据阵列赋值给data变量
fs=1000; %采样频率
dt=1/fs; %采样间隔/周期
N=length(data);
t=[0:N-1]*dt; %时间轴序列
y1=data(:,2); %数据第一通道(第2列)
y2=data(:,6);
y3=data(:,3);
y4=data(:,7);
figure(1);
plot(t,y1,’r’);
grid on;
hold on;
plot(t,y3,’b’);
grid on;
legend(‘信号y1’,’信号y3’);
xlabel(‘Time/s’);
title(‘原始声信号y1和y3’)
hold on;
figure(2);
plot(t,y2,’r’);
grid on;
hold on;
plot(t,y4,’b’);
grid on;
legend(‘信号y2’,’信号y4’);
xlabel(‘Time/s’);
title(‘原始振动信号y2和y4’)

Y = fft(y1,256); %256点快速傅立叶变换FFT
Pyy = Y.*conj(Y)/256; %复共轭(CONJ)计算功率谱密度
f = fs/256*(0:127); %为前128数据点(后128点与支对称)建立频率轴。
figure(3); %开第3个图形窗，并在该窗口画出第1个图形
plot(f,Pyy(1:128)); %画频谱