哈工大实验方法二次大作业.docVIP

哈尔滨工业大学 试验方法及数字信号处理分析 ————第二次大作业 通过位移算速度 指导老师：包钢 学生姓名：陈方鑫 学生学号：15S008043 第一部分 作业题目 设计题目 条件：（1）x为110000个采集的位移点，已给出具体值。 （2）采用周期。 （3）有用信号频率为：0Hz。 要求： （1）绘出杂波信号波形。 （2）根据位移曲线求速度曲线。 （3）对比对位移滤波和不滤波，速度曲线形状差别。 第二部分 具体设计内容 一、设计程序展示 x=[ ]; t1=0.0005; t=t1:t1:55; %对位移滤波 [n1,wn]=buttord(0.001,0.015,0.95,40); [b,a]=butter(n1,wn); y=filter(b,a,x); freqz（b，a）； %计算速度 v=zeros(1,110000); for i=2:110000 v(i)=(y(i)-y(i-1))/t1; end %截取前2000个点（由于滤波过程中，前面滤波数据出现畸变，现在截取，使得之后的绘图更清晰） v1=zeros(1,108000); for i=1:108000 v1(i)=v(i+2000); end %为滤波数据计算的速度 v0=zeros(1,110000); for i=2:110000 v0(i)=(x(i)-x(i-1))/t1; end T=1+t1:t1:55; %绘图 subplot(211); plot(t,v0); title(滤波前的速度图像); subplot(212); plot(T,v1); title(滤波后的速度图像); %计算最大值 for i=2:108000 if v1(i)v(i-1) vm=v1(i); else vm=v1(i-1); end end 二、程序分析 1、对于语句： [n1,wn]=buttord(0.001,0.015,0.95,40); 由于有效频率为0，所以对于截止频率Wn希望尽量小一点。对于Wmax的估计是认定：Wmax=Ws/2=1000Hz。通过多次试验，发现Wn=0.001，Ws=0.015相对比较合适。 2、对于语句 v1=zeros(1,108000); for i=1:108000 v1(i)=v(i+2000); end 由于滤波没有考虑前面数据畸变的问题，所以直接绘出的速度图像前2000个点是有畸变的，这样导致图像整体不美观，最主要的是导致图像畸变部分占刻度多，有用部分反而占刻度少。为了图像更直观（真正科研时不能如此处理，此处只是为了作业美观），对速度信号做了时移处理，丢掉前2000个畸变点。 程序用的是MATLAB现成的滤波程序，滤波参数是经过多次试验得到的相对较好的参数。 求速度公式： for i=2:110000 v0(i)=(x(i)-x(i-1))/t1; end 用的是公式，也可认为是差分方程的形式给出。 结果展示 频率响应图 局部放大图 分析：从程序可看出，我们取的截止频率Wn和过渡频率Ws分别为0.001和0.0015，是非常小的，所以频率响应曲线的低通部分非常小，数量级为10-5，这也是前期通过多次测试，发现这个数量级的滤波效果相对较好。 2.位移图像 图1 图2 分析： 通过比较可以看出，滤波前，信号虽然保持了大致的加速位移的信息，但是具体到每段曲线，毛刺很严重。滤波后，图像个光滑程度有了明显的增加，基本可以确定认为干扰已经被滤出。用滤波后的位移求速度是可行的。 3.速度图像 图3 分析： 第一张图示滤波前的速度图像，可以看出，由于位移数据有很多干扰，而且速度和位移之间是导数关系，相当于对数据做了一次微分，这样就将干扰严重放大，使得真正的速度信息在图像中完全没有体现。 而第二张图示滤波之后的结果，可以发现，虽然速度图像任然不光滑，但是速度的基本信息已经包含在图像中，而且求出最大速度为4.2mm/s。同时可以看出该运动是一个加速度越来越大的加速运动。同时从这里也能看出，速度图像之所以有毛刺，是因为位移滤波是没有滤“干净”的。如果要想速度曲线进一步光滑，就要对速度曲线进一步高质量滤波。 第三部分 实验总结 总的来说，本次实验室很成功的。而且通过这次大作业，我也收获了很多课堂上没学到的知识。 首先，通过两次实验对比，用MATLAB的工具滤波效果好，效率高，程序简单，使用方便。但是在使用的过程中，由于有一些参数没有实践经验，所以不知道如何适当选取。比如对于1、2不知道多少合适，所以选择参