Formation of another matrix from where I can Plot Data

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

Answer 1

[Cris LaPierre]

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.

[Sumit Saha]

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]

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]

@Cris This command is available for Matlab 2017a?

files = dir("**/*.txt");

This shows

Error using dir

Function is not defined for 'string' inputs.

Error in Test_check (line 27)

files = dir("**/*.txt");

[Cris LaPierre]

dir has been around for a long time. MATLAB introduced strings (enclosed in double quotes) starting in R2016b. Each releases expands the functionality to more functions. However, the error suggests that, in 17a at least, dir had not been updated to support it. In that case, use a character array instead (single quotes).

files = dir('**/*.txt');

[Sumit Saha]

@Cris kindly look into this matter...

Undefined function or variable 'readmatrix'.

Error in Test_check (line 34)

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

[Cris LaPierre]

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.