HVD (Simulink)

Fs=1; % Sampling Frequency
N=200;% Remez filter length
n=2048;
t=((0:n-1)/Fs)’;

%% Example 1 (Denoising)
%u=0.2*sin(0.1*t)+ 0.02*randn(n,1);

%% Example 2 (Two harmonics)
u=sin(0.05*t)+ 0.5*sin(0.3*(t));

[t,u1,yout]=sim(‘s16’);
y=yout(:,1);
A=yout(:,5);
w=yout(:,6);

yi=yout(:,2);
Ai=yout(:,3);
wi=yout(:,4);
yi_2=yout(:,7);

figure(1)
subplot(311)
plot([y w A])
title(‘Composition’)
legend(‘Signal’,’Frequency’, ‘Envelope’)

subplot(312)
plot([yi wi Ai])
title(‘Largest component’)

legend(‘Signal component’,’Frequency’, ‘Envelope’)

subplot(313)
plot([yi_2])
title(‘Residual’)
xlabel(‘Time’)
legend(‘Next component’)