繪制Matlab函數的問題

Question

我是Matlab的初學者，我想在一定時間間隔繪制系統濃度與時間曲線，以下是我寫的代碼

％9個樣本的輸入函數，活動和時間校准好，％計數器值約為：1.856，來自3個患者的所有9個輸入值

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

我收到了以下錯誤

操作數到|| 和&&運算符必須可轉換為邏輯標量值。

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

請。讓我知道我的邏輯是否錯誤

Answer 1

您的t不是與0比較的單個值，因此無法評估為true或false。

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