  Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620   Was thanked: 1382 time(s) in 806 post(s)
|
SADEL    
Functions list: manzhuk(5|6|7)DAE (ODE) solvers: manzhuk(init, x1, x2, intvls, D) Arguments: - init is either a vector of n real initial values, where n is the number of unknowns (or a single scalar initial value, in the case of a single ODE). - x1 and x2 are real, scalar endpoints of the interval over which the solution to the ODE(s) is evaluated. Initial values in init are the values of the ODE function(s) evaluated at x1. - intvls is the integer number of discretization intervals used to interpolate the solution function. The number of solution points is the number of intervals + 1. - D is a vector function of the form D(x,y) specifying the right-hand side of the system manzhuk(init, x1, x2, intvls, G, J1) manzhuk(init, x1, x2, intvls, G, J1, J2)      manzhuk. Kinetic curves.sm (14kb) downloaded 68 time(s). manzhuk. Van der Pol oscillator 1.sm (11kb) downloaded 59 time(s). manzhuk. Van der Pol oscillator 2.sm (14kb) downloaded 51 time(s). manzhuk. Lorenz-Attractor.sm (17kb) downloaded 52 time(s). manzhuk. Kinetic curves.pdf (39kb) downloaded 51 time(s). manzhuk. Van der Pol oscillator 1.pdf (47kb) downloaded 41 time(s). manzhuk. Van der Pol oscillator 2.pdf (63kb) downloaded 43 time(s). manzhuk. Lorenz-Attractor.pdf (66kb) downloaded 44 time(s).Manual (rus): DAE-solver-manzhuk.pdf (3,499kb) downloaded 76 time(s).Links: 1. SADEL. Edited by user 15 August 2019 13:57:03(UTC)
| Reason: Not specified |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
 5 users thanked uni for this useful post.
|
on 21/07/2018(UTC), on 21/07/2018(UTC), on 21/07/2018(UTC), on 23/07/2018(UTC), on 07/12/2018(UTC)
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
|
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
2 users thanked uni for this useful post.
|
on 22/07/2018(UTC), on 22/07/2018(UTC)
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
Updated. manzhuk() now more stable. Don't use hmx parameter. By default, its value is dt / 2. |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
1 user thanked uni for this useful post.
|
|
|
|
 Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866  Was thanked: 983 time(s) in 811 post(s)
|
Originally Posted by: uni  Updated. Short form added. Thanks Uni. The complete Smath code might be educative. mwode45 solves as well. Jean ODE manzhuk Arenstorf Orbit.sm (22kb) downloaded 55 time(s).
|
|
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 983 time(s) in 811 post(s)
|
|
|
|
|
|
|
  Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740  Was thanked: 318 time(s) in 268 post(s)
|
Hello uni,
It seems this is a quite interesting ODE-DAE solver. I did not understand what the matrices J1 and J2 are. The manual is quite complicated for me. I just supposed that J1,J2 have something to do with DAE (differential algebraic equations). It seems to me that J1 is the Jacobian matrix regarding the vector of parameters (third argument in the function G()). I suppose J2 is again some kind of Jacobian matrix.
If one of your motivation was to have a DAE solver in SMath, I appreciate that very much.
Looking forward for some more examples from you about using this solver - especially DAE systems.
Regards,
Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
I'll add the examples later.  Edited by user 23 July 2018 16:43:25(UTC)
| Reason: Not specified |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
3 users thanked uni for this useful post.
|
on 23/07/2018(UTC), on 23/07/2018(UTC), on 23/07/2018(UTC)
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 983 time(s) in 811 post(s)
|
Originally Posted by: omorr If one of your motivation was to have a DAE solver in SMath, I appreciate that very much. DAE are each specific, I don't think a generalized solver can be designed, only for odd examples ... I may be wrong , and wish the best. 
|
|
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
I have more than 12 examples of how to use it, but, unfortunately, symbolic interpreter slows down the computation of the algorithm. So I decided to show an alternative version of its use, using the Script component. I hope that I will succeed. |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
2 users thanked uni for this useful post.
|
on 23/07/2018(UTC), on 23/07/2018(UTC)
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
|
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
1 user thanked uni for this useful post.
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
Updated. Changes: - zero point issue partially fixed; - added support for the simple D(t,x) definition. Edited by user 06 December 2018 22:24:09(UTC)
| Reason: Not specified |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
|
|
|
|
|
 Rank: Newbie
Groups: Registered
Joined: 12/03/2020(UTC) Posts: 4 Location: krasnoyarsk
|
Good afternoon! I use sadel to calculate the parameters of the transient in the RLC-circuit. Why doesn't time go by in calculations? file attached model2_sadel_p.sm (26kb) downloaded 35 time(s).
|
|
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
Можно посмотреть на описание задачи в изначальной форме? |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
|
|
|
|
|
Rank: Newbie
Groups: Registered
Joined: 12/03/2020(UTC) Posts: 4 Location: krasnoyarsk
|
|
|
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
Для начала нужно в выражении G(t,x, xp) исправить xp на px. Смущает выражение для второй производной. Кажется, что первая и вторая производные dx/dt должны быть как-то связаны. Я не составлял подобных систем и с ходу не могу сказать что не так. Хорошо бы ещё посмотреть на результат моделирования схемы в специализированной программе, чтобы убедиться, что начальные условия и параметры схемы действительно приводят к какому-то не нулевому результату. Ещё, мне кажется, что не обязательно использовать библиотеку SADEL для решения этой системы, а можно попробовать использовать любые другие решатели, чтобы избавиться от px-параметра. Edited by user 13 March 2020 14:04:44(UTC)
| Reason: Not specified |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
|
|
|
|
|
Rank: Newbie
Groups: Registered
Joined: 12/03/2020(UTC) Posts: 4 Location: krasnoyarsk
|
Originally Posted by: uni Для начала нужно в выражении G(t,x,xp) исправить xp на px. действительно! не заметила. Но конечно это исправление не помогло, программа видимо умнее меня, сама догадалась что имелось ввиду Originally Posted by: uni Смущает выражение для второй производной. Кажется, что первая и вторая производные dx/dt должны быть как-то связаны. Я не составлял подобных систем и с ходу не могу сказать что не так. dx/dt=v, dv/dt=a. Когда составлялась система дифуров вместо второй производной использовалась первая Originally Posted by: uni Хорошо бы ещё посмотреть на результат моделирования схемы в специализированной программе, чтобы убедиться, что начальные условия и параметры схемы действительно приводят к какому-то не нулевому результату. это хорошая идея, попробую. Просто я опираюсь на источники, где описаны именно такая схема, именно такая система дифуров и именно такие входные параметры. Я пытаюсь пока что просто повторить. Это дает определенную уверенность, что должно работать. Хотя не исключено, что где-то упущена ошибка. То, что в результатах расчета время не движется - это такая реакция программы на то, что она не может посчитать? Не могу понять о чем это говорит
|
|
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
Нужно, чтобы ier был равен нулю, тогда можно делать интерпретации результата. Этого можно достичь, если исправить указанную ошибку и отключить задание минимального шага интегрирования (hmn) или сделать его меньше. Дальше нужно проверить опечатки в системе. Я посмотрю под отладчиком почему время не меняется. Edited by user 13 March 2020 14:31:53(UTC)
| Reason: Not specified |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
|
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
Можно писать вот так, чтобы было нагляднее, но промежуточные вычисления в этом случае замедлят работу функции.  |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
|
|
|
|
|
Rank: Newbie
Groups: Registered
Joined: 12/03/2020(UTC) Posts: 4 Location: krasnoyarsk
|
Originally Posted by: uni Можно писать вот так, чтобы было нагляднее, но промежуточные вычисления в этом случае замедлят работу функции. вы имеете ввиду нагляднее для того, кто пишет/читает функцию? т.е. другая форма записи? это здорово - я попробовала так написать, так тоже запускается расчет, правда результат тот же плачевный :'(
|
|
|
|
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,620 Was thanked: 1382 time(s) in 806 post(s)
|
По сути это та же форма записи, так как результат работы функции определяется последним выражением блока программы. Определение переменных из вектора k можно вынести из программного блока на лист, т.к. это константы. Ещё нужно всем выражениям указать численный тип оптимизации в свойствах контекстного меню. Что касается численного моделирования, то можно попробовать посчитать в Proteus данную схемку. Там только нужно подумать как имитировать ключ. А что за книжка-то? Edited by user 13 March 2020 17:27:09(UTC)
| Reason: Not specified |
Russia ☭ forever (git.io/vOZo0)
Viacheslav N. Mezentsev |
|
|
|
|
|
Forum Jump
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.
