Using ode15s to solve daes and lsqcurvefit to fit a curve

Question

Dursman Mchabe 2019-1-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/439460-using-ode15s-to-solve-daes-and-lsqcurvefit-to-fit-a-curve

Hello everyone,

I am trying to use ode15s to solve 4 daes,and using lsqcurvefit to fit the curve. The attached code gives an error message:

Index exceeds the number of array elements (1).
Error in ODE15sMethod>f (line 67)
daes = [(c(1)./c(2)).*((1930.65./1000000)- (x(1,:) .*c(3) .*c(4) ./ 875000 )) - (x(1,:).* c(3).* c(4) - c(6).* ((x(2,:).* x(4,:).^2)./(x(4,:).^2 + c(7).* x(4,:) +
c(7).* c(8))))./((1./c(12)) + c(6)./((1 + c(9).* ((c(7).* c(6).* x(1,:).* c(3).* c(4) ./ x(3,:)) - ((x(2,:).* c(7).* x(4,:))./(x(4,:).^2 + c(7) .* x(4,:) + c(7)
.*c(8))))./2.85E-09.* ((c(6).* x(1,:).* c(3).* c(4)) - ((x(2,:).* x(4,:).^2)./(x(4,:).^2 + c(7).* x(4,:) + c(7).* c(8)))) + c(11).* ((c(8).* c(7).* c(6).* x(1,:).*
c(3).* c(4) ./ (x(3,:)).^2) - ((x(2,:).* c(7).* c(8))./(x(4,:).^2 + c(7) .* x(4,:) + c(7) .*c(8))))./ 2.85E-09.* ((c(6).* x(1,:).* c(3).* c(4)) - ((x(2,:).*
x(4,:).^2)./(x(4,:).^2 + c(7).* x(4,:) + c(7).* c(8))))).* c(13)))
Error in odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
Error in ode15s (line 150)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
Error in ODE15sMethod (line 17)
[t,y] = ode15s(@f,tspan,x0,options);