emd_separation.m ( File view )

  • By bole 2014-10-30
  • View(s):92
  • Download(s):2
  • Point(s): 1
			% EMD_SEPARATION.M
%
% P. Flandrin, Mar. 13, 2003 - modified Mar. 2, 2006
%
% computes and displays an error measure in the EMD
% separation of two tones
%
% produces Figure 4 in
%
% G. Rilling, P. Flandrin and P. Gonalvs
% "On Empirical Mode Decomposition and its algorithms"
% IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing
% NSIP-03, Grado (I), June 2003

N = 256;% # of data samples
t = 1:N;
tt = N/4:3*N/4;

Nf = 129;% # of tested fequencies
f = linspace(0,.5,Nf);

rapp = 8;% amplitude ratio between modes
a1 = sqrt(rapp);
a2 = 1/a1;

x1 = a1*cos(2*pi*f'*t);
x2 = a2*cos(2*pi*f'*t);

se = zeros(Nf);

for k1 = 5:Nf-1
	
	for k2 = 2:k1-1;
		
		y1 = x1(k1,:);
		y2 = x2(k2,:);
		
		sy1 = sum((y1(tt)).^2);
		sy2 = sum((y2(tt)).^2);
		sy = sum((y1(tt)).^2+(y2(tt)).^2);
	
		imf = emd(y1+y2);
        if size(imf,1) == 1
            imf(2,:) = 0;
        end
		se(k1,k2) = sqrt((sy1*sum((imf(1,tt)-y1(tt)).^2) + sy2*sum((imf(2,tt)-y2(tt)).^2))/(sy1+sy2)/sy);	

		[k1 k2 size(imf)]
				
	end

end

imagesc(f,f,flipud(se'))
axis('square')
axis([0 .5 0 .5])
hold on
plot([f(2) f(Nf-1)],[f(Nf-1) f(Nf-1)])
plot([f(Nf-1) f(Nf-1)],[f(2) f(Nf-1)])
plot([f(2) f(Nf-1)],[f(Nf-1) f(2)])
set(gca,'YTick',[]);set(gca,'XTick',[])
xlabel('f_1 > f_2')
ylabel('f_2')
colormap(jet)
%colormap(flipud(gray))
colorbar('vert')
hold off
			
...
Expand> <Close

Want complete source code? Download it here

Point(s): 1

Download
0 lines left, continue to read
Sponsored links

File list

Tips: You can preview the content of files by clicking file names^_^
Name Size Date
bugfix.sh216.00 B27-06-07 22:12
cemdc.m2.30 kB27-06-07 22:12
cemdc2.m2.31 kB27-06-07 22:12
cemdc2_fix.m2.26 kB27-06-07 22:12
cemdc_fix.m2.25 kB27-06-07 22:12
emd.m21.75 kB27-06-07 22:12
emdc.m2.23 kB27-06-07 22:12
emdc_fix.m2.09 kB27-06-07 22:12
emd_local.m9.67 kB27-06-07 22:12
emd_online.m25.70 kB27-06-07 22:12
make_emdc.m1.33 kB27-06-07 22:12
cemdc.c4.85 kB27-06-07 22:12
cemdc2.c4.87 kB27-06-07 22:12
cemdc2_fix.c4.11 kB27-06-07 22:12
cemdc_fix.c4.13 kB27-06-07 22:12
cextr.c14.37 kB27-06-07 22:12
cextr.h770.00 B27-06-07 22:12
cio.c7.38 kB27-06-07 22:12
cio.h1.16 kB27-06-07 22:12
cio_fix.c6.77 kB27-06-07 22:12
cio_fix.h1.09 kB27-06-07 22:12
clocal_mean.c4.97 kB27-06-07 22:12
clocal_mean.h790.00 B27-06-07 22:12
clocal_mean2.c4.79 kB27-06-07 22:12
clocal_mean2.h752.00 B27-06-07 22:12
emdc.c4.60 kB27-06-07 22:12
emdc_fix.c3.43 kB27-06-07 22:12
emd_complex.c508.00 B27-06-07 22:12
emd_complex.h512.00 B27-06-07 22:12
extr.c10.23 kB27-06-07 22:12
extr.h674.00 B27-06-07 22:12
interpolation.c2.06 kB27-06-07 22:12
interpolation.h468.00 B27-06-07 22:12
io.c6.86 kB27-06-07 22:12
io.h1.14 kB27-06-07 22:12
io_fix.c6.04 kB27-06-07 22:12
io_fix.h1.07 kB27-06-07 22:12
local_mean.c2.77 kB27-06-07 22:12
local_mean.h710.00 B27-06-07 22:12
emd_fmsin.m1.82 kB27-06-07 22:12
emd_sampling.m896.00 B27-06-07 22:12
emd_separation.m1.30 kB27-06-07 22:12
emd_triang.m704.00 B27-06-07 22:12
ex_online.m2.02 kB27-06-07 22:12
triangular_signal.m378.00 B27-06-07 22:12
bivariate_EMD_illustration.m1.25 kB27-06-07 22:12
bivariate_EMD_mean_definitions.m2.91 kB27-06-07 22:12
bivariate_EMD_principle.m3.84 kB27-06-07 22:12
dirstretch.m614.00 B27-06-07 22:12
float_position_record.mat8.34 kB27-06-07 22:12
index_emd.m4.38 kB27-06-07 22:12
install_emd.m1.53 kB27-06-07 22:12
ls-R1,004.00 B27-06-07 22:12
revert_bugfix.sh216.00 B27-06-07 22:12
uninstall_emd.m2.08 kB27-06-07 22:12
addtag.m791.00 B27-06-07 22:12
boundary_conditions_emd.m3.46 kB27-06-07 22:12
cemd_disp.m2.49 kB27-06-07 22:12
cemd_visu.m2.76 kB27-06-07 22:12
cenvelope.m2.00 kB27-06-07 22:12
dirstretch.m1.22 kB27-06-07 22:12
disp_hhs.m1.55 kB13-09-14 22:55
emd_visu.m3.14 kB27-06-07 22:12
extr.m1.89 kB27-06-07 22:12
findtag.m1.02 kB27-06-07 22:12
hastag.m659.00 B27-06-07 22:12
hhspectrum.m1.41 kB27-06-07 22:12
io.m504.00 B27-06-07 22:12
plot3c.m798.00 B27-06-07 22:12
plotc.m3.17 kB27-06-07 22:12
rmtag.m820.00 B27-06-07 22:12
toimage.m2.93 kB27-06-07 22:12
<src>0.00 B29-10-14 13:58
<NSIP2003>0.00 B29-10-14 13:58
<SPL2007>0.00 B29-10-14 13:58
<EMDs>0.00 B29-10-14 13:58
<examples>0.00 B29-10-14 13:58
<utils>0.00 B29-10-14 13:58
<package_emd>0.00 B29-10-14 13:58
...
Sponsored links

emd_separation.m (91.07 kB)

Need 1 point
Your Point(s)

Your Point isn't enough.

Get point immediately by PayPal

More(Debit card / Credit card / PayPal Credit / Online Banking)

Submit your source codes. Get more point

LOGIN

Don't have an account? Register now
Need any help?
Mail to: support@codeforge.com

切换到中文版?

CodeForge Chinese Version
CodeForge English Version

Where are you going?

^_^"Oops ...

Sorry!This guy is mysterious, its blog hasn't been opened, try another, please!
OK

Warm tip!

CodeForge to FavoriteFavorite by Ctrl+D