Using dsolve, replace arbitrary constants like C11 with one's own parameter name.
17 次查看(过去 30 天)
显示 更早的评论
when one solves a differential equation using dsolve, without specifying terminal conditions, matlab outputs the solution with an arbitrary parameter name. Here's an example.
syms c T lambda0 x(tau) x0 xT a tau
xDot = @(tau) a*x - exp(-a*tau)/lambda0;
diffEqSolve = dsolve(diff(x,1) == xDot(tau));
diffEqSolve is now:
C3*exp(a*tau) + exp(-a*tau)/(2*a*lambda0)
Since I known that the solution to this problem will have just one parameter, I'd like to pass my own symbol to dsolve. E.g., pass the symbolic variable 'c' to dsolve, and have it return
c*exp(a*tau) + exp(-a*tau)/(2*a*lambda0)
If this isn't possible (I suspect it isn't), could somebody please explain how matlab decides what number to append to C? I've had the same string of commands return C3 or C11, depending on what commands have be processed before dsolve is invoked, and have been unable to see any pattern. If I could compute the number beforehand, then I could do a subs command and get the result I want, e.g.,
diffEqSolve = subs(dsolve(diff(x,1) == xDot(tau)),'C3','c')
Thanks in advance for any suggestions.
0 个评论
采纳的回答
Walter Roberson
2015-12-6
There is no pattern. Or rather, the symbolic engine uses the next one that it has not used before in the session, including using it internally during computations that are not shown to you. Calculations will not necessarily get done in any particular order, either.
You can use symvar() on the expressions before and after to find the names that have been introduced, to give you a list of constants of integration without context in the expression. You could pull the expression apart if you really wanted to, or you could just rename them without caring where they occur or what they mean.
If you want the new names to have meaning reflecting what they are for, you should probably be using initial conditions; it is completely allowed to use symbolic initial conditions such as x(0)==x0 and the computation engine will use that information. Caution, though: it is not uncommon that when you specify an initial condition, the computation engine ends up discarding a valid computation path and is not able to find a solution when one exists. It is also common for pretty strange computations to show up, like x being a function of tau involving a constant omega19 where omega19 is the root of a integral expression involving a function that itself involves the root of an integral expression...
2 个评论
Walter Roberson
2015-12-7
You should not do that. char() of a symbolic expression is not necessarily going to give you something that uses MATLAB syntax. Symbolic expressions use a language that is like MATLAB but has some subtle differences.
matabConst = setdiff(symvar(diffEqSolve),[a,tau,lambda0]);
diffEqSoln = matlabFunction( subs(diffEqSolve, matlabConst, sym('c')) );
更多回答(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!