differential-equations

    -1熱度

    1回答

    因此,我試圖編寫一個代碼來解開開普勒勢阱中的(我們稱之爲)軌道的微分方程,當你進行數學計算時,你將得到一個數學公式微分方程,看起來像這樣: d^2U/d(FI)^ 2 + U - 米/ M^2 = 0 其中u = 1/R ,我們最終尋找r(fi) 和現在我嘗試使用數值方法來解決它,第一I所述的du/DFI = Y 然後definig的函數(I採取了一些任意的M和M) def func(y,fi):

    0熱度

    1回答

    我不理解如何以下第二階方程轉換爲一階等式: X「」(T)= 686 * -0.5 X'(T^2)+ 10(X(t)+ 150) 我使用的替換方法,但我不知道如何對x(t)的轉換: X'(T)= Y(t)的, (t)= x(t)= y^2(t)/ 2, Y'(T)= 686 -0.5 * Y(T^2)+ 10(Y^2(T)/ 2 + 150)

    0熱度

    1回答

    假設我們有以下功能: function f=lorenz(t,x,a,b,c) % solve differential equation like this %dx/dt=a*(y-x) %dy/dt=-x*z+b*x-y %dz/dt=xy-c*z/3 f=zeros(3,1);% preallocate result f(1)=a*

    1熱度

    2回答

    我有一大組ODE代表生物系統中的化學通量。分子正在反應,隔離和循環。我試圖讓這個功能以這樣一種方式發揮作用,讓我知道在一系列條件下可以生產出多少某些產品。 我使用這些包 import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt import math import pylab

    0熱度

    1回答

    我有如下方程組的 是否有可能利用MATLAB的ODE45來解決這個問題?我知道我需要將二階方程轉換爲兩個一階方程,但是我的困惑來自這個術語,它是的衍生物的產物,和theta。 P.S.貝塔僅僅是一個恆定 我已定義的方程組作爲函數爲: function dy = pend(t,y) beta = 1; dy(1) = y(2); dy(2) = -1/(1+y(3))*sin(y(1))-2/

    -1熱度

    1回答

    誰能幫我做使用MATLAB nemerical分化爲這個圖像的功能: 部分分式使用分化 其中,函數F(伽馬)由式(11)給出在下面的圖片: F(伽馬) 我從來沒有做過小號在Matlab上就像這樣,所以我很感激它,如果你能給我一個例子或者試着給我這個方程的初始步驟。

    1熱度

    1回答

    我寫關於解決在使用拍攝二分法下邊界值問題的程序: y''-y+x=0, y(0)=y(1)=0. 我第一次將它轉換爲一階微分方程的系統,設置 y'=z 然後我讓dydt代表矢量(Y 'Z'),並拿出腳本文件: function dydt=shoot(t,y) dydt=[y(2);y(1)-t] end 有了這個,然後我想出了下面的代碼: clear clc a=0; b=1

    2熱度

    1回答

    我創建了一個阻尼擺的相平面的gif。 爲此,我使用包裝DifferentialEquations.jl進行整合,將Plots.jl與繪圖後端和ImageMagick.jl一起用於繪圖和動畫。要做到這一點,我將函數集成到不同的初始值中,並且每個函數都進行了很多次,因此我得到了一些小片段,它們以gif的形式一起形成最終的繪圖。代碼如下: using DifferentialEquations usi

    2熱度

    2回答

    我正在寫一個簡單的python腳本來解決使用歐拉方法的微分方程,現在我必須改變每次我想要解決一個新的方程時的源代碼。 (安全是因爲這是一個個人項目的那一刻任何關心的問題。) 我的問題: 有沒有可能來輸入一個公式作爲用戶輸入,並使其一個變量?例如,使用e=input("enter equation")您輸入y/x,並使e變量供以後使用? #I would like to be able to typ

    0熱度

    1回答

    假設我有BVP NY ^(N-1)-2 *披* Y'= 0,其中n和phi中給出。 Y(0)= Y(5)= 4 爲了使用BVP4c,我首先轉換成一階微分的系統中,這方程 說Y_1 = y和Y_2 = Y」 然後我們得到系統 y'_1 = Y2, y'_2 =(NX ^(N-1))/(2 * PHI) 我現在寫一個matlab函數來做到這這樣做,即 function sol=bvpsolve(n,p