资源描述：
《基于matlab的信号处理实例》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、基于matlab的信号处理实例[键入作者姓名]10基于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:lenx(i)=a(k);%将a中的第一列数据赋给xy(i)=a(k+1);%将a中
2、的第二列数据赋给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'),xlabel('n'),ylabel('z');函数调用1：signalsplot('C:UsersadminDesktop基于matla
3、b的信号处理实例run100m_TROUSERSPOCKET_1_陈佳_1.txt')实验结果1：102.将读出的3D加速度信号分解为偶序列及奇序列，分别绘出波形；绘出偶序列及奇序列的和信号、差信号及积信号的波形。实验代码2：oddandeven.m%将读出的3D加速度信号分解为偶序列及奇序列，分别绘出波形；%绘出偶序列及奇序列的和信号、差信号及积信号的波形functionoddandeven(filename)fid=fopen(filename,'r');a=fscanf(fid,'%d,%d,%d');%将文件中的数据读进来，并将其赋给afclose(
4、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进行反转z1=fliplr(z);%将z进行反转t=[-len:-11:len];%给出横坐标的范围，使之匹配10X1=1/2.*[x1,x];%求x的偶部Y1=1/2.*[y1,y];%求y的偶部Z1=1/2.*[z1,z];%求z的偶部fig
5、ure(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]');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)subplo
6、t(3,1,1);plot(t,X2,'r');xlabel('t'),ylabel('Oddx[n]');title('3D信号的奇信号');subplot(3,1,2)plot(t,Y2,'g');xlabel('t'),ylabel('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'
7、);xlabel('t'),ylabel('sumx[n]');title('奇信号与偶信号的和信号');subplot(3,1,2)plot(t,a2,'g');xlabel('t'),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');xlabel('t'),ylabe
8、l('de
