How can I resolve Error Using max function
显示 更早的评论
clc
clear all;
close all;
for kk =1:44
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(Roof_Drift_max(kk,k,:),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
采纳的回答
Star Strider
2021-1-8
I am not certain what you want.
One option is to enclose all the arguments within square brackets [], effectively concatenating them:
MIDR(kk,k,:) =max([Roof_Drift_max(kk,k,:),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,:)]);
This will produce an ‘MIDR’ matrix with the same dimensions as the argument matrices.
Experiment with other approaches to get different results.
16 个评论
Sumit Saha
2021-1-8
Star Strider
2021-1-8
As always, my pleasure!
Sumit Saha
2021-1-9
Star Strider
2021-1-9
First, this is not the same as the max function, and
Second, using eval is definitely not recommended.
If you want to load .mat files (or certain text files), use the function implementation, not the command implementation.
Also, always load into a variable, creating a structure in that variable from which it is then relatively easily get the desired components. See Load List of Variables into Structure Array for details.
Sumit Saha
2021-1-9
"now from this situation how can i build my code wihtout using eval?"
The basic approaches are to either generate the file/folder names or to read them from the OS using dir. Although you can use eval if you really want to force yourself into writing slow, complex, obfuscated, buggy code that is hard to debug, you really should just use simpler and more efficient approaches, e.g.:
As Star Strider stated, your code would be a lot simpler if you used function syntax rather than command syntax. Your totally inappropriate usage of command syntax is how you forced yourself into unnecessary eval usage.
Rather than using num2str and string concatenation, you should use sprintf and fullfile.
Sumit Saha
2021-1-9
Star Strider
2021-1-9
Note that if you actually execute this example line (without the eval call):
['load IDAOutput\EQ_',num2str(p),'\Scale_',num2str(q),'\NodeDisplacements\','NodeDisplacementLevel1.txt']
the result is:
'load IDAOutput\EQ_1\Scale_2\NodeDisplacements\NodeDisplacementLevel1.txt'
since the square brackets [] are concatenation operators.
See if this does what you want:
p = 1;
q = 2;
r = 3;
filepath = fullfile(sprintf('IDAOutput\\EQ_%d\\Scale_%d\\',p,q), 'NodeDisplacements\', sprintf('NodeDisplacementLevel%d.txt',r))
LD = load(filename)
with ‘filepath’ (or whatever you choose to name it) now being:
filepath =
'IDAOutput\EQ_1\Scale_2\NodeDisplacements\NodeDisplacementLevel3.txt'
and all the file contents appearing in the ‘LD’ structure for you to extract.
If you are not familiar with structure arrays, see Access Data in Structure Array for details on how to work with them
.
Stephen23
2021-1-9
fn = sprintf('NodeDisplacementLevel%d.txt',r);
f1 = sprintf('EQ_%d',kk);
f2 = sprintf('Scale_%d',k);
ff = fullfile('IDAOutput',f1,f2,'NodeDisplacements',fn);
LD = load(ff)
Sumit Saha
2021-1-9
Sumit Saha
2021-1-9
Sumit Saha
2021-1-9
Sumit Saha
2021-1-9
编辑：Sumit Saha
2021-1-9
Star Strider
2021-1-9
Sumit Saha —
‘Can I select only second column of the NodeDisplacementLevel3.txt during execution of loop?’
I have no idea what ‘NodeDisplacementLevel3.txt’ is, or what it contains. If it only contains one variable (presumably a matrix), then this could work:
LD = load(filename);
SecondColumn = LD.SomeVariableMatrix(:,2);
If the variable is called ‘NodeDisplacementLevel3’, then something like this could work:
LD = load(filename);
SecondColumn = LD.NodeDisplacementLevel3(:,2);
Obviously, name the assigned variables something that makes sense in the context of your code.
That is as close as I can come to a reply, given the information available.
Sumit Saha
2021-1-9
Star Strider
2021-1-9
Try this to extract the variables:
uz = unzip('IDAOutput.zip');
k2 = 0;
for k1 = 1:numel(uz)
if ~isdir(uz(k1))
[~,txtname] = fileparts(uz{k1});
k2 = k2+1;
LD{k2,1} = txtname;
LD{k2,2} = load(uz{k1});
end
end
I have to unzip it, so I included tahat as part of my code. Change my code to work with the .zip file or it was created from. My code extracts the contents of the .zip file to a series of cell arrays, with the first element of the cell array being the name of the file it came from, and the second cell array the contents of that file as a double matrix.
So to get the information for the first 5 elements of ‘LD(1)’:
BldgPart_1 = [LD{1,1}]
BldgPrtData_1 = [LD{1,2}(1:5,:)]
produces:
BldgPart_1 =
'Roof'
BldgPrtData_1 =
0.002 6.9378e-08
0.004 6.9378e-08
0.006 6.9378e-08
0.008 6.9378e-08
0.01 6.9378e-08
and for ‘LD(40)’:
BldgPart_40 = [LD{40,1}]
BldgPrtData_40 = [LD{40,2}(1:5,:)]
produces:
BldgPart_40 =
'Story9'
BldgPrtData_40 =
0.002 5.5276e-07
0.004 5.5276e-07
0.006 5.5276e-07
0.008 5.5276e-07
0.01 5.5276e-07
I leave the rest to you.
更多回答（1 个）
Feng Shui Déco
2021-1-9
0 个投票
Yeah, you have to close your arguments with square bracket.
类别
在 帮助中心 和 File Exchange 中查找有关 Get Started with MATLAB 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
