Appending a calculated value during each iteration of a for loop

33 visualizzazioni (ultimi 30 giorni)
Hi! I am trying to use a for loop in a function code to append a calculated value to a vector for each iteration of the loop. For some context, the code I am writing is for calculating eigen values and eigen vectors for a system of ODEs. I am trying to calculate the coefficent values (called eigen_factors in my code) for each value of time, t, and then store all of those values in a vector, using the for loop. Whenever I call my function in a separate scricpt, it says that the variable eigen_factors is not defined in the function. I'm not sure why my function is not running the for loop, so please let me know if you have any suggestions!
Before the for loop, I have my attempt at preallocation. As of now, I have it with a % in front because when I kept in in the fucntion, the output value for eigen_factors would be that preallocation value. I'm not sure if I did that correctly, but my main problem here is just getting the for loop to run.
I had also tried calculating the coefficent values separtely within the for loop and then storing it in the eigen_factors vector, but I was not able to get that to work either.
Here is my code for the function of interest.
function [x, u, lambda, v, eigen_factors] = eigen_diffuse(sigmoid_ic, tspan, n)
k=5;
x=[0:1/(n+1):1];
x=[x(2:end-1)]';
u0=sigmoid_ic(x,k);
A=createA(n);
dx=1/(n+1);
[v,e_values]=eig(A*(1/dx^2));
lambda=diag(e_values);
c_k=v\u0;
%tspan_size=size(tspan); %part of preallocation attempt
%eigen_factors=zeros(tspan_size); %preallocation (or my attempt at that lol)
%this is where the problem begins to occur
for i=1:tspan
eigen_factors=[c_k.*exp(lambda*tspan(i))];
end
%rest of code (no problems here)
u=v*(c_k.*exp(lambda*tspan(end)));
x=[0;x];
x=[x;1];
u=[0;u;0];
end
Thank you!!!
  2 Commenti
Stephen23
Stephen23 il 17 Mar 2021
eigen_factors=[c_k.*exp(lambda*tspan(i))];
% ^ superfluous brackets ^
Square brackets are a concatenation operator. You are not concatenating anything together here, so those square brackets do nothing (except mislead). Get rid of them.
Catherine Zdunek
Catherine Zdunek il 17 Mar 2021
Thank you! I wasn't sure if those were necessary, so good to know for the future also!

Accedi per commentare.

Risposta accettata

David K.
David K. il 16 Mar 2021
I can't see right away what the dimension of
[c_k.*exp(lambda*tspan(i))]
is expected to be so I am just going to put out answers for various situations. The main problem looks to be that you are not indexing eigen_factors in the loop.
Case: 1x1
eigen_factors(i)=[c_k.*exp(lambda*tspan(i))]; % change in loop
Case: 1xn
eigen_factors = zeros(n,length(tspan)); % Preallocation
eigen_factors(i,:)=[c_k.*exp(lambda*tspan(i))];
Case: nx1
eigen_factors = zeros(length(tspan),n); % Preallocation
eigen_factors(:,i)=[c_k.*exp(lambda*tspan(i))];
Case: nxm
eigen_factors = zeros(n,m,length(tspan)); % Preallocation
eigen_factors(:,:,i)=[c_k.*exp(lambda*tspan(i))];

Più risposte (0)

Categorie

Scopri di più su Loops and Conditional Statements in Help Center e File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by