%PlotSimplotSpe

NP=2048;
SW=63000;

spc1=readmatrix('apostolos_1.txt');
spc2=readmatrix('apostolos_2.txt');
spc3=readmatrix('apostolos_3.txt');

xaxis=linspace(1,NP,NP)*SW/NP - SW/2;

figure(1);
%spc1b= real(  (spc1([1:NP],1)+ i * spc1([1:NP],2) ) *exp(-i*3.4)    );
%plot(xaxis,30*spc1b([1:NP],1));
%hold on
spc2b= real(  (spc2([1:NP],1)+ i * spc2([1:NP],2) ) *exp(i*0.5)    );
plot(xaxis,30*spc2b([1:NP],1));
hold on
plot(xaxis,spc3([1:NP],1));
xlabel('^{14}N Frequency / kHz')
ylabel('Signal intensity / a.u.')
hold off
xlabel('^{14}N Frequency / kHz')
ylabel('Signal intensity / a.u.')

%legend('Glycine Quad constant=1.2 MHz x 30','Surface Nitrade Quad constant=0.5 MHz x 30','Interstitial Nitrade Quad constant=0.006 MHz')
legend('Surface Nitride x 30','Interstitial Nitride')

spcfulljp=sprintf('simpson_nitride_new.jpg')
set(gcf, 'PaperSize', [6 3]);
saveas(gcf,spcfulljp);