Index exceeds the number of array elements. Index must not exceed 1.

Question

Jack Verderber 2021-11-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1596309-index-exceeds-the-number-of-array-elements-index-must-not-exceed-1

回答： Walter Roberson 2021-11-26

function res = digestive_model()
syms x(y)
c = 7; tau = 1; K = 1; c0 = 1; c1 = 1; a = 1; b = 1; v0 = 0.001;
Wt = linspace(0,100); Vt = linspace(0,100);
dx = diff(x);
at = tau*(1-(1/c)*diff(x,y,1));
ode = diff(x,y,2) + (K./Wt).*diff(x,y,1) - at.*((c0+c1.*Vt)./(a+b*x)) == 0;
cond1 = x(0) == 0;
cond2 = dx(0) == v0;
[V] = odeToVectorField(ode);
M = matlabFunction(V,'var',{'Y','X'});
res = ode45(M,[0 20],[cond1 cond2]);
figure
fplot(res,[0 20])
end

Any help would be greatly appreciated, thank you.

2 个评论
显示无隐藏无

Bjorn Gustavsson 2021-11-26

Dont attach the code as an image. Insert it as code in the message, using the editor facility for that. If someone is to help you they will for now write off what can be seen in the image - why making it harder for us to help you?

Jack Verderber 2021-11-26

My apologies, first time asking a question on here so was unsure on the best way to do it. Should be fixed now.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2021-11-26

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1596309-index-exceeds-the-number-of-array-elements-index-must-not-exceed-1#answer_841284

在 MATLAB Online 中打开

digestive_model()

ode(y) =

Error using mupadengine/feval_internal
Unable to convert the initial value problem to an equivalent dynamical system. Either the differential equations cannot be solved for the highest derivatives or inappropriate initial conditions were specified.

Error in odeToVectorField>mupadOdeToVectorField (line 171)
T = feval_internal(symengine,'symobj::odeToVectorField',sys,x,stringInput);

Error in odeToVectorField (line 119)
sol = mupadOdeToVectorField(varargin);

Error in solution>digestive_model (line 20)
[V] = odeToVectorField(ode);

function res = digestive_model()

syms x(y)

c = 7; tau = 1; K = 1; c0 = 1; c1 = 1; a = 1; b = 1; v0 = 0.001;

Wt = linspace(0,100); Vt = linspace(0,100);

dx = diff(x);

at = tau*(1-(1/c)*diff(x,y,1));

ode = diff(x,y,2) + (K./Wt).*diff(x,y,1) - at.*((c0+c1.*Vt)./(a+b*x)) == 0;

cond1 = x(0) == 0;

cond2 = dx(0) == v0;

ode

[V] = odeToVectorField(ode);

M = matlabFunction(V,'var',{'Y','X'});

res = ode45(M,[0 20],[cond1 cond2]);

figure

fplot(res,[0 20])

end

Notice your ode is a vector -- a vector the same length as Wt and Vt since those are both vectors. If it was going to work at all, you would need a vector of initial conditions at least as long as that vector.

Also, the first entry in Wt is 0, so K./Wt includes a division by 0, so the first entry in ode involves an infinity.

If Wt and Vt are parameters of the ode, then use them as symbolic, and use odeFunction() instead of matlabFunction and specify Wt and Vt as parameters of the system. You may need to iterate over all the Wt / Vt combinations to solve all of the cases.