numerical-methods

1熱

1答えて

配列xとyで表される値の配列でサンプリングされた関数y（x）を考えてみましょう。私はすでに係数とXを持っている場合は：私は別のx値x0を持っている場合、私は今、スプライン y0 = spline(x,y,x0); を使用して、Y（X0）を評価することができ、私はまた pp = spline(x,y); y0 = ppval(pp,x0); 私の質問を書くことができますmy_coefs（s

0熱

1答えて

常微分方程式の非線形システム

最終状態が与えられたシステムの初期状態と、システムのさまざまな側面の変化率を表す非線形微分方程式系を解くことは可能ですか？例えば： V_x' = n * cos(b) V_y' = (n * sin(b)) - (g_0 * (r_0/r)^2) b' = ((g_0 * (r_0/r)^2) * cos(b))/||V|| b_final = 0 V_y_final = 0 V_x

1熱

1答えて

Matlabで数値的に統合されています。なぜ部門が増えると正確な出力が低下するのですか？

私は現在、分散ロードビームを与えられたMatlabの問題に取り組んでおり、数値的に（複合台形ルールを使用して）積分してせん断力と曲げモーメントを求めています。そこから、私は最大の瞬時値とこれが起こる位置を見つける。それから私は再び2倍に分けてせん断力と分散荷重を求めます。私は、数値積分と微分が正しく計算されていることを確認し、確認するために分析的に統合された関数を使用しています。私の問題は、計算

-2熱

1答えて

匿名関数を掛ける

私は匿名関数を格納する変数を持っています。 temp_f = @(x)((1/h^3)*(((xi+2*i)-x)^3-4*((xi+i)-x)^3)) 私は、次の段階で、私は temp2_f = @(x) (temp_f(x) * varA) がどのように私はそれに近づくことができる得ることができるように、「バラ」を掛けたいですか？あなたはvarAが、新しい匿名関数への入力引数になりたい

-1熱

1答えて

x = 0、無限大からf（x）exp（-x）を最適化

f(x)exp(-x) x = 0と無限大の間の堅牢な積分アルゴリズムが必要です。f(x)は正の微分可能な関数です。私は配列xを先験的に知りません（これは私のルーチンの中間出力です）。 xアレイは、通常、ログが等間隔ですが、非常に不規則です。現在、私はSimpsonアルゴリズムを使用しています。私の問題は、x配列によってドメインが非常にアンダーサンプリングされていることが多いことです。これは積

-1熱

1答えて

コードを最適化するにはどうすればよいですか？

私は、連続分数を含む方程式を持っています。私は方程式のパラメータであるすべてのRについてこの方程式の特定の根を見つけたいと思っています。この私のために：入力多項式として、この連分数MATLABへは計算に表現のすべての根をこのpolynomiumを簡素化し、最終的に私は「ルートを検索します興味があります。私が問題にしているのは、高いRは高い精度を必要とするためです私のコンピュータはそ

1熱

1答えて

Runge-Kutta of fourthth

私は、4つのODE：dX/dt = F（X）のシステムを考えてみましょう。ここでXはベクトル（4次元）で、F：R^4→R^4です。 FはvectorDE_total_functionと呼ばれ、私はRK-4を使って解を計算しようとしています。 def solvingDES(): previous_vector = np.array ([theta_1, omega_1, theta_2,

0熱

1答えて

ルンゲクッタ四次法

私は2つの方程式のシステムに対してルンゲクッタ四次法を実装しています。 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[-1], x2[-1])

-1熱

1答えて

Runge Kutta 4をC++でLorenzシステム用

誰かが私のコードでエラーを見つけることができますか？特殊な点では正常に動作しますが、許容値はメソッドに適合しません。私のErrorstepperはかなりシンプルなので、私はそれの問題だとは思わない。助けてください。 #include <math.h> #include <stdlib.h> #include <stdio.h> #include <conio.h> #include <ti

0熱

2答えて

ストリームを等しい数のビンに分割する

理想的には、ハードディスクからデータを何度も読み取らずに、次のようにします。データは大きく、メモリは同時にすべてのデータを保持することはできません。入力はハードディスクからのストリームx[t]です。数字のストリームには、N要素が含まれています。 xのヒストグラムをmビンとすることは可能です。 n個のビンは、ビンによって定義されるが、E < E 、...、< E メートル縁。例えば、e i = <