android — Giải pháp ODE phi tuyến cho odeint

In lại Tác giả: Vũ trụ không gian Thời gian cập nhật: 2023-11-03 19:42:05 27 4

từ nhập scipy.integrate odeint
nhập numpy dưới dạng np
nhập matplotlib.pyplot dưới dạng plt

def harmonic_motion(X,t,k,m,tau):
    x,v=X
    F=-(k/m)*x


    F_1=tau/v*(np.gradient(x,t)**2) # Dòng này không hoạt động. 

    dx_dt=v
    dv_dt=F-F_1
    dX_dt=[dx_dt,dv_dt]
    trả lại dX_dt


m=1
k=1
tau=0,1
X0=-3
V0=0

t = np.linspace(0, 15, 250)
sol = odeint(harmonic_motion, [X0,V0], t,args=(k,m,tau))
x=sol[:,0]
v=sol[:,1]

plt.plot(t,x,label='x')
plt.plot(t,v,label='v')
plt.legend()
plt.xlabel('t (s)')
plt.ylabel('x,v (m,m/s)')
plt.show()

m*x'' + tau*abs(x')*x' + k*x = 0

def harmonic_motion(X,t,k,m,tau):
    x, v = X
    return v, -(tau*abs(v)*v + k*x)/m

27 4 0

• Vị trí ODE trong hệ thống ODE tạo ra các kết quả khác nhau

  Tôi có đoạn mã sau: def multirk4(funcs, x0, y0, step, xmax): n = len(funcs) table = [[x0] + y0] f1,

• Việc sử dụng một ODE bên trong một ODE khác với các tham số sẽ làm cho mã rất rất chậm

  Tôi đã viết đoạn mã sau. Nó là một ODE có tham số là một ODE khác. Như chúng ta có thể thấy, M(m0,z,b,c) được sử dụng trong một ODE khác, bản thân nó cũng là một hàm ODE. Mã rất chậm, có ai có thể cho tôi lời khuyên về cách cải thiện nó không? Tôi

• scipy.integrate.ode với hai ODE được ghép nối?

  Tôi hiện đang cố gắng sử dụng gói tích hợp.ode của SciPy để giải một cặp ODE bậc nhất được ghép nối: chẳng hạn như phương trình con mồi-kẻ săn mồi Lotka-Volterra.

• Giải quyết hai ODE tách rời trong một vòng lặp bằng python và scipy.integrate.ode

  Tôi gặp sự cố khi sử dụng scipy.integrate.ode để giải quyết hai ODE tách rời rất đơn giản. Ví dụ: đoạn mã đơn giản sau: from scipy.integrate import ode def f(

• Làm cách nào để giải ODE bằng ngưỡng nội bộ?

  Tôi có hàm sau chứa một số ode: myfunction 20) { # đây là ngưỡng nội bộ Y <- 35000!

• Làm thế nào để giải quyết ODE này?

  Bạn đang gặp rắc rối với phương trình này bằng cách sử dụng odeint hoặc Solve_ivp để giải nó. nhập numpy dưới dạng np từ scipy.integrate nhập giải quyết_ivp nhập

• matlab - Giải và vẽ đồ thị ODE từng phần

  Tôi có hàm dφ/dt = γ - F(φ) (trong đó F(φ) -- a là hàm tuần hoàn 2π) và đồ thị của hàm F(φ) . φ(t) Chương trình cho các giá trị khác nhau của γ (γ = 0.

• Julia - Các vấn đề được giải quyết bằng ODE

  Tôi đang cố gắng giải ví dụ chuẩn từ gói DifferentialEquation, theo hướng dẫn trên trang của họ. Đây là ví dụ: sử dụng DifferentialEquations sử dụng Plots fun

• Sao chép mô hình web thực phẩm ODE từ bản thảo

  Tôi đang cố gắng tái tạo mô hình trang web thức ăn trong hồ được đăng ở đây. Mô hình này thể hiện hai chuỗi thức ăn (ven biển và vùng biển) được kết nối bởi loài săn mồi hàng đầu (cá). Tôi đã mã hóa mô hình nhưng khi chạy sau 2-3 bước thời gian, mô hình sẽ tạo ra NaN. TÔI

• Sao chép mô hình web thực phẩm ODE từ bản thảo

  Tôi đang cố gắng tái tạo mô hình trang web thức ăn trong hồ được đăng ở đây. Mô hình này thể hiện hai chuỗi thức ăn (ven biển và vùng biển) được kết nối bởi loài săn mồi hàng đầu (cá). Tôi đã mã hóa mô hình nhưng khi chạy sau 2-3 bước thời gian, mô hình sẽ tạo ra NaN. TÔI

• Giải và vẽ đồ thị ODE với tham số phụ thuộc

  Tôi đang viết một số mã trong đó tôi có hệ phương trình sau đây. Vấn đề là tôi thực sự muốn giải quyết nhiều giá trị k và vẽ biểu đồ mặt phẳng pha/rung động cho mỗi giá trị k. Có ai có thể giúp tôi được không? Đây là những gì tôi biết về bộ giải cho đến nay: import numpy

• matlab - Làm thế nào để giải ODE mà không cần sử dụng các hàm lồng nhau?

  Tôi có một số phương trình vi phân cần giải bằng bộ giải ODE của MATLAB. Mặc dù bản thân các phương trình vi phân khá đơn giản nhưng chúng phụ thuộc vào nhiều "hằng số". Các hằng số này không phổ biến và cần được người gọi cung cấp. Một ví dụ về ODE như vậy là: dx/dt = -

• MATLAB giải ODE trên đa tạp bất biến

  Tôi có một hệ trông giống như dn/dt=f(n,v) dh/dt=g(h,v) Tôi muốn giải phương trình này trên đa tạp F(v,n,h)=0, đa tạp là một hàm phi tuyến trong v. Tôi đã thử sử dụng một cái gì đó như v=fzero(@(x) F(

• ODE với các mục phức tạp

  Tôi đang cố gắng giải một hệ thống ODE có các mục phức tạp. Từ tài liệu GSL có thể thấy rằng nó chỉ chấp nhận các mục thực. Có cách nào để vượt qua phần phức tạp (cách trực tiếp hơn là phân biệt phần thực và phần ảo) không? Nếu không thể, bạn có thể giới thiệu bất kỳ thư viện tốt nào khác cho mục đích này không? câu trả lời hay nhất

• Julia - Bayesian ODE cho Julia

  Tôi đã cố gắng triển khai ODE Bayesian. Trong ngành dầu khí, chúng tôi sử dụng phương trình sau để khớp dữ liệu sản xuất và sau đó đưa ra dự đoán: Phương trình ODE được mô tả là: trong đó 0 Mã ban đầu của tôi: sử dụng DiffEqFlux, OrdinaryDiff

• Kết quả bất ngờ của Python GEKKO ODE

  Tôi đã cố gắng triển khai mô hình ODE để mô phỏng đường truyền tín hiệu insulin, như được mô tả trong tài liệu bổ sung của bài viết này, sử dụng GEKKO của python.

• android — Cách giải phương trình ODE lồng nhau

  Chúng ta có thể sử dụng gói Phương trình vi phân thông thường (ODE) deSolve trong R, tuy nhiên, tôi không thể tìm ra cách giải hai phương trình ODE lồng nhau giả sử `b'(t) = beta - k*b(t) = alpha -b

• ode - Bộ giải Dymola DAE

  Tôi đang cố gắng tìm ra các bước cần thiết để Dymola giải mã Modelica. Qua việc đọc một số tài liệu tham khảo và sách, tôi đã tìm hiểu về Dymola: Convert Modelica code to hybrid DAE (làm phẳng). Thao tác DAE để chuyển đổi

• Điều này có đúng khi mô hình hóa trọng lực như ODE bậc hai không?

  Đây là câu hỏi đầu tiên của tôi ở đây, vì vậy tôi xin lỗi nếu định dạng bị tắt. Tôi đang cố gắng mô hình hóa định luật hấp dẫn của Newton dưới dạng phương trình vi phân bậc hai trong Python, nhưng biểu đồ kết quả không có ý nghĩa. Để tham khảo, đây là phương trình và [đây là kết quả]

• c# - Sympy và ma trận vi phân ODE

  Làm thế nào để giải phương trình vi phân ma trận bằng Sympy? Tôi có một phương trình có dạng y'(t) = A*y(t) + B, trong đó A là ma trận 3x3, y(t) là vectơ 1x3 và B là vectơ 1x3. hơn