MatCont for Homotopy Method

Question

FBBVC 2020-4-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519783-matcont-for-homotopy-method

评论： Andrew Newell 2020-4-29

Hello,

has anybody managed or can give some advice how to apply the homotopy method within matcont? My goal is to follow the path of solution x (equilibrium) to a linear shifted function H(x, t) = F(x(t)) − (1 − t)F(x0). During that path, I want to encounter possible bifurcatin points. F is a different function, that outputs just some kind of vector to the given x(t).

MatCont seems to be perfect for that application but it seems like it only handles differential equations. How can I adapt it to work with a simple (system of) equations?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Andrew Newell 2020-4-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519783-matcont-for-homotopy-method#answer_428251

编辑：Andrew Newell 2020-4-24

The primary purpose of MatCont (which is a GUI front end for Cl_MatCont) is to analyze solution curves of the form

. These may be equilibrium solutions of an ODE system

and that is how MatCont expects you to enter it. So one approach would just be to pretend it's an ODE and then find the equilibrium solutions.

MatCont creates a curve file and then uses Cl_MatCont to solve for the curve and its critical points. So another approach is to create the curve file yourself and run Cl_MatCont. You can then formulate it explicitly as an equilibrium curve instead of an ODE. This is more work and requires a much greater familiarity with Cl_MatCont, but it's a better approach if you're planning to solve a particular kind of problem many times. If you want to do that, I recommend reading the manual. In the first example it shows you how to solve

and there is no ODE.

8 个评论
显示 6更早的评论隐藏 6更早的评论

FBBVC 2020-4-27

编辑：FBBVC 2020-4-27

First of all, thank you for having this discussion with me, I really appreciate it!

I'm looking for solutions to a system F(x) consisting of nonlinear, polynomial equations that result from applying the harmonic balance method. By increasing the order N of the harmonics, the number of polynomial equations increase likewise. An adapted homotopy method (degree-increasing spectral homotopy, or just DISH) allows me to find solutions by starting at x0 (with x0 is a solution of F_N-1(x0) at lower dimension) until I reach F_N(x) at higher dimension. The method itself isn't very spectacular can be easily implemented. So the answers to your question are:

If I can track some bifurcation points during the homotopy path and follow their emerging paths, I hope to find more solutions than just one by following the main branch. So actually, it's both.
Only real solutions are relevant.
As you can find above, it's a system of polynomial equations.
Since I start homotopy at very low dimension, I can afford to compute solutions out of a large grid of initial values for Newton-Rapson-Method. Having done this, I found a set of different solutions from where I hoped to apply matcont for path continuation and track down bifurcation points.

I think that approach is somehow possible with matcont, because even if t increases up and down unpredictably, in some cases it passes t=1. I inserted the coordinates that matcont computed at that point and they are indeed solutions to F_N. Of course there might be isolated solutions that will never fuse with any branch, but some will eventually do and those are the one I'm trying to find.

Do you think it's possible to give the parameter t a generic gradient to force it "up" the curve to t=1?

Andrew Newell 2020-4-28

That doesn't sound like it's likely to be very efficient or effective. To get much advantage out of multiplying solutions by following bifurcations, you'd have to start with a very symmetric polynomial. And you'd probably have to pay close attention to each path, so you'd better not have many solutions.

It turns out that you can get smooth, monotonic homotopies "with probability 1" (i.e., aside from some isolated exceptional cases) if you are doing it in complex projective space and add some randomness into the problem. Matcont can't do that, but a package called Bertini can. And I have created a Matlab front end for this called BertiniLab that is in the Mathworks File Exchange.

The beauty of Bertini is that it is virtually guaranteed to find all of the solutions to a system of polynomial equations, you don't have to worry about finding the starting points, and it has very sophisticated "end games" that can give you as much accuracy as you like. However, you get all the complex solutions, the vast majority of which are generally not real, so you have to sift through them. In general, the number of solutions goes up rapidly with dimension or degree, so you can't do very large systems or high polynomial degree. Mind you, that's an inherent problem, and I don't know if there is any reason to expect that your system has an unusually low number of solutions.

There are some tricks that greatly improve the efficiency. One is paramatopy, which creates a good starting system from which you can rapidly move to any parameter you want. There is also regeneration, which is not easily explained in a sentence or two. These are implemented in BertiniLab. The creators of Bertini have also created Bertini_real, which is presumably much more efficient for finding real solutions, but I haven't created a Matlab front end for it.

FBBVC 2020-4-28

I will have a look into Bertini. Thank you very much so far!

Andrew Newell 2020-4-29

Glad to help!

请先登录，再进行评论。

MatCont for Homotopy Method

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

8 个评论
显示 6更早的评论隐藏 6更早的评论

另请参阅

类别

标签

Community Treasure Hunt

MatCont for Homotopy Method

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

8 个评论 显示 6更早的评论隐藏 6更早的评论

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

8 个评论
显示 6更早的评论隐藏 6更早的评论