北理工-现代数据处理-李明-实验六车辆信号试验数据分析.docVIP

车辆信号试验数据分析 姓名：査小江 学号：1102130099实验6. 车辆信号试验数据分析 1 6.1 实验目的 1 6.2 实验内容 1 6.3 数据处理 2 图1. 原始声信号y1和y3 2 图2. 原始振动信号y2和y4 3 图3. 声音信号频谱分析 3 图4. 第8通道原始震动信号 4 图5. 第8通道原始信号频谱分析 4  车辆信号试验数据分析 实验目的 熟悉工程数据的简单处理 实验内容 读入data3.txt文件中的数据，第一列数据的时间，作出第1通道和第6通道数据的时间曲线。 load data3.txt; %读入数据 data=data3; %将数据阵列赋值给data变量 fs=10000; %采样频率 dt=1/fs; %采样间隔/周期 N=length(data); t=[0:N-1]*dt; %时间轴序列 y1=data(:,2); %数据第一通道(第2列) y2=data(:,6); y3=data(:,3); y4=data(:,7); y5=data(:,9); 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,1024); %256点快速傅立叶变换FFT Pyy = Y.*conj(Y)/1024; %复共轭(CONJ)计算功率谱密度 f = fs/1024*(0:511); %为前128数据点(后128点与支对称)建立频率轴。 figure(3); %开第3个图形窗，并在该窗口画出第1个图形 plot(f,Pyy(1:512)); figure(4); plot(t,y5,r); title(第8通道原始信号); Y = fft(y4,256); Pyy = Y.*conj(Y)/256; f = fs/256*(0:127); figure(6); plot(f,Pyy(1:128)); title(第8通道原始信号频谱分析); 数据处理 原始声信号y1和y3 原始振动信号y2和y4 声音信号频谱分析 第8通道原始震动信号 第8通道原始信号频谱分析 经数据信号频谱分析可以看到，声音以及振动信号的频率都集中在某一数值附近。 4