runge-kutta

0열

1답변

나는 Runge-Kutta 방법을 연구하고 간단한 진자에 적용하려고합니다. 시간 단계 dt = 0.1 (h = 0.1)을 사용하면 진자는 매주기 1 %의 에너지를 발생시킵니다.이 사이트에서 : link 약 0.01 %의 에너지가 감소합니까? 소프트웨어에는 어디에서 계산 된 값을 비교할 수 있습니까? 인수에 주어진 마지막 샘플에서 파생 모델 상태가 아닌 전역

0열

1답변

룽게 - 쿠타 4 차 방법 누적 오류가 나는 간단한 추적 문제에 대한 수치 솔루션을 얻으려고 노력했다

반복 할 때 모든 내 속도 모듈 감소를 반복 (일정 속도 모듈과 대상 + 로켓을 이동) 조금, 오류를 더하고; 그리고 몇 백 번의 반복 오류가 발생하고 속도가 급격히 떨어지면 그러나 이것은 오일러 방법 (큰 블록 이하의 코드)이 아니며 RK4 방법을 사용할 때만 나타납니다. 나는 오류가 어디 있는지 모르겠습니다하고 무슨 일이 일어나고 왜 당신이 할 수있는,

1열

1답변

제 4 Runge-Kutta를 사용하는 코드가 예상 값을 제공하지 않는 이유는 무엇입니까?

나는 나에게 무슨 문제가 있는지 이해하려고하는데 약간의 문제가있다. 어떤 도움이 도움이 될 것이다. 나는 x가 성장함에 따라 y는 상승으로 그러나, 내 코드가 제공 값이 내 책들 또는 텅스텐의 사람과 일치하지 않는이 간단한 방정식을 해결하고 싶었다. 일반적으로 import matplotlib.pyplot as plt from numpy import exp

3열

1답변

혼란 스 캐터링 시뮬레이션의 값이 기본 경우와 일치하지 않습니다.

스택 오버플로에 대한 첫 번째 게시물은 부드럽게하십시오. I는 사용에 의해 해결된다 움직임 M(dv/dt)=-grad V(r), dr/dt=v, 방정식의 4 차원 시스템에 의해 설명되는 전위 V (R)에 질량 M의 입자의 X 위치, Y 평면을 수행하는 코드를 작성 Runge Kutta 4th order method, 여기서 r = (x, y) 및 v

0열

1답변

ODE 해결사가 함께 실행될 때 느리게 실행 됨

클래스에 대한 ODE 해결사의 사용자 지정 구현이 있습니다. 시간 간격을 dt보다 작게 설정하면 .2보다 작은 문제가 발생합니다. 프로그램이 정지합니다. 그러나 Runge-Kutta 솔버 중 하나를 주석 처리하면 매우 신속하게 실행되며 어느 솔버가 주석 처리되었는지를 전환하여 두 솔버의 솔루션을 얻을 수 있습니다. 나는 이것을 고치는 방법을 궁금해. 한 솔

0열

2답변

ode의 시스템을 해결하기위한 C++ 수치 적 통합 자

최근에 C++을 사용하기 시작했고 사용자 정의 된 ode 's 시스템의 통합을 허용하는 클래스를 만들었습니다. 성능을 비교하기 위해 두 개의 서로 다른 통합자를 사용합니다. 여기서, 코드의 일반 레이아웃은 : class integrators { private: double ti; // initial time double *x

0열

1답변

여러 Timesteps 동안 C#에서 Runge-Kutta를 사용하여 궤도를 계산하는 방법은 무엇입니까?

내가 작업하고있는 프로젝트의 일환으로 4 차 Runge-Kutta 방법을 사용하여 별 주변의 궤도를 시뮬레이션해야합니다. 이 개체의 위치와 속도를 RK4 메서드를 구현할 때마다 .dat 파일에 작성해야합니다. 내 문제는 이후의 시간 간격에 걸쳐 전체 메서드를 반복하는 방법이나 통계 및/또는 인스턴스와 혼동하지 않고 각 단계에서 파일에 쓰는 방법을 해결할 수

0열

1답변

Runge-Kutta 4 차 법

2 방정식의 시스템에 대해 Runge-Kutta 4 차 법을 구현 중입니다. T/H 단계 그래서 H는 세그먼트의 수이다. 그런 다음 def cauchy(f1, f2, x10, x20, T, h): x1 = [x10] x2 = [x20] for i in range(1, h): k11 = f1((i-1)*T/h, x1[-

0열

1답변

Euler/RK4 Orbit 시뮬레이션을위한 python

태양이 (0,0)에있는 8 개의 행성에 대한 orbite 시뮬레이션 코드를 시작했습니다. X와 Y에 가속도가 있습니다. 이제이를 사용하여 X와 Y의 속도, 오일러를 먼저 사용하고, RK4를 사용하면 멈추게됩니다. 어디서부터 시작해야할지 모르겠습니다. numpy에 의해 제공

-1열

1답변

Runge Kutta 4 in C++ for Lorenz system

누군가 내 코드에서 오류를 찾을 수 있습니까? 그것은 특별한 점에서 잘 작동하지만, 허용 오차는 방법에 맞지 않습니다. 내 Errorstepper 꽤 간단합니다, 그래서 나는 그것에 문제가 있다고 생각하지 않습니다. 도와주세요. #include <math.h> #include <stdlib.h> #include <stdio.h> #include <co