Formation of another matrix from where I can Plot Data

Question

Sumit Saha 2021-1-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/712823-formation-of-another-matrix-from-where-i-can-plot-data

评论： Cris LaPierre 2021-1-10

IDAOutput.zip

How can I open txt files folder wise and read 2nd column of txt files and store it in another matrix from where I can Plot Data? I've also attached the txt files for reference. I didn't get any confined answer

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Cris LaPierre 2021-1-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/712823-formation-of-another-matrix-from-where-i-can-plot-data#answer_594413

I'd start with learning how to import a single file.

This page will help you learn how to index your data. I'd also suggest completng Ch 5 of MATLAB Onramp. Ch 9 covers plotting as well. Ch 11 talks about importing data.

Perhaps this answer is helpful in seeing how to import multiple files from multiple folders.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Sumit Saha 2021-1-9

在 MATLAB Online 中打开

I've tried in this way

clc
clear all;
close all;
for kk =1:60
for k=1:50
 % here drift values are in the form of drift ratio
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Roof.txt']);
 Roof_Drift_max(kk,k,:)= max(abs(Roof(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story1.txt']);
 Story1_Drift_max(kk,k,:)= max(abs(Story1(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story2.txt']);
 Story2_Drift_max(kk,k,:)= max(abs(Story2(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story3.txt']);
 Story3_Drift_max(kk,k,:)= max(abs(Story3(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story4.txt']);
 Story4_Drift_max(kk,k,:)= max(abs(Story4(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story5.txt']);
 Story5_Drift_max(kk,k,:)= max(abs(Story5(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story6.txt']);
 Story6_Drift_max(kk,k,:)= max(abs(Story6(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story7.txt']);
 Story7_Drift_max(kk,k,:)= max(abs(Story7(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story8.txt']);
 Story8_Drift_max(kk,k,:)= max(abs(Story8(:,2)));
 eval(['load IDAOutput/EQ_',num2str(kk),'/Scale_',num2str(k),'/StoryDrifts/','Story9.txt']);
 Story9_Drift_max(kk,k,:)= max(abs(Story9(:,2)));
% Maximum Interstorey Drift Ratio
 MIDR(kk,k,:) =max([Story1_Drift_max(kk,k,:),Story2_Drift_max(kk,k,:),Story3_Drift_max(kk,k,:),Story4_Drift_max(kk,k,:),...
                   Story5_Drift_max(kk,k,:),Story6_Drift_max(kk,k,:),Story7_Drift_max(kk,k,:),Story8_Drift_max(kk,k,:),Story9_Drift_max(kk,k,:)]);
end
end
 PGA_Scale= [0.1:0.1:5];
 figure(1)
 for tt = 1:60
  plot(Roof_Drift_max(tt,:),PGA_Scale,'-*r','linewidth',0.5);hold on;
 end
 xlabel('Maximum Storey Drift Ratio (MSDR)')
 ylabel('IM (PGA)')     
 title('Incremental Dynamic Analysis (IDAs)for GMSetAll')
 set(gcf,'position',[ 146 378 858 779])
 set(findall(gcf,'-property','Fontsize'),'Fontsize',12)
 print('DT100','-dpng','-r300');
 figure(2)
 subplot(2,2,1)
 plot(Roof_Drift_max(1,:),PGA_Scale,'-*b',Roof_Drift_max(5,:),PGA_Scale,'-*r',Roof_Drift_max(6,:),PGA_Scale,'-*k',Roof_Drift_max(13,:),PGA_Scale,'-*g',...
     Roof_Drift_max(14,:),PGA_Scale,'-*y',Roof_Drift_max(21,:),PGA_Scale,'-*c',Roof_Drift_max(22,:),PGA_Scale,'-*m',Roof_Drift_max(29,:),PGA_Scale,'--*r',...
     Roof_Drift_max(30,:),PGA_Scale,'k',Roof_Drift_max(37,:),PGA_Scale,'-*b',Roof_Drift_max(38,:),PGA_Scale,'-*b',Roof_Drift_max(45,:),PGA_Scale,'-*b',...
     Roof_Drift_max(46,:),PGA_Scale,'-*b',Roof_Drift_max(53,:),PGA_Scale,'-*b',Roof_Drift_max(54,:),PGA_Scale,'-*b'); hold on
 xlabel('Maximum Storey Drift Ratio (MSDR)')
 ylabel('IM (PGA)')     
 title('GMSet1')
  set(gcf,'position',[ 146 378 858 779])
  set(findall(gcf,'-property','Fontsize'),'Fontsize',12)
  print('DT200','-dpng','-r300');
 subplot(2,2,2)
 plot(Roof_Drift_max(2,:),PGA_Scale,'-*b',Roof_Drift_max(7,:),PGA_Scale,'-*r',Roof_Drift_max(8,:),PGA_Scale,'-*k',Roof_Drift_max(15,:),PGA_Scale,'-*g',...
     Roof_Drift_max(16,:),PGA_Scale,'-*y',Roof_Drift_max(23,:),PGA_Scale,'-*c',Roof_Drift_max(24,:),PGA_Scale,'-*m',Roof_Drift_max(31,:),PGA_Scale,'--*r',...
     Roof_Drift_max(32,:),PGA_Scale,'r',Roof_Drift_max(39,:),PGA_Scale,'r',Roof_Drift_max(40,:),PGA_Scale,'r',Roof_Drift_max(47,:),PGA_Scale,'-*b',...
     Roof_Drift_max(48,:),PGA_Scale,'-*b',Roof_Drift_max(55,:),PGA_Scale,'-*b',Roof_Drift_max(56,:),PGA_Scale,'-*b'); hold on
 xlabel('Maximum Storey Drift Ratio (MSDR)')
 ylabel('IM (PGA)')     
 title('GMSet2')
 set(gcf,'position',[ 146 378 858 779])
 set(findall(gcf,'-property','Fontsize'),'Fontsize',12)
 print('DT200','-dpng','-r300');
 subplot(2,2,3)
 plot(Roof_Drift_max(3,:),PGA_Scale,'-*b',Roof_Drift_max(9,:),PGA_Scale,'-*r',Roof_Drift_max(10,:),PGA_Scale,'-*k',Roof_Drift_max(17,:),PGA_Scale,'-*g',...
     Roof_Drift_max(18,:),PGA_Scale,'-*y',Roof_Drift_max(25,:),PGA_Scale,'-*c',Roof_Drift_max(26,:),PGA_Scale,'-*m',Roof_Drift_max(33,:),PGA_Scale,'--*r',...
     Roof_Drift_max(34,:),PGA_Scale,'k',Roof_Drift_max(41,:),PGA_Scale,'r',Roof_Drift_max(42,:),PGA_Scale,'b',Roof_Drift_max(49,:),PGA_Scale,'-*b',...
     Roof_Drift_max(50,:),PGA_Scale,'-*b',Roof_Drift_max(57,:),PGA_Scale,'-*b',Roof_Drift_max(58,:),PGA_Scale,'-*b'); hold on
 xlabel('Maximum Storey Drift Ratio (MSDR)')
 ylabel('IM (PGA)')     
 title('GMSet3')
 set(gcf,'position',[ 146 378 858 779])
 set(findall(gcf,'-property','Fontsize'),'Fontsize',12)
 print('DT200','-dpng','-r300');
 subplot(2,2,4)
 plot(Roof_Drift_max(4,:),PGA_Scale,'-*b',Roof_Drift_max(11,:),PGA_Scale,'-*r',Roof_Drift_max(12,:),PGA_Scale,'-*k',Roof_Drift_max(19,:),PGA_Scale,'-*g',...
     Roof_Drift_max(20,:),PGA_Scale,'-*y',Roof_Drift_max(27,:),PGA_Scale,'-*c',Roof_Drift_max(28,:),PGA_Scale,'-*m',Roof_Drift_max(35,:),PGA_Scale,'--*r',...
     Roof_Drift_max(36,:),PGA_Scale,'k',Roof_Drift_max(43,:),PGA_Scale,'r',Roof_Drift_max(44,:),PGA_Scale,'k',Roof_Drift_max(51,:),PGA_Scale,'-*b',...
     Roof_Drift_max(52,:),PGA_Scale,'-*b',Roof_Drift_max(59,:),PGA_Scale,'-*b',Roof_Drift_max(60,:),PGA_Scale,'-*b'); hold on
 xlabel('Maximum Storey Drift Ratio (MSDR)')
 ylabel('IM (PGA)')     
 title('GMSet4')
 set(gcf,'position',[ 146 378 858 779])
 set(findall(gcf,'-property','Fontsize'),'Fontsize',12)
 print('DT200','-dpng','-r300');
 figure(21)
 for tt = 1:60
  plot(MIDR(tt,:),PGA_Scale,'-*b','linewidth',0.5);hold on;
 end
 xlabel('Maximum Inter-Storey Drift Ratio(MIDR)')
 ylabel('IM (PGA)')     
 title('Incremental Dynamic Analysis (IDAs)for GMSetAll')
 set(gcf,'position',[ 146 378 858 779])
 set(findall(gcf,'-property','Fontsize'),'Fontsize',12)
 print('DT2100','-dpng','-r300');
 

Cris LaPierre 2021-1-9

在 MATLAB Online 中打开

I would load the files without using eval statements. See this answer for a way to get all text files in a folder and all its subfolders.

https://www.mathworks.com/matlabcentral/answers/708923-how-can-i-add-an-extension-to-all-files-in-a-folder

Adapted to your data, this code loads all the files. I'll leave it to you to getting the data organized how you want.

Roof_Drift_max = [];
Story1_Drift_max= [];
Story2_Drift_max= [];
Story3_Drift_max= [];
Story4_Drift_max= [];
Story5_Drift_max= [];
Story6_Drift_max= [];
Story7_Drift_max= [];
Story8_Drift_max= [];
Story9_Drift_max= [];
% Find all text files in the current folder and all its subfolders
files = dir("**/*.txt");
% Use a loop to load all the files
for k = 1:numel(files)
    f = fullfile(files(k).folder,files(k).name);
    [p,nm,ext] = fileparts(f);
    
    d = readmatrix(f,'Range','B1');
    
    switch lower(nm)
        case "roof"
            Roof_Drift_max= [Roof_Drift_max; max(abs(d))];
        case "story1"
            Story1_Drift_max= [Story1_Drift_max; max(abs(d))];
        case "story2"
            Story2_Drift_max= [Story2_Drift_max; max(abs(d))];
        case "story3"
            Story3_Drift_max= [Story3_Drift_max; max(abs(d))];
        case "story4"
            Story4_Drift_max= [Story4_Drift_max; max(abs(d))];
        case "story5"
            Story5_Drift_max= [Story5_Drift_max; max(abs(d))];
        case "story6"
            Story6_Drift_max= [Story6_Drift_max; max(abs(d))];
        case "story7"
            Story7_Drift_max= [Story7_Drift_max; max(abs(d))];
        case "story8"
            Story8_Drift_max= [Story8_Drift_max; max(abs(d))];
        case "story9"
            Story9_Drift_max= [Story9_Drift_max; max(abs(d))];
            
    end
end

Sumit Saha 2021-1-10

@Cris kindly look into this matter...

Undefined function or variable 'readmatrix'.

Error in Test_check (line 34)

d = readmatrix(f,'Range','B1');

Cris LaPierre 2021-1-10

This is the function for loading the data. It was introduced in R2019a, so is not available to you in R2017a. However, just replace this with whatever you normally use to load your files.

请先登录，再进行评论。