为什么我只能在Octave的ODE45模拟中得到NaN och Inf？
在Octave的ODE45模拟中得到NaN和Inf的原因可能有以下几个方面：
- 初始条件问题：ODE45模拟需要提供合适的初始条件，包括函数的初始值和导数的初始值。如果初始条件设置不当，可能会导致数值计算过程中出现NaN和Inf。
- 函数定义问题：ODE45模拟需要提供一个函数来描述系统的微分方程。如果函数定义有误，比如除以零或者计算结果超出了数值范围，就可能导致NaN和Inf的出现。
- 数值稳定性问题：ODE45模拟使用了一种自适应步长的数值积分算法，但在某些情况下，系统的微分方程可能具有数值不稳定性，导致数值计算过程中出现NaN和Inf。
针对这个问题，可以尝试以下方法来解决：
- 检查初始条件：确保初始条件设置合理，并且满足数值计算的要求。
- 检查函数定义：仔细检查函数定义，确保没有除以零或者计算结果超出数值范围的情况。
- 调整数值积分算法：尝试使用其他数值积分算法，例如ODE23或者ODE15s，看是否能够避免NaN和Inf的出现。
- 调整数值参数：尝试调整数值积分算法的参数，例如步长大小、容错限制等，看是否能够改善数值计算的稳定性。
需要注意的是，以上方法仅供参考，具体解决方法可能需要根据具体情况进行调试和优化。
关于ODE45模拟和Octave的更多信息，您可以参考腾讯云的数学建模服务（https://cloud.tencent.com/product/mms）和Octave的官方文档（https://www.gnu.org/software/octave/doc/interpreter/Ordinary-Differential-Equations.html）。
相关·内容
我在网上读到这篇文章，其中有人提到，使用"if语句“和"abs()”函数可能会对MATLAB的可变步长常微分方程求解器(如ODE45)产生负面影响。根据OP，它会显著影响时间步长(对时间步长的要求太低)，并在微分方程最终积分时产生较差的结果。我想知道这是不是真的，如果是真的，为什么。此外，如何在不求助于固定步骤解算器的情况下缓解此问题。我给出了下面的示例代码来说明我的意思：
function [Z,Y] = sauters(We,Re,rhos,nu_G,Uinj,Dinj,theta,ts,SMDs0,Uzs0,...
Uts0,Vzs0,zspan,K)
Y0 = [SM
浏览 2提问于2013-05-16得票数 4
回答已采纳
我正在用ode15在MATLAB中求解一个刚性偏微分方程，它通常会根据初始条件而冻结。我实际上从来没有收到错误，只是它不会结束，即使在10个小时后，它应该需要大约30秒才能运行。我正在尝试不同的空间和时间节点间隔，但这很难，因为我得不到反馈。
对于fsolve，是否存在某种等价于诊断的功能？stats没有什么用处，因为它只在fsolve完成后才显示输出。
浏览 0提问于2013-06-13得票数 0
我在matlab中用ode45和ode15s写了一个解决系统的代码。我想知道我是否可以在我的脚本中使用多核(或并行代码)来提高代码的速度。
有人试过这个吗?？
谢谢
浏览 3提问于2013-05-14得票数 5
回答已采纳
1回答
我有一个关于一阶微分方程的问题。我正在尝试用ode23和ode23s来解决它们。微分方程是：
y'+3y+z=0
z'-y+z=0
初始值: y(0)=1和z(0)=1
我还想将其与解析解进行比较：
y=exp(-2x)(1-2x)
z=exp(-2x)(1+2x)
我之所以这样做，是因为我需要进行比较，以便选择更好的求解器：ode23或ode23s，哪个更接近解析解。
我的代码是：
function dy=projectb1(t,y)
%y'=-4y
%z'= 0
%y(1)=y'
%y(2)=z'
dy = [-4*y(2); 0*y(
浏览 3提问于2015-05-28得票数 0
在我的模型中，每个智能体在每个节拍上求解一个ODE系统。我采用了欧拉法(类似于NetLogo中的系统动力学模型)来求解这些一阶常微分方程。然而，为了获得稳定的解决方案，我被迫使用非常小的时间步长(dt)，这意味着使用这种方法进行模拟的速度非常慢。我很好奇，有没有人有建议可以更快地解决这些常微分方程？我正在考虑实现Runge-Kutta (具有更大的时间步长？)正如这里所做的()。我也会考虑使用R扩展和在R中使用ODE求解器，但同样，ODE是由每个代理求解的，所以我不知道这是不是一种有效的方法。
我希望有人对这些方法的性能有所了解，并能提供一些建议。如果没有，我会尝试分享我的发现。
浏览 0提问于2015-04-10得票数 1
3回答
我这里有一个齿轮系统，我需要用MATLAB和Simulink进行仿真。
要做到这一点，首先我需要得到这个旋转(齿轮)系统的微分方程。我认为成功地做到了，如下所示
iT1 = J * theta2''+D * theta2'+K * theta2,
其中J = i^2 * Jn1 + Jn2 + J2 (同时，i是齿轮比；i = theta1/theta2 = N2/N1 = Tn2 / Tn1 = T2/T1)。
现在我需要用MATLAB中的ode23来求解这个二阶微分方程。问题是，我不知道任何初始条件，也不知道如何从这个特定的系统中得到它们。例如，y(0)
浏览 6提问于2021-02-01得票数 0
1回答
我正在使用以下命令来求解ODE函数"DynamicFunc"：
x0_sim = zeros(12,1);
time_sim = [0 100];
options=odeset('AbsTol',1e-9,'RelTol',1e-6);
[tp,xp]=ode45(@DynamicFunc,time_sim,x0_sim,options);
DynamicFunc可以是任何函数，并且代码运行得很好。模拟时间为0- 100秒。决赛xp的大小为3017x12 (即有12个州和3017个试点数)。如果我想在一个持续时间为1秒的FOR循环中运行这段代码，
浏览 0提问于2013-05-22得票数 0
回答已采纳
1回答
我有一个关于integrate()的编码问题。它产生于几个步骤：
First，我有一个名为Like的初始函数。第二，得到Like的积分，并调用结果norm.cons。第三代，我将Like除以norm.cons并调用结果Like.2。终于得到了，得到了Like.2的积分。(所有R代码如下.所示)
问题：
根据定义，我以上最后一步的答案应该是"1“。但是为什么我会得到以下答案呢？4.573253e-12
这是我的R码：
Like = function(x) sapply(lapply(x, dnorm, x = seq(1, 30), 2), prod) # Initial functio
浏览 1提问于2017-05-01得票数 0
回答已采纳
2回答
我是MATLAB的新手用户。我有一些代码试图找到状态空间模型的时间历史。我想用ode45同时求解四种一阶赋值。待解方程的实质如下：
x1_dot = x2
x2_dot = -[M] * [K] * x1 - [M] * [C] * x2 + constant*[M] * [P3] * x3 + constant*[M] * [P4] * x4
x3_dot = x2 - constant*x3
x4_dot = x2 - constant*x4
其中[M]、[K]、[C]、[P3]和[P4]是3x3矩阵；x1、x2、x3、x4都是3x1向量；x1_dot等表示时间导数( 3x1向量)。
浏览 4提问于2017-01-04得票数 1
我正在尝试使用Wolfram Mathematica在本网站上给出的微分方程来模拟疾病的传播。
我进入：
NDSolve[{i'[t]== 1/2s[t]i[t]-1/3i[t], s[t]==-1/2s[t]i[t],r[t]==1/3i[t], r[0] ==0, s[0]==1, i[0]==1.27*10^-6,s'[0]==0} i, {t, 0, 100}]
并收到使用2个参数调用的错误NDSolve；需要3个或更多参数。
我也试过了
NDSolve[{i'[t]== 1/2s[t]i[t]-1/3i[t], s[t]==-1/2s[t]i[t],r[t]
浏览 2提问于2016-09-18得票数 2
4回答
微分方程解曲线的绘制
、、、
我有几个微分方程，我想为各种起始值N_0画出解决方案
以下是方程式：
dN\dt= bN^2 - aN
dN\dt = bN^2 (1 - N\K) - aN
我该怎么做呢？
我真的不关心所使用的语言。在专门的数学方面，我的电脑上有mathematica和matlab。我有权进入梅普。我必须做更多这样的事情，我希望有来自任何语言的例子，因为这将帮助我找出我想要使用的语言并学习它。
浏览 1提问于2011-07-23得票数 1
我正在设法使这段代码更快。
Nagumo1是计算时间t的两个导数的值的函数。
function x = nagumo(t, y, f)
Iapp = f(t);
e = 0.1;
F = 2/(1+exp(-5*y(1)));
n0 = 0;
x = zeros(2, 1);
z(1) = y(1) - (y(1).^3)/3 - y(2).^2 + Iapp; %z(1) = dV/dt
z(2) = e.*(F + n0 - y(2)); %z(2) = dn/dt
x = [z(1);z(2)];
end
它是一个微分方程组，它代表了一个大
浏览 0提问于2018-03-13得票数 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
回答已采纳
1回答
弹簧/阻尼器计算与绘图
、、、、
给出了两个带阻尼/弹簧的系统：
第一个系统的simulink模型，步长为2，终值为0.5：
具有相同输入的第二个系统的Simulink：
我必须使用dsolve和ode45找到代码，用Simulink生成相同的图。给定的价值是：
m1 = 500
c1 = 1200
k1 = 25000
k2 = 15000
m2 = 50
我试图找到dsolve，但它解决不了它。所以我不得不使用ode45，我完全迷路了。
第一个系统的微分方程：
syms x(t) y(t)
Dy = diff(y,t);
Dx = diff(x,t);
D2
浏览 7提问于2018-01-03得票数 0
我正在尝试在Matlab中模拟Kuramoto的隐含。我试着用ode45来解决这个问题。我还看到其他人使用Runge-kutta方法。我知道ode45使用龙格-库塔方法，然而，我从每个方法中获得的值有可疑的不同。
kuramoto= @(x,K,N,Omega)Omega+(K/N)*sum(sin(x*ones(1,N)-(ones(N,1)*x')))'
%Kuramoto is a model of N coupled ocilators (such as multiple radiowaves)
%The solution to the model is the pha
浏览 1提问于2018-08-29得票数 1
2回答
我正在将MATLAB代码翻译成Python，但在担心翻译之前，我想先了解MATLAB，特别是它的ODE15s求解器是如何解释方程的。
我有一个函数脚本，它在主脚本中被调用，这个函数脚本包含以下等式：
function testFun=testFunction(t,f,dmat,releasevec)
testFun=(dmat*f)+(releasevec.');
在testFunction中，t表示时间，f表示我要求解的值，dmat表示我感兴趣的常量矩阵，refers表示附加常量的向量。
主脚本中的ODE15s解算器使用以下行实现其魔力：
for i=1:1461
浏览 2提问于2015-07-10得票数 0
1回答
我试图用八度法求解一个微分方程，用欧拉法。
给出了Euler方法(并且是正确的)，它适用于给定的初值问题，
y*y'' + (y')^2 + 1 = 0; y(1) = 1;
初值问题在以下八度函数中定义：
function [YDOT] = f(t, Y)
YDOT(1) = Y(2);
YDOT(2) = -(1 + Y(2)^2)/Y(1);
我的问题是关于这个函数的定义。为什么YDOT(1) = 1？什么是Y(2)？我没有找到任何关于使用函数YDOT来定义函数的文档，而不是简单的函数YDOT，我希望对Octave代码所做的事情有任何澄清。
浏览 2提问于2014-03-21得票数 2
回答已采纳
我正在使用ode15s来模拟/求解一组代码。我想实现一个特性，在模拟过程中达到给定条件时，模型中的一个数字在固定的时间内按程序(例如，指示常数)变化，然后返回。
例如，可以使用方程：
dx/dt = alpha_x - beta_x*y
dy/dt = (delta+indicator)_x_y - gamma_y + epsilon_indicator
指示器从0开始。假设x达到10，我想把10倍单位的指示器切换到1，然后把它翻转回0。
这可以通过使用全局变量以一种肮脏的方式完成，但是，这是我想避免的事情(不可能并行+一般避免全局变量)。在使用ode15s时，是否有一种简洁的替代方法(即，我
浏览 0提问于2018-11-18得票数 1
回答已采纳
我使用下面的np.lib.stride_tricks.as_strided生成滑动窗口
wsize=4
overlap=0
vector=np.array(range(31))
fillval=np.nan
part_to_fill=np.full(wsize - (vector.shape[0] - 1) % wsize - 1,fillval)
a_ext = np.concatenate(( vector,part_to_fill))
n = a_ext.strides[0]
strided = np.lib.stride_tricks.as_strided
res=stride
浏览 1提问于2018-01-26得票数 0
回答已采纳
我有一个非常小的数字除以非常小的数字，MATLAB给出了Inf (当x太小时)或NaN (当x为零)。但是正确的限值应该是2。如何解决这个问题？为太小的j或x = 0获得2的x = 0
X= 0.000001
>> x = 0.000001
x =
1.0000e-06
>> j = 2*x*(1+(-1)*exp(-x))^(-1)*exp(-x)
j =
2.0000
X= 0.01
>> x = 0.01
x =
0.0100
>> j = 2*x*(1+(-1)*exp(-x))^(-1)*exp(-x)
j =
浏览 5提问于2016-07-30得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
