matlab处理信息代码

clearvars; clc; % plotCMZ.m
inter_raw_data_CMZ=importdata("rawData\t_vs_cm_inter.csv");
tmp_1=zeros(size(inter_raw_data_CMZ,1),1);
for ii=1:size(inter_raw_data_CMZ,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,'\t');
    tmp_2=str2double(tmp_2(1:1));
    tmp_1(ii,1)=tmp_2;
end
tmp_1=unique(tmp_1);
inter_CMZ=zeros(size(tmp_1,1),4);
for ii=1:size(tmp_1,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,' ');
    tmp_2=tmp_2(1,2:size(tmp_2,2));
    tmp_2=str2double(tmp_2);
    inter_CMZ(ii,1)=tmp_1(ii,1);
    inter_CMZ(ii,2)=tmp_2(1,1);
    inter_CMZ(ii,3)=tmp_2(1,2);
    inter_CMZ(ii,4)=tmp_2(1,3);
end
clear inter_raw_data_CMZ tmp_2 tmp_1 ii;
inter_raw_data_CMZ=importdata("rawData\t_vs_cm_over.csv");
tmp_1=zeros(size(inter_raw_data_CMZ,1),1);
for ii=1:size(inter_raw_data_CMZ,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,'\t');
    tmp_2=str2double(tmp_2(1:1));
    tmp_1(ii,1)=tmp_2;
end
tmp_1=unique(tmp_1);
over_CMZ=zeros(size(tmp_1,1),4);
for ii=1:size(tmp_1,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,' ');
    tmp_2=tmp_2(1,2:size(tmp_2,2));
    tmp_2=str2double(tmp_2);
    over_CMZ(ii,1)=tmp_1(ii,1);
    over_CMZ(ii,2)=tmp_2(1,1);
    over_CMZ(ii,3)=tmp_2(1,2);
    over_CMZ(ii,4)=tmp_2(1,3);
end
clear inter_raw_data_CMZ tmp_2 tmp_1 ii;
inter_raw_data_CMZ=importdata("rawData\t_vs_cm_thinc.csv");
tmp_1=zeros(size(inter_raw_data_CMZ,1),1);
for ii=1:size(inter_raw_data_CMZ,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,'\t');
    tmp_2=str2double(tmp_2(1:1));
    tmp_1(ii,1)=tmp_2;
end
tmp_1=unique(tmp_1);
thincp_CMZ=zeros(size(tmp_1,1),4);
for ii=1:size(tmp_1,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,' ');
    tmp_2=tmp_2(1,2:size(tmp_2,2));
    tmp_2=str2double(tmp_2);
    thincp_CMZ(ii,1)=tmp_1(ii,1);
    thincp_CMZ(ii,2)=tmp_2(1,1);
    thincp_CMZ(ii,3)=tmp_2(1,2);
    thincp_CMZ(ii,4)=tmp_2(1,3);
end
clear inter_raw_data_CMZ tmp_2 tmp_1 ii;
inter_raw_data_CMZ=importdata("rawData\t_vs_cm_overThinc.csv");
tmp_1=zeros(size(inter_raw_data_CMZ,1),1);
for ii=1:size(inter_raw_data_CMZ,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,'\t');
    tmp_2=str2double(tmp_2(1:1));
    tmp_1(ii,1)=tmp_2;
end
tmp_1=unique(tmp_1);
overThinc_CMZ=zeros(size(tmp_1,1),4);
for ii=1:size(tmp_1,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=cell2mat(tmp_2);
    tmp_2=convertCharsToStrings(tmp_2);
    tmp_2=strsplit(tmp_2,' ');
    tmp_2=tmp_2(1,2:size(tmp_2,2));
    tmp_2=str2double(tmp_2);
    overThinc_CMZ(ii,1)=tmp_1(ii,1);
    overThinc_CMZ(ii,2)=tmp_2(1,1);
    overThinc_CMZ(ii,3)=tmp_2(1,2);
    overThinc_CMZ(ii,4)=tmp_2(1,3);
end
clear inter_raw_data_CMZ tmp_2 tmp_1 ii;
inter_raw_data_CMZ=importdata("rawData\t_vs_cm_thincRK2.csv");
tmp_1=zeros(size(inter_raw_data_CMZ,1),1);
for ii=1:size(inter_raw_data_CMZ,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=tmp_2(1:1);
    tmp_1(ii,1)=tmp_2;
end
tmp_1=unique(tmp_1);
thincr_CMZ=zeros(size(tmp_1,1),4);
for ii=1:size(tmp_1,1)
    tmp_2=inter_raw_data_CMZ(ii,:);
    tmp_2=tmp_2(1,2:size(tmp_2,2));
    thincr_CMZ(ii,1)=tmp_1(ii,1);
    thincr_CMZ(ii,2)=tmp_2(1,1);
    thincr_CMZ(ii,3)=tmp_2(1,2);
    thincr_CMZ(ii,4)=tmp_2(1,3);
end
clear inter_raw_data_CMZ tmp_2 tmp_1 ii;
% 
% inter_raw_data_CMZ=importdata("rawData\t_vs_cm_thincRK2.csv");
% tmp_1=zeros(size(inter_raw_data_CMZ,1),1);
% tmp_3=zeros(size(inter_raw_data_CMZ,1),1);
% for ii=1:size(inter_raw_data_CMZ,1)
%     tmp_2=inter_raw_data_CMZ(ii,:);
%     tmp_2=cell2mat(tmp_2);
%     tmp_2=convertCharsToStrings(tmp_2);
%     tmp_2=strsplit(tmp_2,'\t');
%     tmp_2=str2double(tmp_2(1:1));
%     tmp_1(ii,1)=tmp_2;
%     tmp_4=inter_raw_data_CMZ(ii,:);
%     tmp_4=cell2mat(tmp_4);
%     tmp_4=convertCharsToStrings(tmp_4);
%     tmp_4=strsplit(tmp_4,' ');
%     if (size(tmp_4,2)<2)
%         break;
%     else
%         tmp_4=tmp_4(1,size(tmp_4,2));
%     end
%     tmp_3(ii,1)=str2double(tmp_4);
% end
% tmp_1=unique(tmp_1);
% tmp_1=tmp_1(2:size(tmp_1,1),:);
% indexJJ=0;
% for ii=1:size(tmp_3,1)
%     if (tmp_3(ii,1)<0.01)
%         break;
%     else
%         indexJJ=indexJJ+1;
%     end
% end
% tmp_3=tmp_3(1:indexJJ,:);
% tmp_5=zeros(round(size(tmp_3,1)/2),1);
% for ii=1:size(tmp_5,1)
%     tmp_5(ii,1)=(tmp_3(2*ii,1)+tmp_3(2*ii-1,1))/2;
% end
% thincr_CMZ=zeros(size(tmp_1,1),4);
% for ii=1:size(tmp_1,1)
%     tmp_2=inter_raw_data_CMZ(ii,:);
%     tmp_2=cell2mat(tmp_2);
%     tmp_2=convertCharsToStrings(tmp_2);
%     tmp_2=strsplit(tmp_2,' ');
%     tmp_2=tmp_2(1,2:size(tmp_2,2));
%     tmp_2=str2double(tmp_2);
%     thincr_CMZ(ii,1)=tmp_1(ii,1);
%     thincr_CMZ(ii,2)=tmp_2(1,1);
%     thincr_CMZ(ii,3)=tmp_2(1,2);
%     thincr_CMZ(ii,4)=tmp_5(ii,1);
% end
% clear inter_raw_data_CMZ tmp_2 tmp_1 ii tmp_3 tmp_4 tmp_5;
fig1=figure(1);
% x_to_plot=inter_CMZ(:,1); y_to_plot=inter_CMZ(:,4);
% plt1=plot(x_to_plot,y_to_plot);plt1.LineWidth=1.5;plt1.Color='blue';hold on;
% x_to_plot=over_CMZ(:,1); y_to_plot=over_CMZ(:,4);
% plt2=plot(x_to_plot,y_to_plot);plt2.LineWidth=1.5;plt2.Color='red';hold on;
x_to_plot=thincp_CMZ(:,1); y_to_plot=thincp_CMZ(:,4);
plt3=plot(x_to_plot,y_to_plot);plt3.LineWidth=1.5;plt3.Color='magenta';hold on;
x_to_plot=thincr_CMZ(:,1); y_to_plot=thincr_CMZ(:,4);
plt4=plot(x_to_plot,y_to_plot);plt4.LineWidth=1.5;plt4.Color='black';hold on;
x_to_plot=overThinc_CMZ(:,1); y_to_plot=overThinc_CMZ(:,4);
plt5=plot(x_to_plot,y_to_plot);plt5.LineWidth=1.5;plt5.Color='green';hold on;
xL1=xlabel('time (s)'); xL1.FontSize=15; yL1=ylabel('Height (m)'); yL1.FontSize=15;
% leg1=legend(); leg1.FontSize=17; leg1.String=["inter","overInter","thinc-p","thinc-RK2","overThinc"];
leg1=legend([plt3,plt4,plt5]); leg1.FontSize=17; leg1.String=["thinc-CN","thinc-RK2","overThinc"];
saveas(gcf,"image\poltCMZ.png");