1回答
我试图在区间0-3中近似并绘制u"(x) = exp(x)的解，边界条件为x(0)=1，x(3)=3。我能够绘制近似解与精确解的关系图，但图看起来有点不对劲：a=0;n=10;
alpha=1;%grid size这是我正在遵循的方案：
浏览 4提问于2020-02-22得票数 0
回答已采纳
我有一个代数方程的模型，偏微分方程和常微分方程。我用有限差分法做了偏微分方程的离散化，就像Tiller的例子6.9，用Modelica介绍PhysicaL建模一书。然后，我的模型的初始条件是在时间零点离散化的所有不同点的数组。在我的模型中有一个输入：问题是我需要在模型外部设置这个数组，并出现以下错误
浏览 5提问于2020-09-21得票数 1
3回答
我是提亚诺的新手，我试图实现一个反应扩散系统的数值积分器--这个版本的：
目前，我的表达方式是：import theano.tensor as'e')a1 = T.dscalar('a1')
dvdt = e*(u - a1*v - a0)
所以我还没有实现有限差分我的</em
浏览 7提问于2016-02-21得票数 8
回答已采纳
1回答
线性回归卷积核
、、
换句话说，具有可变大小的移动窗口的线性回归。return np.convolve(aSignal, aKernel, 'same')例如，对于一维信号，一阶导数核是[-1,0,1]。可以计算窗口大小为5 ([a,b,c,d,e])时的斜率吗？
浏览 0提问于2015-05-26得票数 0
1回答
我正在尝试设计一种用于有限差分法的算法，但我有点困惑。所讨论的常微分方程是y''-5y'+10y = 10x，其中y(0)=0，y(1)=100。然后将得到的系数存储到一个矩阵中，这将是我将通过高斯-乔丹求解的系统的矩阵。问题归结为如何获得这些系数并将它们移动到矩阵中。我想过手动计算系数，然后输入矩阵，但我需要对大小为0.1，0.001和0.001的步长执行此操作，所以
浏览 16提问于2019-04-07得票数 0
1回答
假设存在一个具有圆区域的泊松方程，且边界条件为零。我想用fininte差分法进行数值求解。然而，据我所知，有限差分在角区域有着广泛的应用。那么，如果我仍然使用有限差分在圆区域，我应该做什么？
浏览 7提问于2022-04-13得票数 0
我很好奇MATLAB将如何解决这个特殊问题的有限差分法。我知道，如果我们有一个线性常微分方程，例如y'' + (e^x)y = 0，具有相同的边界条件，那么程序就相当简单。假设我们使用间隔0,1的分区为20个相等的子间隔，则以下代码将工作：h = 1/N;A(1,1) = 1;
F(1) = 0;= 2/h^2+e
浏览 6提问于2017-01-13得票数 0
我在excel电子表格中得到了这些数据。所以导入之后，我只需要做velocity = cumtrapz(t,y)，然后做position = cumtrapz(velocity)
浏览 0提问于2014-04-25得票数 0
回答已采纳
F(i) = integral(f,(i-1)*dt,i*dt);end
在for循环中，这两种计算都给出了f或a的粗略逼近在数学堆栈交换中，我找到了一个，它在一点上导出了一个像积分一样的有限差分法。然而，当我在Matlab中进行计算时，它输出了一个缩小版本的f，这一点在绘图后就很明显了(请参阅上面我所说的缩小的意思)。我认为这是因为对于较小的</e
浏览 2提问于2018-09-12得票数 1
回答已采纳
4回答
用C#求解偏微分方程
、、、
我正在做一个项目(C#和.NET框架)，它需要我解决一些偏微分方程。有没有基于.NET框架的特定库可以让我看到，让我的工作变得更简单？
我曾经使用过MATLAb，在那里求解偏微分方程非常简单。
浏览 4提问于2009-09-20得票数 5
回答已采纳
我已经开始编写一个程序，用有限差分法计算二维Schrödingers方程的解。我想要使用等高线图或一些其他图形显示，通过接受用户输入的尺寸和网格点数量，以图形方式显示解决方案。我将hbar^2/2m设置为1，并将电位(V)设置为0，从而简化了Schrödinger方程，从而得出： -(dψ^2/dx^2 + dψ^2/dy^2) = E*ψ 使用有限差分法，方程
浏览 16提问于2020-04-27得票数 0
我用MATLAB的函数'pdepe‘来解一个偏微分方程的问题，一个抛物方程。谢谢你帮我。
浏览 4提问于2017-04-02得票数 0
回答已采纳
2回答
我编写了一个简单的程序来计算函数的一阶和二阶导数，使用函数指针。我的程序计算正确的答案(或多或少)，但对于某些功能，准确性低于我想要的。这就是我要区分的功能： return (x * x);这些是导数函数，使用中心有限差分法：
浏览 3提问于2019-07-16得票数 2
回答已采纳
1回答
如果R是平面(0,1) x (0,2)的区域，设L是2-d Laplace算子，并考虑R上的Poisson方程Lu =4，其中一个解是函数v(x，y) =(X)^2。设g是v对R边界的限制，设h=k= 1/2，m=3 =5。
这是我的代码:但是它不起作用。
浏览 0提问于2017-04-19得票数 0
回答已采纳
1回答
在我的课程作业中，我被告知要用Python编写一个算法来求解具有Dirichlet边界条件的欧拉-拉格朗日方程。非常感谢!
浏览 0提问于2018-11-17得票数 0
1回答
python中三对角矩阵的求解
、、、、
我一直在研究求解化学反应微分方程的数值方法。通常用有限差分法将微分方程放入三对角矩阵，然后用列向量表示边界条件。一旦我有了矩阵和向量，我就使用了枕木的直线。
浏览 4提问于2022-11-10得票数 0
回答已采纳
我试图让maxima执行一些“抽象的”Taylor级数展开，我遇到了一个简化的问题。这个问题的原型可能是梯度的有限差分模拟，taylor(g(x1,dx1), [dx1但是现在试试二阶导数的有限差分模拟(Hessian)，
h(x1,dx1) := (f(x1+dx1)
浏览 1提问于2018-06-04得票数 3
回答已采纳
1回答
Solve_ivp是一个初值问题的求解函数。，method=‘RK45 45’，t_eval=None，dense_output=False，events=None，vectorized=False，args=None，**选项)参数有趣，系统右侧可调用.呼叫签名很有趣。或者，它可以具有形状(n，k)；然后，fun必须返回一个形状为(n，k)<
浏览 2提问于2020-02-18得票数 3
1回答
我把布莱克·斯科尔斯方程转化为热火方程。我尝试用显式有限差分法来求解这个PDE，得到看涨期权的价格。我也用“解析”的黑色schols方程来解决这个问题。import numpy as npsigma
浏览 1提问于2016-04-09得票数 1
回答已采纳
1回答
我必须解Euler Bernoulli微分光束方程，它是：(x是梁轴点的坐标)u(0)=0, u’(0)=0, u’’(1)=0,u’’’(1)=a U’k = (1/2*h)(Uk+1 - Uk-1)
U’’’’k = (1/h4)(U
浏览 4提问于2014-11-16得票数 3
云服务器
ICP备案
云直播
对象存储
实时音视频
腾讯云开发者
