differential-equations

23热度

2回答

我有一组定义3D轮廓的3D点。我想要做的是获得对应于这个轮廓的最小表面表示(see Minimal Surfaces in Wikipedia)。基本上这需要求解一个非线性偏微分方程。在Matlab中，使用pdenonlin功能(see Matlab's documentation)几乎是直接的。可以在这里找到其用于解决最小表面问题的示例：Minimal Surface Problem on t

-1热度

2回答

ODE函数的连续绘图

我有一个ODE求解器，工程很好，很流畅，但我需要在一个图中绘制所有图。连接图（1）+（3）和图（2）+（4），我必须设置启动和停止条件，但它不适用于我，我是在死胡同。我试图用x_m设置结束条件而没有结果。 options = odeset('Events',@events); [t,y] = ode45(@ph1,[0,w_max],[0,0], options); figure(1),pl

0热度

1回答

的Scilab Xcos，插入在Xcos图中的差分方程

我怎样才能插入Xcos图中的差分方程，如： Y（K + 1）= Y（k）的[一个 SQRT（Y（K -1））] + b * y（k-1）; ？感谢，问候编辑查找互联网上我认为要解决我的问题，最好的办法就是使用： scifunck block。

0热度

1回答

解决微分方程绝对值

我想解决这种形式的方程： X” = -A.x + B. | SIN（100 * PI * T）| 我用ODE45这样的： function find_x t = 0:0.001:10; x0 = 0; R1 = 90000; R2 = 1000; C = 0.001; [t,x]=ode45(@rhs, t , x0); plot(t,x); function dxdt

-1热度

1回答

整合轨道轨迹2

原来的二阶微分方程是 x'' - 2 * omega * y' - omega ** 2 * x = - mue * (x + pi2 * r12)/np.sqrt((x + pi2 * r12) ** 2 + y ** 2) ** 3 - mum * (x - pi1 * r12)/np.sqrt((x - pi1 * r12) ** 2 + y ** 2) y'' + 2 * omega *

0热度

2回答

如何数值测试振荡曲线的稳定性？

我现在有一个振荡曲线，它是一组非线性常微分方程解的一部分。随着时间的流逝，我需要测试这条曲线的稳定性/收敛性。如何用Matlab做到这一点？这个数字看起来是这样的：

4热度

1回答

如何在Python中使用dorpi5或dop853

我已经浏览了scipy.integrate.ode，但我无法找到如何实际使用这些集成方法，dorpi5和dop853。我想尝试整合ode integration python versus mathematica我的Python代码与这两种方法，看看它是如何影响结果，但不知道如何。

1热度

1回答

在Matlab中求解微分方程，ode45

我想用Matlab中的函数ode45来求解一个带有三个微分方程的系统。我真的不明白我得到的错误，我可以用一些帮助来理解我做错了什么。微分方程如下： F1 = -k1y1+k2(y2-y1) F2 = -k2(y2-y1)+k3(y3-y2) F3 = -k3(y3-y2) 而且我在Matlab代码是这样的： function dz = kopplad(t, z) global m1 m

-5热度

1回答

计算ODE与MATLAB

微分方程： α '（T）= S（β-βα+α-Qα^ 2） β'（T）=（S^-1 ）（ - β-αβ+γ） γ'（T）= W（α-γ）初使值 α（0）= 30.00 β（0）= 1.000 γ（0）= 30.00 计算我要解决从T_0 = 0的问题至t = 10，而使用数值• = 1, q = 1和w = 0.1610 我不知道如何编写ODE的函数，并非常感谢帮助！

2热度

2回答

使用ODE45查找函数的最大值

我试图在MATLAB中的微分方程系统中找到其中一个方程的位置。我试图使用odeset的事件属性。我如何挑选我的函数中的特定方程？ options = odeset('Events',@event); [t x tm xm ie] = ode45(@Lorenz,[0 200],I,options); function X = Lorenz(t,x) r = 15; sigma = 10