Solving dependent functions using ODE45
显示 更早的评论
Hi All,
I am trying to solve a set of equations that are dependent on each other. dx_dt(1) is dependent on dx_dt(2). dx_dt(2) and dx_dt(3) are dependent on dx_dt(1). I have not been able to solve this. Any suggestions are greatly appreciated.
采纳的回答
Star Strider
2016-6-24
0 个投票
The ode45 and other functions (choosing the appropriate solver is important) are able to solve such systems of differential equations relatively easily. Please post your code for your differential equation funciton, and original equations if necessary.
5 个评论
Star Strider
2016-6-27
My pleasure.
Do not use globals! They cause the problems you’re seeing, because you cannot control them and they make debugging your code almost impossible. I don’t know what the ‘params’ argument refers to, since you never use it in your code.
I would do this instead:
function dx_dt = dyn(t,x, Vs, D, fsw, L, C, C2, Rload, rL, Lp1, Lp2, Cj1, Cj2)
and delete the global call.
I don’t understand what you’re doing in the if block, but be aware that the ODE solvers do not handle integrating across discontinuities at all well (No numerical solvers do.) The same caution applies to:
dx_dt(7,:) = 2*fsw*sign(cos(2*pi*fsw*t));
although if the magnitude of ‘dx_dt(7,:)’ isn’t large, especially with respect to the other values, the effect of the discontinuity may not be significant.
Star Strider
2016-6-28
The use of global variables is a remnant of earlier computer languages (like FORTRAN II that I learned programming with back in 1968) when that was the only way to pass variables and avoid some restrictions on subroutine (MATLAB function) usage. They are vestigial, and will eventually disappear. They are definitely to be avoided in MATLAB code.
One other problem you mentioned in a subsequent Question was that ‘dxdt(1)’ requires ‘dxdt(2)’. I would eliminate the double subscripts (that could be part of the problem), and instead insert a line defining ‘dxdt’ as a column vector from the outset:
...
iLoad = vC/Rload;
dxdt = zeros(7,1);
dxdt(1) = (x*(Vs - vC) + y*(-vj1 - Lp1*dxdt(2) - vC))/L; % diL/dt
...
and so for the others. I’m not certain that will solve your problem, since I didn’t run your code, but it should at least eliminate the problem you mentioned.
Torsten
2016-6-29
@Josh:
Take the first and the second equation (dxdt(1,:) = ... and
dxdt(2,:) = ...) and solve for dxdt(1) and dxdt(2) (2 linear equations in two unknowns).
Alternatively, use the mass-matrix option of the ODE integrators.
Best wishes
Torsten.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
