1回答
我正在尝试使用Assimulo包来求解一组微分代数方程(DAE)。我正在尝试使用类似于该shown here的算法。但是，似乎没有传递稀疏矩阵的选项。您是否知道是否有一种方法可以在计算上更高效地解决我的DAE？
浏览 41提问于2020-09-02得票数 0
我的问题是：
对于如何在Modelica模型中求解DAE方程，有什么解释吗?这样我就可以理解Jacobian矩阵的使用了吗？
浏览 3提问于2020-07-12得票数 4
回答已采纳
每当我尝试运行使用.dae扩展的代码时，都会收到此错误消息。 我在代码中遇到了这个错误，这是我一直在解决的问题。最后，我运行了我在网上找到的针对不同问题的代码，并得到了相同的错误。这使我认为这是我运行求解器(ipopt)而不是代码本身的一些错误 下面是指向我收到错误消息的GitHub代码的链接：https://github.com/Pyomo/pyomo/blob/master/examples/dae/Optimal_Control
浏览 20提问于2019-06-12得票数 0
回答已采纳
1回答
具有大矩阵的渐近编译函数
、、、
我一直在使用sympy来处理微分方程组。我象征性地编写方程，使用autowrap通过cython编译它们，然后将得到的函数传递给compile求解器。这样做的主要好处之一是，我可以象征性地使用渐近jacobian函数为jacobian求解，编译它，并将其编译到ODE求解器中。
对于大约30个变量的系统来说，这是非常有效的。最近我试着用150个变量来做这件事，结果是我在编译jacobian函数时内存用完了。这是在和microsoft C
浏览 0提问于2017-09-18得票数 3
回答已采纳
我使用MATLAB的函数，并试图通过选项JacobPattern设置用户定义的Jacboian模式。我设置了要使用的trust-region-reflective算法的首选项，来自lsqnonlin的output表示这确实是求解器使用的算法(使用JacobPattern选项所必需的)。我发现的问题是，如果我的JacobPattern太稀疏(例如，在500x500Jacobian中只有几行)，它将被求解器忽略，而整个Jacobian则被计算出来。我希望能
浏览 2提问于2012-07-26得票数 4
我试图用SciPy中的精确牛顿法求解一个大型非线性系统.在我的应用程序中，Jacobian易于组装(并分解)为稀疏矩阵。似乎在scipy.optimize.root中可用的所有方法都以某种方式近似于Jacobian，而且我也找不到使用中讨论的API的方法。尽管如此，使用内部API，我还是成功地使用了牛顿的方法，并编写了以下代码：x, info = nonlin_solve(f
浏览 3提问于2021-07-08得票数 2
我试图在python中求解下面的多项式。其中一些常量(2.11和.125567481)是动态的，并将在程序中更改。我需要一种方法来计算这个方程在python中的解。我试过使用Sympy，但它一直在执行，但始终没有想出解决方案。Wolfram给出了4种解决方案，其中3种是复杂的，1种是真实的。我只需要真正的那个。在python中有什么方法可以有效地求解这个方程吗？在mathematica中很容易做到，但我无法找到一种在python3中高效</em
浏览 4提问于2017-10-23得票数 2
回答已采纳
Easy() = default;
else if( f_jacobian.empty() ) f_jacobian.create(1, 2, CV_64F);
Mat jacobian = f_jacobian</e
浏览 32提问于2022-04-12得票数 0
回答已采纳
问题是如何使用该函数。#print jacobian(f,np.array[0.])
/drake/lib/python2.7/site-packages/pydrake&
浏览 39提问于2019-12-22得票数 0
回答已采纳
我正在使用本征：：LevenbergMarquardt求解器作为模型拟合应用程序。我提供给它的函子包括计算误差向量和雅可比的函数。这些函数包含大量类似的代码，一些昂贵的计算重复进行。用于计算误差向量的()操作符的原型包含一个指向Jacobian矩阵的可选指针。如果特征：：LevenbergMarquardt求解器可以同时计算误差向量和Jacobian，那么在相关情况下，它将确实加快了我的算法。operator()(const Eigen::VectorXf& z, Eigen::VectorXf&
浏览 4提问于2022-06-23得票数 0
我试图用牛顿法求解一些非线性系统，解的精度对我的问题很重要。
- an algorithm suitable for general nonlinear systemcases
浏览 4提问于2013-11-12得票数 0
回答已采纳
我对将旧的Fortran代码重写成Python很感兴趣。该代码用于求解任何一般的场变量，称为F (速度、温度、压力等)。但要求解每个变量，我们必须将该变量的EQUIVALENCE定义为F。例如，如下所示：上面的概念有没有Python版本？
浏览 27提问于2019-11-27得票数 1
我试图通过package scikits.odes ()，使用Sundials (scikits.odes)的求解器IDA来求解DAE系统(2 ODE和1代数方程)。我使用的是scikits.odes 2.4.0、Sundials 3.1.1和Python3.664位。以下是代码：from scikits.odes.dae import dae
extra_op
浏览 5提问于2019-06-24得票数 0
回答已采纳
1回答
我以前唯一接触DAE的时间是整行为0 --但在这种情况下，我的质量矩阵并不总是奇异的。
在数字上实现这一点的最佳方法是什么？到目前为止，我已经尝试使用Python，在M的主要对角线上添加一个小的数(0.0001)，并将其反演，求解方程y‘= M^{-1}(t，y) F(t，y)。Python没有任何内置函数来处理质量矩阵，所以我也尝试用Julia编写这个函数。然而，DifferentialEquations.jl 认为“不直接支持非恒定质量矩阵:建议用户通过用常数质量矩阵替换DAE来将问
浏览 8提问于2021-02-08得票数 2
回答已采纳
1回答
下面是引用dll中的函数stl_example_model和stl_example_jacobian的R代码。= getNativeSymbolInfo(name = "example_jacobian_stl",我已经使用Rcpp::sourceCpp获得了另一个函数(test.cpp)。如何使用Rcpp引用此函数(类似
浏览 2提问于2019-11-15得票数 2
回答已采纳
我正在使用python ODE来求解Van der Pol方程。我想要做的是找出在一次调用integrate(t)方法期间，ODE的派生函数被调用的次数。我使用的求解器是isoda。我试图找出在没有jacbian函数的一次调用期间，导数被调用的次数，以及当我包含jacobian函数时，jacobian函数被调用的次数。callsydot)
r2 = ode(f2.ydot, f2.jac)
r2.set
浏览 0提问于2015-03-19得票数 0
我正在生成smt2文件，并使用parseSMTLIB2File和solver.check()来解析和求解约束。做其他C++/C/Python执行与parse
浏览 0提问于2019-01-22得票数 1
我想在python中实现伴随灵敏度分析，以便确定我的目标函数相对于某些参数的梯度。具体地说，目标函数依赖于微分方程的解，而微分方程的解又取决于这些参数，我正在寻找这些参数的最佳值。为了做到这一点，Julia (see here)中有许多很好的包，以及来自SUNDIALS的CVODES，然而后者显然有一个为python制作的包装器，根据这个link不包括敏感性分析功能。因此，我的问题是，python中是否存在具有敏感度分析功能的CVODES版本，或者是否有其他包可以用于执行伴随敏感度分析？
浏览 48提问于2020-12-15得票数 1
回答已采纳
1回答
我想求解系统A.b=x，其中A是，几乎是，是python中的三对角矩阵：a b 0 0 .... 0 0 b0 b a b .... 0 0 0我可以使用numpy解决程序来解决和集成我的系统：这是可行的，但非常慢，因为我的矩阵是巨大的，我不认为它利用了稀疏和接近三对角的A阵列。如果它是一个纯三对角系统，我知道如何使用经典的正反向替换算法快速高效地求解，但
浏览 3提问于2020-06-10得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
