Question

我是Matlab的初学者，我想在一定的时间间隔绘制系统浓度与时间曲线，以下是我编写的代码

％9个样品的输入函数，其活性和时间用Well校准％计数器值约为：来自3名患者的所有9个输入值的1.856

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

我收到以下错误

操作数到||和＆amp;＆amp;运算符必须可转换为逻辑标量值。

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

请帮助。让我知道我的逻辑是否不正确

Answer 1

您的t不是与0比较的单个值，因此无法评估为真或假。

Answer 2

t被写为数字数组。因此，它无法与标量值ex进行比较。 0。在for循环中尝试它

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

Answer 3

您希望使用逻辑索引

执行此操作

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)));

通过将逻辑索引分配给变量，您还可以使它更漂亮（更容易阅读）：

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

然后，例如，第二个赋值变得更具可读性：

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

绘制Matlab函数的问题

3 个答案: