ode45 no enough input argument

Question

0 个投票

hi guys, i really need help.I have been trying to solve 3 differential equations with ODE 45, but it keeps giving the error that there are not enough input arguments.

This is my function:

function dxdt=funsystest(t,x) 
%given conditions
theta1=1;
theta2T =[0.5,2];
theta2=transpose(theta2T);
%a=abs(b);
phi1=x(1)^2;
n1=norm(phi1);
phi2=[x(1),sin x(2)];
n2=norm(phi2);
%umin=0.02;
c1=5;
c2=5;
etha1=1;etha2=1;
gamma=4;
sigma=0.05;
n=0.1; kd=1;
s=0.6;
%update law
 theta_cap=integral(dxdt(3));
%derivations
   z1=x1;  %first error
   alpha1=-(c1+1)*z1-(1/4*etha)*theta_cap*n1*n1*z1;%virtual control input
   alpha1dot=-(c1+1)-(1/4*etha)*theta_cap*n1*n1;
   z2=x2-alpha1;%second error
   v=(c2+0.5+kd)*z2+(1/2*etha2)*theta_cap*n2*n2*z2-alpha1dot*(z2+alpha1)+(1/2*etha2)*theta_cap*alpha1dot^2*n1^2*z2-alphacap*z2;
   %control law
   u=-(z2*v^2)/((1-s)*sqrt(z2^2*v^2+n^2));
  % udot=diff(u);
   %update law
   %tuning function
   tau1=(gamma/4*etha1)*n^2*z1^2-gamma*sigma*theta_cap;
   dxdt=zeros(2,1);
  dxdt(1) = x(2)+ phi1*theta1;
  dxdt(2) = u + phi2*theta2;
  dxdt(3)=tau1+(gamma/2*etha2)*n2^2*z2^2+alpha1dot^2*n1^2*z2^2;
  dxdt = dxdt';
end

and this is my odefile

%odefile%
clc ;
close all;
timerange = [0 3];
IC        = [0.5;0;0];
[t,x]=ode45('funsystest',timerange,IC);
subplot(2,2,1);
plot(t,x(:,1));
subplot(2,2,2);
plot(t,x(:,2));
subplot(2,2,3);
plot(t,x(:,3));

And this is the error i get.

Error using sin
Not enough input arguments.
Error in funsystest (line 12)
phi2=[x(1),sin x(2)];
Error in odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
Error in ode45 (line 115)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
Error in test (line 9)
[t,x]=ode45('funsystest',timerange,IC);
>>