How i modifed this for loop so it get the last value Kf_LMax of previous iteration not from the computed values before the loop?

Question

Ehtisham 2024-10-1，7:16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2156635-how-i-modifed-this-for-loop-so-it-get-the-last-value-kf_lmax-of-previous-iteration-not-from-the-com

编辑： Ehtisham about 21 hours 前

function [Kf_LMax] = Kf_Cal(Kf_Max, RLC, TauKf_ON, TauKf_OFF, PhaseTimes, t)

%compute the maximum Kf_L based on RLC values

Kf_LMax_values = Kf_Max * (RLC ./ (1 + RLC));

% Determine which phase each time step falls into

num_phases = numel(RLC);

phase_idx = zeros(size(t));

for j = 1:num_phases

T_start = PhaseTimes(j);

if j < num_phases

T_end = PhaseTimes(j + 1);

else

T_end = inf;

end

phase_idx(t >= T_start & t < T_end) = j;

end

% Calculate Kf_L for each time step based on the phase

Kf_LMax = zeros(size(t));

for i = 1:length(t)

j = phase_idx(i);

if j == 1

Kf_LMax(i) = Kf_LMax_values(j);

else

prev_end = PhaseTimes(j - 1);

if RLC(j - 1) < RLC(j)

Kf_LMax(i) = Kf_LMax_values(j) - (Kf_LMax_values(j) - Kf_LMax_values(j - 1)) * exp(TauKf_ON * (t(i) - prev_end));

else

Kf_LMax(i) = Kf_LMax_values(j) + (Kf_LMax_values(j) - Kf_LMax_values(j - 1)) * exp(TauKf_OFF * (t(i) - prev_end));

end

How i modifed this for loop so it get the last value Kf_LMax of previous iteration not from the computed values in bold before the loop?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

nick 2024-10-1，7:54

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2156635-how-i-modifed-this-for-loop-so-it-get-the-last-value-kf_lmax-of-previous-iteration-not-from-the-com#answer_1525145

在 MATLAB Online 中打开

Hi Ehtisham,

I understand that you want to modify the given function so that each iteration of the loop uses the last computed value of 'Kf_LMax' from the previous iteration instead of 'Kf_LMax_values(j - 1)'.

However, it would be really helpful to have more background information regarding the objective of the given function and the variables and any related files or data involved to accurately reproduce the results. Nevertheless, here's a suggestion based on my understanding on how you can modify the loop:

 Kf_LMax = zeros(size(t));
    Kf_LMax(1) = Kf_LMax_values(1); % Start with the first phase value
    % Calculate Kf_L for each time step based on the phase
    for i = 2:length(t) 
        j = phase_idx(i);
        if j == 1
            Kf_LMax(i) = Kf_LMax_values(j);
        else
            prev_end = PhaseTimes(j - 1);
            if RLC(j - 1) < RLC(j)
                Kf_LMax(i) = Kf_LMax(i-1) - (Kf_LMax_values(j) - Kf_LMax(i-1)) * exp(TauKf_ON * (t(i) - prev_end));
            else
                Kf_LMax(i) = Kf_LMax(i-1) + (Kf_LMax_values(j) - Kf_LMax(i-1)) * exp(TauKf_OFF * (t(i) - prev_end));
            end
        end
    end

You can refer to the following documentation to learn more about 'Array Indexing' :

https://www.mathworks.com/help/matlab/math/array-indexing.html

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

Ehtisham about 23 hours 前

编辑：Ehtisham about 21 hours 前

@nick @Star Strider thanks for your response but my question was that Kf_LMax(i) = Kf_LMax_values(j) - (Kf_LMax_values(j) - Kf_LMax_values(j - 1)) * exp(TauKf_ON * (t(i) - prev_end));

else

Kf_LMax(i) = Kf_LMax_values(j) + (Kf_LMax_values(j) - Kf_LMax_values(j - 1)) * exp(TauKf_OFF * (t(i) - prev_end));

end

in this loop for each iteration as is using the time of the prev_end so for Kf_LMax values is used the last value of the array of Kf_LMax_value for the next iteration.

It given the exact calculation of the Kf_LMax values given in plot