MATLAB - 如何從矢量讀取單個值？ （ode45參與）

Question

我正在嘗試使用一個涉及ode45的代碼，在將要解決的方程式中，我必須從矢量中引入一個單值，這個值將根據解決方程式的時間而改變，我的意思是如果x = x（0），則u1（1）和u2（t），當x = x（1）時，則u1（2），u2（2）...這是我的代碼：

function [sal] = scobe1 
      clear 
      clc 
    global a; 
      a=1; 
      assignin('base','a',a) 
      tspan = [0; 180];    
      x0 = [80; 0]; 
      [t,x] = ode45(@f,tspan,x0);   
      sal=[t,x] 
      figure 
      subplot(2,1,1) 
      plot(t,x(:,1),'k-'); 
      subplot(2,1,2) 
      plot(t,x(:,2),'b-'); 
    function dxdt = f(t,x) 
     global a s1 u1 u2 newData; 
newData=importdata('datosnoembarazo.mat'); 
assignin('base','newData',newData); 
u1=getfield(newData,'glucosa4dias'); 
assignin('base','u1',u1); 
u2=getfield(newData,'insulina4dias'); 
assignin('base','u2',u2); 
u1=getfield(newData,'glucosa4dias'); 
assignin('base','u1',u1); 
u2=getfield(newData,'insulina4dias'); 
assignin('base','u2',u2); 
      if (a<=s1) 
      dxdt = [  
       (-4.9e-2-x(2))*x(1) + (4.42 + u1(a)) %(1) IN HERE u1 IS THE VECTOR'S NAME 
       -9.1e-2*x(2) + u2(a) %(2) IN HERE u2 IS THE VECTOR'S NAME 
       ]; 
      a=a+1 
      assignin('base','a',a); 
      else return 
      end 

問題是，它給我發出了這個錯誤：我不知道代碼出了什麼問題，或者我還能做些什麼才能讀取它，你能幫我嗎？由於

Error in scobe1>f (line 36) 
     global a s1 s2 u1 u2 newData; 
Error using feval 
Output argument "dxdt" (and maybe others) not assigned during call to 
"C:\Users\AnnieA\Dropbox\Tesis (A. Olay)\GUI Examples\Resumen de Modelos\scobe1.m>f". 
Error in odearguments (line 87) 
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. 
Error in ode45 (line 113) 
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 
Error in scobe1 (line 25) 
     [t,x] = ode45(@f,tspan,x0); 

Answer 1

你的功能function dxdt = f(t,x)已在if聲明兩個分支：如果a<=s1，然後dxdt得到正確分配。如果a>s1，dxdt根本沒有分配。在return之前，您可以通過在if陳述的else分支中爲dxdt分配值來修復錯誤。我不確定在這種情況下什麼是合適的價值。也許0？