alumdt=1E-8;
alumt0=0+alumdt;
alumFs=1/alumdt;
L=9999;
alumt=(0:L-1)*alumdt;
alumAmplP=alumdata (9:10008,5);
NFFT=2^nextpow2(L);
alumfP=alumFs/2*linspace (0,1,NFFT/2+1);
alumAmplP_2=alumAmplP (1:L);
alumYP=fft (alumAmplP_2,NFFT)/L;
absalumYP=abs(alumYP);
AlumSpectrum = [transpose(alumfP),2*absalumYP(1:NFFT/2+1)];
absoluteRange = (0:1500:1.25*10^8);
alumYYP = interp1(transpose(AlumSpectrum(1:8193,1)),AlumSpectrum(1:8193,2),absoluteRange);
alumPick = find(alumAmplP > .025, 1, 'first') - 200;
alumPickEnd = alumPick + 600;
alumTimePick = (alumPick*alumdt:alumdt:alumPickEnd*alumdt);
ClippedAlumWave = alumAmplP(alumPick:alumPickEnd);
ClippedAlumT = (alumTimePick);
TaperClippedAlumWave = window(@hamming,numel(ClippedAlumT)).*transpose(ClippedAlumWave)';
alumL2 = length(alumTimePick);
NFFT=2^nextpow2(L);
alumfP=alumFs/2*linspace (0,1,NFFT/2+1);
ClippedalumAmplP_2=TaperClippedAlumWave;
ClippedalumYP=fft (ClippedalumAmplP_2,NFFT)/alumL2;
absClippedalumYP = abs(ClippedalumYP);
ClippedAlumSpectrum = [transpose(alumfP),2*absClippedalumYP(1:NFFT/2+1)];
absoluteRange = (0:1500:1.25*10^8);
ClippedalumYYP = interp1(transpose(ClippedAlumSpectrum(1:8193,1)),ClippedAlumSpectrum(1:8193,2),absoluteRange);
figure(1)
subplot(2,1,1)
hold on
plot(alumt, alumAmplP(1:9999))
title('Aluminum Wave-Form')
xlabel('Time (s)')
ylabel('Amplitude')
hold on
line([alumPick*alumdt alumPick*alumdt], [-.15 .15])
line([alumPickEnd*alumdt alumPickEnd*alumdt], [-.15 .15])
subplot(2,1,2)
hold on
plot(absoluteRange,alumYYP)
title('Aluminum Spectrum')
xlabel('Frequency (Hz)')
ylabel('Amplitude')
xlim([0 1.5E6])
figure(2)
subplot(2,1,1)
hold on
plot(ClippedAlumT,TaperClippedAlumWave)
title('Aluminum Wave-Form Initial Peak')
xlabel('Time (s)')
ylabel('Amplitude')
subplot(2,1,2)
hold on
plot(absoluteRange, ClippedalumYYP)
title('(Clip) Aluminum Spectrum')
xlabel('Frequency (Hz)')
ylabel('Amplitude')
xlim([0 1.5E6])
