ode

1熱

1答えて

私はODEと協力しており、rk(x, times1, func1, parms)とrk(x, times2, func2, parms)を20回積み重ねる（更新する）ので、新しい実行は前の実行の結果から始まります。前もって感謝します！ここで require(deSolve) func1 <- function(t, x, parms) { with(as.list(c(parms, x

0熱

1答えて

メープルにODEをプロットする

私はメープルの初心者です。このようなODEシステムをどのようにプロットするのか分かりません。ここでのコードは、事前に g := 9.81; m := 0.25e-2; h := 0.73e-3; xi := -.5; yi := -.5; zi := .3; vxi := 8; vyi := 2; vzi := 3; wx := (0*2)*Pi; wy := (13*2)*

0熱

2答えて

関数が別の関数内で実行されるたびに特定のパラメータの値を格納する

Rパッケージの別の関数によって何度も実行されている関数を書きましたが、内部パラメータのうちのいくつか関数が使用されているたびに機能しますが、Rパッケージの関数によって関数が使用されているときにはそれを行うことは不可能です。 library(pracma) model <- function(t,x,params){ y0 <- x[1] y1 <- x[2]

3熱

2答えて

MATLABのオードソルバーの質量行列は何ですか？

質量行列を使用しない場合、ode45のようなオードソルバーは、y '= f（t、y）を解くことができます。しかし、「質量」行列M（t、y）y '= f（t、y）を含む問題については、オードソルバーのオプションがあります。正確に「質量」マトリックスとは何ですか？この用語は、質量 - ばね系振動の質量から来るのだろうか？このドキュメントのサンプルコードは見つかりません。また、y '= f（t、y）

0熱

1答えて

暗黙のDAE質量行列Python

Pythonで解くことになるM * y '= f（y）問題があります。ここで、Mは質量行列、y'は微分、yはベクトルです.y1 、y2等はrの異なる点を指す。 Pythonで同様の問題に質量行列を使用した人はいますか？問題は、rおよびz方向の2D問題です。 r方向は離散化され、問題を1次元問題に還元する。質量行列は、対角線上に1と0の対角行列です。

1熱

1答えて

ODE Iは、式の次のシステムを有するソルバ

する機能を一般： #dY0=/dt = k1*S - k2*Y1*Y0 #dY1/dt = k3*S - k4*Y1 S（t）がtに一の単位によってように= 4、8、12とを増加させる階段関数である。私のコードは次のとおりです： import numpy as np import matplotlib.pyplot as plt from scipy.integrate import

1熱

1答えて

MATLAB - ode45は

に正弦波強制関数を渡す私はMatlabのに新たなんだと、本当に基本にグリップを得るためにも、苦労しています。 Iは、位置と減衰との質量/スプリングシステムの速度及び駆動力を解決する機能、myspringを持っています。私はode45を実行する前に、コマンドウィンドウに係数（c）、および質量（m）を減衰、バネ剛性（k）の値を指定することができます。私は何をすることができませんこと（g = sin(t

0熱

1答えて

Matlabで二次ODEを解く

Matlabの初期条件y（0）= 1とy '（0）= 0を使って方程式y "+ 4y = 0をモデル化しようとしています。上記ODEのための時間の関数としての位置の作業近似： clear syms y(t) %Differential Equation Conversion P = odeToVectorField(diff(y, 2) == -4*y); M = matlabFunc

1熱

1答えて

RパッケージdeSolve：checkFuncのエラー

私はRでかなり新しく、deSolveパッケージで微分方程式を解くことを試みています。私は合計10000のタイムステップと100のタイムステップごとに変数の値を変更する必要があります。私はエラーがどこから来るのか、非常にわからないんだけど Error in checkFunc(Func2, times, y, rho) : The number of derivatives return