問題タブ [finite-difference]

diff(x, lag = L, differences = K)RではラグL K次数の前方差分が得られることを知っています。

L Kラグ- 次の後方差分と中心差分も計算する同等の関数はありますか (理想的xには、ベクトルまたは行列のいずれかにすることができます)。私は CRAN を検索しましたが、何も見つからないように見えますが、関数が既にどこかに存在しているはずのような些細なことのように思えます。誰か考えはありますか？

有限差分法を使用した 2D および 3D の速度/圧力解析データがあります。したがって、メッシュの各ノードの座標と、メッシュの各ノードの解 (速度と圧力) の値があります。

通常、有限要素法には paraview を使用しますが、三角形があり、メッシュの各要素の数値解がわかっているため、paraview を使用するのが自然です。

しかし、有限差分では、点の解しか知らないので、メッシュの各正方形の解を補間する方が良い (より正確な) かどうか (たとえば、各正方形の区分線形関数)、またはメッシュの各要素の数値解、プロットは区分定数解です。

私は有限差分で十分な経験を持っていないので、質問したいのですが、有限差分法の数値結果を示すのに適した方法はどれですか? おそらく、paraview はこれに適したソフトウェア オプションではありませんが、私が知っているソフトウェアです。

アルゴリズムを読んでも、非線形有限差分法のコーディングに問題があります。これは、本Numerical analysis (Burden RL, Faires JD, 2010)、セクション 11.4、694 ページの例です。

以下は本から得たアルゴリズムです。これをコーディングする必要があります：

非線形有限差分

私は線形問題に対して同じものをコーディングすることに成功しました（上記の本のセクション11.3）：

非線形でも同じだと思っていたのですが、そうではありません。助けてください、本当に感謝しています。どうもありがとう！