資源描述:
《基于matlab的信號(hào)處理實(shí)例》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫。
1、基于matlab的信號(hào)處理實(shí)例[鍵入作者姓名]10基于matlab的信號(hào)處理實(shí)例一.實(shí)驗(yàn)過程1.讀取給定的3D加速度信號(hào)文件,繪出信號(hào)波形.實(shí)驗(yàn)代碼1:signalsplot.mfunction[x,y,z]=signalsplot(filename)fid=fopen(filename,'r');a=fscanf(fid,'%d,%d,%d');%將文件中的數(shù)據(jù)讀進(jìn)來,并將其賦給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中
2、的第二列數(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加速度信號(hào)文件并繪出信號(hào)波形');subplot(3,1,2);plot(y,'g'),xlabel('n'),ylabel('y');subplot(3,1,3);plot(z,'b'),xlabel('n'),ylabel('z');函數(shù)調(diào)用1:signalsplot('C:UsersadminDesktop基于matla
3、b的信號(hào)處理實(shí)例run100m_TROUSERSPOCKET_1_陳佳_1.txt')實(shí)驗(yàn)結(jié)果1:102.將讀出的3D加速度信號(hào)分解為偶序列及奇序列,分別繪出波形;繪出偶序列及奇序列的和信號(hào)、差信號(hào)及積信號(hào)的波形。實(shí)驗(yàn)代碼2:oddandeven.m%將讀出的3D加速度信號(hào)分解為偶序列及奇序列,分別繪出波形;%繪出偶序列及奇序列的和信號(hào)、差信號(hào)及積信號(hào)的波形functionoddandeven(filename)fid=fopen(filename,'r');a=fscanf(fid,'%d,%d,%d');%將文件中的數(shù)據(jù)讀進(jìn)來,并將其賦給afclose(
4、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進(jìn)行反轉(zhuǎn)y1=fliplr(y);%將y進(jìn)行反轉(zhuǎn)z1=fliplr(z);%將z進(jìn)行反轉(zhuǎn)t=[-len:-11:len];%給出橫坐標(biāo)的范圍,使之匹配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信號(hào)的偶信號(hào)');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信號(hào)的奇信號(hào)');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('奇信號(hào)與偶信號(hào)的和信號(hào)');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