diff --git a/MAS_testing.m b/MAS_testing.m deleted file mode 100644 index 934111b..0000000 --- a/MAS_testing.m +++ /dev/null @@ -1,105 +0,0 @@ -clear all -close all -clc - - -signals = loadAllSignals('Supelec_2012_SIR_Spectral_Analysis_EA_v001.mat'); -Fe = 1000000.0; -signal = signals(:,1); -t = (0: 1: length(signal)-1)/Fe; - -%filtering smooth -nb_coeff = 10; -coeff = ones(1, nb_coeff)/nb_coeff; -signalm = filter(coeff, 1, signal); - -%envelope time -signalenv = filter(coeff, 1, abs(hilbert(sgolayfilt(signal, 1, 3)))); -plot_param = {'Color', [0.6 0.1 0.2],'Linewidth',1}; - -% savitzky-golay -signalsg = sgolayfilt(signal, 1, 3); - -%fft -t0pad=2^12;%0-padding -f=((1: t0pad)-1)*Fe/t0pad-Fe/2;%echelle des frequences centrée en 0 -TFx0pad=filter(coeff, 1, fft(signal, t0pad));%fft avec 0-padding -ModulTFx0pad=abs(TFx0pad);%module -CenterModulTFx0pad=fftshift(ModulTFx0pad);%shift zero-frequency component to centrer of spectrum -signalF=20*log10(CenterModulTFx0pad);%log -%fft smooth -TFx0pad=filter(coeff, 1, fft(signalm, t0pad));%fft avec 0-padding -ModulTFx0pad=abs(TFx0pad);%module -CenterModulTFx0pad=fftshift(ModulTFx0pad);%shift zero-frequency component to centrer of spectrum -signalmF=20*log10(CenterModulTFx0pad);%log -%fft hilbert -TFx0pad=filter(coeff, 1, fft(signalenv, t0pad));%fft avec 0-padding -ModulTFx0pad=abs(TFx0pad);%module -CenterModulTFx0pad=fftshift(ModulTFx0pad);%shift zero-frequency component to centrer of spectrum -signalenvF=20*log10(CenterModulTFx0pad);%log -%savitzky-golay -TFx0pad=filter(coeff, 1, fft(signalsg, t0pad));%fft avec 0-padding -ModulTFx0pad=abs(TFx0pad);%module -CenterModulTFx0pad=fftshift(ModulTFx0pad);%shift zero-frequency component to centrer of spectrum -signalsgF=20*log10(CenterModulTFx0pad);%log - -%plot -figure -subplot(221)%%%%%%%%%%%%%%%%%%%%%%% -plot(t, signal) -xlabel('Time (s)') -ylabel('Amplitude (a.u.)') -hold on -plot(t, signalm) -xlim([0 length(signal)/Fe]) -title('smooth filtering') -legend('signal', 'avg signal') -subplot(223)%%%%%%%%%%%%%%%%%%%%%% -plot(t, signal) -xlabel('Time (s)') -ylabel('Amplitude (a.u.)') -hold on -plot(t,[-1,1].*signalenv,plot_param{:}) -plot(t, signalsg) -xlim([0 length(signal)/Fe]) -title('Hilbert Envelope') -legend('signal', 'Hilbert Envelope', 'savitzky-golay') -subplot(222)%%%%%%%%%%%%%%%%%%%%%% -plot(f, signalF) -xlabel('f(Hz)'); -ylabel('Module de la TF en dB'); -grid on -hold on -plot(f, signalmF) -legend('signal', 'avg signal') -title('FFT') -xlim([0 Fe/2]) -subplot(224)%%%%%%%%%%%%%%%%%%%%%% -plot(f, signalF) -xlabel('f(Hz)'); -ylabel('Module de la TF en dB'); -grid on -hold on -plot(f, signalenvF) -plot(f, signalsgF) -legend('signal', 'signal Hilbert', 'savitzky-golay') -title('FFT') -xlim([0 Fe/2]) - -figure -ax1 = subplot(3,1,1); -plot(t,signal) -ax2 = subplot(3,1,2); -pspectrum(signal,Fe,'spectrogram','OverlapPercent',99, 'Leakage',1,'MinThreshold',-150) -colorbar(ax2,'off') -ax3 = subplot(3,1,3); -pspectrum(signal,Fe,'spectrogram','OverlapPercent',99, 'Leakage',1,'MinThreshold',-150,'TimeResolution', 10e-6) -colorbar(ax3,'off') -linkaxes([ax1,ax2,ax3],'x') - -figure -[sp,fp,tp] = pspectrum(signal,Fe,'spectrogram','OverlapPercent',99, 'Leakage',1); -mesh(tp,fp,sp) -view(-15,60) -xlabel('Time (s)') -ylabel('Frequency (Hz)')