% 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);
- 1
- 2
前往页