plots2com.m

算荧光衰减时间或者寿命的软件.matlab语言编写.

function [ExportData, RedChib,DWb]= plots2com(FittedDecayParams,Time,pp,InstrumentResponse,Data,fid2)

%%%%%%%%%%%%%
Amp1 = FittedDecayParams(1); Amp2 = FittedDecayParams(2);
Life1 = FittedDecayParams(3); Life2 = FittedDecayParams(4);
iback = FittedDecayParams(5);
Dback = FittedDecayParams(6);
TimeShift = FittedDecayParams(7);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

TotalDecay = lifeFunc2com(FittedDecayParams,Time,pp);

Residb = (1./sqrt(Data)).*(TotalDecay-Data);
ResidSqb = Residb.^2;
RedChib = (sum(ResidSqb))./length(Data); RedChib = RedChib(1,1);
ResDifSqb = (diff(Residb)).^2;
DWb = sum(ResDifSqb)/sum(ResidSqb);

figure;
subplot('Position',[0.1 0.40 0.8 0.55]);
semilogy(Time,TotalDecay,'r');
hold on;
semilogy(Time,InstrumentResponse);
semilogy(Time,Data,'g');
ylabel('Counts');axis tight;
title(['Life1: ' num2str(Life1) '   Amp1: ' num2str(Amp1) '     Life2: ' num2str(Life2) '   Amp2: ' num2str(Amp2)], 'BackgroundColor',[.7 .9 .7]);


subplot('Position',[0.1 0.1 0.8 0.2]);
plot(Time,Residb);
xlabel('Time (ns)');
ylabel('Residual');axis tight;
title(['Red Chi: ' num2str(RedChib) '       DW: ' num2str(DWb)], 'BackgroundColor',[.7 .9 .7]);

fprintf(fid2,'\n');
fprintf(fid2,'COMPONENT 1 \n');
fprintf(fid2,'  Lifetime: %d',Life1);
fprintf(fid2,'  Amplitude: %d',Amp1);
yield1b = 100*Amp1*Life1/(Amp1*Life1+Amp2*Life2);
fprintf(fid2,'  yield: %d\n',yield1b);

fprintf(fid2,'\n');
fprintf(fid2,'COMPONENT 2 \n');
fprintf(fid2,'  Lifetime: %d',Life2);
fprintf(fid2,'  Amplitude: %d',Amp2);
yield2b = 100*Amp2*Life2/(Amp1*Life1+Amp2*Life2);
fprintf(fid2,'  yield: %d\n',yield2b);

fprintf(fid2,'\n');
AverageLifetimeb = yield1b*Life1/100 + yield2b*Life2/100;
fprintf(fid2,'Average Lifetime: %d\n',AverageLifetimeb);
fprintf(fid2,'\n');


fprintf(fid2,'\n');
fprintf(fid2,'STATS \n');
fprintf(fid2,'Background Correction Intensity: %d\n',Dback);
fprintf(fid2,'Instrument Response Correction Intensity: %d\n',iback);
fprintf(fid2,'Shift in Instrument Response: %d\n',TimeShift);
fprintf(fid2,'Reduced Chi Sq: %d\n',RedChib);
fprintf(fid2,'Derb Wat: %d\n',DWb);

ExportData = [Time TotalDecay InstrumentResponse Data Residb];