Python: TypeError：'function‘对象没有属性'__getitem__’。无法将函数的元素传递到python中的ode求解器
在Python中，当我们尝试将函数作为参数传递给ode求解器时，可能会遇到TypeError: 'function'对象没有属性'getitem'的错误。这个错误通常发生在我们错误地将函数名作为参数传递给求解器的某个方法时。
要解决这个问题，我们需要确保正确地将函数作为参数传递给ode求解器。下面是一些可能导致此错误的常见原因和解决方法：
- 检查函数名是否正确：首先，确保函数名正确无误。检查是否拼写错误或者函数名是否与定义时一致。
- 确保函数被正确定义：确保函数被正确定义，并且具有正确的参数和返回值。如果函数定义有误，可能会导致无法将其作为参数传递给ode求解器。
- 使用lambda函数：如果函数的定义没有问题，但仍然遇到此错误，可以尝试使用lambda函数来包装函数。lambda函数是一种匿名函数，可以将其作为参数传递给ode求解器。
下面是一个示例，展示了如何使用lambda函数来解决此错误：
from scipy.integrate import ode
def my_function(t, y):
# 函数定义
return y
# 使用lambda函数将my_function作为参数传递给ode求解器
solver = ode(lambda t, y: my_function(t, y))在这个示例中，我们使用lambda函数将my_function作为参数传递给ode求解器。lambda函数接受t和y作为参数，并调用my_function来计算结果。
需要注意的是，以上解决方法是基于Python中的ode求解器的常见用法。如果你在使用其他求解器或特定的库，可能需要查阅相关文档以了解正确的函数传递方式。
推荐的腾讯云相关产品：腾讯云函数（云原生应用开发工具），腾讯云云服务器（云服务器产品），腾讯云数据库（云数据库产品）。
腾讯云函数介绍链接：https://cloud.tencent.com/product/scf
腾讯云云服务器介绍链接：https://cloud.tencent.com/product/cvm
腾讯云数据库介绍链接：https://cloud.tencent.com/product/cdb
相关·内容
我在python的complex_ode解决程序上遇到了麻烦。
我试图解决以下问题：
dy/dt = -iAy - icos(Omegat)_B_y
其中A和B是NxN数组，未知y是Nx1数组，I是虚单位，欧米茄是参数。
这是我的密码：
import numpy as np
from scipy.integrate import ode,complex_ode
N = 3 #linear matrix dim
Omega = 1.0 #parameter
# define symmetric matrices A and B
A = np.random.ranf((N,N))
A =
浏览 8提问于2016-01-21得票数 3
回答已采纳
我正在尝试编写一个ODE求解器(类似于MATLAB中的ODE45 )。我希望这个求解器是“独立的”，即我可以传递给它我想要它求解的颂歌的名称，然后它运行那个特定微分方程的代码。
到目前为止，我有我的主类(程序的主要部分)和一个名为ODEsolver的类( ODE求解器所在的位置)。我希望能够从我的主类调用ODEsolver，并让ODE查看由字符串参数引用的方法(我想要求解的ODE的名称)
我设想的方式是：
public class Main {
double y[];
double x0,xf,y0;
x0 = 0;
xf = 10;
y0 = 1;
ODEsolve
浏览 1提问于2011-10-08得票数 0
回答已采纳
1回答
我正在运行长时间的模拟，这涉及到一个ode的集成。由于运行时间长，我决定产生后援，允许从检查点继续模拟。对于ode集成，我决定选择实例。这种方法适用于可重入的求解器“dopri5”和“vode 853”，但对于“vode”和“lsoda”则失败。对于我的特殊问题，'vode‘algo似乎是最快的，所以我感兴趣的是，是否可以绕过这个问题。我希望这是可能的，因为在重新启动时也只需要一个求解器实例。
运行以下两个代码片段将演示我的问题，请查看此
initialize.py
启动集成，并通过酸洗(使用dill，也可以对lambdas进行腌制)来保存ode解决程序实例。
import scipy
浏览 1提问于2018-01-04得票数 2
下面是一个list子类，它将其条目委托给。
from itertools import compress
class WeirdList(list):
def __getitem__(self, item):
return compress(self, item)
l = WeirdList([1, 2, 3, 4])
print(*l[0, 1, 0, 1]) # 2 4
虽然我的印象是可以直接将compress分配给__getitem__，但是上面的工作很好。
class WeirdList(list):
__getitem__ = compress
浏览 0提问于2019-02-13得票数 2
回答已采纳
我试图通过修改找到的代码来解决问题。
import sympy as sp
t = sp.symbols('t', real=True)
x = sp.Function('x', real=True)
diffeq = sp.Eq(x(t).diff(t), sp.sqrt(1 - x(t)**2))
res = sp.dsolve(diffeq, ics={x(0): 0})
math.se文章建议，分段函数满足方程，但也许我的代码没有添加所有必要的假设。我得到以下回溯：
Traceback (most recent call last):
File
浏览 11提问于2022-10-22得票数 0
回答已采纳
1回答
我在使用时遇到了问题。
我希望从Python中获得ODEs的近似解(使用类似于MATLAB中的ode45函数)，但问题是，ODE近似需要一个我似乎无法从MATLAB创建的ODE函数规范。
它可以很好地从Python调用MATLAB函数，比如isprime，但是似乎没有办法在Python中指定MATLAB函数。
因此，我的问题是:是否有任何方法从Python生成MATLAB函数代码，还是从Python中指定MATLAB函数？
浏览 11提问于2016-09-18得票数 3
回答已采纳
3回答
我正在尝试为python中的类实现__getitem__方法。我想试试类似..。
def __getitem__(self, key):
if isinstance(key, int):
return foo(key)
elif isinstance(key, slice):
return [foo(k) for k in key]
else:
raise TypeError("Invalid index.")
...but --每当我试图将它与片一起使用时，我就会得到一个'slice' object is not itera
浏览 2提问于2015-03-01得票数 0
回答已采纳
我使用数值积分一个常微分方程组(ODEs)。请看我的笔记本电脑。
与API文档中的示例代码一样，函数ode_fn(t, y, theta)定义了待解决的ODE系统。我能够接受ode_fn wrt theta的梯度，并将ODEs与tfp.math.ode.BDF集成。
但是，当我尝试接受ODE解决方案的梯度结果wrt theta时，我会得到以下错误。当我用一组简单的ODEs替换ode_fn时，代码运行时不会出现任何问题。是否应该调整解决程序设置以避免此错误？
InvalidArgumentError Traceback (most recent call
浏览 1提问于2021-11-04得票数 0
回答已采纳
我知道ode45求解器具有由Matlab程序本身控制的自适应步长。下面的描述由给出
使用两个以上的元素指定tspan不影响>求解器用于遍历从tspan(1)到tspan(end)的间隔的内部时间步骤。ODE >套件中的所有求解者都通过对基本公式的连续扩展获得输出值。尽管>一个求解器不一定精确到tspan中指定的时间点，但在指定时间点产生的>解与在内部时间点上计算的>解具有相同的精度。
但是，如果我在very_small_step中指定了tspan=[to:very_small_step:tf]，这会影响程序控制的步骤大小吗？这个力的步长会小于very_sm
浏览 0提问于2012-09-06得票数 0
回答已采纳
我目前正在使用Python来用rk方法求解函数。r8_rkf45是一个文件，它帮助绘制函数()。
import numpy as np
import matplotlib.pyplot as plt
from numpy import zeros, linspace, exp, sqrt
from rkf45 import *
from r8_rkf45 import *
def rungekutta(ode2, x0, t, n):
n = 200
z = zeros(n)
a = zeros(n)
f = ode2
neqn = 1
ab
浏览 4提问于2015-05-12得票数 0
回答已采纳
有人能解释一下这里发生了什么吗？
class Test(object):
__getitem__ = getattr
t = Test()
t['foo']
提供错误(在Python 2.7和3.1中)：
TypeError: getattr expected at least 2 arguments, got 1
鉴于：
def f(*params):
print params # or print(params) in 3.1
class Test(object):
__getitem__ = f
打印我想要的两个参数。
浏览 2提问于2010-09-13得票数 3
回答已采纳
1回答
如何使用Sympy解矩阵微分方程？
我有一个方程y'(t) = A*y(t) + B，其中A是3x3矩阵，y(t)是1x3向量，B是1x3向量。
更具体地说，我正在研究一个使用微分方程在3D中移动点的计算机图形学问题。我在3D空间y(t)中有一个点，一个3x3旋转矩阵和一个平移向量。y'=Ay+B方程是我正在研究的实际问题的简化，但我甚至无法使用Sympy来求解y'=Ay+B。我在Sympy中寻找闭合形式的解(而不是数值解)。
我知道如何解等式y'=Ay+B，但我想使用Sympy找到相同的解，然后使代码适应我试图解决的更复杂的问题。
我当前的代码是：
from
浏览 3提问于2018-08-25得票数 4
我知道如何在MATLAB中使用ode15s或其他ode求解器，但我不确定的是，这段代码(来自)对我来说似乎很模糊：
[VOI, STATES] = ode15s(@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS), tspan, INIT_STATES, options);
更具体地说，以下(？)是什么意思：
@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS)
该函数的头文件"computeRates“如下：
function [RATES, ALGEBRAIC] = computeRa
浏览 1提问于2011-03-29得票数 3
回答已采纳
我有一个从S-函数调用热力学性质计算(CoolProp)库的模型.如果它的输入没有意义，这个库会产生一个错误并突然停止仿真，这有时发生在Simulink干扰状态时，特别是当模型中存在一个Step块时，即使它的输出是浮动的！我认为这是一种预测系统状态变化的尝试，我注意到减小步长会减少问题，但我无法在10 us运行20,000 s的模拟(是的，微秒！)时间步骤，当它已经花了15分钟与5s步长！
有没有办法捕捉这个错误(我想是通过调用error()引起的)，并告诉Simulink“嗯，您使用的步长太大了”？
编辑：我使用的是ode45求解器，但是除了模拟时间急剧增加之外，即使是ode23tb或ode
浏览 3提问于2015-03-26得票数 0
现在，除了我之前的两篇文章之外，我还试图重构我的代码并修复一些问题。我决定，在逻辑上创建这样的类：Solver、Problem。
因此，ODE_Solver和FE类的代码最终将编码并工作。
# ODS.py
import numpy as np
class ODE_Solver(object):
def __init__(self, f):
if not callable(f):
raise TypeError('f is not %s, function' % type(f))
self.f =
浏览 8提问于2022-05-23得票数 0
回答已采纳
1回答
时变输入ODE求解
、、、、
我试图用ode45来求解一组ODE方程。我的一些参数已经是时间的函数，但我总是会出错。
function odo
dx(1,1) = (vl+vr)/2*cos(x(3));
dx(2,1) = (vl+vr)/2*sin(x(3));
dx(3,1) = obz
其中obz、vr和vl都是向量，例如：
obz = sin(t), t = 0:dt:tf;
我使用以下语法：
[t, x1] = ode45(@(t,x) odo(t,x,b,obz,vr,vl), 0:dt:tf, [0;0;0]);
使用R15，但继续获取错误：
赋值的非单例rhs维度比非单例rhs的下标要多。
浏览 0提问于2019-04-13得票数 0
问题如上。在阅读了文档之后，我可以更改积分器本身(RK45、RK23、DOP853等)，但是我无法找到关于这些积分器的顺序的信息，也无法找到将积分器限制为1阶的方法的信息。
这是如何做到的呢？我是否必须使用默认为1阶的特定ODE求解器方法，还是可以编辑任何一阶的方法？
浏览 3提问于2022-01-04得票数 0
回答已采纳
来自
x[i]大致相当于type(x).__getitem__(x, i)。
与看似简单的x.__getitem__(i)相比，上述方法的好处是什么？
编辑:为什么Python会这样做？
作为标准行为的一个缺点，让我展示这个示例代码，在这里我惊讶地发现最后一个断言失败了，而第二个断言(直接调用__getitem__ )通过了。
def poww_bar(base):
class Bar():
def __getitem__(self, x):
return lambda: base**x
return Bar()
def pow
浏览 4提问于2021-07-12得票数 4
我试图在ode45中使用MAtlab，并希望修复MAtlab使用的点数(时间步骤数)。在'refine'中使用ode45选项似乎没有帮助。例如，如果我将'refine'设置为10，Matlab将返回一个101数组。改变'RelTol'和'AbsTol'也没有帮助。我知道将tspan写成[0,t1,t2,t3,...,tn]是可能的，这解决了这个问题，但是我想通过'refine'选项来修正点数。
浏览 2提问于2014-01-03得票数 2
大家好，我是python的初学者，我似乎无法让函数__getitem__和__setitem__在自定义类中使用2D数组。
我一直在尝试为OOP战舰项目实现一个2D数组类。我希望能够使用坐标输入改变网格上的随机元素。
class Grid:
def __init__(self, width=10, height=10):
self.grid = [["~" for i in range(0, width)] for i in range(0, height)]
def __getitem__(self, point):
retu
浏览 0提问于2019-09-03得票数 0
1回答
我正在使用ode15s函数在MATLAB中求解一个DAE系统，在仿真中的某个时刻，求解器无法达到所需的公差，执行将终止。如果我将最大步长减小到一个小值(~1e-5)，程序就能够继续并到达指定的端点。
基于缩短时间步长允许程序继续运行的事实，我怀疑ode15s脚本在某个时候采取了太大的时间步骤，无法从它中恢复过来。然而，保持这个小的最大步骤意味着我有数千个不必要的数据点来减缓模拟，占用空间。
首先，是否有可能在ode15s中有一个可调整的最大步长？到目前为止，我发现的一切都表明这是不可能的，但我可能遗漏了一些东西。
如果这是不可能的，我可以通过等待函数通过临界时间(即仿真将在没有时间步长约束的情
浏览 1提问于2015-07-15得票数 0
回答已采纳
我目前正在使用Python和R探索Lorenz系统，并注意到ode包中的细微差别。Python和的都表示，他们使用lsoda来计算衍生品。但是，对这两个命令使用lsoda命令似乎会产生截然不同的结果。我对R中的ode函数尝试了ode45，以获得更类似于Python的东西，但我想知道为什么我不能获得完全相同的结果：
from scipy.integrate import odeint
def lorenz(x, t):
return [
10 * (x[1] - x[0]),
x[0] * (28 - x[2]) - x[1],
x[0]
浏览 49提问于2021-05-21得票数 2
回答已采纳
在中，假设progress_message的缺省值是ODE_DEFAULT_PROG_MESSAGE (第58行)，它缺省设置为显示dt、t和最大值u (第21-22行)。然而，我在TerminalLoggers中使用它来运行Jupyter Lab和/或终端中的代码，而不是，我得到消息"ETA：(某个时间号)“。
using Logging: global_logger
using TerminalLoggers: TerminalLogger
global_logger(TerminalLogger())
using OrdinaryDiffEq
function lorenz(
浏览 0提问于2021-05-20得票数 0
2回答
我正试图用Python解决一阶ODE：
其中Gamma和u是方阵。我并不是在任何时候都明确地知道u(t)，但在离散的时间步骤中，我知道它是从先前的计算中得到的。
我在网上找到的Python求解器的每一个例子(例如，表示scipy.integrate.odeint和scipy.integrate.ode)都知道作为时间函数解析的导数的表达式。
有没有一种方法叫这些(或其他微分方程解)，而不知道导数的解析表达式？
现在，我已经编写了自己的Runge解决程序，并将其与numba混合使用。
浏览 3提问于2017-11-15得票数 2
回答已采纳
1回答
首先，我知道还有其他类似的问题，事实上，其中一个有着相同的标题，但我无法从其他问题中得到任何帮助。
我正在使用功能模拟界面。我有一个功能模拟单元(FMU)，它基本上是python中一个类的实例。
我想把它传递给一个函数来得到Jacobian矩阵，但是我不能将类的一个实例传递给一个模块中定义的函数。如果我在模块中定义了这样的函数
# Function to take an FMU instance and do something using it's methods
def do_something_to_fmu(fmu):
N=fmu.get_ode_sizes() # ge
浏览 2提问于2014-06-19得票数 0
1回答
我们可以将R中的deSolve包用于常微分方程( ODE )，但是，我无法找到一种方法来求解两个嵌套的ODE方程。
b'(t) = beta - k*b(t);
a'(t) = alpha -b(t)*gamma;
其中'的意思是分化。那么，我们如何解决a和b呢？由于a'是b的一个函数，我们必须同时求解a和b。
我有个错误：
lsoda中的错误(y，time，func，parms，.)：使用的求解器组合不能嵌套。
当我试图将b的解决方案添加到a的ode解决方案中时。
浏览 8提问于2011-04-12得票数 2
回答已采纳
我可以运行：
import chef
chef.autoconfigure()
for node in chef.Node.list():
if "auto" in node.lower():
print "deleting node " + node
nodeObj = chef.Node(node)
nodeObj.delete()
直接在控制台中，但当我尝试将其作为脚本运行时：python2.7 test.py，我得到以下错误：
Traceback (most recent call last):
浏览 0提问于2013-05-31得票数 1
回答已采纳
2回答
我的问题是关于当前的赋值颂歌求解者。在中，它们的用法是：
# A problem to integrate and the corresponding jacobian:
from scipy.integrate import ode
y0, t0 = [1.0j, 2.0], 0
def f(t, y, arg1):
return [1j*arg1*y[0] + y[1], -arg1*y[1]**2]
def jac(t, y, arg1):
return [[1j*arg1, 1], [0, -arg1*2*y[1]]]
# The integration:
r =
浏览 5提问于2014-01-21得票数 2
回答已采纳
我正在使用ode45来求解位置变量对时间的影响，但是在初始条件之后，我只得到了所有4个变量的NaN响应。下面是我的函数代码：
function dxdt = ode(t,x)
global m b a x_teta k_lin k_tor Ip U Pinf gama
m = 5; % Mass
b = 1.5; % Semi-chord
a = 0; % Reference point to displacement -1 <= a <= 1
x_teta = 0; % Static unbalanced parameter -&g
浏览 1提问于2018-09-18得票数 0
我有一个python代码()，它使用scipy.integrate.ode来解决ODE系统。代码工作良好，结果是合理的。然而，我注意到了一些对我来说没有意义的关于诗歌解译器的东西。
我有一个put打印函数(print("t inside ODE function", t)) --计算导数向量(__call__(self, t, y))的函数，以及在while循环(print("t outside ODE function", solver.t);)中该函数之外的函数。我希望当求解者进行时间整合时，必须调用内部打印，然后调用外部打印。换句话说，如果没有"t
浏览 2提问于2021-05-16得票数 0
回答已采纳
我认为这个问题有点长，所以我认为最好先考虑一下简化的版本：
有两个类A和B.B继承自A。B (add)中有一个成员函数需要使用A中的一个成员函数来运行。
class A;
typedef int(A::*operation)(int c);
typedef void (A::*myfunction)(operation, int);
class A
{
public:
int a;
int b;
int do_something(myfunction f, operation op)
{
(this->*f)(op, 1);
}
浏览 0提问于2021-06-13得票数 1
回答已采纳
我正在使用Python进行关于ODE数值集成的练习。我遇到了这个问题。
在每次运行ODE积分器(即，对integrator (T)方法的一次调用)期间，积分器内部将时间间隔划分为多个离散步骤，并使用这些步骤运行指定的方案(例如Runge)。让我们考察在此过程中调用ODE导数函数的次数。
有没有找出次数的方法？谢谢。
浏览 1提问于2015-03-18得票数 1
回答已采纳
我试图用Matlab中的ode45来求解一个微分方程。我正在运行两个脚本：
function xdot=linearproblem(t,x)
global kappa mass F
xdot(1)=-(kappa/mass)*x(2)+(F/mass)*(cos(omega1*t));
xdot(2)=x(1);
xdot=xdot';
end
在第二个脚本中，我
close all
clear
clc
global kappa mass F
kappa=4;
F=2;
mass=0.5;
options=odeset('omega1',[1.4 1.5
浏览 2提问于2013-11-20得票数 0
我使用一个简单的if循环来更改ode脚本中的参数值。下面是我写的一个示例脚本，它展示了同样的问题。因此，首先是工作的版本：
function aah = al(t,x)
if (t>10000 && t<10300)
ab = [0; 150];
else
ab = [150; 0];
end
aah = [ab];
这可以使用以下命令运行
t = [0:1:10400];
x0 = [0,0];
[t,x] = ode23tb(@al, t,x0);
和可视化的
plot(t,x(:,1))
plot(t,x(:,2))
好了，这就是好的版本。
浏览 2提问于2012-11-22得票数 1
回答已采纳
1回答
我理解数值求解的逻辑，即matlab从未知值开始循环，直到方程收敛到一个值为止。然而，我不明白的是如何正确地进入我的方程式。我认为使用ode45函数是最好的方法。我有下面的等式U^(n+1) = U^n - (t'*3250/10)-(t'/2)(.004(v^n)^2/10。我怀疑我的方程需要不同的形式，但是我不确定让matlab求解方程的正确方法。
我试着在matlab中按原样输入方程，但是它抱怨v和n是未知变量，我不确定如何处理它们。求解这个方程的最终目标是找到v的值。
C=@(t,v) u^n-(3250*t'/10)-(t'/2)*((.004*(v^
浏览 1提问于2014-09-16得票数 0
代码： from unittest.mock import MagicMock, call
mm = MagicMock()
mm().foo()['bar']
print(mm.mock_calls)
print()
mm.assert_has_calls([call(), call().foo(), call().foo().__getitem__('bar')]) 输出： [call(), call().foo(), call().foo().__getitem__('bar')]
Traceback (most recent c
浏览 25提问于2019-08-24得票数 5
我目前正在尝试使用solve_ivp命令在Python语言中使用scipy.integrate来解决ODE系统。有没有办法使用decimal数据类型来集成这个系统？我需要16位以上的精度。例如，我需要28位数的输出。
浏览 30提问于2021-03-22得票数 0
我正在尝试使用ode45来解决一个ODE系统：
[X,Y]= ode45(@sys,[0, T],y0);
哪里,
function dy = sys(t,y)
dy(1) = f_1(y)
dy(2) = f_2(y)
dy(3) = f_3(y)
end
问题是，函数ode45要求y0是初始值[y_1(0), y_2(0), y_3(0)]，而在我的系统中，我只有值[y_2(0), y_3(0), y_3(T)]可用。
从数学上讲，这组初始/终止条件应该足以约束系统，但是有什么方法可以通过ode45或MATLAB中的任何其他函数来使用它吗
浏览 3提问于2013-06-14得票数 7
回答已采纳
3回答
我对“物理系统建模”这门学科很陌生。我阅读了一些基本文献，并在Modelica和Simulink/Simscape中做了一些教程。我想问你，我是否正确地理解了以下内容：
符号操作是将一个微分-代数方程组(物理模型: DAE)转化为一个微分方程组(ODE)，可由标准求解者(Runge，Kutta，BDF，.)求解的过程。
也有可以直接解决DAE问题的求解器。但是Modelica (openModelica，Dymola)和Simscape将系统转换为ODE (为什么这种方法比直接DAE求解器更好？)
“平面Modelica代码”是转换的结果(= ODE)。
非常感谢你的回答。
浏览 4提问于2016-05-20得票数 7
我正在使用非matlab的ODE模拟软件来重现一个用matlab中的simbiology工具箱创建的模型。
其中一个问题是重复赋值的表示。有没有可能用一种只支持ODE系统的标准Runge Kutta (或其他迭代方法)来模拟重复赋值？或者，将具有重复赋值的模型重新表达为ODE系统是不可能的？
浏览 0提问于2014-04-17得票数 0
1回答
我有一个切片对象，我想从它中提取第一个int。此对象类型似乎是只读/不可变的。使用python 2.7
foo = slice(22,24)
foo
Out[44]: slice(22, 24, None)
foo[0]
---------------------------------------------------------------------------
TypeError Traceback (most recent call last)
<ipython-input-45-d8976430aff2&g
浏览 1提问于2014-07-11得票数 2
回答已采纳
1回答
我试图积分一组常微分方程，并在给定的点上计算解及其一阶导数。Matlab提供了求解方程的"ode45“和计算解及其一阶导数的"deval”，但是在我的情况下，我无法使用它们。
Matlab给出了以下示例：
sol =ode45(@sol 1，0 20,2 0)；
X=直线空间(0，20,100)；
Y= deval(sol，1)；
地块(x，y)；
我猜"vdp1“是由Matlab定义的，这个例子运行得很好。但是，当我尝试使用以下代码时：
sol =ode45(‘vdp1 1’，0.20，2.0)；
X=直线空间(0，20,100)；
Y= deval(sol
浏览 4提问于2014-10-12得票数 1
回答已采纳
我有一个长度为N的向量x，我想用它的值来解微分方程：dy/dt = x - 4*y。对于每个步骤，我希望ode求解器函数使用向量的一个值，然后在下一步使用矩阵的下一个值。
为此，我尝试将向量声明为全局变量，并在ode求解器中这样使用它：
global x
tspan = 0:0.01:10;
[t,filter_coef] = ode45(@ode_filter,tspan,0);
这样的解题功能：
function dx = ode_filter(t,fil)
global x
dx = x - 4*fil(1);
end
产生跟随错误
ODE_FILTER retur
浏览 0提问于2019-05-19得票数 0
回答已采纳
我有一个带有句柄@func的函数，它的初始条件为y0，我想在t = [0, tspan]上进行测试。确切地说，我需要做什么来增加ode45使用的时间步骤的数量，而不改变tspan
我看到了ode45的MATLAB文档，并看到我需要以某种方式将options输入更改为ode45。但是，我并不真正理解如何做到这一点，因为没有提供示例代码。
浏览 11提问于2015-11-04得票数 2
回答已采纳
1回答
我想使用切换功能来集成一个简单的ODE。切换函数将确保t的某些值的梯度为零，并使梯度成为t的任何其他值的数值常数。
我可以通过使用虚拟状态来获得所需的结果，但是当我重复对单个状态的计算时，来自solve_ivp的输出是不同的。我想知道为何会这样。
下面是再现结果的代码：
#==============================================================================
# import modules
#==================================================================
浏览 3提问于2022-01-06得票数 1
回答已采纳
1回答
是否可以在不加载整个文件的情况下从hdf5文件中读取给定的行集？我有相当大的带有大量数据集的hdf5文件，下面是我为减少时间和内存使用而想到的一个例子：
#! /usr/bin/env python
import numpy as np
import h5py
infile = 'field1.87.hdf5'
f = h5py.File(infile,'r')
group = f['Data']
mdisk = group['mdisk'].value
val = 2.*pow(10.,10.)
ind = np.whe
浏览 3提问于2015-02-09得票数 4
回答已采纳
我正在使用Matlab来尝试解决一个由三个一阶常微分方程组成的系统，但是我得到的错误信息是'syms‘需要Symbolic Toolbox。爬行器中出现错误(第1行) syms f(t) s(t) v(t) r W c h Q a k b H K e 在以前的情况下，我收到一个错误，说这个ODE系统不能显式求解(即以封闭形式)。我认为数值积分是唯一的方法。r、W、c、h等都是参数。有人能告诉我如何模拟/求解和绘制下面的ODEs吗？ syms f(t) s(t) v(t) r W c h q a k b H K e
r = 1;
W = 0.5;
c = 0.4;
h = 0.9;
q
浏览 34提问于2020-02-02得票数 0
回答已采纳
2回答
我是Python初学者。我正在尝试切换matlab中的一些程序。我需要解一个刚性常微分方程，它的输入都是矩阵。在matlab中，我使用
[ttT,uT] = ode23s('SST',t,fT);
浏览 0提问于2012-01-05得票数 4
回答已采纳
2回答
我试图用MatLab和ode45函数来求解微分方程。问题是如何更改最小步长？我想把它做大一点，避免太小的步子。我试着这样做：
tspan = [3, 4]; (boundaries of time line)
[t, q] = ode45('dqdt2', tspan, q0);
T-时间
Q-解查找
dqdt2 -我的文件有函数
q0 -初始条件
浏览 5提问于2017-05-03得票数 1
回答已采纳
我想使用这样的东西：
class Board():
...
def __getitem__(self, y, x):
return self.board[y][x]
但不幸的是，当我打电话给：
board[x][y]
我明白了：TypeError: __getitem__() takes exactly 3 arguments (2 given)
浏览 0提问于2013-07-05得票数 19
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
