2 views (last 30 days)

Show older comments

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!!!

David K.
on 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))];

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

Start Hunting!