www.gusucode.com > 局部均值分解源代码 难得的matlab程序代码源码 > lmd/lmd2_example.m

    fs=1000;
N=1000;
%ts=1/fs;
n=0:N-1;
t=n/fs;
x=@(t) (1+0.5*cos(9*pi*t)).*cos(200*pi*t+2*cos(10*pi*t))+sin(pi*t).*sin(30*pi*t);
%t=0:0.001:1-0.001;
y=x(t);
%figure;
%plot(t,y);
%--------------------------
[pf,a,si]=lmd2(y);
subplot(311);plot(t,pf(1,:));
subplot(312);plot(t,pf(2,:));
subplot(313);plot(t,pf(3,:));
% subplot(514);plot(t,pf(4,:));
% subplot(515);plot(t,pf(5,:));
figure;
subplot(211);plot(t,si(1,:));
subplot(212);plot(t,si(2,:));
% subplot(313);plot(t,si(3,:));
% subplot(414);plot(t,si(4,:));
figure;
subplot(211);plot(t,a(1,:));
subplot(212);plot(t,a(2,:));
figure;
ff=shunpin(si).*fs;%实际频率要乘以采样率
subplot(211);plot(t,ff(1,:));
subplot(212);plot(t,ff(2,:));
%--------------------------------------------
fs=1000;
N=1000;
%ts=1/fs;
n=0:N-1;
t=n/fs;
x=@(t) (1+0.5*cos(2*pi*7*t)).*cos(2*pi*100*t+1.5*sin(2*pi*9*t))+0.5*sin(2*pi*30*t);
%t=0:0.001:1-0.001;
y=x(t);
figure;
plot(t,y);
%--------------------------
[pf,a,si]=lmd2(y);
ff=shunpin(si).*fs;%实际频率要乘以采样率
subplot(311);plot(t,pf(1,:));
subplot(312);plot(t,pf(2,:));
subplot(313);plot(t,pf(3,:));
% subplot(514);plot(t,pf(4,:));
% subplot(515);plot(t,pf(5,:));
figure;
subplot(211);plot(t,si(1,:));
subplot(212);plot(t,si(2,:));
% subplot(313);plot(t,si(3,:));
% subplot(414);plot(t,si(4,:));
figure;
subplot(211);plot(t,a(1,:));
a1=a(1,:)-mean(a(1,:));
%subplot(212);hua_fft1(an,fs);
subplot(212);hua_xihua(a1,fs,0,100,256,2);
figure;
%a1_f=a1_f-mean(a1_f);
f1=ff(1,:)-mean(ff(1,:));
subplot(211);plot(t,ff(1,:));
%subplot(212);hua_fft1(a1_f,fs);
subplot(212);hua_xihua(f1,fs,0,100,256,2);
figure;
pinduan=300;
subplot(211);hua_baoluo(pf(1,:),fs,1);
subplot(212);hua_baoluo(pf(2,:),fs,1);
%------------------------------------------