Hilbert+Huang 源代码

%产生非平稳信号 sig=real(amgauss(128).*fmlin(128)); %%时域波形 figure(1) plot(sig,'LineWidth',2); xlabel('时间 t'); ylabel('幅值 A'); %设置窗…

%产生非平稳信号 sig=real(amgauss(128).*fmlin(128)); %%时域波形 figure(1) plot(sig,'LineWidth',2); xlabel('时间 t'); ylabel('幅值 A'); %设置窗…

function Mf3=Mf3_func(s,fs,Ac) h=hilbert(s); f=dif(h,fs); f0=f0_func(f); Mf3=M_func(A_func(A_func(f0)));

function Mf2=Mf2_func(s,fs,Ac) h=hilbert(s); f=dif(h,fs); f0=f0_func(f); Mf2=M_func(A_func(f0));

function Mp1=Mp1_func(s,fs,fc) h=hilbert(s); fai=fai0_func(h,fs,fc); p1=A_func(fai); Mp1=M_func(p1);

function Mp2=Mp2_func(s,fs,fc) h=hilbert(s); fai=fai0_func(h,fs,fc); p1=A_func(fai); Mp2=M_func(A_func(p1));

%单脊提取 clc; clear all; a(1)=25; Te=0.001; t=0:0.001:1; s=cos(2*pi*(30*(t.^2)+40*t)); %改了一下原信号 %s=hilbert(s); %hudu=4*pi; …

ts = 1000;%采样率 t = 0:1/ts:1; x = sin(2*pi*2*t)+2;%原始信号，频率2Hz x2 = sin(2*pi*50*t);%载波，频率50Hz y = x.*x2;%调制 plot(y);%调制信号 y = …

function u=lssb_mod(m,ts,fc) % u=lssb_mod(m,ts,fc) %LSSB_MOD takes signal m sampled at ts and carrier % freq. fc as input a…

function z=analytic(x) % z=analytic(x) %ANALYTIC Returns the analytic signal corresponding to signal x. % z=hilbert(x);