Matlab的时频分析程序(共1页).doc

clcclear allfname,fpath=uigetfile(*.xls,选择文件);du=fopen(strcat(fpath,fname),r);data=xlsread(strcat(fpath,fname);fclose(du);for n=1:8000t(n)=data(n,1);endfor n=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;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;