资源描述:
《毕业论文-基于matlab的信号处理实例说明书.docx》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、2011级信息工程基于matlab的信号处理实例[键入作者姓名]2013/5/2323基于matlab的信号处理实例一.实验过程1.读取给定的3D加速度信号文件,绘出信号波形.实验代码1:signalsplot.mfunction[x,y,z]=signalsplot(filename)fid=fopen(filename,'r');a=fscanf(fid,'%d,%d,%d');%将文件中的数据读进来,并将其赋给afclose(fid);len=length(a)/3;%确定循环的总次数k=1;fori=1
2、:lenx(i)=a(k);%将a中的第一列数据赋给xy(i)=a(k+1);%将a中的第二列数据赋给yz(i)=a(k+2);%将a中的第三列数据赋给zk=k+3;endfigure(1)subplot(3,1,1);plot(x,'r'),xlabel('n'),ylabel('x');title('读取给定的3D加速度信号文件并绘出信号波形');subplot(3,1,2);plot(y,'g'),xlabel('n'),ylabel('y');subplot(3,1,3);plot(z,'b'),xlabe
3、l('n'),ylabel('z');函数调用1:signalsplot('C:UsersadminDesktop基于matlab的信号处理实例run100m_TROUSERSPOCKET_1_陈佳_1.txt')实验结果1:232.将读出的3D加速度信号分解为偶序列及奇序列,分别绘出波形;绘出偶序列及奇序列的和信号、差信号及积信号的波形。实验代码2:oddandeven.m%将读出的3D加速度信号分解为偶序列及奇序列,分别绘出波形;%绘出偶序列及奇序列的和信号、差信号及积信号的波形functionodd
4、andeven(filename)fid=fopen(filename,'r');a=fscanf(fid,'%d,%d,%d');%将文件中的数据读进来,并将其赋给afclose(fid);len=length(a)/3;%确定循环的总次数k=1;fori=1:lenx(i)=a(k);%将a中的第一列数据赋给xy(i)=a(k+1);%将a中的第二列数据赋给yz(i)=a(k+2);%将a中的第三列数据赋给zk=k+3;endx1=fliplr(x);%将x进行反转y1=fliplr(y);%将y进行反转z
5、1=fliplr(z);%将z进行反转t=[-len:-11:len];%给出横坐标的范围,使之匹配23X1=1/2.*[x1,x];%求x的偶部Y1=1/2.*[y1,y];%求y的偶部Z1=1/2.*[z1,z];%求z的偶部figure(2)subplot(3,1,1);plot(t,X1,'r');xlabel('t'),ylabel('Evex[n]');title('3D信号的偶信号');subplot(3,1,2);plot(t,Y1,'g');xlabel('t'),ylabel('Evey[n]'
6、);subplot(3,1,3);plot(t,Z1,'b');xlabel('t'),ylabel('Evez[n]');X2=1/2.*[-x1,x];%求x的奇部Y2=1/2.*[-y1,y];%求y的奇部Z2=1/2.*[-z1,z];%求z的奇部figure(3)subplot(3,1,1);plot(t,X2,'r');xlabel('t'),ylabel('Oddx[n]');title('3D信号的奇信号');subplot(3,1,2)plot(t,Y2,'g');xlabel('t'),ylab
7、el('Oddy[n]');subplot(3,1,3);plot(t,Z2,'b');xlabel('t'),ylabel('Oddz[n]');a1=X1+X2;%x的奇部与偶部之和a2=Y1+Y2;%y的奇部与偶部之和a3=Z1+Z2;%z的奇部与偶部之和figure(4)subplot(3,1,1);plot(t,a1,'r');xlabel('t'),ylabel('sumx[n]');title('奇信号与偶信号的和信号');subplot(3,1,2)plot(t,a2,'g');xlabel('t'
8、),ylabel('sumy[n]');subplot(3,1,3);plot(t,a3,'b');xlabel('t'),ylabel('sumz[n]');s1=X1-X2;%x的奇部与偶部之差s2=Y1-Y2;%y的奇部与偶部之差s3=Z1-Z2;%z的奇部与偶部之差figure(5)subplot(3,1,1);plot(t,s1,'r');xlab