differential-equations

    0熱度

    1回答

    ode的主要功能如下所示。 function dxdt = state(t,x,vgth,vgval1,vgval2) vgval=vgval1+vgval2; p=1; k=10^0.7; window1=1-((2*x)-1).^(2*p); dxdt=k*(vgval-vgth+1.2)*window1; end 腳本如下。 step=0.01; t = 0:step:10

    1熱度

    2回答

    我試圖在Sage筆記本上繪製斜坡場和幾個獨特的解決方案,但是我的斜坡場與解決方案不匹配。 c = var('c') y = (ln(x)+c)/x show(y) p1 = plot_slope_field(y, (x,0,5), (c,-5,5)) for i in range(-5,5): p1 = p1 + plot(y(c=i),x,0,5, ymin=-5, ymax

    3熱度

    1回答

    美好的一天。我感興趣的解決形式的問題: x_dot = A + F, 用Matlab。數值解算器(ode23/ode45)的用法看起來很簡單,但在我的情況下,矩陣A和矢量F是狀態相關的。因此,我需要在每個迭代步驟後使用新派生的狀態來更新它們。 實際上可以使用ode23/ode45來完成嗎?我需要遵循另一條路徑嗎? 在此先感謝,任何洞察讚賞。

    1熱度

    1回答

    XMDS是一個軟件包,它允許通過將它們編譯爲C並以本地機器體系結構允許的最大速度運行它們來解決普通,部分和隨機微分方程。 通常情況下,XMDS支持帶有GOMP OpenMP庫的GCC編譯器。 我試圖安裝XMDS,但每當我運行它,我得到以下錯誤: Checking for program mpic++ : /Users/dg813/Downloads/XMDS.app/Conten

    1熱度

    1回答

    我第一定義函數爲DY/DT = Y和dy/dt的= T: function dy=d(y): dy=y end function ddy=dd(t): ddy=t end 然後我分別使用ODE45,: [t,y]=ode45('d',[1 10],1) [t,y]=ode45('dd',[1 10],1) 返回以下錯誤:使用d 太多輸入參數出錯。 我的問題是:

    0熱度

    1回答

    我試圖解決的4個偏微分方程耦合非線性系統,下面的代碼: import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint from odeintw import odeintw # constants global rH, L, q, p, Q, mu rH = 1.0 ; L

    0熱度

    1回答

    我想運行求解&將我的ode解決方案繪製爲以下三個不同的時間,每次只改變一個初始條件。我想解決q.0 = 16,q.0 = 20,q.0 = 18 (這將是很好,如果這都繪製在矩陣..!在此先感謝) rm(list = ls()) par(mfrow = c(1,1)) D = 6 ; A = pi* D^2/4 ; # ODE solver library(deSolve) p

    0熱度

    1回答

    好吧,所以我想寫一個實現四階runge kutta方法的數值近似的微分方程,作爲我的數學課程的一部分但也要學習一些編程,但問題是它在近似的每一步中使用這些係數。所以我們從x值和y值開始,它想要爲後面的x值找到y值,我給它一個步長,通常爲0.1,它在x中向上移動0.1,並給出了一個新的近似值y值,並在這些步驟中的每一個步驟中執行4個近似值,我稱之爲k1,k2,k3和k4。那麼它需要這4個近似值的加權

    0熱度

    1回答

    我對任何編碼都非常新,並且有很多python問題(python 3)。 我正試圖在一篇有關艾滋病的科學論文中重現數字。他們使用了我必須重新創建的模型的計算機模擬。我不斷變化的方面,並不斷收到不同的錯誤信息 - 目前是 TypeError: odeint() missing 1 required positional argument: 't' 這是我使用的代碼,基於我以前的任務。 #Import

    -1熱度

    1回答

    明確解決常微分方程 我的兩個第一階微分如下 y1' = -sin(y0) + (gamma)*cos(y0)sin(beta * x) 和 y0' = y1 其中 (theta)'' = y, (theta)' = y1, theta = y0 我原來的方程爲 (((d^2)*theta)/dt^2)=-sin(theta)+(gamma)cos(theta)sin(Bx) 從t