这样算是完成了MDSC的数据处理吗? 跟Q100结果对比...

%function y=simulation_MDSC;
close all
clear
clc

load MDSCdata1
N=120;
L=length(data(:,1));
t=data(:,1)*60; % Convert time unit from minute to second
for i=1:15,
    [num2str(i),' : ',textdata{i}]
end
AmplitudeHF=amplitude(data(:,13));
figure, set(gcf,'color','w')
plot(AmplitudeHF),hold on, plot(data(N/2:L,10),'r:')
legend('Result by Matlab','Q100 Result'), legend boxoff
 title('Amplitude of Heat Flow')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')

ModulatedT=data(:,12)-smooth(data(:,12),N);

figure, set(gcf,'color','w')
Tm_Q100=data(:,12)-data(:,2);
plot(ModulatedT),hold on, plot(Tm_Q100,'r:')
legend('Result by Matlab','Q100 Result'), legend boxoff
 title('Modulation of Temperature')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')
 
dTmdt=mmspder(t,ModulatedT,t);

dTmdt_Q100=mmspder(t,data(:,12),t)-mmspder(t,data(:,2),t);

figure, set(gcf,'color','w')
plot(dTmdt),hold on, plot(dTmdt_Q100,'r:')
legend('Result by Matlab','Q100 Result'),legend boxoff
title('Heating Rate Modulation')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')
 
AmplitudeHR=amplitude(dTmdt,0);%amplitude(mmspder(t,data(:,12),t));%
AmplitudeT=amplitude(data(:,12)).';
OmgB=2*pi*AmplitudeT./(120*diff(t(N/2:L-N/2+1)));

figure, set(gcf,'color','w')
plot(AmplitudeHR),hold on, plot(OmgB,'r:')
legend('Result by Amplitude heating Rate modulation','Result by Omiga * Amplitude of Temperature Modulation')
legend boxoff, title('Amplitude of Heating Rate Modulation')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')
%break
RevCp=AmplitudeHF./(AmplitudeHR+eps);
RevCp=RevCp.';
figure, set(gcf,'color','w')
plot(RevCp),hold on, plot(data(:,4),'r:')
legend('Result by Matlab','Q100 Result'),legend boxoff
title('Reversible Heat Capability')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')
 
Temperature=smooth(data(:,12),N);
%Temperature=Temperature(N/2:L-N/2);
%t=t(N/2:L-N/2);
dTdt=mmspder(t,Temperature,t);
RevHF=-smooth(RevCp.*dTdt(N/2:L-N/2),N);

figure, set(gcf,'color','w')
plot(RevHF),hold on,plot(data(N/2:L-N/2,8),'r:')
legend('Rev Heat Flow by Matlab','Rev Heat Flow by Q100'), legend boxoff
title('Reversible Heat Flow')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')
 
figure, set(gcf,'color','w')
TotalHF=smooth(data(:,13),N);
TotalHF=smooth(TotalHF,N);
plot(TotalHF),hold on, plot(data(:,6),'r:')
legend('Total Heat Flow by Matlab','Total Heat Flow by Q100'),legend boxoff
title('Total Heat Flow')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')
 
figure, set(gcf,'color','w')
NonRevHF=TotalHF(N/2:L-N/2)-RevHF;
plot(NonRevHF),hold on, plot(data(N/2:L-N/2,8),'r:')
legend('Non-Reversible Heat Flow by Matlab','Non-Reversible Heat Flow by Q100')
legend boxoff, title('Non-Reversible Heat Flow')
 set(gcf,'name',get(get(gca,'title'),'string'),'Numbertitle','off')
 
figure, set(gcf,'color','w')
index=N/2:L-N/2;
subplot(1,2,1)
plot(TotalHF(index)),hold on,
plot(RevHF,'m'), plot(NonRevHF,'g')
legend('Total Heat Flow','Reversible Heat Flow','Non-Reversible Heat Flow')
legend boxoff, title('MDSC Heat Flow Deconvolution by Matlab')

%figure, set(gcf,'color','w')
subplot(1,2,2)
plot(data(:,6)),hold on, plot(data(:,7)), plot(data(:,8))
legend('Total Heat Flow','Reversible Heat Flow','Non-Reversible Heat Flow')
legend boxoff, title('MDSC Heat Flow Deconvolution by Q100')

 set(gcf,'name',[get(get(gca,'title'),'string'),' and Matlab'],'numbertitle','off')
%end