含lsode的倍频程微分方程组
是指一类微分方程组，其中包含了lsode函数，该函数是用于求解常微分方程组的数值解的一种方法。倍频程微分方程组是指方程组中包含了多个频率的项，即方程组中的各个方程中包含了不同频率的函数。
这类微分方程组的优势在于可以描述多个频率的振动或波动现象，适用于许多领域，如物理学、工程学、生物学等。通过求解倍频程微分方程组，可以得到系统在不同频率下的响应和行为。
在云计算领域，倍频程微分方程组的应用场景较为有限。然而，在科学计算和工程仿真领域，倍频程微分方程组的求解是一个重要的问题。例如，在电路设计中，可以使用倍频程微分方程组来描述电路中的振荡现象，通过求解方程组可以得到电路的稳定性和频率响应。
腾讯云提供了一系列的云计算产品和服务，其中包括计算、存储、网络、人工智能等方面的解决方案。然而，与具体的倍频程微分方程组相关的产品和服务在腾讯云中并没有明确的提及。因此，无法给出具体的腾讯云产品和产品介绍链接地址。
总结：含lsode的倍频程微分方程组是一类包含lsode函数的微分方程组，用于描述多个频率的振动或波动现象。在云计算领域的应用相对较少，但在科学计算和工程仿真领域具有重要意义。腾讯云作为云计算服务提供商，没有明确提供与倍频程微分方程组相关的产品和服务。
相关·内容
我有，我想解决X(t)，Y(t)，Z(t)，希望用我熟悉的Octave/Matlab，但我不介意用任何其他必要的方法来解决它。
现在，Fsolve对于正则方程组是有用的，而Ode45，Lsode对于微分方程是有用的。但是，这个特殊的系统呢？请注意，底部的微分方程不仅包含Y，还包含X和Z，它们依赖于上面的两个非微分方程。
老实说，我不太确定如何使用基本的代码来解决这个系统，在花了一些时间思考之后，我决定寻求帮助。如果可能的话，我真的很感谢任何能有效地解决这个问题的指导。现在几乎所有的回复对我来说都是有用的。
浏览 9提问于2020-02-15得票数 0
回答已采纳
我试图用octave4 .It来解一个常微分方程，.It对sin很好，而余弦函数.But用于tan，它显示了误差
# Define the right-hand side of the equation:
function ret=f(y,x);ret=tan(x) ;end;
# x will be on the interval [0,1]; y(0)=0.5
# x will be the set of moments of time:
# ywill be the values of the function at these moments of time.
y=lsode (&#
浏览 2提问于2015-09-11得票数 0
2回答
我很高兴，建模工具(在我的例子中是OpenModelica和Dymola -建模语言Modelica)是如何解决方程组(线性和/或非线性)的。这些工具是为求解微分代数方程而设计的。我对将微分代数方程系统转换为ODE (关键字“索引-约简”)的理论有一点了解。我的问题：
这些工具是如何求解没有微分方程的方程组的？然而，系统是否将(索引缩减)转换为ODE？
如果我有一个有几个代数方程和几个ODE的模型，但它们不是耦合的呢？
非常感谢。
浏览 4提问于2017-03-31得票数 3
4回答
我正在用c++编写程序来数值求解微分方程(常微分方程和偏微分方程组)。原因很简单--对于大规模计算，c++比matlab快得多。然而，当涉及到绘图时，matlab显然是赢家。有没有什么方法可以让matlab像在c++中绘图一样？
浏览 1提问于2014-05-26得票数 4
1回答
我刚刚开始使用Python绘制微分方程的数值解。我知道如何使用scipy.odeint求解和绘制单个微分方程，但我对微分方程组一无所知。如何绘制下面的耦合系统？
N' = a * N - (C/(1+C)) * b * N
C' = (C/(1+C)) * N - C + 1
a = 4
b = 7
N(0) = 100
C(0) = 5
浏览 3提问于2017-03-18得票数 5
回答已采纳
引言
利用ODE45对二阶常微分方程组进行数值求解，用Matlab对一些动态系统进行了数值模拟。我从Mathworks那里找到了一个很好的教程(最后是教程链接)，介绍了如何做到这一点。
在本教程中，方程组在x和y中是显式的，如下所示：
x''=-D(y) * x' * sqrt(x'^2 + y'^2)
y''=-D(y) * y' * sqrt(x'^2 + y'^2) + g(y)
上述两个方程均为y'‘= f(x，x'，y，y')
问题
然而，我遇到的方程组中，变量不能显式地求解，如
浏览 4提问于2017-01-24得票数 1
3回答
我对“物理系统建模”这门学科很陌生。我阅读了一些基本文献，并在Modelica和Simulink/Simscape中做了一些教程。我想问你，我是否正确地理解了以下内容：
符号操作是将一个微分-代数方程组(物理模型: DAE)转化为一个微分方程组(ODE)，可由标准求解者(Runge，Kutta，BDF，.)求解的过程。
也有可以直接解决DAE问题的求解器。但是Modelica (openModelica，Dymola)和Simscape将系统转换为ODE (为什么这种方法比直接DAE求解器更好？)
“平面Modelica代码”是转换的结果(= ODE)。
非常感谢你的回答。
浏览 4提问于2016-05-20得票数 7
我必须解Euler Bernoulli微分方程，即：
w’’’’(x) = q(x)
和边界条件：
w(0) = w(l) = 0
和
w′′(0) = w′′(l) = 0
光束如下图所示：
连续力q为2N/mm。
我必须使用射击方法和scipy.integrate.odeint()功能。
我甚至无法开始，因为我不知道如何将微分方程写成一个方程组
有谁能理解用python中的边界条件来求解微分方程，请帮助！
谢谢:)
浏览 6提问于2018-01-02得票数 2
回答已采纳
1回答
假设我有几个如下形式的微分方程 ? 其中尖括号中的变量是复数，例如 ? 我的问题是，在极大值中，是否可以先写下上面提到的形式的微分方程组，而不计算导数，然后对尖括号中的变量进行替换，然后计算导数并分离实部和虚部。
浏览 13提问于2021-07-25得票数 3
回答已采纳
1回答
我有一个微分方程组的目标解，它有一些未知的参数。我想要找到这些参数的值，它的解更接近于target.Can，我用神经网络来做这件事?如果是，怎么做？
我之所以问这个问题，是因为我正在读的一篇论文(不幸的是希腊文)意味着要做这件事。下面是下列微分方程组
想要的输出是
其中控制输入u有一些未知的非线性，用神经网络进行逼近。由于没有数据来训练一个网络，我无法理解它是如何实现的？有什么想法吗？
浏览 2提问于2014-07-13得票数 0
回答已采纳
我正在通过伴随法建立范德波尔模型的相位响应曲线。为此，有必要对原微分方程组进行求解，并将其解包含在伴随系统中。如下所示。
原系统对应于：
def vanderpol(t, v, mu):
x, y = v
dxdt = mu*(x-(x**3)/3-y)
dydt = (1/mu)*x
return (dxdt, dydt)
参数如下：
mu = 1.5
T1 = 7.09637359
omega1 = (2*np.pi)/T1
tend = 25*T1
y0 = [1, 1]
和scipy.integrate.solve_ivp函数：
sol =
浏览 13提问于2022-07-18得票数 0
如何解非线性微分方程组并得到其解的曲线图？系统是没有初始条件的。例如,
x'= (x + y)^2 -1 y'= -y^2 -x+1
浏览 1提问于2016-05-21得票数 1
1回答
我知道我可以用八度内置的ODE求解器来解决一个ODE系统。现在我想要解相同结构但系数不同的多个ODE系统。例如，我有以下的颂歌：
函数xdot=f(x，t)
a = 0.1;
xdot(1) = a * x(1) ^ 2 + x(2);
xdot(2) = x(1) - x(2);
端函数
初始条件x(1)=0，x(2)=0在t=0。我想在t=1点得到x(1)和x(2)，我知道用lsode命令可以很容易地解决这个问题。
现在我想用步骤0.1将"a“的值从0.1改为20，即a=0.1:0.1:20，因此我可以从不同的"a”中得到x(1)和x(2)。我如何在循环中做到这一点？我需
浏览 0提问于2014-04-19得票数 0
回答已采纳
1回答
我需要解下列微分方程组：
$\frac{dx_1}{dt} = -k_1x_1+k_2x_2-(K_R)x_1y_1$
$\frac{dx_2}{dt} = k_1x_1-k_2x_2-k_3x_2-(K_R)x_2y_2$
$\frac{dx_3}{dt} = k_3x_3$
$\frac{dy_1}{dt} = -k_1y_1+k_2y_2-(K_R)x_1y_1$
$\frac{dy_2}{dt} = k_1y_1-k_2y_2-k_3y_2-(K_R)x_2y_2$
$\frac{dy_3}{dt} = k_3y_3$
$\frac{dz_1}{dt} = -k_1z_1+k_
浏览 4提问于2020-03-30得票数 0
回答已采纳
1回答
嗯，我想用这样的形式来解决一个复杂的ODE系统：
$ i\hbar \dfrac{部分t} =\左H，\rho \右$
其中$\rho$和$H$是$n x n$矩阵。
我试图使用以下方法定义一个通用函数：
function xdot = f(x,t)
i*xdot(1)=
i*xdot(2)=
endfunction
x0=[0;0];
t=linspace(0,20,200);
y=lsode("f",x0,t).
但我收到了一条关于我的错误信息。那么，我的问题是:如何在八倍频程中求解复杂的微分方程？
浏览 0提问于2012-11-20得票数 1
所以我想问，是否有任何方法可以用矩阵符号来定义和求解R中的微分方程组。
我知道你通常会做这样的事情
lotka-volterra <- function(t,a,b,c,d,x,y){
dx <- ax + bxy
dy <- dxy - cy
return(list(c(dx,dy)))
}
但我想做的是
lotka-volterra <- function(t,M,v,x){
dx <- x * M%*% x + v * x
return(list(dx))
}
其中x是长度为2的向量，M是2*2矩阵，v是长度为2的向量
浏览 0提问于2018-05-01得票数 3
1回答
在了解如何在八度/matlab中实现事件时，在微分方程的解算方面，我遇到了一些困难。
例如，考虑求解微分方程y‘=-y的简单代码：
function dy = odefun(t,y)
dy = -y;
endfunction
options = odeset('RelTol',1e-4,'AbsTol',[1e-4]);
[T,Y] = ode45(@odefun,[0 12],[1],options);
现在我想介绍一个随机事件。例如，在某个固定的时间步长上，我想随机改变y的值，然后根据微分方程继续演化。我该怎么做？
浏览 2提问于2014-02-01得票数 2
回答已采纳
2回答
ODE，多根和事件，R
、、、
我喜欢解一个包含多个阈值的耦合微分方程组。通过R信息，这将导致我将ODE与根函数和事件函数结合使用。
通过各种例子，例如温度模型，第14页 --下面粘贴的代码--我可以让我的模型在一个阈值上工作，即为一个变量找到根/到达阈值触发一个事件。
library(deSolve)
yini <- c(temp = 18, heating_on=1)
temp <- function(t,y, parms) {
dy1 <- ifelse(y[2] == 1, 1.0, -0.5)
dy2 <- 0
list(c(dy1, dy2))
}
rootfunc <
浏览 2提问于2013-09-18得票数 4
我想用Mathematica 7解决这个微分方程组，但我发现了一个错误，说明函数是在没有依赖于所有自变量的情况下指定的。这些方程式是：
感谢大家的帮助
浏览 1提问于2011-12-10得票数 1
回答已采纳
我在试着画出微分方程组的解，但我的解包括初始条件的时间导数，我不知道如何把点放在变量x.Does的上面，有人知道怎么做吗？谢谢。
浏览 0提问于2019-03-26得票数 1
我用MATLAB的函数'pdepe‘来解一个偏微分方程的问题，一个抛物方程。
我需要知道函数使用的是哪种数值方法，因为我必须在报告中通知这一点。
MathWorks中函数的描述是“在一个空间变量和时间内求解抛物型和椭圆型偏微分方程组的初边值问题”。是有限差分法吗？
谢谢你帮我。
浏览 4提问于2017-04-02得票数 0
回答已采纳
这里有点新手。我试图解决一个由6个常微分方程组成的系统，但是有一个条件，如果(x(1)<=0.1)，那么x(6)=x(6)/2.0。我如何将其放入代码中？谢谢!
浏览 0提问于2015-07-04得票数 0
这就是我的问题。多亏了下面的两个函数，我正在尝试解决一个由两个微分方程组成的系统。代码中给我带来麻烦的部分是变量"rho“。" rho“是一个函数，它的值是从一个文件中给出的，我试着把它放在变量rho中。 function [xdot]=f2(x,t)
# Parameters of the equations
t=[1:1:35926];
x = dlmread('data_txt.txt');
rho=x(:,4);
beta = 0.68*10^-2;
g = 1.5;
l = 1.6;
# Definit
浏览 11提问于2019-02-05得票数 0
1回答
问题
我现在有一个由四个方程组成的系统。两个是二阶微分方程，两个是一阶微分方程：
最初的条件是：
x = 0 |
y = 0.3 |
f(x) = 2.05 |
f(y) = 0.55 |
所有的角度都是度的。
我已经尝试过的
我曾尝试使用，并与SciPy和NumPy合作。不幸的是，我不知道如何编程，因为这些方程是非线性的。有人能给我一些其他模块的提示吗？
浏览 4提问于2020-10-15得票数 1
回答已采纳
我正在尝试使用遗传算法来优化大约10个参数。我需要将main函数中参数的中间值传递给另一个函数，该函数基本上是线性常微分方程组的ODE求解器，参数是微分方程中的系数。然后，我会将ODE的解返回到我试图最小化的main函数中。我尝试将变量声明为全局变量，但似乎没有帮助。请建议如何最好地在两个函数之间传递变量。
浏览 1提问于2012-04-28得票数 0
2回答
我目前正在用deSolve来解一个常微分方程组，我想知道是否有任何方法可以防止微分变量值降到零以下。我看过一些关于在向量、数据框架等中将负值设置为零的帖子，但由于这是一个生物学模型(T细胞计数变为负数是没有意义的)，我需要从一开始就阻止它发生，这样这些值就不会扭曲结果，而不仅仅是替换最终输出中的负数。
浏览 5提问于2017-01-14得票数 8
回答已采纳
1回答
如何用和NumPy求解线性方程组：
我们有一个方程组，在等号后面有一个值的右边。我们将所有系数写入矩阵matrix = np.array(...),，将右侧写入向量vector = np.array(...)，然后使用命令np.linalg.solve(matrix, vector)查找变量。
但是，如果我在等号之后有导数，并且我想对微分方程组做同样的事情，我该如何实现呢？
( lambda已知值在哪里，我需要找到A
我从库y = odeint(f, y0, t)中看到了这个命令scipy的用法，但是我不知道如何设置我自己的函数f如果我在那里有一个矩阵，初始值y0和什么t。
浏览 6提问于2020-10-23得票数 0
回答已采纳
目前，我确实使用来求解一个形式为0= F(t，y(t)，y'(t))的隐式微分方程。我想使用一个解决器，它附带了一个本机枕安装，并查看了的文档。文档中说，求解者解决了表单的一个系统：
求解一个系统的初值问题。该函数数值积分了给定初值dy / dt = f(t, y) y(t0) = y0的常微分方程组。
我确实怀念解决表单dy / dt = f(t, y, y')系统的能力。这与scipy.integrate.solve_ivp是可能的吗？
浏览 1提问于2018-06-11得票数 3
回答已采纳
1回答
我想用boost::numeric::odeint::runge_kutta4.求解一个耦合微分方程组这是一个三维的晶格系统，所以我很自然(也很方便)使用3D阵列。runge_kutta4是否有一种方法可以使用用户定义的数据结构或增强multi_array？
浏览 2提问于2015-10-03得票数 3
1回答
我正在尝试绘制一个常微分方程组的解。代码是：
tspan = [0 10];
z0 = [0.01 0.01 0.01 0.01];
[t,z] = ode45(@(t,z) odefun3(t,z), tspan, z0);
plot(z(:,3))
为什么输出绘制在区间0,60上，而不是代码中的0,10？
浏览 3提问于2018-06-10得票数 0
我想把一个n阶常微分方程简化为一个一阶方程组。这是为数值分析做准备。我使用Sympy和Sagemath for Computer Algebra，但我没有在它们中找到任何函数来执行这种类型的约简。我不确定是否有其他人可以指出此功能是否存在于Sympy或Sagemath中。
这方面的一个例子是减少方程式：
x''' - 2x'' + x' = 0
一阶方程组：
[0 1 0
0 0 1
0 -1 2]
浏览 0提问于2014-08-31得票数 2
3回答
迭代列表以创建函数
、、、
我正试图用solve_ivp来求解一个大的微分方程组。
from scipy import integrate
sol = integrate.solve_ivp(func_system, (0,100), initial_value_array, t_eval)
func_system是一个依赖于列表X长度的微分方程组。
X = ['a','b','c','d',...]
K1, K2, K3, K4 = [1, 2, 2, 1]
def dadt(X,t):
return -K1*X[0]+K2*X[1]
def
浏览 4提问于2020-04-06得票数 1
回答已采纳
1回答
我是个新手，还在学习中。我浏览了堆栈交换中关于使用渐近符号求解具有初始条件的微分方程系统的文档和问题。
我有一个简单的ODE-s系统。
( dV/dt ) = -( 1 / RC ) * ( V(t) ) + I(t)/C
( dI/dt ) = -( R1 / L ) * ( I(t) ) - ( 1 / L) * V(t) + Vs/L
具有初始条件V(0) = V0和I(0) = I0
我在堆栈交换中浏览了很多问题，但没有成功地找到合适的答案。如果有人能告诉我一个语法来输入一个有初始条件的耦合微分方程组，那将是很有帮助的。
浏览 1提问于2014-10-03得票数 6
1回答
如果我有一个非线性常微分方程组，M(t，y) y‘= F(t，y)，那么当我的质量矩阵M有时是奇异的时候，最好的解法是什么？
我正在处理下列方程组：
如果是t=0，这可以归结为一个微分代数方程。然而，即使我们限制t>0，无论何时y4=0，这都变成了一个微分代数方程，我不能设置一个域限制来避免(并且是我试图建模的系统的一个组成部分)。我以前唯一接触DAE的时间是整行为0 --但在这种情况下，我的质量矩阵并不总是奇异的。
在数字上实现这一点的最佳方法是什么？到目前为止，我已经尝试使用Python，在M的主要对角线上添加一个小的数(0.0001)，并将其反演，求解方程y‘= M^{-
浏览 8提问于2021-02-08得票数 2
回答已采纳
1回答
求解两个耦合二阶边值问题
、、、
用模solve_bvp求解了一个具有两个边界条件的单二阶微分方程。然而，现在我正试图求解二阶微分方程组；
U'‘+ a*B’=0
B'‘+ b*U’=0
在边界条件U(+/-0.5) = +/-0.01和B(+/-0.5) =0。我把它分解成一个第一常微分方程组，我试图用solve_bvp数值解它们。但是，我的解决方案只是得到满是零的数组。我相信我执行边界条件是错误的。我不清楚如何处理文档中的两个以上方程。我的尝试如下
import numpy as np
from scipy.integrate import solve_bvp
import matplotlib.pyp
浏览 5提问于2020-02-26得票数 1
回答已采纳
1回答
我有一个由16个耦合微分方程组成的系统，我使用的是Scipy.integrate.odeint软件包。
代码必须在0到5千兆的范围内运行，然后绘制数据，微分方程的本质是它们是非常振荡的。最重要的是，我使用的时间步骤是1e+8，或者我得到了这个常量错误：
‘在这个调用上完成的进程工作(可能是错误的Dfun类型)。运行full_output =1来获取数量信息。
我只是在想，还有什么其他方法可以让我的积分器更快地解决呢？对于大型耦合ODEs系统，是否有更好的集成方法？代码已经运行了20多个小时，但没有效果。
浏览 0提问于2018-07-18得票数 0
回答已采纳
Introduction
注:代码和指示：
单d是一阶导数双d是二阶导数
通过数值求解控制LaGrange方程，用Matlab对一些动态系统进行了数值模拟。基本上是一组二阶常微分方程。我正在使用ODE45。我从Mathworks (下面的教程链接)上找到了一个关于如何求解基本的二阶常微分方程组的很好的教程。
根据本教程，我通过获得两个二阶常微分方程(一个用于角θ，另一个用于弹簧伸长)来模拟弹性弹簧摆的运动，如下所示：
theta双素数方程：
M*thetadd*(L + del)^2 + M*g*sin(theta)*(L + del) + M*deld*thetad*(2*L + 2*de
浏览 4提问于2017-01-25得票数 0
我有一个具有外部确定性输入(控制)和随机分量的常微分方程组。除了初始状态之外，我如何安全地(执行良好的代码风格)通过tf.contrib.integrate.odeint()将这些额外的输入参数传递给等式函数？如果有这样的方法的话。或者在外部作用域中定义它们，并从方程函数中引用它们是迄今为止唯一的方法？
浏览 8提问于2017-08-17得票数 0
1回答
求解Fipy中的多个PDE
、、、、
我试图用来解一个偏微分方程组。我有至少3份政策制定文件的系统。
我想知道哪种方法最适合解决这样的系统？Fipy支持两种用于PDE系统的方法。我想知道，首先，如果Fipy甚至支持两个以上的耦合方程，如果不支持(或者如果支持)，那么解决这样一个系统的最佳方法是什么？
非常感谢。
浏览 2提问于2016-06-27得票数 2
回答已采纳
我用这个八度码来解微分方程。
# Define the right-hand side of the equation:
xvall= -11 ;#xvall
xvalu= 10 ;#xvalu
range=5000;
function ret=f(x,t);ret= t ;end;
# ywill be the values of the function at these moments of time.
t=linspace(xvall,xvalu,range);
y=lsode ('f', 2, linspace(xvall,xvalu,range));
y
plo
浏览 3提问于2015-10-19得票数 0
回答已采纳
我想解决一个由4个耦合的微分方程组成的系统，其中包含python (sympy)：
eqs = [Eq(cP1(t).diff(t), k1*cE1(t)**3), Eq(cE1(t).diff(t), -k1 * cE1(t)**3 + k6 * cE3(t)**2), Eq(cE2(t).diff(t), -k8 * cE2(t)), Eq(cE3(t).diff(t), k8 * cE2(t) - k6 * cE3(t)**2)]
当我尝试用“dsolve_system”来解决系统问题时：
solution = dsolve_system(eqs, ics={cP1(0): 0, cE1
浏览 7提问于2020-12-21得票数 0
1回答
python中的微分方程求解系统
、、、、
我在试着用python解一个微分方程组。我有一个由两个方程组成的系统，其中我有两个变量，A和B，初始条件是A0=1e17和B0=0，它们同时变化。我使用ODEINT编写了以下代码： import numpy as np
from scipy.integrate import odeint
def dmdt(m,t):
A, B = m
dAdt = A-B
dBdt = (A-B)*A
return [dAdt, dBdt]
# Create time domain
t = np.linspace(0, 100, 1)
# Initial condi
浏览 48提问于2021-02-25得票数 -1
回答已采纳
我正在尝试用DSolve来求解数学中的线性双调和方程。我认为这个问题不仅限于双调和方程，而且当我试图求解它时，MATHEMATICA只是吐出了方程。
我试过解决其他偏微分方程，没有问题。
双调和方程就是：
Laplacian^2[f]=0
这是我的方程式：
DSolve[
D[f[x, y], {x, 4}] + 2 D[D[f[x, y], {x, 2}, {y, 2}]] +
D[f[x, y], {y, 4}] == 0,
f,
{x, y}]
解决方案如下所示
DSolve[(f^(0,4))[x,y]+2 (f^(2,2))[x,y]+(f^(4,0))[x,y]==0
浏览 0提问于2011-10-18得票数 3
回答已采纳
警告，我不确定这个问题是否属于这里或数学堆栈交换。这主要是一个性能问题，但很大程度上涉及到一个微分方程求解器。
无论如何，在Matlab中，我目前正在使用以下代码语法来求解一个微分方程：
f0 = [tot;zeros(size-1,1)];
options=odeset('RelTol',1e-7, 'AbsTol', 1e-10);
[T,S] = ode15s(@myfunc, [t0,tF], f0, options, k1, k2, k3, n, size, r);
%Rest of code, including a function myfunc
浏览 2提问于2015-07-07得票数 1
回答已采纳
0回答
用RK4实现的二阶ODE
、、、
我需要绘制一个显示两个变量的图形，并使用RK4的二阶ODE，到目前为止，我已经完成了这项工作
from numpy import arange
from pylab import plot,xlabel,ylabel,show
Qger = 400
K = 20
T1 = 150
T2 = 60
N = 1000
h = (T2-T1)/N
rpoints = arange(6.0,8.0,h)
xpoints = []
x = 423
def df(s,t):
dTdt = -Qger*t/(2*K) + 172.8/t
return dTdt
for r in rpoin
浏览 5提问于2017-06-22得票数 0
我正在做一些关于解微分方程组的练习。我想在同一张图上绘制许多解决方案曲线，所以我做了一个for循环：
for i=0:6
[t, x] = ode45(@exc_345, [-4 4], [-i i]);
hold on
plot(x(:,1),x(:,2));
hold off
end
它起作用了，但现在我发现它不是“正确的”：hold on和hold off命令必须在循环之外。为什么它会起作用呢？谢谢。
浏览 1提问于2014-12-30得票数 0
2回答
我正在使用MATLABs套件来解决一个由13个微分方程组成的系统，这些微分方程决定了神经元的行为。现在我想添加第二个神经元，它由同一组微分方程运行，但受到第一个神经元的影响。更重要的是，第二个神经元也会影响第一个神经元。(这两个单元之间的前馈和反馈。)
有没有一种方便的方法可以做到这一点？我可以将微分方程分布在两个函数文件中吗?或者我是否必须将它们复制到原始文件下面，以便在同一文件中有一个更长的方程列表？我希望每个单元格都有一个文件，并以某种方式将其组织起来。(另外，以防我想再次将其扩展到三个或四个神经元。)
如果我的问题有任何不清楚或不够具体的地方，请指出。我会试着解释我在做什么/试得更好。
浏览 1提问于2014-05-20得票数 1
2回答
我正在实现的任务是将ECG信号作为一维数组文件输入，做一些处理以检测心率。
第一步是用5点差分方程求微分值，去掉低频值，我搜索过八度差分，但都是关于多项式的。那么如何在octave/mat-lab命令中实现这一点呢？
谢谢
浏览 2提问于2014-11-30得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
