Substitution of initial numerical values before odefunction
2 次查看(过去 30 天)
显示 更早的评论
Hi guys,
I am using ode45 to solve a system of ODEs. According to the syntax of ode45, the system is formatted into the form of
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/230083/image.png)
where x(t) is a 32*1 symbolic vector in our project whose elements are to be solved in ode45.
Then odefunction is applied to build the handle for mass matrix like
Mass_handle = odefunction( M( t, x(t) ) )
However, the size of M( t, x(t) ) is very huge and the elements inside it are really complicated, then it is quite time-consuming to return Mass_handle: by now it has taken over 96 hours and consumes over 45GB in RAM in our PC and it is still running at this moment.
Here I have an idea: since we have the initial numerical values of x(t), are we able to substitute the 32 symbolic variables in x(t) by their initial numerical values before running the code
Mass_handle = odefunction( M( t, x(t) ) )
Because we notice that in ode45, it is necessary to input the initial values of the symbolic variables, according to its syntax like ode45(odefun, tspan, y0), where y0 is the initial value vector.
If we substitute the initial numerical values of x(t) before running the code odefunction( M( t, x(t) ) ), then M( t, x(t) ) will be numerical rather than symbolic. It will just take a second to return Mass_handle.
Here we want to know: is this way feasible ? if not, are there any methods to deal with this problem?
Many thanks!
Cheers
0 个评论
采纳的回答
Torsten
2019-7-17
Use "matlabFunction" to convert every symbolic input to ODE45 into numerical function handles before calling ODE45.
2 个评论
Torsten
2019-7-19
Yes, that's what I meant.
Before calling ode45 with the function handle, you should output the handle into a file to see the list of numerical parameters and the order in which they have to be supplied.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!