Geek asks and answers

Matlab function build issues

I start in Matlab, I would like to calculate the concentration of the system and the timeline in a certain period of time following the code that I wrote

% Input function of 9 samples with activity and time calibrated using a well% counter value approximately: 1.856 of all 9 input values ​​of 3 patients

function c_o = Sample_function(td,t_max,A,B)

   t   =(0 : 100  :5000); % time of the sample post injection in mins
   c   =(0 : 2275.3 :113765);

   A_max= max(c);   %Max value of Concentration (Peak of the curve)

   if (t >=0 && t <= td)
      c_o(t)=0;
   else if(td <=t && t<=t_max)
      c_o(t)= A_max*(t-td);
   else if(t >= t_max)
      c_o(t)=(A(1)*exp(-B(1)*(t-t_max)))+(A(2)*exp(-B(2)*(t- t_max)))+...
             (A(3)*exp(-B(3)*(t-t_max)));
   end

   fprintf('plotting Data ...\n');
   hold on;
   figure;
   plot(c_o);
   xlabel('Activity of the sample Ba/ml ');
   ylabel('time of the sample in minutes');
   title (' Input function: Activity sample VS time ');
   pause;
   end

I get the following error

Operands to || and && operators must be converted to logical scalar values.

Error in Sample_function (line 18)
if (t >=0 && t <= td)

Request. Let me know if my logic is wrong

+3
source share
3 answers

You want to do this with logical indexing

c_o = zeros(size(t));

c_o(t>=0 & t<=td) = 0; % this line is actually redundant and unnecessary since we initialized the vector to zeros
c_o(t>td & t<=t_max) = A_max*(t(t>td & t<=t_max)-td);
c_o(t>t_max) = (A(1)*exp(-B(1)*(t(t>t_max)-t_max)))+(A(2)*exp(-B(2)*(t(t>t_max)- t_max)))...
      +  (A(3)*exp(-B(3)*(t(t>t_max)-t_max)));

You can also make it a little prettier (and easier to read) by assigning logical indexes to variables:

reg1 = (t>=0 & t<=td);
reg2 = (t>td & t<=t_max);
reg3 = (t>t_max);

, , :

c_o(reg2) = A_max*(t(reg2)-td);
+1

t 0, true false.

+2

t is written as an array of numbers. Thus, it cannot be compared with the scalar value ex. 0. Try in a for loop

for i=1:length(t)
   if (t(i) >=0 && t(i) <= td)
      c_o(t(i))=0;
   else if(td <=t(i) && t(i)<=t_max)
      c_o(t(i)))= A_max*(t(i)-td);
   else if(t(i) >= t_max)
      c_o(t)=(A(1)*exp(-B(1)*(t(i)-t_max)))+(A(2)*exp(-B(2)*(t(i)- t_max)))...
      +  (A(3)*exp(-B(3)*(t(i)-t_max)));

   end
end
+1
source

More articles:


All Articles