differential-equations

1热度

1回答

我有一个简单的线性系统，其二阶ODE。系统的解析解欧拉方法现在我想解决系统和比较近似值与真正的价值。这是我在Matlab代码 clear all; clc; t = 0; dt = 0.2; tsim = 5.0; n = round((tsim-t)/dt); A = [ -3 0; 0 -5]; B = [2;3]; XE = [1; 2]; u = 1;

1热度

1回答

如何从文本文件中读取微分方程组来解决scipy.odeint的系统问题？

我有一个庞大（> 2000方程）的ODE系统，我想用python scipy的odeint解决。我有三个问题，我想解决（也许我会问3个不同的问题？）。为了简单起见，我会在这里用玩具模型来解释它们，但请记住我的系统很大。假设我有ODE的以下系统： dS/dt = -beta*S dI/dt = beta*S - gamma*I dR/dt = gamma*I 使用β= C p我，其中

3热度

1回答

我如何使用MATLAB来解决这个PDE

我有一个考试的做法以下问题：我需要使用MATLAB来解决它。问题是，我之前没有看到过这样的问题，而且我正在努力开始。我有我的1x1网格，分成10x10。我知道我可以使用1/10 * x * 2来计算除角落以外的整个底部行。我也知道我可以使用（1/10）（1 + t）^ 2来计算整个右行。但是，我无法弄清楚如何获得足够的点数来填充整个网格的值。我知道它必定与问题中的偏导数有关，但我不太确定它们在

2热度

2回答

使用Runge-Kutta求解耦合微分方程

我有一个耦合方程系统：流体静力平衡方程，质量连续性方程和理想气体状态方程。这些是，在数学语法， \frac{dP}{dr}=- \rho*g，其中\rho是密度和g是重力加速度。 \frac{dM}{dr}=4*pi* r^2*\rho 和 p=\rho* k_B* T/(\mu *m_p)，其中k_B是玻尔兹曼常数，\mu是平均分子量，而质子质量是m_p。我要解决使用龙格 - 库塔数值技术

0热度

1回答

scipy.odeint返回二阶非线性微分方程不正确的值

我一直在试图解决牛顿万有引力定律（平方反比定律）二阶非线性微分方程： x(t)'' = -GM/(x**2) 的卫星接近地球的在一名维使用numpy.odeint用一系列第一ORD的运动（在这种情况下，点质量） er微分方程，但与Mathematica或简化形式的法则（Δx=（1/2）在^ 2）相比，该操作产生了不正确的结果。这是程序中的代码： import numpy as np from

0热度

1回答

两个变量的图形在x轴上同时变化

我的下面的代码生成一个带有x轴循环变量的图形。具体而言，eta_22从0到1变化，循环迭代大小为0.01。该行下面的代码是源函数文件。我的问题是：我该如何生成一个图表，其中eta_1从0到1变化，循环迭代大小为0.01？（我想在y轴上绘制AA图，并且eta_1，eta_2从0到1变化）。我尝试了：我试图创建嵌套的“for”循环，但绘图本身循环。我试图将情节线放在“for”循环之外，但那不起

1热度

1回答

在Matlab

我如何解决这个耦合随机微分方程我有单一的Hindmarsh - 罗斯（HR）的神经元模型如下 x' = y - a*x^3 + b*x^2 -z + I0 + I1*cos*w*t + D*Zyi(t); y' = c - d*x^2 -y; z' = r[s(x - x0) - z ]; 其中a，b，I1，I0，d，C，d，R， s，x0是参数。并且包括具有D强度的白高斯噪声（Zyi

3热度

1回答

在Mathematica中求解一个微分方程

我有一个语法问题在Mathematica（第10版）中求解一个微分方程。的公式，我需要解决的输入如下： solv = DSolve[{ a*u''[y] - b*u[y] == d, u'[0] == 0, u[1] == 0}, u, {y, -1, 1}] 其中使用ExpToTrig和FullSimplify后，我得到我要的答案： (d (-1 + Cosh[(Sqrt[b] y)/Sq

-2热度

1回答

ODE45系统二阶

我有写这个程序，但它不工作 function xprime = exh(t,z) xprime = [z(2); 30*z(3) - 40*z(1); z(4); 30*z(1)+30*z(5)-60*z(3); z(6); 30*z(3)-40*z(5)]; tspan=[0,200]; Z0=[0.05 0.04 0

1热度

1回答

fipy-level set 1D

我试图运行下面的代码来设置1D级别的问题（例如在fipy网页中）。我收到此错误： var.calcDistanceFunction 20行（回溯（最近通话最后一个）：文件 “/用户/ sgowda /文档/ PDE求解器代码/级别设置1D.py C” ）文件“C：\ Users \ sgowda \ AppData \ Local \ Continuum \ Anaconda \ lib \