How to make dynamic variable names (A1, A2, ..., An) with "for" loop? Using "eval" "num2str"?

Question

Italo 2013-11-13

4
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str

评论： João 2022-5-6

采纳的回答： Stephen23

在 MATLAB Online 中打开

Hello community,

my knowledge of Matlab is limited, I admit it. So, I am sorry if I am going to make mistakes.

I have to create a series of variable using a "for" loop.

I make this:

for i=1:3
    eval(['A' num2str(i) '= i'])
end

and it works well, it makes 3 variables A1, A2, A3.

But I need to use this variables to make other variables B1, B2, B3 where Bi=Ai*i. So I should have B1=A1*1=1, B2=A2*2=2*2=4, B3=A3*3=3*3=9

I tried something like this:

for i=1:3
    eval(['A' num2str(i) '= i'])
    eval(['B' num2str(i) '= 'A' num2str(i) '*i])
end

but it gives me error and it doesn't work. Of course I simplified things, the real code is a bt more complicated but I need to know the basics of how this could work.

Thank you

11 个评论
显示 9更早的评论隐藏 9更早的评论

Stephen23 2020-7-6

编辑：Stephen23 2020-7-7

"...unbalanced or unexpected parenthesis or bracket error has occurred in second eval. Why??"

Looking at the character vector shows at least one obvious syntax error (which looks like an implicit multiplication, but without any code comments, this is just a guess), which is quite possibly the cause of that error. There might be more syntax errors.

"Please explain me"

If this code had been written better** using basic indexing then the location of any syntax errors would get highlighted by the MATLAB Editor, and the code would be much easier to debug. So whoever wrote that code chose an approach which forces the user into writing slow, complex, buggy code which is hard to debug. Now you are finding that (otherwise very simple) code is buggy and hard to debug, which will come as no surprise to any readers of this forum (or for that matter to anyone who reads the MATLAB documentation).

** better in the sense simpler, neater, faster, more efficient, less buggy, easier to debug, e.g. by using basic indexing.

Stephen23 2022-1-19

编辑：Stephen23 2022-1-19

"With dynamic naming you can call them pic1.png, pic2.png, pic3.png, pic4.png, etc and have them automatically saved in the directory of your choosing."

No, you are mixing up filenames with the names of variables in the MATLAB workspace. Two totally different things.

Mixing up unrelated topics will not make it easy to understand or use MATLAB. Or any programming language, for that matter.

Also note that MATLAB variable names cannot have the dot character in them, so your example variable names are invalid and would not work, even though you incorrectly state that "you can call them pic1.png, pic2.png, pic3.png, pic4.png". This should also give you a hint as to one reason why your approach of dynamically naming variables based on a filenames is a really bad idea (hint: consider all of the Unicode characters that are valid in filenames. Are they also valid in variable names? hint: no).

"Using indexing, how can you call them pic1, pic2, pic3, and pic4 in order to automatically increment and save each photo?"

Very easily: the two main approaches are:

read a list of existing files and alter the names (e.g. append a suffix, change the extension).
generate the filenames (e.g. using SPRINTF).

Both of these methods generally use indexing, as shown in the MATLAB documentation here:

https://www.mathworks.com/help/matlab/import_export/process-a-sequence-of-files.html

and also on thousands and thousands and thousands of threads on this forum (none of which require dynamic variable names, even if their filenames increment numbers just as you show), because, surprisingly, importing and exporting file data is a very common topic that people want to achieve. Using indexing. Because that is how MATLAB works.

Steven Lord 2022-1-19

在 MATLAB Online 中打开

I think of dynamic file names as somewhat different from dynamic variable names, as the file names are usually themselves data (passed into a function that reads the file contents then returns a variable whose name is not necessarily tied to the file name, or performs some other operation on the file.)

t = fullfile(tempdir, 'writeASeriesOfFiles');
if ~isfolder(t)
    mkdir(t);
end
dir(t) % Nothing in this directory yet
.   ..  
for k = 1:5
    fn = fullfile(t, "file" + k + ".txt");
    % Now use fn as data in several function calls
    fclose(fopen(fn, 'wt')); % Open and immediately close the file
    doesItExist = isfile(fn);
    if doesItExist
        fprintf("File %s does exist!\n", fn)
    else
        fprintf("File %s does NOT exist!\n", fn)
    end
end
File /tmp/writeASeriesOfFiles/file1.txt does exist!
File /tmp/writeASeriesOfFiles/file2.txt does exist!
File /tmp/writeASeriesOfFiles/file3.txt does exist!
File /tmp/writeASeriesOfFiles/file4.txt does exist!
File /tmp/writeASeriesOfFiles/file5.txt does exist!
dir(t)
.          ..         file1.txt  file2.txt  file3.txt  file4.txt  file5.txt  

请先登录，再进行评论。

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

Answer 1

Stephen23 2016-1-12

8
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str#answer_205854

编辑：Stephen23 2019-6-19

https://www.mathworks.com/help/matlab/matlab_prog/string-evaluation.html

https://www.mathworks.com/matlabcentral/answers/304528-tutorial-why-variables-should-not-be-named-dynamically-eval

2 个评论
显示无隐藏无

Jan 2018-4-30

+1. Another time: Thanks Stephen.

Giovanni Barbarossa 2018-6-27

Hello, I have a similar need deriving from what it seems to be a limitation of the Classification Learner app. I am generating a cell array of tables with my program and want to try each of the tables as training features in the Classification Learner, but the app does not see cell arrays in the workspace. I do understand it's a bad idea to generate dynamic variable names, but what's the solution to my need? Thanks

请先登录，再进行评论。

Answer 2

Amit 2014-1-22

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str#answer_121865

编辑：Amit 2014-1-22

在 MATLAB Online 中打开

Italo, I am glad you're trying out things on matlab.

Matlab is very good with matrixes. What you're doing in these can be done much easily. For example:

A = 1:3;
B = A.*A; %

You can find many tutorials online on how to use MAtlab. For example: http://www.mathworks.com/academia/student_center/tutorials/launchpad.html

2 个评论
显示无隐藏无

Amit 2014-1-22

在 MATLAB Online 中打开

In you code:

for i=1:3
eval(['A' num2str(i) '= i'])
eval(['B' num2str(i) '= A' num2str(i) '*i']) % Issue was Here. Notice the positions of '
end

Bruno Pop-Stefanov 2014-1-22

I would like to add that eval is #1 in this list: Top 10 MATLAB code practices that make me cry

请先登录，再进行评论。

Answer 3

Gabor 2021-9-20

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str#answer_791504

在 MATLAB Online 中打开

T=table;
Date=datetime(2014,12,31);
eval(['Dynamic_var_name_' datestr(Date,'mm_dd_yyyy') '=T;']);

This is how you name dynamically a variable or a table or a cell or anything regardless if it is recommended or not.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Stephen23 2021-9-20

编辑：Stephen23 2021-9-20

Whether it is recommended or not is irrelevant.

The fact is that it forces beginners into writing slow, complex, inefficient code, as this thread already demonsrates:

https://www.mathworks.com/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str#comment_598847

请先登录，再进行评论。

Answer 4

Jerry Olup 2018-4-28

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str#answer_317694

编辑：Jerry Olup 2018-4-28

在 MATLAB Online 中打开

I did find dynamically named variables useful for a specific need. I have several files in a directory. Each file has a few columns that are read into a structure. It's convenient to read and simultaneously name the variables per a filename. For instance I read to a temp var then use the ill-advised exec capability to assign this to a variable specific to the file I read in. Looping, I read all files in the dir efficiently. 'files' is a cell array of file names discovered in that directory....

    tmp = table(dataArray{1:end-1}, 'VariableNames', {'max_mag','max_noise','max_cnt'});
    eval([files{i} '=tmp'])

I do see where this is not good for general variable naming, but I believe there's a niche usage in this case. Opinions? Best, Jerry

3 个评论
显示 1更早的评论隐藏 1更早的评论

Stephen23 2018-4-28

编辑：Stephen23 2018-5-2

在 MATLAB Online 中打开

" It's convenient to read and simultaneously name the variables per a filename."

A structure would be simpler, as well as be more efficient and easier to debug:

S.(files{i}) = table(...);

Or simply store the tables themselves in a cell array, and keep the corresponding files data in another cell array: this would be trivial to achieve using indexing, and also has the significant advantage that then the files data are not restricted to valid field/variable names.

Or just use one table, and add a files column to distinguish between the different cases.

"I do see where this is not good for general variable naming"

How do you then access those dynamically defined variables? Your "niche case" does not see to be any different to any of the other cases that have been described on this forum: can you explain why writing slow, complex code is better than using a simple structure, like I showed? Or better than using very simple (and very efficient) indexing?

"but I believe there's a niche usage in this case"

So far you have not actually described anywhere in your answer why it is necessary to dynamically create/access variable names: other methods would also store your data more efficiently than what you have shown, and you have not given any reason why none of those methods would work.

Jan 2018-4-30

编辑：Stephen23 2018-5-2

I cannot emphasize it enough how important Stephen's suggestions are. Creating the variables dynamically impedes the JIT acceleration massively, such that the complete function can runs ways slower, not just the line containing these variables. The JIT must be able to handle the case, that one of the "files{i}" is called "plot". Afterwards plot is not a command anymore, but a table object. To catch such problems, all function calls must be checked manually if they have been overwritten dynamically.

Using dynamic field names instead by S.(files{i})=table(...) is safe and efficient. In addition reading the code allows to identify these variables directly: If I find "S.function" in the source code, I know where it is coming from, while for a simple "function" I cannot be sure, if this is a typo or a dynamically created variable. The editor cannot distinguish this also and will assume an illegal definition of a function. Such strange effect can never happen with dynamic field names.

Summary: No, there is no niche for dynamically created variables. eval is a shot in your knee. Sometimes it does not hurt immediately, but the later it starts to bite you, the harder it will be to clean up the code.

This has been discussed exhaustively in the links, Stephen has posted in his answer above. Believe him and the bunch of corresponding threads.

Jerry Olup 2018-8-12

编辑：Stephen23 2018-8-12

Hi Stephen, I'll give the idea of provide parallel constructs where the file-specific variable references are in one cell while data is collected in structure. I've shifted in the last month or so to use python for preprocessing - so this problem has largely gone away and been replaces with xlsread, csvread etc.

>>files data are not restricted to valid field/variable names

My file names are conventionally, ascii with underscores and timestamps yyyymmdd_hhmmss type, so I do see where this is problem. They are pretty specific so I thought I was giving a concrete counter example in contrast to @Italo and @Muna, who started this thread naming generic matrices.

>>you have not given any reason why none of those methods would work.

You're absolutely correct. I latched onto that solution as it was elegant for what I needed.

Thanks for the interaction, best, Jerry

请先登录，再进行评论。

Answer 5

César Silva Proaño 2018-8-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str#answer_332454

在 MATLAB Online 中打开

Hi, look I was doing something similar in order to create a function with the Krammer Method, for solving any size of matrix-equation problem. Look I created a function in order to do this, therefore, I needed to do it automatically, creating variables dinamically. Look at this fucntion.

function [ X ] = funcionMetodoKrammer( A )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
% ex: A = [-3 2 1 10;5 -8 3 15;-8 -9 10 20];
sizeA = size(A);
sizeFilasA = sizeA(1,1);
sizeColumnasA = sizeA(1,2);
matrizASinResultado = A(1:sizeFilasA,1:sizeFilasA);
matrizResultados = A(:,sizeColumnasA);
if (sizeFilasA == sizeColumnasA - 1)
    for (i = 1:sizeColumnasA)
        genvarname('matrizColumna', num2str(i));
        eval(['matrizColumna' num2str(i) '=A(:,i)']);
    end
      for (i = 1:sizeFilasA)
          genvarname('matrizTemporalColumna', num2str(i));
          eval(['matrizTemporalColumna' num2str(i) '=matrizASinResultado(:,i)']);
      end
      matrizIncognita1 = [matrizResultados,matrizASinResultado(:,2:sizeFilasA)];
      for (i = 2:sizeFilasA)
          genvarname('matrizIncognita', num2str(i));
          eval(['matrizIncognita' num2str(i) '=[matrizASinResultado(:,1:i- 
 1),matrizResultados,matrizASinResultado(:,i+1:sizeFilasA)]']);
      end
else
    msgbox('La matriz ingresada para el método Krammer es incorrecta');
end
for (i = 1:sizeFilasA)
    genvarname('X', num2str(i));
    genvarname('matrizIncognita', num2str(i));
    eval(['X' num2str(i) '= det(matrizIncognita' num2str(i) ')/det(matrizASinResultado)']);
end
X = zeros(1,sizeFilasA);
for (i = 1:sizeFilasA)
    genvarname('X', num2str(i));
    eval(['X' '(1,' num2str(i) ')=X' num2str(i)]);
end
end
%{
Asuminedo que existen variables B1, B2 y B3 se pueden acceder a las
variables de la siguiente manera
B1 = 10; B2 = 11; B3 = 12;
for (i = 1:3)
    genvarname('Prueba', num2str(i));
    genvarname('B', num2str(i));
    eval(['Prueba' num2str(i) '=' 'B' num2str(i) '+10']);
end
%}

I strongly disagree with the idea that you shouldn't use this because, if it exists I mean why not? It is not prohibited, and additionally it is very helpful in some cases.

Regards,

Sebastián

5 个评论
显示 3更早的评论隐藏 3更早的评论

Stephen23 2018-8-12

编辑：Stephen23 2020-7-6

在 MATLAB Online 中打开

@César Silva Proaño: thank you for showing us your function. It is a good example of what we warn against: eval causing pointlessly slow, pointlessly complex, obfuscated, inefficient code. One bad design decision that ruins what should be a simple function.

"I strongly disagree with the idea that you shouldn't use this because, if it exists I mean why not? It is not prohibited, and additionally it is very helpful in some cases."

You can "strongly disagree" all you want, but in this case it "helped" you to write pointlessly slow, pointlessly complex, pointlessly obfuscated code. Which is exactly "why not" !

"Look I created a function in order to do this, therefore, I needed to do it automatically, creating variables dinamically."

There is absolutely nothing in your algorithm that you "needed to" use dynamically named variables for. I proved this with my own simpler version of your function (see betterKrammer below), which uses basic indexing to give a much faster, much simpler, less buggy, and easier-to-debug function. Actually you used exactly the same indexing, you just obfuscated it behind lots of dynamically named variables.

OBFUSCATION: using eval made it easy to hide pointless "features". For example:

eval(['X' '(1,' num2str(i) ') = ...

Why do you need eval to simply use basic indexing? Answer: you don't. If you had read the eval documentation you would know that in that case it is recommended to put the eval only on the RHS of the assignment.

None of the genvarname calls do anything useful (just waste time), because their outputs are not assigned to any variable and so are immediately discarded:

genvarname('X', num2str(i)); <- discarded output.
genvarname('matrizTemporalColumna', num2str(i)); <- discarded output.
... etc

So they just slow down your code and produce nothing.

What do you use the ['matrizColumna'...] and ['matrizTemporalColumna'...] variables for? Answer: also nothing. So why do you bother to generate them? If you had written better code (without eval) then the MATLAB editor would simply warn you that all these variables are unused. But by choosing to use eval you turned off all of these warnings, and so have lots of pointless calculations in your code: literally pointless because they are not used for anything. That is the "help" that eval gave you: it turned off those useful warnings!

Your function also inefficiently and pointlessly duplicates parts of the input array into their own variables: obviously all of that data duplication escaped your attention, but really it is not necessary to copy the same data into lots of (dynamically named) variables, when you can access it very efficiently using indexing, and without any data duplication whatsoever.

SLOW and COMPLEX: I rewrote your function without any dynamically named variables, by using basic MATLAB indexing. This is what I came up with:

function X = betterKrammer(A)
% ex: A = [-3 2 1 10;5 -8 3 15;-8 -9 10 20];
S = size(A);
assert(numel(S)==2&&diff(S)==1,'La matriz ingresada para el método Krammer es incorrecta')
X = nan(1,S(1)); % preallocate output.
for k = 1:S(1)
    X(k) = det(A(:,[1:k-1,S(2),k+1:S(1)])) ./ det(A(:,1:S(1)));
end
end

Notice that your function has sixty lines, whereas my function has just nine lines (including the same comments and better input checking). Of course it returns exactly the same answer, it just does so much more efficiently. After I stopped your function from printing to the command window I did some timing tests using tic and toc, here are the results (1e4 iterations):

Elapsed time is 109.192323 seconds.  % your slow function using EVAL.
Elapsed time is 0.842797 seconds.    % my fast function sans EVAL.

And we can compare the outputs:

>> funcionMetodoKrammer([-3 2 1 10;5 -8 3 15;-8 -9 10 20]) % your slow function
ans =
       6.9898       8.0612       14.847
>> betterKrammer([-3 2 1 10;5 -8 3 15;-8 -9 10 20]) % my fast function
ans =
       6.9898       8.0612       14.847

So by using basic indexing and avoiding eval I wrote a simpler function that is over one hundred times faster than yours. This is what happens when beginners use eval: they force themselves into writing inefficient, complex, obfuscated code. Which, no matter how much you "strongly disagree" with it, is the reality that we warn against.

So there you have it: your code is a very good example of exactly what we are warning about: pointless eval usage forcing beginners into writing really slow, complex code. Thank you for the example!

César Silva Proaño 2018-8-12

在 MATLAB Online 中打开

Thanks for the feedback. But I can see that in every post, in which this question about the generation of dinamic variables with eval and genvarname is lauched, you just attack this method and you do not present an formal solution. Thanks for the critics, yes in fact, that code is not the final one. This is the one:

function [X] = funcionMetodoKrammer(A)
%{
Esta función resuelve un sistema se ecuaciones lineales con el método 
de
Krammer, se debe ingresar una matriz de tamaño n x n+1, es decir la 
parte
n x n contiene los valores de x1, x2, etc, y la columna n+1 tiene los
valores resultado de cada ecuación.
%}
% ex: A = [-3 2 1 10;5 -8 3 15;-8 -9 10 20];
sizeA = size(A);
sizeFilasA = sizeA(1,1);
sizeColumnasA = sizeA(1,2);
matrizASinResultado = A(1:sizeFilasA,1:sizeFilasA);
matrizResultados = A(:,sizeColumnasA);
if (sizeFilasA == sizeColumnasA - 1)
    % Aquí se forman las matrices para cada respuesta:
    matrizIncognita1 = 
[matrizResultados,matrizASinResultado(:,2:sizeFilasA)];
      for (i = 2:sizeFilasA)
          genvarname('matrizIncognita', num2str(i));
          eval(['matrizIncognita' num2str(i) '=[matrizASinResultado(:,1:i-1),matrizResultados,matrizASinResultado(:,i+1:sizeFilasA)]']);
      end
      % Se determinan las respuestas
      for (i = 1:sizeFilasA)
      genvarname('X', num2str(i));
      genvarname('matrizIncognita', num2str(i));
      eval(['X' num2str(i) '= det(matrizIncognita' num2str(i) ')/det(matrizASinResultado)']);
      end
      % Ahora se determina el vector con las respuestas:
      X = zeros(1,sizeFilasA);
      for (i = 1:sizeFilasA)
          genvarname('X', num2str(i));
          eval(['X' '(1,' num2str(i) ')=X' num2str(i)]);
      end
% En caso de que la matriz ingresada sea incorrecta:
else
    msgbox('La matriz ingresada para el método Krammer es incorrecta');
end
end

It seems that you have a personal issue with this form of approaching the dinamic generation of variables. Perhaps you are an excellent Matlab user, so why don't you give us an enlightening solution for this problem, instead of criticizing?

The question should be, how do you generate dinamic variables?

But please make us all a favor and present a code, not just words or links, because as my previous code, that is useless for people who want to learn matlab in a good manner. And yes, I am a beginner in this program, that is why perhaps I can commit mistakes, but I beg my pardon, everyone can make mistakes, specially me and others who are not so smart as you ;) Regards,

Walter Roberson 2018-8-13

César Silva Proaño:

Imagine that you were building a bird house by using a beer bottle to hammer in nails. You get something made and present it to us as a model for how to make bird houses.

Imagine that we notice that you used the wrong size of nails for the materials, that you used the wrong kind of nails for the materials, that using a beer bottle as a hammer is dangerous, and that our considered experience is that it is a lot better to build bird houses with screws instead of with nails.

In those circumstances, should we advise you to use screws and a screwdriver instead of nails, or should we refrain from saying anything, seeing as you created a birdhouse that looks good enough to you? Should we just shrug our shoulders about the dangers of bottles shattering in hands when used as hammers, figuring that people will eventually figure it out for themselves?

jonas 2018-9-11

编辑：Stephen23 2021-9-20

"you just attack this method and you do not present an formal solution"

This is simply untrue. On the contrary, from what I've seen Stephen always provides a concrete solution when possible. Here's an example from 5 minutes ago ( link ) which is nothing out of the ordinary.

"But please make us all a favor and present a code, not just words or links..."

This is the mindset I see in a lot of younger students today - avoid instructions like the plague and feel entitled to a solution. Remember that no one is getting paid to read about your dynamic variables, so you should appreciate any replies including links and words (whatever that means).

请先登录，再进行评论。

Answer 6

João 2022-5-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105936-how-to-make-dynamic-variable-names-a1-a2-an-with-for-loop-using-eval-num2str#answer_957280

编辑：João 2022-5-4

I too have a potential reason to want to do what the original OP posted. My example is that I want to numerically solve a system of nonlinear functions, in which the number of functions and variables depend on the input. I don't think any of the suggested alternatives solve this case. Am I wrong?

14 个评论
显示 12更早的评论隐藏 12更早的评论

João 2022-5-5

在 MATLAB Online 中打开

Dear @Stephen Just not to ask more of your time I'll go straight to the point: I'm a python user, but the case I'm trying to model involves solving a system of mix linear and nonlinear equations that change as the program progresses. To detail, the equations consist in equations of motion of objects during their impact. Over time some objects which were at rest become affected by the pressure wave, so no longer in equilibrium, so the equations of motion will change over time as the net force also changes. The knowns are which objects are in motion or at rest. I have more 2 variables than equations of motion, but I can use the law of conservation of momentum and the law of conservation of total energy (plus some assumptions) to get a n equations and n variables system. The equations of motion are linear (depend on displacement and velocity), as the aw of conservation of momentum. The law of conservation of total energy is nonlinear since it depends on the square of the velocity. I am using very small time increments so that the velocity and displacement can be assumed constant and modelled not as functions of time but as time independent variables at each increment. I want to use a loop to define the equations for each set of objects at rest and in motion, and skip having to write them one by one. Also the number of objects is an initial input from the user. So inside each loop, I will define a set of linear an nonlinear equations, in which the number of equations and the name of variables depend on the range of the loop. So for a 2 objects system and a specific time increment I might have (in Python):

    def equations(vars):
        g = 9.8
        c = 7.
        d = 0
        k = 1000
        m = 100
        dt = 0.1
        x, y, z, j = vars # variables
        eq1 = (m*g - k*(y-x)) * dt - m * (j - c) # where m, k, g and c are known values
        eq2 = (k*(y-x) + m*g - k*x) * dt - m * (z - d) # where m, k, g and d are known values
        eq3 = m * (c + d) - m * (z + j)
        eq4 = m*g*y - 1/2*k*(y-x)**2 + m*g*x - 1/2*k*x**2 - 1/2*m*(j**2 - c**2) - 1/2*m*(z**2 - d**2)
        return [eq1, eq2, eq3, eq4]
    from scipy.optimize import fsolve
    sol =  fsolve(equations, [0.0, 0.0, 0.0, 0.0])

The issue is that for another time increment I might have slightly different equations, and also I want to make the number of objects a user defined variable. So having to define functions for all increments and all possibilities is not feasible.

Stephen23 2022-5-5

编辑：Stephen23 2022-5-6

"So having to define functions for all increments and all possibilities is not feasible."

Sure.

"I want to use a loop to define the equations for each set of objects at rest and in motion, and skip having to write them one by one"

So presumably there is some kind of pattern to these things: those equations must be abe to be defined systematically depending on the number of objects/whatever (otherwise your task is moot), and that is really what you need to utilize if you want a general solution to this (regardless of Python or MATLAB).

Please show one (or all) of the equations for three objects, to give me an idea of how they change.

"Over time some objects which were at rest become affected by the pressure wave, so no longer in equilibrium, so the equations of motion will change over time as the net force also changes"

That sounds like something that could be handled using indexing. Or at worst some loops and IFs.

The more you explain, the more your approach that "the name of variables depend on the range of the loop" seems ... unhelpful in either language. Numpy and MATLAB both work with arrays, so should you.

João 2022-5-6

在 MATLAB Online 中打开

@Stephen23 Yes, equations must be abe to be defined systematically depending on the number of objects. And yes, for each time increment I know which equations to use and which x, y, z, ... to define.

For two objects I would have something like:

x, y, z, j = vars # variables
eq1 = (m*g - k*(y - x)) * dt - m * (j - c) # where m, k, g and c are known values
eq2 = (k*(y - x) + m*g - k*x * dt - m * (z - d) # where m, k, g and d are known values
eq3 = m * (c + d) - m * (z + j)
eq4 = m*g*y - 1/2*k*(y-x)**2 + m*g*x - 1/2*k*(x-l)**2 + m*g*l - 1/2*k*(l)**2 
- 1/2*m*(j**2 - c**2) - 1/2*m*(z**2 - d**2)

For three objects I would have something like:

x, y, z, j, l, f = vars # variables
eq1 = (m*g - k*(y - x)) * dt - m * (j - c) # where m, k, g and c are known values
eq2 = (k*(y - x) + m*g - k*(x - l) * dt - m * (z - d) # where m, k, g and d are known values
eq3 = (k*(x - l) + m*g - k*l) * dt - m * (f - e) # where m, k, g and e are known values
eq4 = m * (c + d + e) - m * (z + j + f)
eq5 = m*g*y - 1/2*k*(y-x)**2 + m*g*x - 1/2*k*(x-l)**2 + m*g*l - 1/2*k*(l)**2 
- 1/2*m*(j**2 - c**2) - 1/2*m*(z**2 - d**2) - 1/2*m*(f**2 - e**2)
eq6 = j - c # assuming that the velocities of object 1 and object 2 are the same during impact, i.e. when object 1 hits objects 2 they have the same velocity

So this made me realise I have to keep making assumptions as the number of objects grows.

So I need to solve the equations explicitly by Runge-Kutta or Central Difference methods...

I need to restructure my code. Thank you for the constructive input!

请先登录，再进行评论。

How to make dynamic variable names (A1, A2, ..., An) with "for" loop? Using "eval" "num2str"?

11 个评论
显示 9更早的评论隐藏 9更早的评论

采纳的回答

2 个评论
显示无隐藏无

更多回答（5 个）

2 个评论
显示无隐藏无

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论

14 个评论
显示 12更早的评论隐藏 12更早的评论

另请参阅

类别

标签

Community Treasure Hunt

How to make dynamic variable names (A1, A2, ..., An) with "for" loop? Using "eval" "num2str"?

11 个评论 显示 9更早的评论隐藏 9更早的评论

采纳的回答

2 个评论 显示 无隐藏 无

更多回答（5 个）

2 个评论 显示 无隐藏 无

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

3 个评论 显示 1更早的评论隐藏 1更早的评论

5 个评论 显示 3更早的评论隐藏 3更早的评论

14 个评论 显示 12更早的评论隐藏 12更早的评论

另请参阅

类别

标签

Community Treasure Hunt

11 个评论
显示 9更早的评论隐藏 9更早的评论

2 个评论
显示无隐藏无

2 个评论
显示无隐藏无

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论

14 个评论
显示 12更早的评论隐藏 12更早的评论