|
From: Emery C. <ec...@vt...> - 2005-09-12 15:37:59
|
<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Vijay,<br>
<br>
Now that limit cycle continuation works, that would be your preferred
method to get period/min/max done (although it may take longer). You
first must have a good 1 parameter bifurcation diagram; then continue
from a HB using limit cycle continuation. The output is the same as
that from sample period/min/max, but is more intelligent (see attached
png). BEWARE that limit cycle continuation can take a VERY long time if
you leave the settings ds, dsmin and dsmax as they are, but you will
have nice results. If you set ds, dsmin too large, you will get
spurious results (they should be obviously bad...).<br>
<br>
BUT, if you want a quick sampling, this can be very quick IF YOU SET
THE PARAMETERS CORRECTLY. This requires that you study the wave forms
and changes in period for some small perturbations BEFORE RUNNING THE
SAMPLING! If you want to use the sample period/min/max, you need to
ensure that dp is small enough so that you get several samples within
the region of oscillation (i.e., look at the 1p bd first), but make
sure that x_match_tol and period_match_tol is large enough (I changed
period_max = 10000, x_match_tol = 60, period_match_tol = 10 to get the
second attached graph). Sometimes you will get inaccuracies or
anomalies, but you can cut these out of the data set...<br>
<br>
It all takes work, but you just must be patient and analyze the system
locally to get the basic behavior, and then set your parameters.<br>
<br>
FYI, limit cycle continuation works MUCH BETTER, but is slower; once
you have a basic picture with the "sample period/min/max" run a
lengthy, accurate limit cycle continuation to get nicer graphs.<br>
<br>
[BTW, I added you and Herbert to the oscill8-help list; this is where
all the hints and help will go to from now on]<br>
<br>
Emery<br>
<br>
<br>
<br>
Vijay Chickarmane wrote:
<blockquote
cite="mid...@de..."
type="cite">
<blockquote style="margin-right: 0px;" dir="ltr">
<span class="104151705-10092005"><font face="Tahoma" size="2">Hi
Emory,</font></span>
<blockquote style="margin-right: 0px;" dir="ltr">
<div><span class="104151705-10092005"><font face="Tahoma" size="2">I am
having trouble with the min-max of oscillations for a couple of models,
so I would be</font></span></div>
<div><span class="104151705-10092005"><font face="Tahoma" size="2">grateful
for your help. I tried to get the max-min for the foll</font></span></div>
<div><span class="104151705-10092005"><font face="Tahoma" size="2">two
models, but was unable to. Could you try and generate the plots, </font></span><span
class="104151705-10092005"><font face="Tahoma" size="2">and send them
as well</font></span></div>
<div><span class="104151705-10092005"><font face="Tahoma" size="2">as
the parameter settings you used to generate them.</font></span></div>
<div><span class="104151705-10092005"></span> </div>
<div><span class="104151705-10092005"></span><span
class="104151705-10092005"><font face="Tahoma" size="2">model 1:
mapkfinal.xml ------ bifurcation parameter J5_k (currently set at
0.0001)</font></span></div>
<div><span class="104151705-10092005"></span> </div>
<div><span class="104151705-10092005"><font face="Tahoma" size="2">model
2: mapknew.xml ------- bifurcation parameter kI (currently set at
100)</font></span></div>
<div><span class="104151705-10092005"></span> </div>
<div><span class="104151705-10092005"><font face="Tahoma" size="2">Thanks
in advance,</font></span></div>
<div><span class="104151705-10092005"><font face="Tahoma" size="2">best,</font></span></div>
<div><span class="104151705-10092005"><font face="Tahoma" size="2">vijay</font></span></div>
</blockquote>
</blockquote>
</blockquote>
</body>
</html>
|