Skip to content
This repository has been archived by the owner on May 3, 2020. It is now read-only.

Commit

Permalink
Input plots
Browse files Browse the repository at this point in the history
  • Loading branch information
Jesper Dramsch committed Apr 11, 2013
1 parent d6332e9 commit 1e2e6e1
Showing 1 changed file with 25 additions and 2 deletions.
27 changes: 25 additions & 2 deletions COmig.m
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@
%}

clear all % Clear workspace
close all % Close figures
close all % Close figures
clc % Clear command line window

format long % Double precision
Expand Down Expand Up @@ -106,7 +106,7 @@
set(gca,'Fontsize',24)
set(gca,'XTickMode','manual')
set(gca,'XTick',[0;2000;4000;6000;8000;10000])
set(gca,'XTickLabel',[' 0 ';'2 / 0';'2 / 0';'2 / 0';'2 / 0';' 2 ']) % rescaling x-axis
set(gca,'XTickLabel',[' 0 ';'2 / 0';'2 / 0';'2 / 0';'2 / 0';' 2 ']) % rescaling x-axis
print('-dpng',sprintf('v%g.png',v));

if v == vfinal % If loop reaches the correct velocity (estimated with constant velocity scan)
Expand Down Expand Up @@ -152,6 +152,29 @@

fprintf('Verbesserung der Signal-to-Noise ratio von %f2 auf %f2\n',SNRin,SNRout)

%Input signal normalized
figure
plot(((1:nt)-1)*dt,filtdata(:,51,1)/max(filtdata(:,51,1)),'r')
hold on
plot(((1:nt)-1)*dt,data(:,51,1)/max(data(:,51,1)),'k')
ylabel('Normalisierte Amplitude','Fontsize',24)
xlabel('Zeit [s]','Fontsize',24)
legend('Filtered data','Original data','Location','best')
set(gca,'Fontsize',24)
print('-dpng','wavelet.png');

%Input signal not normalized
figure
plot(((1:nt)-1)*dt,filtdata(:,51,1),'r')
hold on
plot(((1:nt)-1)*dt,data(:,51,1),'k')
ylabel('Amplitude','Fontsize',24)
xlabel('Zeit [s]','Fontsize',24)
legend('Filtered data','Original data','Location','best')
set(gca,'Fontsize',24)
print('-dpng','wavelet_nN.png');


% Fileoutput of datamatrices
dlmwrite('mig.dat',mig)
dlmwrite('COGatherh0.dat',Kirchhoffdepth(:,:,1));
Expand Down

0 comments on commit 1e2e6e1

Please sign in to comment.