振动信号EMD分析matlab算法

振动信号EMD分析matlab算法

振动信号EMD分析matlab算法，介绍了EMD分析处理的基本原理，可在此基础上改进算法应用到具体问题上

程序如下:

clc

clear

%z=[4319.0,4358.5,4392.2,4418.9,4438.5,4453.1,4465.1,4473.0,4476.5,4477.0,4477.5,4478.8,44 77.7,4472.1,4462.2,4449.9,4435.9,4419.0,4399.5,4380.5,4364.3,4351.5,4341.5,4334.7,4332.7,433 5.2,4338.7,4339.5,4336.7,4332.5,4327.9,4322.4,4316.7,4312.3,4310.9,4313.0,4317.4,4324.3,4335 .3,4351.0,4369.4,4386.5,4398.6,4405.0,4408.6,4411.6,4415.5,4421.3,4430.1,4441.2,4450.6,4456. 0,4458.0,4460.1,4464.2,4467.7,4468.2,4465.7,4461.9,4458.3,4454.2,4448.6,4442.0,4434.9,4427.0 ,4417.2,4405.9,4396.0,4388.8,4383.3,4377.6,4370.3,4363.2,4357.3,4350.8,4341.8,4330.9,4321.2, 4314.9,4311.8,4311.2,4313.7,4320.7,4330.7,4340.5,4347.3,4351.9,4354.9,4355.6,4352.8,4346.9,4 341.4,4338.9,4338.9,4340.1,4341.7,4345.0,4350.3,4355.0,4357.3,4356.6,4355.0,4353.5,4351.0,43 48.0,4345.9,4347.3,4352.1,4357.5,4361.9,4365.6,4370.0,4374.3,4377.3,4378.8,4380.3,4383.3,438 6.1,4387.3,4386.9,4386.8,4388.6,4390.8,4392.2,4391.6,4389.7,4387.0,4382.0,4375.8,4370.6,4368 .2,4369.2,4371.0,4371.0,4369.4,4367.2,4365.8,4364.9,4365.2,4368.3,4374.3,4380.6,4383.5,4381. 1,4375.9,4371.0,4366.8,4364.1,4363.6,4366.5,4371.9,4376.6,4379.0,4380.4,4382.7,4385.9,4387.8 ,4387.2,4386.0,4385.8,4387.4,4389.7,4392.4,4396.9,4403.7,4410.8,4416.1,4419.2,4421.6,4424.7, 4427.2,4428.5,4429.3,4431.3,4435.1,4437.9,4438.4,4437.0,4435.2,4433.2,4429.2,4422.9,4416.7,4 412.6,4410.5,4408.5,4405.4,4402.1,4399.1,4394.5,4385.7,4372.2,4354.5,4332.6,4303.7];

%z=[5210.1,5291.7,5362.0,5424.2,5481.8,5534.9,5582.4,5624.1,5662.4,5698.9,5732.0,5759.4,57 79.6,5793.0,5800.1,5799.7,5792.1,5780.1,5766.4,5751.9,5734.5,5712.3,5685.5,5653.8,5614.0,556 2.3,5496.4,5418.3,5331.7,5238.9];

%z=[5585.5,5650.5,5699.2,5737.4,5765.4,5782.2,5787.8,5785.3,5778.3,5767.8,5755.5,5743.7,57 35.3,5731.1,5728.3,5723.9,5716.2,5703.8,5684.7,5655.8,5616.0,5566.8,5509.7,5444.9,5372.3,529 4.4,5214.8,5136.8,5061.1,4985.6,4907.9,4828.8,4749.6,4671.7,4597.6,4530.5,4473.1,4426.7,4388 .7,4356.6,4330.5,4311.6,4298.8,4289.3,4281.1,4273.9,4269.0,4266.8,4265.9,4264.9,4264.6,4266. 3,4270.1];

z=[0.03 -0.11 0.18 -0.17 -0.05 -0.36 0.14 -0.18 0.22 0.03 0.17 -0.15 -0.32 0.01 0.12 -0.02 -0.19 -0.02 -0.08 0.06 -0.12 0.15 0.27 -0.01 0.24 0.00 -0.21 0.08 0.04 -0.05 0.19 0.02 0.29 -0.05 0.08 -0.04 0.06 -0.03 -0.12 -0.02 -0.20 0.15 0.19 0.16 0.05 -0.04 0.31 0.03 -0.09 0.07 0.29 0.08 -0.12 -0.22 0.04 0.20 -0.10

];

figure(1)

t=1951:1:2007;

set(gcf,'color','w')

plot(t,z,'k')

set(gca,'fontname','times New Roman')

set(gca,'fontsize',14.0)

axis([1951,2007,-Inf,Inf])