問題タブ [orthogonal]

http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.linalg.orth.html

どうすれば調整できますか？予想よりもベースが少なくなっています。Google は、最初の数ページで明らかに有用な結果を提供してくれません。

scipy.special.legendre()最近、 ( scipy のドキュメント)に関する興味深い問題に出くわしました。ルジャンドル多項式はペアごとに直交する必要があります。ただし、ある範囲でそれらを計算しx=[-1,1]、異なる次数の 2 つの多項式のスカラー積を作成すると、常にゼロまたはゼロに近い値が得られるとは限りません。関数の動作を誤解していませんか? 以下に、ルジャンドル多項式の特定のペアのスカラー積を生成する短い例を書きました。

単一の多項式のプロットは、実際にはうまく見えます。

しかし、スカラー積を手動で計算すると、異なる次数の 2 つの凡例多項式を要素ごとに乗算して合計します (500 は正規化用です)...

...出力として次の値を取得します。

最初の多項式とそれ自体のスカラー積は (予想どおり) 1 に等しく、他の結果の半分はほぼゼロですが、次の値がいくつかあり、その10e-3理由はわかりません。私もscipy.special.eval_legendre(n, x)関数を試しました – 同じ結果:-\

scipy.special.legendre()これは機能のバグですか？それとも私は何か間違ったことをしていますか？私は建設的な回答を探しています:-)

乾杯、マーカス

350 のドキュメント スコアがあり、それらをプロットすると、次のような形状になります。

ここに完全な配列を投稿しました（以下のコードのリストにpastebin対応しています）。dataPoints

さて、もともとelbow pointこの曲線の を見つける必要がありましたが、この記事L-shapeのおかげで見つかりました。

次のプロットでは、赤いベクトルpはエルボ ポイントを表します。の正射影に対応するベクトルx=(?,?)上の点 (黄色い星)を見つけたいと思います。bpb

プロットの赤い点は私が得たものです (明らかに間違っています)。次のようにして取得します。

pここで、への射影がb開始点と終了点、つまりsとx(黄色の星)によって定義される場合、 proj_p_onto_b = x - s、したがってx = proj_p_onto_b + s?

私はここで間違いを犯しましたか？

編集: @cxw への回答として、エルボ ポイントを計算するためのコードを次に示します。

編集:プロットのコードは次のとおりです:

八量体のセット (8 文字のセット) に対して直交エンコード関数を実行し、nx160 の数値の行列として返す必要があります (n はデータの八量体の数です)。

直交符号化関数は次のとおりです。

一部の人が尋ねたように、これが機能していない部分ではない場合でも、例を次に示します。

この関数は個々のオクタマーで動作していますが、lapply を使用しようとすると、結果は 160 個の数値ベクトルになり、今回はコードが変更されています (そして無意味です)。

結果は次のようになります。

オルソコード機能は実際に機能しています。私が知る必要があるのは、データフレームからオクタマーを取得し、それらに対して機能を実行する方法です。その結果、次のようなマトリックスになります。

出力は n 行 160 列の行列でなければなりません。実行する必要があるデータでは、結果のマトリックスは 947x160 になるはずです。

何か案は？