differential-equations

1열

1답변

나는이 방정식의 오른쪽 측면 용어를 정의하는 클래스를 가지고 있습니다 : 이 클래스는 rhs 함수와 그 파생어를 계산하는 방법을 제공합니다. 함수는 이런 방식으로 벡터 컨테이너에 저장됩니다 이 클래스는 미분 방정식의 시스템에도 적합합니다. 여기 내가 변경하는 것이 방법이 정의 된 인터페이스 template <typename Type = double> cl

0열

1답변

Simulink 적분기 루프의 초기 조건 t = 0 설정 방법

나는 미분 방정식이 dx/dt = a * x입니다. Matlab Simulink를 사용하여이 방정식을 풀고 Scope 블록을 사용하여 출력해야합니다. 문제는 초기 조건 값 t = 0을 지정하는 방법을 모르겠다는 것입니다. 지금까지 나는이처럼 보이는 솔루션을 만들 처리했다 : 내가 내부 통합이 설정 "초기 조건"하는 possiblity가 있다는 것을 알고을하

-1열

1답변

견고한 방정식을 풀기 위해 Microsoft Oslo Library를 어떻게 사용할 수 있습니까?

Microsoft 오슬로 라이브러리의 어떤 구가 단단한 상미 분 방정식을 해결할 수 있습니까? 전혀 해결할 다른 방법이 있습니까?

2열

1답변

math.js로 Lotka-Volterra 방정식을 풀다

math.js 홈페이지의 미분 방정식을 풀 때 한 가지 예가 있지만 상당히 복잡하며 개인적으로 수학을 적용하기에 충분한 정보를 제공하지 못합니다. js와 비슷한 문제가 있습니다. 그래서 내가하려는 것은 Lotka–Volterra equations for predator-prey simulation을 해결하는 것입니다. 이 방정식이 시스템에 있습니다, 그러나

1열

1답변

int 또는 분수를 예상하고 sympy에서 %와 %를 얻습니다.

미분 방정식 y '= 6e^(2x-y)를 풀 필요가 있습니다. dsolve()를 사용하여 sympy에서이를 수행하려고합니다. sol = dsolve(Derivative(f(x), x) - 6 *(e**(2*x-f(x))), f(x)) 그러나 항상 오류 기대의 int 또는 분수를 얻을 수는 7.38905609893065022723042746058 6 문제

0열

1답변

0 값에 대한 deSolve 오류

"deSolve"패키지를 사용하여 R의 진자를 미분 방정식으로 풀려고했습니다. 진자는 2 차원으로 움직이며, 가장 중요한 힘과 코리올리 스 (coriolisforce) 및 측면으로부터의 바람의 이동을 포함합니다. 이 스크립트입니다 : 내가 원하는대로 install.packages("deSolve") library("deSolve") #parameter

0열

1답변

푸앵카레 단면도 그리기 방법? (Duffing Oscillator)

Duffing 방정식에 대한 간단한 한계 사이클을 성공적으로 보여주는 프로그램을 작성했습니다. 그러나 이제는이 경우 Poincaré 섹션을 그려야합니다. t*omega = 2*pi*n과 같은 일정한 시간 간격으로 위상 공간 다이어그램의 스냅 샷을 찍어서이 작업을 수행해야합니다. 이 경우에 오메가가 1로 설정되어 있기 때문에 이것은 단지 t = 2*pi*n

1열

1답변

왜 이러한 미분 방정식에서 비슷한 결과가 발생합니까?

2 차 미분 방정식을 푸는 데있어 MATLAB 코드에 문제가 있는지 궁금합니다. ,가 >> tspan = [0 40]; y0 = [1,2]; [t,y]=ode45(@order2,tspan,y0); plot(t,y(:,1)) 문제는이다. function [ dydt ] = order2(t,y) dydt = zeros(size(y)); c=2.5

0열

1답변

DAE의 암시 적/역방향 오일러

암울한 오일러에 관한 질문이 있습니다. 내 암시 적 오일러 메서드를 계산하는 방법을 알고 있지만 내 문제는 DAE (차동 대수 방정식)에서이를 사용하는 방법입니다. 원본 DAE에서 인덱스 축소를 적용한 후 올바른 해결책을 얻었으므로 ODE를 얻은 다음 내 암시 적 오일러를 적용했습니다. 그러나 DAE에 암시 적 오일러를 배포하는 작업이있었습니다. 누구든지

0열

1답변

odeint를 사용하여 분기를 결정합니다.

분기 다이어그램의 데이터, 즉 [N]의 정상 상태 (nullcline) 값을 제공하는 trigger [A] (Aspace의 linspace로 지정) 값을 얻으려고합니다. 미분 방정식의 세트에서 : 방정식의 집합을 다시 배열로 nullclines을 얻기 위해 시도하는 데 시간이 너무 오래있는 그대로 , 나는 주어진 시간이 지남에 통합 할 수 odeint를 사