How to sort table more than one variable with specified condition ?

Question

Arif 2024 年 4 月 21 日

0
リンク

この質問への直接リンク

https://jp.mathworks.com/matlabcentral/answers/2109761-how-to-sort-table-more-than-one-variable-with-specified-condition

回答済み: Steven Lord 2024 年 4 月 21 日

MATLAB Online で開く

hi guys, can u tell me a different way to sort table by 2 variable columns which more simple than my method?

I want to sort jointTOS by its X and Z from minimum to maximum but X is the priority.

Here the code :

clear;

clc;

H = 5;

L = 15;

Ns = 2;

Nb = 1;

PortalHeight = H;

Bayspan = L;

NumberStory = Ns;

NumberBay = Nb;

% Variable1

h = [0.8,1.00,1.20,1.40,1.6,1.80];

% Variable2

Lss = [2.00,2.40,2.80,3.20,3.60,4.00];

% Variable3

Nss = [2,3,4];

% Definition

Var1HeightofTrussBeam = h;

Var2LengthSpecialSegment = Lss;

Var3NumberSpecialSegment = Nss;

% % Xbase = zeros(1,length(Totalvariant));

% % Zbase = zeros(1,length(Totalvariant));

% % Xcol = zeros(1,length(Totalvariant));

% % Zcol = zeros(1,length(Totalvariant));

% % Xos = zeros(1,length(Totalvariant));

% % Zos = zeros(1,length(Totalvariant));

% % Xss = zeros(1,length(Totalvariant));

% % Zss = zeros(1,length(Totalvariant));

% % X = zeros(1,length(Totalvariant));

% % Z = zeros(1,length(Totalvariant));

% % jointbase = zeros(1,length(Totalvariant));

% % jointcolumn = zeros(1,length(Totalvariant));

% % jointordinarysegment = zeros(1,length(Totalvariant));

% % jointtopordinarysegment = zeros(1,length(Totalvariant));

% % jointbottomordinarysegment = zeros(1,length(Totalvariant));

% % jointspecialsegment = zeros(1,length(Totalvariant));

% % jointtopspecialsegment = zeros(1,length(Totalvariant));

% % jointbottomspecialsegment = zeros(1,length(Totalvariant));

ijk = 0;

for i = 1:numel(h)

for j = 1:numel(Lss)

for k = 1:numel(Nss)

ijk = ijk+1;

Variantmodel{ijk,:} = table(ijk,h(i),Lss(j),Nss(k),'VariableNames',{'VariantID','Height','Length','Number'});

%% Cooridinate Point Base(base)

Xbase{ijk,:} = L*(0:NumberBay);

Zbase{ijk,:} = repelem(0,1,length(Xbase{ijk,:}));

%% Column(col)

Xcolumn = L*(0:NumberBay);

Zcolumn = repelem(H*(1:NumberStory),1,length(Xcolumn));

Xbotbeam = L*(0:NumberBay);

Zbotbeam = repelem((H*(1:NumberStory)-h(i)),1,length(Xbotbeam));

% replicate [X] to be the same size as [Z] for plotting:

Xcol{ijk,:} = repmat([Xcolumn,Xbotbeam],1,NumberStory);

Zcol{ijk,:} = [Zcolumn,Zbotbeam];

%% Cooridinate Point Ordinary Segment (os)

% calculate span length of ordinary segment

Los = (L-Lss(j))/2;

Nos = Los/h(i);

Nos = round(Nos);

los = Los/Nos;

totalspan{ijk,:} = los*Nos*2+Lss(j);

if totalspan{ijk,:} == 15.000

checkspan{ijk,:} = "correct span";

else

checkspan{ijk,:} = "span must be 15m";

end

% top and bottom joint of ordinary segment(os)

XTopOrdinarySegment = reshape([los*(1:Nos),(L/2+Lss(j)/2)+los*(0:Nos-1)].'+L*(0:(NumberBay-1)),1,[]);

XBottomOrdinarySegment = XTopOrdinarySegment;

ZTopOrdinarySegment = repelem(H*(1:NumberStory),1,2*Nos*NumberBay);

ZBottomOrdinarySegment = repelem(H*(1:NumberStory)-h(i),1,2*Nos*NumberBay);

% replicate [X] to be the same size as [Z] for plotting:

Xos{ijk,:} = repmat([XTopOrdinarySegment,XBottomOrdinarySegment],1,NumberStory);

Zos{ijk,:} = [ZTopOrdinarySegment,ZBottomOrdinarySegment];

%% Cooridinate Point Special Segment (ss)

% calculate span length of special segment

lss = Lss(j)/Nss(k);

% top and bottom joint of special segment(ss)

XTopSpecialSegment = reshape((lss*(0:Nss(k))+L/2-Lss(j)/2).'+L*(0:(NumberBay-1)),1,[]);

XBottomSpecialSegment = XTopSpecialSegment;

ZTopSpecialSegment = repelem(H*(1:NumberStory),1,length(XTopSpecialSegment));

ZBottomSpecialSegment = repelem(H*(1:NumberStory)-h(i),1,length(XBottomSpecialSegment));

% replicate [X] to be the same size as [Z] for plotting:

Xss{ijk,:} = repmat([XTopSpecialSegment,XBottomSpecialSegment],1,NumberStory);

Zss{ijk,:} = [ZTopSpecialSegment,ZBottomSpecialSegment];

%% Plot Cooridinate Point (X,Z)

X{ijk,:} = [Xbase{ijk,:},Xcol{ijk,:},Xos{ijk,:},Xss{ijk,:}];

Z{ijk,:} = [Zbase{ijk,:},Zcol{ijk,:},Zos{ijk,:},Zss{ijk,:}];

scatter(X{ijk,:},Z{ijk,:},'filled'),grid on

%check angle

tangent = (Lss(j)/Nss(k))/h(i);

anglespecialsegment(ijk) = atand(tangent);

%% Name the Joint

jointbase{ijk} = table((1:numel(Xbase{ijk,:}))',Xbase{ijk,:}',Zbase{ijk,:}','VariableNames',{'ID','X','Z'});

jointcolumn{ijk} = table((max(jointbase{ijk}.ID)+1:numel(Xcol{ijk,:})+max(jointbase{ijk}.ID))',Xcol{ijk,:}',Zcol{ijk,:}','VariableNames',{'ID','X','Z'});

jointordinarysegment{ijk} = table((max(jointcolumn{ijk}.ID)+1:numel(Xos{ijk,:})+max(jointcolumn{ijk}.ID))',Xos{ijk,:}',Zos{ijk,:}','VariableNames',{'ID','X','Z'});

jointtopordinarysegment{ijk} = jointordinarysegment{ijk}(ismember(jointordinarysegment{ijk}.Z,Zcolumn),:);

jointbottomordinarysegment{ijk} = jointordinarysegment{ijk}(ismember(jointordinarysegment{ijk}.Z,Zbotbeam),:);

jointspecialsegment{ijk} = table((max(jointbottomordinarysegment{ijk}.ID)+1:numel(Xss{ijk,:})+max(jointbottomordinarysegment{ijk}.ID))',Xss{ijk,:}',Zss{ijk,:}','VariableNames',{'ID','X','Z'});

jointtopspecialsegment{ijk} = jointspecialsegment{ijk}(ismember(jointspecialsegment{ijk}.Z,Zcolumn),:);

jointbottomspecialsegment{ijk} = jointspecialsegment{ijk}(ismember(jointspecialsegment{ijk}.Z,Zbotbeam),:);

end

% Verify Variants Which Passed the Angle Requirement

tablecheck = table((1:numel(anglespecialsegment))',totalspan,anglespecialsegment','VariableNames',{'VariantID','Total Span','Angle'});

VariantIDPassed = tablecheck.VariantID(tablecheck.Angle>=30 & tablecheck.Angle<=60);

jointbase = table((1:numel(jointbase))',jointbase','VariableNames',{'VariantID','Joint Base'});

jointbase = jointbase(ismember(jointbase.VariantID,VariantIDPassed),:);

jointcolumn = table((1:numel(jointcolumn))',jointcolumn','VariableNames',{'VariantID','Joint Column'});

jointcolumn = jointcolumn(ismember(jointcolumn.VariantID,VariantIDPassed),:);

jointtopordinarysegment = table((1:numel(jointtopordinarysegment))',jointtopordinarysegment','VariableNames',{'VariantID','Joint Top Ordinary Segment'});

jointtopordinarysegment = jointtopordinarysegment(ismember(jointtopordinarysegment.VariantID,VariantIDPassed),:);

jointbottomordinarysegment = table((1:numel(jointbottomordinarysegment))',jointbottomordinarysegment','VariableNames',{'VariantID','Joint Bottom Ordinary Segment'});

jointbottomordinarysegment = jointbottomordinarysegment(ismember(jointbottomordinarysegment.VariantID,VariantIDPassed),:);

jointtopspecialsegment = table((1:numel(jointtopspecialsegment))',jointtopspecialsegment','VariableNames',{'VariantID','Joint Top Special Segment'});

jointtopspecialsegment = jointtopspecialsegment(ismember(jointtopspecialsegment.VariantID,VariantIDPassed),:);

jointbottomspecialsegment = table((1:numel(jointbottomspecialsegment))',jointbottomspecialsegment','VariableNames',{'VariantID','Joint Bottom Special Segment'});

jointbottomspecialsegment = jointbottomspecialsegment(ismember(jointbottomspecialsegment.VariantID,VariantIDPassed),:);

%% Name the Frame

% column

for i = 1:height(VariantIDPassed)

jointIcolumn = 0;

jointJcolumn = 0;

jointcolumnnew{i,:} = jointcolumn.("Joint Column"){i}(ismember(jointcolumn.("Joint Column"){i}.Z,Zcolumn),:);

jointcolumnnew{i,:} = vertcat(jointbase.("Joint Base"){i},jointcolumnnew{i});

for j = 1:height(jointcolumnnew{i,1})-2

jointIcolumn(j) = jointcolumnnew{i,1}.ID(j);

jointJcolumn(j) = jointcolumnnew{i,1}.ID(j+2);

end

framecolumn{i,:} = table((1:j)',jointIcolumn',jointJcolumn','VariableNames',{'Frame ID','JointI','JoinJ'});

end

% top ordinary segment (TOS)

for i = 1:height(VariantIDPassed)

endTOS{i,:} = jointcolumn.("Joint Column"){i}(ismember(jointcolumn.("Joint Column"){i}.Z,Zcolumn),:);

jointTOS{i,:} = vertcat(jointtopordinarysegment.("Joint Top Ordinary Segment"){i},endTOS{i});

jointTOSnew1{i,:} = sortrows(jointTOS{i}(ismember(jointTOS{i}.Z,5),:),2);

jointTOSnew2{i,:} = sortrows(jointTOS{i}(ismember(jointTOS{i}.Z,10),:),2);

jointTOSnew{i,:} = vertcat(jointTOSnew1{i},jointTOSnew2{i});

end

here the resut that i want :

jointTOSnew{1}
ans = 36x3 table
    ID      X       Z
    __    ______    _

       0    5
  0.8125    5
   1.625    5
  2.4375    5
    3.25    5
  4.0625    5
   4.875    5
  5.6875    5
     6.5    5
     8.5    5
  9.3125    5
  10.125    5
  10.938    5
   11.75    5
  12.562    5
  13.375    5