ODE函数中的非负选项
是指在求解常微分方程（Ordinary Differential Equation，简称ODE）时，限制解的取值范围为非负数。常微分方程是描述自然现象中变化规律的数学模型，常见于物理学、工程学、生物学等领域。
在某些实际问题中，解的非负性是非常重要的，例如人口增长模型、化学反应速率方程等。为了保证解的物理意义和合理性，需要将解限制在非负范围内。
在求解ODE函数时，可以通过引入非负选项来约束解的取值范围。常见的非负选项有以下几种：
- 非负初值条件：在ODE问题中，给定初值条件时，可以要求解在初始时刻的取值为非负数。这可以通过在求解过程中添加约束条件或者选择适当的数值方法来实现。
- 非负约束条件：在某些问题中，除了初值条件外，还需要满足一定的非负约束条件。例如，化学反应速率方程中的物质浓度必须为非负数。这可以通过在ODE函数中引入非负约束条件来实现。
- 非负解的物理意义：在某些问题中，解的非负性是问题本身的物理意义所要求的。例如，人口增长模型中的人口数量必须为非负数。这可以通过选择合适的ODE模型和参数来保证解的非负性。
腾讯云相关产品和产品介绍链接地址：
腾讯云提供了丰富的云计算服务和解决方案，包括计算、存储、数据库、人工智能等领域。以下是一些与ODE函数求解相关的腾讯云产品：
- 云服务器（Elastic Compute Cloud，简称CVM）：提供弹性计算能力，可用于部署和运行求解ODE函数的计算任务。产品介绍链接：https://cloud.tencent.com/product/cvm
- 云数据库MySQL版（TencentDB for MySQL）：提供高性能、可扩展的关系型数据库服务，可用于存储和管理与ODE函数求解相关的数据。产品介绍链接：https://cloud.tencent.com/product/cdb_mysql
- 人工智能平台（AI Platform）：提供丰富的人工智能服务和工具，可用于开发和部署与ODE函数求解相关的机器学习和深度学习模型。产品介绍链接：https://cloud.tencent.com/product/ai
请注意，以上仅为示例，腾讯云还提供其他与云计算和云服务相关的产品和解决方案，具体可根据实际需求进行选择和使用。
相关·内容
我有一个带有句柄@func的函数，它的初始条件为y0，我想在t = [0, tspan]上进行测试。确切地说，我需要做什么来增加ode45使用的时间步骤的数量，而不改变tspan
我看到了ode45的MATLAB文档，并看到我需要以某种方式将options输入更改为ode45。但是，我并不真正理解如何做到这一点，因为没有提供示例代码。
浏览 11提问于2015-11-04得票数 2
回答已采纳
我有一个长度为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
回答已采纳
我试图在每个空间网格单元上应用一个常微分方程组(ODE)。因此，每个景观单元都有一个相关的ODE模型。易感蚊子(Sv)、暴露蚊子(Se)和感染蚊子(St)的数量在每一时间步骤上都被更新，ODE模型与一个基于离散agent的动物运动模型相耦合。下面是一个运行ODE模型的示例：
library(deSolve)
mod1 <- function(out_tab, time_step, var){
Sv <- out_tab[time_step,c("Sv")]
Ev <- out_tab[time_step,c("Ev")]
Iv &l
浏览 0提问于2018-04-17得票数 6
回答已采纳
2回答
这首歌怎么解？
、、、
无论我是用odeint还是solve_ivp来解决这个方程，我都遇到了一个问题。
import numpy as np
from scipy.integrate import solve_ivp
import matplotlib.pyplot as plt
def ODE(E, p):
u, v = p
n = 3.25
dudE = v
dvdE = -(u**n)-(2*v/E)
return [dudE, dvdE]
P0 = [1,0]
solve = solve_ivp(ODE, (0.001,10), P0, t_eval=np.lin
浏览 3提问于2020-03-23得票数 2
回答已采纳
我试图解决一个基本的RLC电路(串行电路)的阶跃响应，基于电路中电流i(t)的ODE。电路由Heaviside(t)函数驱动(电压节流阀为1V)，因此，在ODE中，第二个成员是Dirac函数(diff(Heaviside(t),t) )。
因此，ODE如下(L是电感，R是电阻，C是电容)：
L diff(i(t), t, t) + R diff(i(t), t) + i(t)/C = Dirac(t)因为diff(Heaviside(t), t) = Dirac(t)
初始条件如下：init_cond := D(i)(0) = 1, i(0) = 0.
组件值为R = 2, L = 1, C
浏览 3提问于2021-02-08得票数 1
1回答
不幸的是，我注意到BOOST的odeint不能及时地解决一个ODEs系统，也就是说，当我改变条件时
typedef std::vector< double > state_type;
void ode_function(const state_type &x, state_type &dxdt, const double
{
dxdt[0] = x[0];
}
using namespace std;
using namespace boost::numeric::odeint;
state_type x(1);
x[0] = std::exp(1)
浏览 2提问于2016-05-26得票数 0
回答已采纳
1回答
我感兴趣的方程系统如下所示：
通过执行以下操作，我可以通过修改某人发布的一个示例来绘制它们：
import scipy as sp
import pylab as plt
import numpy as np
import scipy.integrate as spi
#Constants
c13 = 4.2
c14 = 4.2
c21 = 4.3
c32 = 4.4
c34 = 4.4
c42 = 4.4
c43 = 4.4
e12 = 1.9
e23 = 2.5
e24 = 2.2
e31 = 2.0
e41 = 2.0
#Time
t_end = 700
t_start
浏览 0提问于2015-12-11得票数 0
1回答
R中带函数工具的非负常数解？
、、、、
我正在尝试实现一个处理非负ODE系统的R算法。我需要像MATLAB中的ode45这样的东西来定义必须为非负的状态。 我在3年前就已经讨论过这个问题，但没有真正的解决方案。deSolve仍然不是我们要走的路。我发现了一些python代码，它们看起来很有前途。也许这在R中也是可能的。最后，我必须定义一个函数包装器，就像python中的functools一样。它的作用非常简单。以下是python包装器的代码： def wrap(f):
@wraps(f)
def wrapper(t, y, *args, **kwargs):
low =
浏览 7提问于2019-06-16得票数 0
我在regex方面不是很多才多艺，尤其是多行，所以我希望有人能在这里帮我。
基于下面的示例，我试图找到没有"TableRelation"-property集的代码类型的所有字段定义。
因此，在本例中，这将是"Holding Name“字段。
table 123 "MyTable"
{
fields
{
field(1000; "Created on"; Date)
{
Caption = 'Created on';
DataCl
浏览 13提问于2022-01-10得票数 2
回答已采纳
1回答
我正在尝试用ode45模拟神经元的Morris-Lecar模型。
我在初始化ode45调用时遇到了问题，文档也无法帮助我。我知道我必须通过一个函数调用ode45，然后从我的主脚本调用该函数。
一般来说，我对ODE的掌握有限，而且似乎难以理解初始化ODE45-call所需的语法。
此外，我被指示为变量'pulse‘使用时间范围，但函数中没有时间范围的输入(似乎是一个变量，而不是固定的)，该函数从主脚本接收输入，并将其与其他函数一起发送到ode45函数。提供给ode45的函数也有时间输入，但是我还是不知道如何输入时间范围。说明非常清楚，主脚本中使用的函数不接受任何时间变量。
如果您能指出我
浏览 0提问于2012-10-04得票数 0
回答已采纳
我需要对以下系统进行数值集成：
dA/dR = f(R，A，B) dB/dR = g(R，A，B)
我在为一个初值稳定问题解ODEs。在这个问题中，系统最初是稳定的，但在一定半径下却不稳定。然而，虽然稳定，但我不希望振幅从起始值(例如O(10^-5) )衰减，因为这是非物理的，因为系统的稳定性仅限于背景噪声幅度。振幅应保持在起始值1，直到系统不稳定。因此，当导数估计为负值时，我想将其改写为零。
我已经编写了一些4阶Runge代码来实现这一点，但我更愿意简单地传递ODE45 (或任何内置的求解器)一个参数，让它在它为负值时覆盖导数。这个是可能的吗？
浏览 2提问于2014-05-17得票数 2
回答已采纳
1回答
势能函数同零
、、、
我用ode45数值解太阳和某个行星的引力系统。界定了以下职能：
function eqs = prob4_4(t,in2,in3)
%PROB4_4
% EQS = PROB4_4(T,IN2,IN3)
% This function was generated by the Symbolic Math Toolbox version 5.11.
% 05-Mar-2014 15:56:04
G = in3(:,2);
M = in3(:,1);
vx = in2(3,:);
vy = in2(4,:);
x = in2(1,:);
y = in2(2,:);
t2 =
浏览 6提问于2014-03-07得票数 0
回答已采纳
我试图解决一个基于ODE的仿真问题。但我摆脱了记忆问题，有人能帮我一个小提示吗？
警告:矩阵接近单数或严重缩放。结果可能是
不准确。RCOND = 1.453788e-16。在ode23s中为207，在tr_ode为77
警告:矩阵在工作精度上是单数的。在ode23s中为207，在tr_ode为77
从记忆中消失。键入选项的帮助内存.
odenumjac (第127行) ydel = y(:，1，ny)+diag(Del)的误差；
ode23s (第224行) dfdy，Joptions.fac，nF = odenumjac(odeFcn，{t，y，odeArgs{：}，f0，Joptions
浏览 7提问于2014-01-28得票数 1
1回答
在FUSE ESB中，我尝试使用命令(install ode)安装ODE特性，但得到如下错误。任何帮助我们都将不胜感激
FuseESB:karaf@root> features:install ode
Error executing command: No feature named 'ode' with version '0.0.0' available
浏览 1提问于2013-03-18得票数 1
回答已采纳
当使用MATLAB的ode45求解器求解ODE时，当满足(一般)布尔条件时，我试图触发一个事件来终止该解决方案。例如，在下面的代码中，条件是x > 1.5 (因为事件在值为零时触发，也就是说，当x <= 1.5为false时)。从满足时间为零但随后不满足的点开始解决方案，将永远不会触发该事件。到底怎么回事？这是越界检测中的漏洞吗？我知道，当使用连续值时，过零比较容易，但我希望在这种情况下，求解器仍然工作，尽管速度较慢。
f = @(t, x) 1;
tspan = [0, 10];
x0 = 1.5;
[t, x] = ode45(f, tspan, x0, odeset('
浏览 0提问于2021-08-08得票数 0
回答已采纳
在中，假设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回答
在不使用质量矩阵的情况下，像ode45这样的ode求解器可以求解y'=f(t，y)。
但是，对于涉及“质量”矩阵，M(t，y)y'=f(t，y)的问题，在ode求解器中有一个质量矩阵的选择。
“质量”矩阵究竟是什么？这个术语来自于质量弹簧系统的振荡吗？我在文档中找不到这方面的示例代码。而且，我似乎可以在y'=f(t，y)的方程中编码关于f(t，y)中的t和y的信息。在什么情况下，M(t，y)y'=f(t，y)在需要M(t，y)的情况下出现？
浏览 0提问于2016-06-26得票数 7
回答已采纳
1回答
我在使用时遇到了问题。
我希望从Python中获得ODEs的近似解(使用类似于MATLAB中的ode45函数)，但问题是，ODE近似需要一个我似乎无法从MATLAB创建的ODE函数规范。
它可以很好地从Python调用MATLAB函数，比如isprime，但是似乎没有办法在Python中指定MATLAB函数。
因此，我的问题是:是否有任何方法从Python生成MATLAB函数代码，还是从Python中指定MATLAB函数？
浏览 11提问于2016-09-18得票数 3
回答已采纳
所以我需要解x''(t) = -x(t)^p，初始条件是x(0)= 0，v(0) = x'(0) = v_o = 1，参数p的值是1。
这就是我所拥有的：
function [t, velocity, x] = ode_oscilation(p)
y=[0;0;0];
% transform system to the canonical form
function y = oscilation_equation(x,p)
y=zeros(2,1);
y(1)=y(2);
y(2)=-(x)^p;
浏览 5提问于2012-11-05得票数 0
回答已采纳
Scikit-learn提供了大量有用的线性模型。然而，我正在处理一个具有非负性约束的线性问题(即解决方案变量应该是非负的)。
我想使用scikit learn，但我看到的唯一可以支持非负性的函数是Lasso()，通过'positive=true‘标志。LassoCV()或lasso_path()不支持非负性。
在scikit-learn中还有其他非负线性模型吗？是否有计划在不久的将来包括这些内容？
浏览 0提问于2013-01-15得票数 3
回答已采纳
1回答
我在上遇到了一个奇怪的问题。下面是一个最低限度的工作示例：
import sys
import time
import numpy as np
import matplotlib.pyplot as plt
import matplotlib
from scipy.integrate import ode, complex_ode
def fun(t):
return np.exp( -t**2 / 2. )
def ode_fun(t, y):
a, b = y
f = fun(t)
c = np.conjugate(b)
dt_a = -2j*f
浏览 2提问于2021-03-13得票数 0
回答已采纳
1回答
我必须在整个区间上画出二阶常微分方程y''(x)+f(x)y(x)+y(x)=0的解，初始条件是y(0)=0.1，y'(0)=0.1，其中f : [0,500] -> [0,1]， ? 我使用了函数 function dwdt=equat(t,w)
f = @(x) zeros(size(x));
for i= 2:500
Bounds_1 = [i - (1/i^2),i];
Bounds_2 = [i,i + (1/i^2)];
f1 = @(x) (i^2.*x - i^3 + 1).*(Bounds_1(1) < x
浏览 28提问于2021-03-25得票数 0
我试图用use来数值解以下微分方程
x''+x=\sum_{k=1}^{20}\delta(t-k\pi), y(0)=y'(0)=0.
这是代码
from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt
from sympy import DiracDelta
def f(t):
sum = 0
for i in range(20):
sum = sum + 1.0*DiracDelta(t-(i+1)*np.pi)
浏览 2提问于2020-08-03得票数 0
回答已采纳
2回答
我正试图解决以下几个问题：
function [eta, sol] = compressible_similarity_wo
global Gamm Ma Pr omega;
Gamm = 1.4;
Ma = 2;
Pr = 0.7;
omega=0.76;
global eta_max_ode;
eta_max_ode = 20;
opt = optimset('Display','off','TolFun',1E-20);
F = fsolve(@(F) eva
浏览 4提问于2020-02-06得票数 1
回答已采纳
我对ode45有个问题。我已经定义了一个函数，并试图用ode来解决它，但当我运行它时，它花费了这么长的时间。我试着在我的函数中显示"t“输入，它显示时间步长是10^-8！我没有从ode45得到任何错误
所以我在我的函数的末尾放了一个断点，在我一步之后，它转到sym.m文件并调用函数delet(h)
function dxr=Dynfun(t,x)
...
dxr=[A;B]
在步骤之后，它将转到
function delete(h)
if builtin('numel',h)==1 && inmem('-isloaded',
浏览 3提问于2012-09-18得票数 1
1回答
是否有可能使用ODE求解器，例如ode45，并且仍然能够“更改”被调用函数中的参数值？
例如，如果我使用以下函数：
function y = thisode(t, Ic)
% example derivative function
% parameters
a = .05;
b = .005;
c = .0005;
d = .00005;
% state variables
R = Ic(1);
T = Ic(2);
y = [(R*a)-(T/b)+d; (d*R)-(c*T)];
使用此脚本：
clear all
% Initial conditions
It = [0.06 0.0
浏览 5提问于2012-10-16得票数 0
回答已采纳
1回答
我正在学习sympy，并希望验证ODE的解决方案。我还不太明白同情的命名规则。
我不想执行在顶部加载所有包的标准方法，我只想使用import sympy，然后使用显式的长名称来引用所有其他名称。关于最新的conda python
Python 3.7.3 (default, Mar 27 2019, 22:11:17)
[GCC 7.3.0] :: Anaconda, Inc. on linux
打字时
import sympy
x = sympy.symbols('x')
y = sympy.Function('y')
ode = sympy.Eq(
浏览 0提问于2019-07-15得票数 2
回答已采纳
我想在R中解一个微分方程(用deSolve?)对此，我没有初始条件，而只有状态变量的最后条件。这是如何做到的呢？
典型的代码是：ode(times, y, parameters, function ...)，其中y是初始条件，function定义了微分方程。
浏览 9提问于2016-11-14得票数 2
回答已采纳
定义像下面这样的微分方程在MATLAB中没有问题，我可以用ODE45函数来解决它。
# example.m
x = pi / 2;
x_span=[0 pi/2];
ic=[0 1];
[X OUT] = ode45(@fun,x_span,ic)
#fun.m
function out=fun(x,s)
y1 = s(1)
y2 = s(2)
out=[y2;x*y1];
但是，如果我添加基于x的新函数，则ODE45无法工作。
# example.m
x = pi / 2;
A1 = sin(x);
A2 = sin(x)+1;
x_span=[0 pi/2];
ic=[0 1];
[X
浏览 2提问于2015-07-10得票数 1
回答已采纳
让我们假设，对于一个一阶ODEs系统，我们有以下求解器：
% func.m
function dydt = func(t,y)
dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];
主要代码是：
% solver.m
tspan=0:1:10;
[t,y] = ode45(@func,tspan,[2; 0]);
如何在不等待整个代码完成的情况下，实时显示y1(t)和y2(t)对ode45所做的每个时间步骤t (t=0,1,2，...,10)的结果？
浏览 0提问于2016-10-20得票数 4
回答已采纳
我试着用ODE45来求解一个非线性微分方程组，我成功地做了很多次，但是这一次我得到了下面的错误，我真的不知道是什么错了，我完全迷惑了。以下是代码。
%%这是错误：
Subscript indices must either be real positive integers or logicals.
Error in non_L_ss (line 6)
(-Fk*(ds0+x(3)-x(1))+Fk*ds0-Fc(x(4)-x(2)))/ms +Fa/ms ] ;
Error in odearguments (line 87)
f0 = feval(ode,t0,y0,ar
浏览 8提问于2017-07-07得票数 0
在将ode与and集成时，ode接受比t和y具有更多参数的函数。
def fun(t, y, param1, param2):
并且可以使用set_f_params方法设置这些参数的值。
但是，当使用set_solout方法并试图在此函数中使用set_f_params更新params时，集成仍然保持不变，就好像没有修改params一样。
您将如何使用sol_out修改参数？--我希望从dopri5密集输出中获益，但我需要在每个时间步骤更新非齐次项。
下面是一个最小的例子。
import numpy as np
import matplotlib.pyplot as plt
from scipy
浏览 2提问于2017-08-15得票数 2
回答已采纳
2回答
我有一个通用函数，它定义了我计划使用集成的一种ODE形式，例如：
def my_ode(K, tau, y, u):
return K*u/tau - y/tau # dydt
我的代码中有几个对象，它们都具有在my_ode中定义的表单的动态性，但是具有唯一的参数K和tau。我希望能够向my_ode传递一个惟一的句柄，在初始化对象时已经设置了这些参数，所以当我更新对象时，我所要做的就是在一些模拟时间( t )中做一些类似soln = odeint(my_ode, t, y, u)的事情。
例如，如果我定义了一个类：
class MyThing:
def __init__(self,
浏览 3提问于2016-01-14得票数 7
回答已采纳
函数dy = g2(x，y)
dy = -0.1 * y；
ym = ode45('g2'，0，5，4)
结束
我收到以下信息：
g2(0.5,4)错误使用odearguments (第83行) tspan中的最后一个条目必须与第一个条目不同。ode45 (第115行)中的错误(FcnHandlesUsed，solver_name，ode，tspan，y0，options，varargin)；g2中的错误(第9行) ym = ode45('g2'，0，5，4)；
我可以补充说，这个方法很好用：
函数dy = g1(x，y)
dy =3*x^ 2；
ym =
浏览 3提问于2022-03-07得票数 1
首先，我对BPEL非常陌生。我有一个托管在Tomcat7.0中的webservice和一个部署在Apache ODE中的BPEL。我已经使用Eclipse JUNO完成了所有这些工作。现在，当我无法从eclipse调用web服务时，因为Apache ODE无法在Tomcat实例运行时启动。此外，我不能在Apache ODE中托管web服务，也不能在Tomcat中托管BPEL流程。现在如何调用Web服务？
浏览 6提问于2012-10-03得票数 1
回答已采纳
1回答
我正在使用R包解决deSolve中的ODE问题。为了加快计算速度，我想使用编译后的代码，使用指令
下面是一个ODE系统的例子，下面是我使用Rcpp编写的代码。ODE系统的详细信息取自一个示例MATLAB代码(可以找到)。我想模拟一组非平凡的ODE，以便只看到R和编译代码之间的差异。下面是我的驱动器文件，其中我以两种不同的方式计算质量平衡(只有R和编译的代码)。
library(deSolve)
library(ggplot2)
library(microbenchmark)
source('parameters_gprotein.R')
p <- parameter
浏览 5提问于2016-04-06得票数 4
回答已采纳
1回答
我编写了以下代码，以查看t my "exponential_decay“跨越零行的代码。我是Brent Method。
odr, hr, dr, cr, m = np.genfromtxt('data.txt',unpack=True)
n=0
with open('RDE_nob_trans.txt', 'w') as d:
for i in range(len(dr)):
c = cr[i]
initp = dr[i]
exponential_decay = lambda t, y: -(1/(
浏览 1提问于2019-05-09得票数 0
回答已采纳
我目前正在尝试使用SciPy的integrate.ode包来解决一对耦合的一阶ODE:比方说，。然而，这意味着在集成循环期间，我必须在每次迭代时更新发送给方法的参数，而简单地跟踪前一个值并在每次迭代中调用set_f_params()似乎并不管用。
hprev = Ho
pprev = Po
yh = np.zeros(0)
yp = np.zeros(0)
while dh.successful() and dp.successful() and dp.t < endtime and dh.t < endtime:
hparams = [alpha, beta, pprev
浏览 0提问于2011-04-20得票数 6
回答已采纳
我们有一个10个ode的耦合系统。耦合出现在最后一个方程中。 我考虑使用一个10乘2的矩阵作为初始条件。我在这里也使用了一个具有相同标题的类似问题，但我仍然得到了相同的错误(“输入参数太多”)。 time = [0 5];
x1_0 = [0 0 0 0 0 0 0 0 0 0];
x2_0 = [0 0 0 0 0 0 0 0 0 0];
initial = [x1_0;x2_0];
x = NaN(length(initial),2*length(time));
[t,x] = ode45(@ode,time,initial);
function [dxdt] = ode(x)
浏览 27提问于2019-09-09得票数 0
与MATLAB类似，ODE求解器具有Outputfcn，其中ODE求解器在每个成功的时间步之后调用函数
options=odeset('OutputFcn',@odeprint)
[T,Y]=ode15s(@(T,X)sys(T,X,vin),[t0 .0005],X(:,1),options)
在每个时间步之后调用的scipy的solve_ivp中是否有任何等效的输出函数来存储解向量？或者，有没有办法在scipy的solve_ivp中创建这个输出函数？
对于我的代码:将A，R定义为4*4，将B，S定义为4*1 Z as4*1`
def func(t,Z):
浏览 0提问于2020-04-28得票数 0
1回答
我在两个.m文件中用两个函数创建了ode。我正在尝试将来自这两个函数的ode组合起来，并同时运行它。
以下是功能，
主要功能：运行玩具-使用ODE45两次解决ode的存在，分别在两个函数中分别命名为toy和yot。
function runtoy()
%% toy
Nnodes = 5;
vc = 50;
vz0 = [280 100 150 180 110];
vz = vz0(2:4)';
tspan = [0 2];
[tz,dvzdt] = ode45(@(t,y) toy(t,y,vc), tspan,vz)
plot(tz,dvzdt)
hold all
plot(tz,
浏览 0提问于2018-10-30得票数 0
在中，
据说如果要传递参数u，需要使用匿名函数：
u = 1.2;
[t y] = ode45(@(t, y) ypdiff(t, y, u), [to tf], yo);
最初，在不传递参数u的情况下，ode行为：
[t y] = ode45(@ypdiff, [to tf], yo);，其中@ypdiff只创建一个函数句柄。
如果我们只想传递u，为什么我们还需要在匿名函数@(t, y) ypdiff(t, y, u)的创建中包含t和y，而不是像@ypdiff(u)这样的东西
浏览 0提问于2016-10-02得票数 0
我是Matlab的新手。我希望你能帮助我。我必须解决一个使用ODE45函数的ODE系统。这是描述我的等式的函数。
function dNdt = rateEquations(t, y)
%populations of corresponding state
Ng = y(1);
Ns = y(2);
Nt = y(3);
%All constants used are dropped for the sake of easy reading.
请注意参数F。
%rate equations
dNs = s0 * Ng * F - Ns/ t_S1;
dNt
浏览 2提问于2011-10-07得票数 6
回答已采纳
我正在尝试访问ode45提供的当前解决方案。我见过的所有示例，它们都求解ode，然后存储解决方案以供进一步操作。虽然这是另一种解决方案，但我需要在每次迭代时访问该解决方案，以便为动画目的绘制它们。这是简单摆的代码。 clear
clc
theta0 = [pi/2; 0];
tspan = [0 10];
[t, theta] = ode45(@odeFun, tspan, theta0);
plot(t, theta(:,1));
function dydx = odeFun(t, theta)
g = 9.8;
l = 1;
dydx = zeros(2,
浏览 19提问于2021-07-29得票数 0
回答已采纳
我试着做这样的事。
class ODESolver
{
public:
// Fourth-order Runge-Kutta ODE solver.
static void RungeKutta4(ODE ode, double stepSize)
{
// Define some convenience variables to make the
// code more readable
int j;
int numEqns = ode.getNumberOfEquations();
浏览 5提问于2014-10-07得票数 0
1回答
我正试图解一个二阶微分方程。我已经设法解决了一个角度(函数中有一个因素取决于一个角度)使用ode45在matlab中。在我的程序脚本中，也有与ode45的输入参数相同的角度。
我想要的是求解不同角度的方程。我用一个for -循环进行了测试，我只需在ode45中运行给定的角度(作为输入参数)，这样就可以绘制出不同的曲线。但是，我不太确定它是否正确，因为我没有改变ode45解决的函数的角度。
所以问题是，对于每个循环，我是如何改变函数内部的角度的。
浏览 3提问于2014-04-18得票数 0
回答已采纳
我有以下颂歌：
x_dot = 3*x.^0.5-2*x.^1.5 % (Equation 1)
我正在使用ode45来解决它。我的解是dim(Kx1)的向量(通常k= 41，由tspan给出)。
另一方面，我从(1)建立了一个近似于模型的模型，但是为了比较第二个模型的精度，我想用ode45来求解它(求解第二个ODE)。我的问题是，第二个ode是离散的：
x_dot = f(x) % (Equation 2)
f是离散的，不是像(1)中那样的连续函数。我为f提供的值是：
0.5644
0.6473
0.7258
0.7999
0.8697
0.9353
0.9967
1.0540
1.107
浏览 3提问于2013-11-01得票数 1
回答已采纳
我只想问一些关于我的code..got的问题--这个来自互联网的代码--这可能是mathworks网站上的一个例子。我在MATLAB中对这个代码进行了仿真，得到了与串行计算相比，使用并行计算(由2名工作人员组成的匹配池)求解ODE所花费的时间更少的结果。
我的问题是，能否有人向我解释一下并行性是如何解决matlab..below中的ODE方程的，这是我在matlab中使用的并行计算的代码。
此代码是编译计算代码和显示代码的主要代码。
% Main Coding
% Initialize the k and b ranges.
%
bVals = 0.1 : 0.05 : 5;
kV
浏览 0提问于2013-12-25得票数 2
回答已采纳
我正在尝试使用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
回答已采纳
我试图使用以下代码：
ode1='D2y1=-1256.4*y1-5*Dy1+255.1*y2+182.781';
ode2='D2y2=-151.5*y2-5*Dy2+255.1*y1-14.0459';
CI='y1(0)=2,y2(0)=-2,Dy1(0)=0,Dy2(0)=0';
sol=dsolve(ode1,ode2,CI,'t');
sol.y1
sol.y2
matlab返回一个错误：
?？?使用==> mupadengine.mupadengine>mupadengine.feval在
浏览 2提问于2015-04-05得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
