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 \