华南理工大学信号及系统实验基于Matlab的信号处理实例.docVIP

第2讲 基于Matlab的信号处理实例实验内容（1）(1) 读取给定的3D加速度信号文件，绘出信号波形；程序源代码： function sy2fid = fopen(run 100m_TROUSERS POCKET_1_陈佳_1.txt,r);a=fscanf(fid,%d,%d,%d\n);fclose(fid);len=length(a)/3;k=1;for i=1:len x(i)=a(k); y(i)=a(k+1); z(i)=a(k+2); k=k+3;endsubplot(3,1,1);plot(x(1:700),g,LineWidth,3);grid on; %标示格子xlabel(n),ylabel(x);%标示坐标legend(x(i)); %标示名称subplot(3,1,2);plot(y(300:700),m,LineWidth,3)grid on; %标示格子xlabel(n),ylabel(y);%标示坐标legend(y(i)); %标示名称subplot(3,1,3);plot(z(300:700),b,LineWidth,3);grid on; %标示格子xlabel(n),ylabel(z);%标示坐标legend(z(i)); %标示名称执行结果如下所示：(2) 将读出的3D加速度信号分解为偶序列及奇序列，分别绘出波形；绘出偶序列及奇序列的和信号、差信号及积信号的波形；a.分解为偶序列及奇序列，分别绘出波形:程序源代码：function sy2oefid = fopen(run 100m_TROUSERS POCKET_1_陈佳_1.txt,r);a=fscanf(fid,%d,%d,%d\n);fclose(fid);len=length(a)/3;k=1;for i=1:len x(i)=a(k); y(i)=a(k+1); z(i)=a(k+2); k=k+3;endx0=[zeros(1,len-1),x];%x前面补0x1=fliplr(x0);%x1即是反转的x0ex=(x1+x0)/2;%x的偶部ox=-(x1-x0)/2;%x的奇部y0=[zeros(1,len-1),y];y1=fliplr(y0);ey=(y1+y0)/2;oy=-(y1-y0)/2;z0=[zeros(1,len-1),z];z1=fliplr(z0);ez=(z1+z0)/2;oz=-(z1-z0)/2;subplot(3,2,1);plot([-len+1:len-1],ex ,g,LineWidth,3);grid on; %标示格子xlabel(n),ylabel(x);%标示坐标legend(ex(i)); %标示名称subplot(3,2,3);plot([-len+1:len-1],ey,m,LineWidth,3)grid on; %标示格子xlabel(n),ylabel(y);%标示坐标legend(ey(i)); %标示名称subplot(3,2,5);plot([-len+1:len-1],ez,b,LineWidth,3);grid on; %标示格子xlabel(n),ylabel(z);%标示坐标legend(ez(i)); %标示名称subplot(3,2,2);plot([-len+1:len-1],ox,g,LineWidth,3);grid on; %标示格子xlabel(n),ylabel(x);%标示坐标legend(ox(i)); %标示名称subplot(3,2,4);plot([-len+1:len-1],oy,m,LineWidth,3)grid on; %标示格子xlabel(n),ylabel(y);%标示坐标legend(oy(i)); %标示名称subplot(3,2,6);plot([-len+1:len-1],oz,b,LineWidth,3);grid on; %标示格子xlabel(n),ylabel(z);%标示坐标legend(oz(i)); %标示名称sumoex=ex+ox;sumoey=ey+oy;sumoez=ez+oz;执行结果如下所示：b.绘出偶序列及奇序列的