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); %------------------------------------------