How to parse a Matrix ODE45. Getting simulation error

3 次查看（过去 30 天）

Shy 2020-7-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/570928-how-to-parse-a-matrix-ode45-getting-simulation-error

评论： Shy 2020-7-27

Hi Dears!

I am simulating a set of Ordinary Differntial Equations ODEs using ODE45 solver. I need to parse a matrix (i.e. zeta=rand(46,16)) into the ODE45, however I am getting this error:

Unable to perform assignment because the left and right sides have a different number of elements.

Error in OscillatorExample_V2>LinearOscillatorFree (line 26)

y(2)=(-2*zeta*y(2))-(w*y(1));

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 OscillatorExample_V2 (line 9)

[T, Y] = ode45(@LinearOscillatorFree, tspan, initcond, [],zeta);

Can you please help me with this to get the right results? Many thanks in advance!

close all; clear all; clc;
zeta=rand(46,16);
tspan=[0 4700];
initcond= [0 1];
for n = 1:length(zeta)
[T, Y] = ode45(@LinearOscillatorFree, tspan, initcond, [],zeta);
figure(1);
plot(T, Y(:,1));
hold on
end
m= length(Y);
for i=1:m
y(i,:)=feval('LinearOscillatorFree',T(i),Y(i,:),zeta);
end
y=y';
function dy = LinearOscillatorFree(t, y,zeta)
w=0.5;
y(1)=y(2);
y(2)=(-2*zeta*y(2))-(w*y(1));
dy= [y(1); y(2)];
end

7 个评论
显示 5更早的评论隐藏 5更早的评论

James Tursa 2020-7-27

编辑：James Tursa 2020-7-27

Your equations don't make sense to me. How can y1 and y2 be scalars and yet zeta be a matrix? The dimensions of the equations don't match. It might make sense for zeta fo be a scalar that is randomly selected from a distribution, and then run Monte-Carlo on that. But what you have written doesn't make sense. Where did these equations come from? Why the 46x16 size of zeta?

Shy 2020-7-27

I am running a similar system like that consisting of four ODEs, so the question is how to parse a matrix into the ode45 solver? Can you please provide me with your email so that I can email you on?

请先登录，再进行评论。

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