Matlab的时频分析程序.doc

《Matlab的时频分析程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、clcclearall[fname,fpath]=uigetfile('*.xls','选择文件');du=fopen(strcat(fpath,fname),'r');data=xlsread(strcat(fpath,fname));fclose(du);forn=1:8000t(n)=data(n,1);endforn=1:8000data1(n)=data(n,2);end%subplot(2,1,1);plot(t,data1)xlabel('时间/ms');ylabel('幅值');z0=fft(data1);Ts=6.64e-7;Ws=1/Ts;Wn=Ws/2;

2、n1=linspace(0,Wn,length(t)/2);z1=abs(z0(1:length(t)/2));subplot(2,2,3);plot(n1,z1);xlabel('频率/Hz');ylabel('幅值');z3=fft(data3);Ts3=6.64e-7;Ws3=1/Ts;Wn3=Ws3/2;n3=linspace(0,Wn3,length(t3)/2);z4=abs(z3(1:length(t3)/2));subplot(2,2,4);plot(n3,z4);xlabel('频率/Hz');ylabel('幅值');