滚动轴承振动信号的包络谱分析.rar

% EMD的Hilbert包络解调 load 2004.02.17.08.02.39; x1=X2004_02_17_08_02(1:20000,2); load 2004.02.19.06.02.39; x2=X2004_02_19_06_02(1:20000,2); y1=x1; y2=x2; figure(1) fs=20000; N=20000; t=1/fs:1/fs:N/fs; plot(t(1:1000),x1(1:1000),'k'); set(gca, 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'XLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'YLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); xlabel('采样时间t/s','fontsize',10.5,'Fontname', '方正宋黑简体'); ylabel('幅值A/g','fontsize',10.5,'Fontname', '方正宋黑简体'); figure(2) plot(t(1:1000),y1(1:1000),'k'); set(gca, 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'XLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'YLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); xlabel('采样时间t/s','fontsize',10.5,'Fontname', '方正宋黑简体'); ylabel('幅值A/g','fontsize',10.5,'Fontname', '方正宋黑简体'); figure(3) envelope_xqy(y1,20000); figure(4) plot(t(1:1000),x2(1:1000),'k'); set(gca, 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'XLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'YLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); xlabel('采样时间t/s','fontsize',10.5,'Fontname', '方正宋黑简体'); ylabel('幅值A/g','fontsize',10.5,'Fontname', '方正宋黑简体'); figure(5); plot(t(1:1000),y2(1:1000),'k'); set(gca, 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'XLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); set(get(gca,'YLabel'), 'Fontname', 'Times New Roman', 'Fontsize', 10.5); xlabel('采样时间t/s','fontsize',10.5,'Fontname', '方正宋黑简体'); ylabel('幅值A/g','fontsize',10.5,'Fontname', '方正宋黑简体'); figure(6); envelope_xqy(y2,20000);