資源描述：
《畢業(yè)論文-基于matlab的信號處理實例說明書.docx》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
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');%將文件中的數(shù)據(jù)讀進來，并將其賦給afclose(fid);len=length(a)/3;%確定循環(huán)的總次數(shù)k=1;fori=1
2、:lenx(i)=a(k);%將a中的第一列數(shù)據(jù)賦給xy(i)=a(k+1);%將a中的第二列數(shù)據(jù)賦給yz(i)=a(k+2);%將a中的第三列數(shù)據(jù)賦給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');函數(shù)調(diào)用1：signalsplot('C:UsersadminDesktop基于matlab的信號處理實例run100m_TROUSERSPOCKET_1_陳佳_1.txt')實驗結(jié)果1：232.將讀出的3D加速度信號分解為偶序列及奇序列，分別繪出波形；繪出偶序列及奇序列的和信號、差信號及積信號的波形。實驗代碼2：oddandeven.m%將讀出的3D加速度信號分解為偶序列及奇序列，分別繪出波形；%繪出偶序列及奇序列的和信號、差信號及積信號的波形functionodd
4、andeven(filename)fid=fopen(filename,'r');a=fscanf(fid,'%d,%d,%d');%將文件中的數(shù)據(jù)讀進來，并將其賦給afclose(fid);len=length(a)/3;%確定循環(huán)的總次數(shù)k=1;fori=1:lenx(i)=a(k);%將a中的第一列數(shù)據(jù)賦給xy(i)=a(k+1);%將a中的第二列數(shù)據(jù)賦給yz(i)=a(k+2);%將a中的第三列數(shù)據(jù)賦給zk=k+3;endx1=fliplr(x);%將x進行反轉(zhuǎn)y1=fliplr(y);%將y進行反轉(zhuǎn)z
5、1=fliplr(z);%將z進行反轉(zhuǎn)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
