save loop data with non integer steps

Question

Jerry 2023-2-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1915975-save-loop-data-with-non-integer-steps

评论： Jerry 2023-2-21

采纳的回答： Torsten

在 MATLAB Online 中打开

I have following script and I am trying to save the results of the loop in an array. I get following error message:

Array indices must be positive integers or logical values.

How can I save the results in an array even though my steps are not integer values?

reduced_scattering_coeff = 10;
absorbing_coeff = 1;
source_detector_distance= 1.5;
refection_parameter= 3.25;
pi = pi;
detectorsignal=zeros(size(0.01:0.01:1));
for absorbing_coeff = 0.01:0.01:1
    
    mu_t = reduced_scattering_coeff+absorbing_coeff;
    mu_eff = sqrt(3*(reduced_scattering_coeff+absorbing_coeff));
    r1 = sqrt(((1/mu_t)^2)+source_detector_distance);
    r2 = sqrt((((4/3*refection_parameter+1)/mu_t)^2)+(source_detector_distance^2));
    detectorsignal= (1/(4*pi*mu_t))*((mu_eff+1/r1)*((exp(mu_eff*r1))/r1^2)+((4/3*refection_parameter)+1))*(mu_eff+1/r2)*((exp(-mu_eff*r2))/(r2^2));
    detectorsignal(absorbing_coeff)= detectorsignal;
end

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

请先登录，再进行评论。

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

Answer 1

Torsten 2023-2-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1915975-save-loop-data-with-non-integer-steps#answer_1176000

在 MATLAB Online 中打开

reduced_scattering_coeff = 10;

absorbing_coeff = 1;

source_detector_distance= 1.5;

refection_parameter= 3.25;

absorbing_coeff = 0.01:0.01:1;

detectorsignal = zeros(size(absorbing_coeff));

for i=1:numel(absorbing_coeff)

mu_t = reduced_scattering_coeff+absorbing_coeff(i);

mu_eff = sqrt(3*(reduced_scattering_coeff+absorbing_coeff(i)));

r1 = sqrt(((1/mu_t)^2)+source_detector_distance);

r2 = sqrt((((4/3*refection_parameter+1)/mu_t)^2)+(source_detector_distance^2));

detectorsignal(i)= (1/(4*pi*mu_t))*((mu_eff+1/r1)*((exp(mu_eff*r1))/r1^2)+((4/3*refection_parameter)+1))*(mu_eff+1/r2)*((exp(-mu_eff*r2))/(r2^2));

end

plot(absorbing_coeff,detectorsignal)

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

Jerry 2023-2-21

this works perfectly, thank you!

请先登录，再进行评论。

Answer 2

John D'Errico 2023-2-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1915975-save-loop-data-with-non-integer-steps#answer_1175995

在 MATLAB Online 中打开

acvalues = 0.01:0.01:1;
detectorsignal=zeros(size(acvalues));
for ind = 1:numel(acvalues)
    absorbing_coeff = acvalues(ind);
    
    % stuff. you should get the idea...
    
    
    
    detectorsignal(ind)= (1/(4*pi*mu_t))*((mu_eff+1/r1)*((exp(mu_eff*r1))/r1^2)+((4/3*refection_parameter)+1))*(mu_eff+1/r2)*((exp(-mu_eff*r2))/(r2^2));
end