Plot_sin3source_wavelet.m 1.67 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79
clear all
%close all

% Definition von Parametern
% Homogener Vollraum
npts=10000;
delta=0.0001;




t=(0:delta:(npts-1)*delta);

% Parameter: Homogener Vollraum
%T_d=0.004;
%T_d=0.0033;
T_d=1/300.0;
%T_d=0.002;
F0=1;

npts_start=1;
npts_end=T_d/delta;

wavelet=zeros(1,npts);


% Definition des Quellwavelets
wavelet(1:npts_end)=F0*(sin((pi*t(1:npts_end))/T_d)).^3;

%figure
%plot(t,wavelet)
%xlabel('Zeit in s')

WAVELET=delta*fft(wavelet);
maximum=max(max(abs(WAVELET)))
WAVELET=WAVELET./maximum;
frequenz=[0:fix(npts/2),-ceil(npts/2)+1:-1]'/(npts*delta);  %Definition des kompletten Frequenzvektors
f=[0:fix(npts/2)]'/(npts*delta);

%{
figure
subplot(211)
plot(f(1:fix(npts/2)+1),abs(WAVELET(1:fix(npts/2)+1)))
xlab=xlabel('Frequenz in Hz')
set(gca,'FontSize',14)
tit=title(['Laenge des Quellsignals: ' num2str(T_d) ' s']);
set(xlab,'FontSize',14)

subplot(212)
plot(f(1:fix(npts/2)+1),abs(WAVELET(1:fix(npts/2)+1)))
xlab=xlabel('Frequenz in Hz')
set(gca,'FontSize',14)
ylim([0 0.01])
%set(gca,'XTick',[0 30 50 70 90 110 150 200 250])
%set(gca,'XTickLabel',{'0';'30';'50';'70';'90';'110';'150';'200';'250'})
set(xlab,'FontSize',14)
%}


figure
plot(t,wavelet,'-k')
xlab=xlabel('Time in s');
ylab=ylabel('s(t) (dimensionless)');
set(gca,'FontSize',18)
set(xlab,'FontSize',20)
set(ylab,'FontSize',20)
xlim([0 1.5*T_d])


figure
plot(f(1:fix(npts/2)+1),abs(WAVELET(1:fix(npts/2)+1)),'-k')
xlab=xlabel('Frequency in Hz');
ylab=ylabel('|S(\omega)| in 1/Hz');
set(gca,'FontSize',18)
set(xlab,'FontSize',20)
set(ylab,'FontSize',20)
%tit=title(['Laenge des Quellsignals: ' num2str(T_d) ' s']);
%set(xlab,'FontSize',14)
xlim([0 600])