-
Notifications
You must be signed in to change notification settings - Fork 304
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
question about second-order ode identification #422
Comments
Can you share the code to create the generalized or customized library? I'm not sure if this is a coding issue or a data/regression issue. Also, FWIW, you called this a PDE, but I think you mean ODE since there are no spatial derivatives in your equation. |
Thank you for your reply, please see my code below. |
Typically, |
Hi Jacob, Thank you for you reply, please let me clarify my question. but for the first equation it should be In terms of the second equation, I want to cancel out the I add the new code related to constrainedSR3 in same repository named as R-P-SINDY-GEN-constrained.py as below. |
But If you want to set certain elements to zero, you need to understand the order that constraint_lhs = np.array([[1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 1, 0])
constraint_rhs = np.array([0, 0]) This ordering is |
Thanks for reply. And if I want to make both |
And I managed to define what I wanted about the customized library but it gave me the "ConvergenceWarning: SR3._reduce did not converge after 30 iterations." what would be the issue? |
We don't have a great way to identify the names of equations before fitting the data. The best way is to fit the data once without constraints, then check |
Try increasing the |
Hi Jacob, My code get stuck when using simulate function and after I manually stop it, it give me "capi_return is NULL Thanks in advance |
The infinity value is because you have discovered an unstable model and the integration blows up. You have to catch it by using a warnings filter to turn Also, see scipy/scipy#15940. |
Hi,
$$(R 1)^{\prime}=1.000 R 2 $$
I would like to identify second-order ode below:
but I whatever I use customised or generalized library, it always comes with R1 and R2 term in second equation which I don't want or that two terms disappear but the first equation is not what I want with$(R 1)^{\prime}=1.000 R 2 $
The text was updated successfully, but these errors were encountered: