Geology Reference
In-Depth Information
armdinthighspline=interp1(strat,arm,dinthigh,&rsq
uo;spline’);
figure;plot(strat,arm); hold all; plot(dinthigh,ar
mdinthighspline);
Figure 4.4
Smoothing options for the Δd = 0.05 m linearly interpolated
ARM stratigraphic series: A Δd = 0.05 m linearly interpolated Arguis ARM
series is smoothed in three different ways over an 80m window using
smooth.m and plotted together.
armdinthigh=interp1(strat,arm,dinthigh);
span=80/0.05;
armmoving=smooth(armdinthigh,span,'moving');
armloess=smooth(armdinthigh,span,'loess');
armrloess=smooth(armdinthigh,span,'rloess');
figure; plot(dinthigh,arminthigh,'red');
hold all; plot(dinthigh,armrloess,'grey');
hold all; plot(dinthigh,armmoving,'green');
hold all; plot(dinthigh,datalowess,'blue');
hold all; plot(dinthigh,dataloess,'black');
Figure 4.6
Gauss and Taner filter bode plots. Create an impulse function of
length 4096 using padarray.m, set dt = 1 and passband to fl = 0.035, fc = 0.050,
fh = 0.065, pass through gaussfilter.m and tanerfilter.m, then apply fft.m and
calculate modulus and phase plots.
impulse=padarray(1,4095,'post');
fc=0.05;fl=fc-0.015;fh=fc+0.015;dt=1.;
gaussbandx=gaussfilter(impulse,dt,fc,fl,fh);
tanerbandx=tanerfilter(impulse,dt,fc,fl,fh);
tanerfft=fft(complex(tanerbandx));
tanermag=abs(tanerfft.^2);
tanerfaz=angle(tanerfft);
f=0.:1/4096:1-1/4096;
figure;plot(f,tanermag);figure;plot(f,tanerfaz);
gaussfft=fft(complex(gaussbandx));
gaussmag=abs(gaussfft.^2);
gaussfaz=angle(gaussfft);
figure;plot(f,gaussmag);figure;plot(f,gaussfaz);
Figure 4.7
Application of Gauss and Taner filters to extract precession
index from insolation: In
Analyseries 2.0.4.2
, with the La2004 model,
create 65° North mean summer insolation from March 21 to September
21, 36,000-40,000 ka with a 1 kyr step; select and save data. Create
and save precession index for the same time interval using the same
procedure.
% Import insolation file as 'data' (use File menu)
time=data(:,1);
summer=data(:,2);
figure;plot(time,summer); % Figure 4.7a
