【实例简介】

【实例截图】

【核心代码】

samplingrate = 0.00005;

[s,t,f]=st(amplitude,0,0.5/samplingrate,samplingrate,1);

T=max(t);

subplot(3,1,1),plot(t,amplitude,'color','k','linewidth',1)

xlabel('时间/s');

ylabel('幅值/V');

xlim([0 T]);

% title('微震波形','FontName','Times New Roman','FontSize',14);

set(gca,'FontName','Times New Roman','FontSize',14);

set(get(gca,'XLabel'),'FontName','Times New Roman','FontSize',14);

set(get(gca,'YLabel'),'FontName','Times New Roman','FontSize',14);

subplot(3,1,2:3),contour(t,f,abs(s),'fill','on');

xlabel('时间/s');

ylabel('频率/Hz');

xlim([0 T]);

ylim([0 3000]);

% title('时频谱','FontName','Times New Roman','FontSize',14);

set(gca,'FontName','Times New Roman','FontSize',14);

set(get(gca,'XLabel'),'FontName','Times New Roman','FontSize',14);

set(get(gca,'YLabel'),'FontName','Times New Roman','FontSize',14);