华工DSP实验2讲述.doc

实验二 加速度信号的处理 ------黄晓荣 13电联 ------201330260151 注：本实验所用的数据为junpData 一、Question one 1.1、问题描述 读出给定的3D加速度信号，并分别画出三个轴的信号波形； 1.2、实验代码 clc tt=fopen(h:\jumpData.txt); var = textscan(tt,%f %f %f,delimiter, ,); data=cell2mat(var); a=size(data); figure(1) subplot(3, 1 ,1); plot(1:a(1), data(:,1)) xlabel(the original signal in x-ray); subplot(3, 1 ,2); plot(1:a(1), data(:,2)) xlabel(the original signal in y-ray); subplot(3, 1 ,3); plot(1:a(1), data(:,3)) xlabel(the original signal in z-ray); 1.3、实验结果 1.4、结果分析 可以看出，三个轴线的采集数据很凌乱，基本看不出什么区别，因此需要后续的数据预处理 二、Question two 2.1、问题描述 用M点（M分别取值3、5、8）滑动平均滤波器对输入的3D加速度信号进行滤波，绘出滤波后信号的波形；观察并分析M值对输出波形的影响；并绘出该滑动滤波器传输函数的零极点图； 2.2、实验代码 %windowSize = 3 figure(2) windowSize = 3; winData3=filter(ones(1,windowSize)/windowSize,1,data); subplot(3, 1 ,1); plot(1:a(1), winData3(:,1)) subplot(3, 1 ,2); plot(1:a(1), winData3(:,2)) subplot(3, 1 ,3); plot(1:a(1), winData3(:,3)) xlabel(the length of window = 3) %windowSize = 5 figure(3) windowSize = 5; winData5 = filter(ones(1,windowSize)/windowSize,1,data); subplot(3, 1 ,1); plot(1:a(1), winData5(:,1)) subplot(3, 1 ,2); plot(1:a(1), winData5(:,2)) subplot(3, 1 ,3); plot(1:a(1), winData5(:,3)) xlabel(the length of window = 5) %windowSize = 8 figure(4) windowSize = 8; winData8 = filter(ones(1,windowSize)/windowSize,1,data); subplot(3, 1 ,1); plot(1:a(1), winData8(:,1)) subplot(3, 1 ,2); plot(1:a(1), winData8(:,2)) subplot(3, 1 ,3); plot(1:a(1), winData8(:,3)) xlabel(the length of window = 8) 2.3、实验结果 M=3的结果 M=5的结果 M=8的结果 2.4、结果分析 可以看出，随着滑动平均滤波器的M的取值增大，信号的平滑性就越好，但是也会损失 的细节信息； 三、Question two 3.1、问题描述 对输入的3D加速度信号进行高斯滤波（窗长取32，参数取2.25）和幅度归一化处理，绘出处理后的信号ax[n],ay[n] 和az[n]的波形； 3.2、实验代码 %Guas Filter figure(5) sigma = 2.25; gausFilter = fspecial(gaussian,32,sigma); guasData = imfilter(data,gausFilter,replicate); subplot(3, 1 ,1); [a1, p1] = mapminmax(guasData(:,1)); plot(1:a(1), a1) subplot(3, 1 ,2); [a2, p2] = mapminmax(guasData(:,2)); plot(1:a(1), a2) subplot(3, 1 ,3); [a3, p3] = mapminmax(guasData(:,3)); plot(1:a