別の分布から測定された分布を削除する


8

粒子ビームを多くの粒子の集合として捉えます。2つの独立した確率変数と仮定δ水平位置まで追加X粒子のを:XβδX

X=Xβ+Dxδ

は単純な数値で、ビームダイナミクスの「分散」関数です。)Dx

Iは、ビームプロファイルの水平方向の測定値を有する、縦運動量分布の他の測定、F δ。両方を単一領域に正規化し、それらをXδの確率密度関数の測定値と見なします。fXfδXδ

水平($ f_X $)および縦($ f_ \ delta $)平面で測定されたビームプロファイル

今、私はの分布/プロファイルを決定したいと思いXβ

どうすればよいですか?

最初に考えたデコンボリューションすることであったF D X δ Iは、位置の同じセットに両方のデータセットを補間した後、。残念ながら、私は失敗しました...私はスペクトルに等しいエラー量で終わります、すなわち、どこにも行きません。fXfDxδscipy.signal.deconvolve

2つを畳み込むと、、に拡張がられます。fXfDバツδ

水平および長手方向プロファイル測定の畳み込み

numpy.convolve(f_x, f_Dxdelta, 'same')両方の配列が同じ長さで、同じ位置にある場合)

今は反対を行い、分散部分を「追加」する代わりに「削除」したいと思います。または私は完全に間違った方向に行っていますか?

もう1つの重要な可能性のある情報:はではなく正規分布を持っていると思います。私は、対応する標準偏差を抽出したいから。バツβδバツβfバツ

助けてくれてありがとう、エイドリアン

PS:物理スタック交換フォーラムで同じ質問をしたので、コミュニティに質問することを提案しました:-)/physics/224671/remove-measured-distribution-from-別の配布

回答:


2

デコンボリューションで岩だらけの道を進む代わりに、可能なアプローチは、仮定されたガウス分布を PDF 、上記のコンボリューションにプラグインすることです。結果の曲線は、仮定されたガウス分布の求められた標準偏差である変化させる反復アルゴリズムを使用して、測定されたプロファイルにフィットするように作成できます。バツβDバツδfDバツδσバツβ

この方法で妥当な結果が得られました。それにもかかわらず、私はあなたの提案と他の、おそらくより良いアプローチのためにオープンです... :-)ありがとう。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.