Filter löschen
Filter löschen

ODE45 error must return column vector

1 Ansicht (letzte 30 Tage)
Brianna Biondo
Brianna Biondo am 11 Mär. 2023
Kommentiert: Walter Roberson am 11 Mär. 2023
No matter what I do I keep receiving the error of a column vector. I have tried making it all zeros first befort completing the function to create a column vector but nothing is working. Any insight would be helpful.
CO(2,1) = 10;
tRange(2,1) = 10;
[tSol,CSol] = ode45(@ConvFunction,tRange,CO);
Error using odearguments
CONVFUNCTION must return a column vector.

Error in ode45 (line 107)
odearguments(odeIsFuncHandle,odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin);
plot(tSol,CSol(:,1))
function F = ConvFunction(t,~)
CaO=10;
k1=.45;
k2=.05;
Ca=CaO*exp(-t*(k1+k2));
dCadt=-(k1*(Ca)-(k2*(Ca)));
dCbdt=k1*(Ca);
dCcdt=k2*(Ca);
dFdt = zeros(3,1);
dFdt = [dCadt dCbdt dCcdt];
F= dFdt;
end

Akzeptierte Antwort

Paul
Paul am 11 Mär. 2023
Try doing what the error message says and make sure that F is a column vector. Here's one option
dFdt = [dCadt ; dCbdt ; dCcdt];
Also, the code is only specifying two initial conditions. It needs to specify three.
  2 Kommentare
Brianna Biondo
Brianna Biondo am 11 Mär. 2023
Bearbeitet: Brianna Biondo am 11 Mär. 2023
Thank you for the help, im getting this as an answer now. Sorry I am newer to matlab and am getting really confused.
CO(2,1) = 10;
tRange(2,1) = 10;
[tSol,CSol,~] = ode45(@ConvFunction,tRange,CO);
Error using odearguments
CONVFUNCTION returns a vector of length 3, but the length of initial conditions vector is 2. The vector returned by CONVFUNCTION and the initial conditions vector must have the same number of
elements.

Error in ode45 (line 107)
odearguments(odeIsFuncHandle,odeTreatAsMFile, solver_name, ode, tspan, y0, options, varargin);
plot(tSol,CSol(:,1))
function F = ConvFunction(t,~,~)
CaO=10;
k1=.45;
k2=.05;
Ca=CaO*exp(-t*(k1+k2));
dCadt=-(k1*(Ca)-(k2*(Ca)));
dCbdt=k1*(Ca);
dCcdt=k2*(Ca);
dFdt = zeros(3,1);
dFdt = [dCadt;dCbdt;dCcdt];
F= dFdt;
end
Walter Roberson
Walter Roberson am 11 Mär. 2023
Unless you have assigned something to CO before this,
CO(2,1) = 10;
creates CO as a column vector with exactly two values. You are passing that vector of length 2 as the initial state. You ignore the state inside ConvFunction and return a vector of length 3 . The length of the vector you return must be the same as the number of input values.

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

Mehr zu Programming finden Sie in Help Center und File Exchange

Produkte


Version

R2022b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by