mathematica中方程组的NDsolve问题
在Mathematica中，NDsolve函数用于求解常微分方程组。它是Mathematica内置的数值求解器，可用于解决一阶或高阶常微分方程组。
NDsolve函数的语法如下： NDsolve[{eqn1, eqn2, ...}, {y1, y2, ...}, {x, xmin, xmax}] 其中，eqn1, eqn2, ...是方程组的方程，y1, y2, ...是待解函数，x是自变量，xmin和xmax是解的区间。
NDsolve函数的返回值是一个规则列表，其中包含了待解函数的数值解。
方程组的分类：
- 一阶常微分方程组：方程组中的未知函数只出现了一阶导数项。
- 高阶常微分方程组：方程组中的未知函数出现了高阶导数项。
NDsolve函数的优势：
- 高精度：NDsolve使用数值方法求解，可以得到高精度的数值解。
- 强大的功能：NDsolve支持求解各种类型的常微分方程组，包括线性和非线性方程组。
- 自动选取合适的数值算法：NDsolve会根据方程组的性质和求解要求，自动选择适合的数值算法进行求解。
- 可视化功能：NDsolve可以将求解结果以图形的形式展示出来，帮助用户更直观地理解解的特性。
方程组的应用场景：
- 物理建模：方程组可以用于描述各种物理现象，如运动学、电路分析等。
- 工程问题：方程组可以应用于工程领域中的建模和仿真，如控制系统、信号处理等。
- 生物学模型：方程组可以用于建立生物学系统的数学模型，如生态系统、遗传网络等。
推荐的腾讯云相关产品： 腾讯云提供了多种云计算相关的产品，可以满足不同用户的需求，以下是一些推荐的产品：
- 云服务器（CVM）：提供可扩展的计算资源，适用于搭建和部署各种应用程序和服务。 产品介绍链接：https://cloud.tencent.com/product/cvm
- 云数据库MySQL：提供稳定可靠的云数据库服务，支持高并发访问和弹性扩展。 产品介绍链接：https://cloud.tencent.com/product/cdb_mysql
- 弹性伸缩（Auto Scaling）：根据负载情况自动扩展或缩减云服务器的数量，实现按需分配计算资源。 产品介绍链接：https://cloud.tencent.com/product/as
请注意，以上推荐的腾讯云产品仅作为示例，并不代表其他云计算品牌商的产品。
相关·内容
在尝试在两个单独的绘图命令中使用来自NDSolve的答案时，我遇到了一个问题。为了说明这个问题，我将使用一个简单的微分方程和一个绘图命令。如果我写这样的东西：
{Plot[x[t], {t, 0, 10}], x[4]}
/. NDSolve[{x'[s] == - x[s], x[0] == 1}, x, {s, 0, 10}]
它解决了方程，没有问题地计算x4，但是图变空了，我不知道为什么。
在我的实际问题中，我的方程对于多个函数来说是一个相当复杂的系统，而不是x4，而是画出了求解函数的参数图。我最终打算在一个操作语句中包含所有这些内容，所以我不希望NDSolve语句出现多次(花
浏览 6提问于2009-08-25得票数 6
回答已采纳
用C语言编写的常微分方程求解器，可能使用GSL库，与Mathematica8.0 NDSolve相比有显着的速度优势吗？在准确性方面，它如何公平？
我的理解是，原则上编译的代码可以更快，但现在NDSolve已经以某种方式使用了大量的编译代码本身？
另外，有没有什么方法可以使用MathLink或Mathematica的编译函数来加快解算常微分方程的速度？
浏览 4提问于2012-10-19得票数 3
回答已采纳
我有一个关于Mathematica NDSolve中的约束的问题
如果我有两个线性微分方程
NDSolve[x'[t]= A x[t]+B y[t],
y'[t]= C x[t]+D y[t],
x[0]==0.002, y[0]==0.005, {x,y}, {t,0.10000}]
我如何施加约束？
x[t]+y[t] == 1
对于任何时间t
谢谢
浏览 0提问于2013-06-07得票数 1
1回答
我使用NMaximize从NDSolve函数中获取值：
Flatten[NDSolve[{x''[t] == (F Cos[\[CapitalOmega] t] -
c x'[t] - (k + \[Delta]kb) x[t] + \[Delta]kb y[t])/m,
y''[t] == (-c y'[t] - (k + \[Delta]kb) y[t] + \[Delta]kb x[t])/m,
x'[0] == 0, y'[0] == 0, x[0] == 0, y[0] == 0}, {x[t], y[t]}, {
浏览 4提问于2013-09-27得票数 1
回答已采纳
我在Mathematica中用NDSolve解决常微分方程问题时发现了一个令人困惑的问题，代码如下:我尝试了'Clear[Derivative]'并重启内核，但这些方法都不起作用。代码： Clear[Derivative];
ClearSystemCache;
r = 0.3; a = 3; delta = 0.45; M0 = 0.975; T = 20;
u[t] = 0.5*l[t]*delta*M[t];
eql1 = M'[t] == r*M[t] Log[1/M[t]] - u[t]*delta*M[t];
e
浏览 391提问于2019-05-07得票数 0
出现这些问题是因为我想用一个PDE的解作为初始条件来解决另一个问题。因为在NDSolve中，解决方案是由InterpolatingFunction给出的，所以我必须在第二个PDE中使用InterpolatingFunction。这个是可能的吗？为什么我的电脑要永远运行？我的代码：
(*********************Parameters********************)
\[Gamma] = 10^(-5);
T = 500;
tm = -250;
\[Sigma] = 100;
L = 1;
L0 = -(L/2);
L1 = L0 + L;
c = 29.97
浏览 0提问于2011-10-28得票数 3
我正在尝试使用NDsolve函数来求解一个PDE集。
我对数学很陌生，这是我输入的代码。
NDSolve[{D[Cm[t, x], t] == Dm*D[Cm[t, x], x, x] + Kg*Cs[t, x] - Ka*Cm[t, x],
D[Cs[t, x], t] == Ds*D[Cs[t, x], x, x] + Ka*Cm[t, x] - Kg*Cs[t, x],
Cm[0, x] == Cm0,
Cs[0, x] == Cs0,
Dm*ND[Cm[t, 0]] == 0.5*FT,
Ds*ND[Cs[t, 0]] == 0.5*FT,
浏览 3提问于2013-10-03得票数 0
1回答
我在用Mathematica 8编程时遇到了问题，下面是我的代码：
f[t_, y_] := {y, y};
RungeKutta3[a_, b_, Alpha_, n_, f_] :=
Module[{h, j, k1, k2, k3},
h = (b - a)/n;
Y = T = Table[0, {100 + 1}];
Y[[1]] = Alpha;
T[[1]] = a;
For[j = 1, j <= n, ++j,
k1 = f[T[[j]], Y[[j]]];
k2 = f[T[[j]]
浏览 1提问于2011-12-31得票数 0
1回答
我一般不使用Mathematica，我需要将它与其他程序进行比较。我想要解一个由三个微分和非线性方程组组成的系统。为此，我使用了解决方案。当我把非线性项(指数)加进去时，一切都会出错。
这是我的代码：
equa = {x'[t] == z[t] - Exp[y[t]],
y'[t] == z[t] - y[t],
z'[t] == x[t] + y[t] - z[t],
x[0] == 0,
y[0] == 0,
z[0] == 0};
slt = DSolve[equa, {x, y, z}, t]
Plot[{x[t] /. slt}, {t, 0, 10
浏览 1提问于2012-03-05得票数 0
我一直在学母猪/收割。它们是很酷的构造。但我需要帮助，看看我是否可以用它们来做我将解释的事情。
我想做的是:在NDSolve运行时绘制它的解决方案。我在想，在NDSolve使用EvaluationMonitor运行时，我可以使用EvaluationMonitor收集解决方案(x，yx)。但是我不想等到最后，Reap它然后绘制解决方案，而是想在它运行时执行它。
我将展示基本的设置示例
max = 30;
sol1 = y /.
First@NDSolve[{y'[x] == y[x] Cos[x + y[x]], y[0] == 1},
y, {x, 0, max}]
浏览 1提问于2011-12-28得票数 10
回答已采纳
1回答
在Matlab中能象征性地求解下面的方程组吗？例如，a和b没有在下面定义：
x1 - (5-a)x2 =2
2x1 + 5x2 =b
x1+x2+x3=0
是否可以使用Matlab (没有先进的库)来解决这个系统？我知道Mathematica提供了这种可能性，但我被要求在Matlab中做到这一点。
浏览 0提问于2013-02-01得票数 0
回答已采纳
7回答
示例：
在……里面
CT =表[预备[10^4*xrange2 2 / ]。NDSolve{...series of equations .}，{t，range1，range2}，MaxSteps -> 10000，PrecisionGoal -> 11，delay]，{ delaymin，delaymin，delaymax，0.1}；//定时
分号之后的// Timing是什么意思？
在……里面
D-I*=(H0+ V).ρ-ρ.(H0 +V)-Γ*ρ//简化；
这里的// Simplify呢？
我哪儿也找不到这个解释！
谢谢你，蒂亚戈
浏览 4提问于2011-05-04得票数 1
1回答
当我试图处理我的实验数据时，在Mathematica中遇到了一个内存问题。我正在使用Mathematica找出一个由三个偏微分方程组成的系统的最优参数。
当e参数大于0.4时，Mathematica会消耗大量内存。对于e < 0.4，程序运行正常。
我尝试过使用$HistoryLength = 0，减少AccuracyGoal和WorkingPrecision，但都没有成功。
我正在尝试了解我犯了什么错误，以及我如何限制内存使用。
Clear[T, L, e, v, q, C0, R, data, model];
T = 13200;
L = 0.085;
e = 0.41;
v
浏览 0提问于2011-11-29得票数 2
回答已采纳
我们有这样的系统：
x'[t] == x[t] - 5 y[t] + z[t]
y'[t] == 3 x[t] - 3 y[t] - 3 z[t]
z'[t] == -2 x[t] + 10 y[t] + 4 z[t]
和初始条件：
x[0] == .01
y[0] == 3
z[0] == 0
我制作了具体的情节：
eqn = {x'[t] == x[t] - 5 y[t] + z[t], y'[t] == 3 x[t] - 3 y[t] - 3 z[t],
z'[t] == -2 x[t] + 10 y[t] + 4 z[t]};
s
浏览 0提问于2012-12-18得票数 5
回答已采纳
7回答
我必须解决多项式方程组，因为它有无限的解，并且我只需要很少的解(任意2或3)，那么我如何才能得到它们呢？，我可以在值在1到10之间的解上指定条件，这样我就可以得到很少的值吗？方程实际上很长很复杂，但无穷多的解是由于根部的"sin(0)“。
浏览 1提问于2010-08-11得票数 0
1回答
我需要解一个形式为w'=g(t,w(t))的微分方程，其中g的定义如下
g[t_, w_] := {f1[t, {w[[3]], w[[4]]}], f2[t, {w[[3]], w[[4]]}], w[[1]],w[[2]]};
和f1，f2是
f1[t_, y_] := Sum[\[Mu][[i]] (s[[i]] - y)/Norm[s[[i]] - y]^2, {i, 1, 5}][[1]];
f2[t_, y_] := Sum[\[Mu][[i]] (s[[i]] - y)/Norm[s[[i]] - y]^2, {i, 1, 5}][[2]];
其他一切都是正确定义的，不
浏览 1提问于2017-12-15得票数 1
1回答
我想用Mathematica和复制一个使用创建的ODE求解器。
下面是使用NDSolve的数学代码
result[r_] := NDSolve[{
s'[t] == theta - (mu*s[t]) - ((betaA1*IA1[t] + betaA2*IA2[t] + betaB1*IB1[t] + betaB2*IB2[t]) +
(betaA1T*TA1[t] + betaA2T*TA2[t] + betaB1T*TB1[t] + betaB2T*TB2[t])) * s[t] -
浏览 0提问于2013-10-06得票数 1
回答已采纳
我是mathematica的新手，所以请容忍我！
我正在尝试用NDSolve来解决mma中的一个非线性偏微分方程。由于奇点发生在模拟时间耗尽之前很久，求解过程被缩短。我意识到，具有这种奇点的僵硬系统可以通过减小步长来处理(至少通过蛮力)。
然而，"MaxSteps“或"MaxStepSize”似乎对我的代码没有明显的影响。
怎么回事？还有其他我可能遗漏的方法吗？
**
代码：
**
Needs["VectorAnalysis`"]
Needs["DifferentialEquations`InterpolatingFunctionAnatomy`"
浏览 0提问于2011-09-23得票数 2
回答已采纳
我需要使用NDSolve，这反过来又将来自另一个ODE的解决方案作为函数从另一个NDSolve输出。
如果我使用NDSolve内部第一个微分方程的精确解，就可以了。但是，当我以函数的形式(使用InterpolatingFunction)使用相同的解决方案时，它不起作用。
我相信，这与NDSolve输出的结构有关。有人能告诉我这件事吗。会有很大帮助的！
守则是：
feq = 2 V alpha fip F''[fi] - (V^2 - (V^2 + sigma - 2 fi) (F'[fi])^2 + (F'[fi])^4
Frange[lo_, hi_] :=
浏览 2提问于2014-05-05得票数 1
回答已采纳
2回答
我不确定这是不是一个合适的地方。其中N、L、H、p和d是参数。我需要解这个方程组。具体地说，我需要求解b(t)和e(t)。
Variables | t=1 | t>1
----------|--------|------------------------
n(t) | N |N(1-p)^(t-1)
s(t) | 1 |((1-p+dp)/(1-p))^(t-1)
b(t) | L |b(t-1)+p(H-b(t-1))
e(t) |(H-L)/2 |e(t-1)+(p(H-b(t-1)))/2
浏览 54提问于2016-08-12得票数 1
我试图解决一个方程组(5个未知变量，5个方程)，但Solve[]函数只是挂起，我不得不中止评估。我能理解其中的一些方程是相当混乱的--至少在我看来(我不是数学家)。
我用“已知/真实”的模拟值代替了Solve[]中使用的方程，它们都算出来了。
因此，我的问题是:是否可以“帮助”Solve[]说，例如.
Solve[{eq1, eq2, eq3, eq4, eq5},{var1, var2, var3, var4, var5}, (*code here along the lines of { 0 < var1 < 10, var2 < 25, ...}*)]
如果有帮助，我可
浏览 5提问于2012-09-07得票数 1
1回答
我是第一次接触Mathematica，正在尝试弄清楚如何动态生成ODE系统。例如，我有一个由100个方程组成的系统，其中每10个方程本质上是相同的，但参数略有不同，可以从一个长度为10的向量中读取。我想写出10个方程式，然后循环一些迭代器来生成所有100个方程式。有没有一种标准的方法来做到这一点？
例如，下面是一个由30个方程组成的系统(对于1:10中的i)：
dX_i/dt = -\beta*X_i*Y_i + \delta_{i-1}*X_{i-1} - \delta_i*X_{i}
dY_i/dt = \beta*X_i*Y_i - \gamma_i*Y_i + \delta_{i-
浏览 0提问于2012-09-18得票数 1
回答已采纳
我遇到了一个涉及NDSolve的问题，我用不同的参数值多次运行它。对于其中的一些值，解决方案会产生奇点，NDSolve使用NDSolve::ndsz或其他相关警告进行警告。
我只想捕捉这些警告，抑制它们的显示，并跟踪这些参数的特定值出现问题的事实。我想到了以下选项(这两个选项都不起作用)：
我知道我可以通过使用Check来确定命令是否导致了警告或错误。然而，这仍将显示警告。如果我用Off关闭它，Check也无法报告警告。
可以使用NDSolve方法停止EventLocator，因此我可以检查函数或其导数的非常大的值，并在这种情况下停止计算。然而，在实践中，这仍然不时地产生警告，大概是因
浏览 1提问于2011-03-15得票数 12
回答已采纳
在数值求解微分方程并绘制结果后，我想确定绘制范围内的单个最大值，但不知道如何确定。
下面的代码用于数值求解微分方程并绘制结果。
s = NDSolve[{x''[t] + x[t] - 0.167 x[t]^3 == 0.005 Cos[t + -0.0000977162*t^2/2], x[0] == 0, x'[0] == 0}, x, {t, 0, 1000}]
Plot[Evaluate[x[t] /. s], {t, 0, 1000},
Frame -> {True, True, False, False}, FrameLabel -> {&#
浏览 1提问于2010-12-24得票数 3
回答已采纳
我想用Mathematica 7解决这个微分方程组，但我发现了一个错误，说明函数是在没有依赖于所有自变量的情况下指定的。这些方程式是：
感谢大家的帮助
浏览 1提问于2011-12-10得票数 1
回答已采纳
m = 10; c = 2; k = 5; F = 12;
NDSolve[{m*x''[t] + c*x'[t] + (k*Sin[2*Pi*f*t])*x[t] == F*Sin[2*Pi*f*t],
x[0] == 0, x'[0] == 0}, x[t], {t, 0, 30}]
{f，0，5} ( 0=< f <= 5)
如何绘制三维图像：
X= u(t，f)
.
如果f= 0.1,0.2，. 5，我们可以解方程：
NDSolve[{m*x''[t] + c*x'[t] + (
浏览 4提问于2011-11-18得票数 3
回答已采纳
2回答
我有一个由8个具有约束(最小和最大可接受值)的8个变量的7阶8个多项式方程组成的系统。
在你的经验和知识中，什么是解决这类问题的最有效的“包”？数学，Matlab，C++/ Java数学库，...
我尝试了一些天真的方法，但效果很差……所以现在我尽量做到稳妥(尽可能)
非常感谢
浏览 0提问于2013-02-21得票数 2
我正在寻找一种在Python 3.2中转换方程的通用方法。我最近才开始使用它，偶然发现了一些我以前的MATLAB作业。我可以在MATLAB中计算这个，但是pylab对我来说仍然是一个谜。
所以，我有一个方程式的文本文件，我试图将它转换成相同形式的A=b，然后在PYLAB中解决一些与之相关的线性代数问题。
文本文件"equations.txt“包含以下格式的线性方程集合：
-38 y1 + 35 y2 + 31 y3 = -3047
11 y1 + -13 y2 + -34 y3 = 784
34 y1 + -21 y2 + 19 y3 = 2949
等。
该
浏览 0提问于2013-04-13得票数 2
回答已采纳
我正在尝试使用NDSolve来求解波动方程，以检查使用它是否更容易和/或更快，而不是我以前的特征方程。方法实现。
我得到了很多不稳定的东西，而用特征方法我不知道，因为这些是简单的方程，我想知道是怎么回事.(希望不是问题的物理方面.)
ans = Flatten@NDSolve[{
u[t, x]*D[d[t, x], x] + d[t, x]*D[u[t, x], x] + D[d[t, x], t] == 0,
D[d[t, x], x] + u[t, x]/9.8*D[u[t, x], x] +
1/9.8*D[u[t, x], t] + 0.0001 u[t, x]*Abs[u[t,
浏览 1提问于2011-09-08得票数 7
2回答
我有一个关于数学中的NDSolve函数的问题。我有一个由这两个方程定义的振荡器：
x' = v
v' = -x - u*v^3
其中u是某个常数。
如何创建一个NDSolve来解决这个问题？我尝试了下面的代码(它必须依赖于时间)，但它不起作用：
eq1 = x'[t] == v;
eq2 = v' == -x[t] - u*v^3;
eq3 = x[0] == 2;
(初始位移为2m)。
s = NDSolve[{eq1, eq2, eq3}, x, {t, 0, 30}]
非常感谢..。
浏览 4提问于2012-04-30得票数 0
回答已采纳
问题是：
我试图解决这个微分方程：
K[x_, x1_] := 1;
NDSolve[{A''[x] == Integrate[K[x, x1] A[x1], {x1, 0, 1}],
A[0] == 0, A'[1] == 1}, A[x], x]
我得到了错误(Function::slotn和NDSolve::ndnum)
(它应该返回一个等于3/16 x^2 + 5/8 x的数值函数)
我正在寻找一种解决这个微分方程的方法:是否有一种更好的方式来写出它，以便NDSolve能够理解它？还有其他功能或包可以帮助吗？
注1:在我的全部问题中，K[x,
浏览 5提问于2011-08-07得票数 4
回答已采纳
1回答
我有以下两个微分方程 y'[t] == k1*y[t] - k3*y2[t]
y2'[t] == k2*y2[t] - k4*y[t] 我可以使用数学DSolve函数来求解这些方程： DSolve[{y'[t] == k1*y[t] - k3*y2[t], y2'[t] == k2*y2[t] - k4*y[t],
y[0] == y10, y2[0] == y20}, {y[t], y2[t]}, t] Mathematica给出了这个微分方程组的分析结果，形式如下： y[t] = f{t} 和 y2[t] f{t} 但是我想要的是一个yt的函数，它
浏览 15提问于2019-03-18得票数 0
在尝试重新创建由NDSolve生成的NDSolve时，NDSolve的InterpolationOrder选项遇到了非常奇怪的问题。考虑下面的InterpolationFunction (一个示例函数)：
ifun = First[
x /. NDSolve[{x'[t] == Exp[x[t]] - x[t], x[0] == 1}, x, {t, 0, 10}]]
现在让我们试着重建它。以下是数据：
Needs["DifferentialEquations`InterpolatingFunctionAnatomy`"]
data = Transpose@{Int
浏览 3提问于2011-10-10得票数 4
回答已采纳
2回答
创建一个使用克雷默法则 (带有矩阵的行列式)来求解线性方程组解的算法。代码应该对变量的"n“数起作用。
您可以使用任何数据结构来保存矩阵并返回结果
注意:另外，考虑到一个可能的系统可能没有解决方案，并且可能有无限多的解决方案:)，所以，您的解决方案应该考虑到这一点.
如果是这样，只需打印出来或返回“无”或“无限多”。
既然这是暗号高尔夫，最小的代码就赢了..。
编辑:为了使这更具有挑战性，您不能使用语言的内置矩阵操作库。
另外，您的算法不需要处理如何获取输入，只需要处理输入并返回正确的输出。如前所述，您可以将此输入存储在任何结构中。
浏览 0提问于2014-07-23得票数 2
2回答
计算四次方差
、、
我正在使用Modelica来求解一个热传递问题的方程组，其中一个方程是辐射，它写为
Ta^4-Tb^4
有没有人可以说，如果方程写成：
(Ta-Tb)(Ta+Tb)(Ta^2+Tb^2)
浏览 2提问于2011-12-15得票数 1
回答已采纳
1回答
在math.js中，在math.js主页上有一个关于求解微分方程的例子，但是它相当复杂，并且没有提供足够的信息让我个人能够将math.js应用于其他类似的问题。所以，我要做的是解决。在这个系统中有两个方程：
dx/dt = ax - bxy
dy/dt = cxy - y
在math.js里，我得到了
math.import({ndsolve:ndsolve});
const sim2 = math.parser();
sim2.eval("dxdt(x, y) = x - x * y");
sim2.eval("dydt(x, y) = x * y - y"
浏览 15提问于2017-11-19得票数 2
回答已采纳
我正在尝试使用NDSolve对一长串常微分方程进行数值求解，如下所示：
sols = NDSolve[deqs, operons, {t, 0, 1000}];
其中deqs保存微分方程列表，operons存储函数列表。
我的错误是"NDSolve::deqn:方程或方程列表应该是第一个参数中的False“
它随后给我的错误代码字符串示例是：
acs'[t]==-0.04 acs[t]+(0.08Sin[0.01t]^4)/(1+Sin[<<1>>]^4), acs[0]==0, adaalkB'[t]==-0.04adaalkB[t]+0.08
浏览 1提问于2013-05-16得票数 1
我正在研究一种解决偏微分方程的方法，确切地说，是菲克扩散第二定律。我能够使用NDSolve和Plot3D函数生成3D绘图。使用的代码：
NDSolve[{D[c[t, h], t] == 1*D[c[t, h], h, h],
c[0, h] == Erfc[h/(2*81.2)],
c[t, 0] == 1,
c[t, 4000] == 3.08*^-18}, c, {t, 0, 900}, {h, 0, 274}]
我希望找到t= 900处的图形的数值点，而不是图形表示。我想知道如何在NDSolve (
浏览 0提问于2013-06-26得票数 1
回答已采纳
2回答
我对python完全陌生，事实上，任何一种基本的编程语言，我都使用Mathematica来进行所有的符号和数值计算。我正在学习如何使用python，并发现它真的很棒！这里有一个问题，我正试图解决，但却毫无头绪！例如，我有一个数据文件。
0. 1.
0.01 0.9998000066665778
0.02 0.9992001066609779
... ..
其中只有{t，Cos2t}。我想从这些数据中定义一个函数，并使用它来求解python中的一个方程。我的数学直觉告诉我，我应该定义如下函数：
iFunc[x_] = Interpol
浏览 1提问于2017-05-07得票数 3
回答已采纳
4回答
微分方程解曲线的绘制
、、、
我有几个微分方程，我想为各种起始值N_0画出解决方案
以下是方程式：
dN\dt= bN^2 - aN
dN\dt = bN^2 (1 - N\K) - aN
我该怎么做呢？
我真的不关心所使用的语言。在专门的数学方面，我的电脑上有mathematica和matlab。我有权进入梅普。我必须做更多这样的事情，我希望有来自任何语言的例子，因为这将帮助我找出我想要使用的语言并学习它。
浏览 1提问于2011-07-23得票数 1
2回答
所以我试图使伽马泊松的似然函数最大化，我把它编程成R，如下所示：
lik<- function(x,t,a,b){
for(i in 1:n){
like[i] =
log(gamma(a + x[i]))-log(gamma(a))
-log(gamma(1+x[i] + x[i]*log(t[i]/b)-(a+x[i])*log(1+t[i]/b)
}
return(sum(like))
}
其中x和t是数据，我有n个数据行。
我需要同时解决a和b问题。R中是否存在内置函数？或者我需要硬编码
浏览 3提问于2016-10-03得票数 0
我需要用数学来解这个方程：
d/dx v(x) = A。诉(X)
这里，v是列向量 {v1(x),v2(x),v3(x),v4(x)}，A是4×4矩阵。我想对函数v1, v2, v3, v4进行任何初始条件的求解。X的范围从0到1000。
如何用NDSolve编写这类微分方程的数学程序
浏览 0提问于2012-01-25得票数 2
回答已采纳
有没有办法在Mathematica中获取表达式求值过程中生成的实际消息？假设我在数值求解一个常微分方程，它爆炸了，就像这样
In[1] := sol = NDSolve[{x'[t] == -15 x[t], x[0] == 1}, x, {t, 0, 1},
Method -> "ExplicitEuler"];
在本例中，我将得到NDSolve::mxst错误，告诉我在t == 0.08671962566152185达到了10000步的最大数量。现在，如果我查看$MessageList变量，我只收到消息名称；特别
浏览 1提问于2010-01-05得票数 3
回答已采纳
1回答
指数函数参数
、、、
我有三点[x0 y0]，[x1 y1]，[x2 y2]与严格的条件x0<x1<x2，y0<y1<y2。所有这些点都依赖于一些指数函数y=ae^(bx)+c。我需要找到a,b,c..。不可能精确地解3个方程组，所以我需要近似它。java中有没有数学库可以帮我解决这个问题？我在mathcad上找到了类似的东西
，但在java中找不到。
另一种方法-如何近似求解3个方程和3个数值的系统。
ae^(bx_0)+c=y_0 ae^(bx_1)+c=y_1 ae^(bx_2)+c=y_2
浏览 0提问于2019-02-21得票数 0
1回答
我在密码学应用程序中遇到了以下问题:我已经给出了一组线性同余
a[1]*x[1]+a[2]*x[2]+a[3]*x[3] == d[1] (mod p)
b[1]*x[1]+b[2]*x[2]+b[3]*x[3] == d[2] (mod p)
c[1]*x[1]+c[2]*x[2]+c[3]*x[3] == d[3] (mod p)
这里，x是未知数，a，b，c，d是给定的
系统很可能是欠确定的，所以我有一个很大的解决方案空间。我需要一个算法，使用伪随机数生成器(或失败)找到该问题的等分布解决方案(即在解空间中等分布)。
据我所知，在我的线性代数课程中，大多数线性方程组的标准算法都不能直接应
浏览 3提问于2011-12-14得票数 2
回答已采纳
2回答
我试图解出下列方程组：
x1 + x2 + x3 = 1
0.5 * x1 + 0.75 * x2 + 0.25 * x3 = 0.25
我首先从均匀分布中随机选择$x3$。然后，我更新我的方程组来求解剩下的x1和x2。然而，我想确保x1和x2也在0到1之间。目前，解决上述方程组的例子会导致负的x1和大于1的x2。有什么方法可以将这个约束结合起来吗？
library(matlib)
set.seed(3)
x3 <- runif(1, 0.01, 0.99)
A <- matrix(c(1, 2/4, 1, 3/4), 2, 2)
b <- c(1 - x3, 0.25 -
浏览 0提问于2019-08-20得票数 0
回答已采纳
I用NDSolve求解一个非线性偏微分方程.I希望其中一个变量(Kvar)是当前正在求解的时间步骤的函数，并因此使用Piecewise。Mathematica生成一条错误消息，上面写着：
设置延迟：：写:标签在0.05t_中是受保护的。>> NDSolve::deqn:第一个参数中所期望的方程或方程列表，而不是$Failed .
ReplaceAll::reps：.
为了便于阅读，我没有包括全部的错误消息。
我的代码如下：
Needs["VectorAnalysis`"]
Needs["DifferentialEquations`Interpolating
浏览 2提问于2011-09-21得票数 2
回答已采纳
当我用下面的代码数值求解一个ode时，出现了名为"evfrf“的警告。
我想知道如何强制微分方程的变量为实数
NDSolve[{y''[t] + .1 y'[t] + Sin[y[t]] == 0, y'[0] == 1,
y[0] == 0}, y, {t, 0, 20},
Method -> {"EventLocator", "Event" -> y[t],
"EventCondition" -> y'[t] > 0,
"EventAction"
浏览 0提问于2011-09-17得票数 1
回答已采纳
对Mathematica语法来说仍然是个新手。当我这样做的时候：
DSolve[{
D[u[x, t], {x, 2}] == (1/(v*v))*D[u[x, t], {t, 2}],
u[0, t] == 0,
u[l, 0] == 0
}, u, {x, t}]
它只返回我输入的内容
DSolve[{(u^(2,0))[x,t]==(u^(0,2))[x,t]/v^2,u[0,t]==0,u[l,0]==0},u,{x,t}]
然而，当我去掉边界条件时，我得到
{{u->Function[{x,t},C[1][t-(Sqrt[v^2] x)/v^2]+C[2][t+
浏览 2提问于2011-09-22得票数 4
回答已采纳
我想要这样的东西
w[w1_] :=
NDSolve[{y''[x] + y[x] == 2, y[0] == w1, y'[0] == 0}, y, {x, 0, 30}]
这看起来效果更好，但我想我错过了smtn
w := NDSolve[{y''[x] + y[x] == 2, y[0] == w1, y'[0] == 0},
y, {x, 0, 30}]
w2 = Table[y[x] /. w, {w1, 0.0, 1.0, 0.5}]
因为当我试着做一张桌子的时候，它不工作：
Table[Evaluate[y[x] /. w
浏览 1提问于2010-12-08得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
云直播活动推荐
腾讯云开发者
