ode

2열

2답변

ode (패키지 deSolve)에서 rk4와 rootfun을 사용할 수 있습니까

LV 모델을 기반으로 미분 방정식을 사용하여 먹이 - 먹이 - 포식자 시스템을 모델링하려고합니다. 정밀도를 위해 나는 runge-kutta4 방법을 사용할 필요가있다. 그러나 방정식이 주어지면 일부 인구는 빠르게 부정적으로 변합니다. : eventFunc <- function(t, y, p){ if (y["N1"] < 0) { y["N1"] = 0

2열

1답변

R : 행렬을 입력으로 사용하여 deSolve 패키지로 ODE를 해결합니다.

R에서 상미 분 방정식을 풀려고합니다. 이미 정의한 초기 값과 매개 변수로 행렬이 있습니다. 그러나 그것을 해결하려고하면 행렬 대신에 단일 값을 입력 할 때 나타나지 않는 아래의 오류가 발생합니다. (베타) %의 *의 %의 S에이 오류 : 필요 숫자/복잡한 행렬/벡터 인수 송시를 해결하기위한 내 코드는 내가 왜 이해가 안 S = matrix(c("S1","

0열

1답변

여러 초기 조건을 테스트하기 위해 루프에서 ODE45 실행

나는 ODE45를 사용하여 해결하려고하는 ODE 시스템을 가지고 있습니다. 나는 테스트하고 싶은 약 80 가지 초기 조건을 가지고있다. 아이디어는 각 반복 후에 종속 변수의 최종 값을 얻고 싶습니다. 사실, ODE45는 초기 조건이 바뀔 때마다 다른 크기의 매트릭스를 제공하며, 각 반복마다 새로운 매트릭스를 생성하려고하지만 온라인 사람들은 이것이 이라고 말

0열

1답변

MATLAB에서 매개 변수 변경

하나의 플롯에 ODE에 대한 근사 솔루션, ODE에 대한 정확한 솔루션 및 상대 오류를 표시하려고합니다. 내 함수는이 그래프 중 하나만 출력합니다. 한 번에 세 가지를 그리는 방법과 k에 대한 매개 변수를 변경하는 방법을 잘 모르므로 .02, .1 및 1.5로 k를 표시합니다. 이 기능은 별도의 파일에 있습니다. 대략적인 해결책은 다음과 같습니다. func

0열

1답변

MATLAB에서 플로팅

내부에 추가 플롯을 추가하려는 코드가 있습니다. 내 기능 중 두 가지의 상대적인 오류를 음모 싶습니다. 루프가 어디에 있는지 알 수 없으므로 오류 메시지가 나타나지 않고 세 가지 경우 모두 단일 그래프로 표시됩니다. 또한, 내 상대 오류 코드 문제가 발생할 수있는 몇 가지 실수가있을 수 있습니다. 이 상대 에러 코드 : rel_error = (y_exact

3열

2답변

파이썬 : 행렬로 작은 뉴럴 네트워크 작성하기

파이썬에서 두 개의 뉴런 네트워크를 시뮬레이트하려고합니다. 그것은 각 뉴런에 대해 별도의 방정식을 작성하는 것만 큼 간단하지만 코드를 좀 더 일반화하여 방정식을 반복해서 다시 작성하지 않고도 뉴런의 수를 늘릴 수 있기를 바랍니다. 기본적으로, I는 전압 방정식에서 마지막 항을 통해 함께 결합 된 두 지킨 - 헉슬리 뉴런을 가지고 두 개의 신경망 방정식은 다

0열

1답변

BVP4c가 알 수없는 경계를 해결합니다.

4 개의 시스템을 해결하기 위해 bvp4c를 사용하려고합니다. 문제는 경계 중 하나를 알 수 없다는 것입니다. bvp4c에서 처리 할 수 있습니까? 내 코드에서 L은 내가 알지 못하는 부분이다. 아래에 오류 메시지가 표시됩니다. function mat4bvp L = 8; solinit = bvpinit(linspace(0,L,100),@mat4ini

0열

1답변

van der pol oscillator의 방향 필드를 그리는 방법은 무엇입니까?

나는이 위키 페이지에 표시된 방향으로 필드 및 위상 초상화를 얻기 위해 노력하고 있어요 : Van der Pol oscillator in wikipedia 내 코드 : options = odeset('MaxStep',0.5); temp = inputdlg('Enter mu value'); mu = str2double(temp{1,1}); [t,y] =

6열

1답변

Modelica 이벤트 및 하이브리드 모델링

하이브리드 모델링 (특히 상태 이벤트)의 일반적인 아이디어를 수치 적 관점에서 이해하고 싶습니다. (저는 수학자는 아니지만 :)). model BouncingBall constant Real g=9.81 Real h(start=1); Real v(start=0); equation der(h)=v; der(v)=

0열

1답변

matlab에 ODE 해결사 인해 작은 한 걸음에 영원히 갇혀가

내가 MATLAB과이 ODE 해결사로 비교적 새로운 해요 크기, 아래에있는 내 코드입니다 : main.m format short; tspan=[0 5]; y0=[0.30;-0.30;0;-0;0]; [t,y]=ode23s(@(t,y) pend(t,y),tspan,y0); figure(1) %subplot(2,1,1); plot(t,y(:,1