粒子ビームを多くの粒子の集合として捉えます。2つの独立した確率変数と仮定とδ水平位置まで追加X粒子のを:
(は単純な数値で、ビームダイナミクスの「分散」関数です。)
Iは、ビームプロファイルの水平方向の測定値を有する、縦運動量分布の他の測定、F δ。両方を単一領域に正規化し、それらをXとδの確率密度関数の測定値と見なします。
今、私はの分布/プロファイルを決定したいと思い。
どうすればよいですか?
最初に考えたデコンボリューションすることであったとF D X δ Iは、位置の同じセットに両方のデータセットを補間した後、。残念ながら、私は失敗しました...私はスペクトルに等しいエラー量で終わります、すなわち、どこにも行きません。scipy.signal.deconvolve
2つを畳み込むと、、に拡張がられます。
(numpy.convolve(f_x, f_Dxdelta, 'same')
両方の配列が同じ長さで、同じ位置にある場合)
今は反対を行い、分散部分を「追加」する代わりに「削除」したいと思います。または私は完全に間違った方向に行っていますか?
もう1つの重要な可能性のある情報:はではなく正規分布を持っていると思います。私は、対応する標準偏差を抽出したいから。
助けてくれてありがとう、エイドリアン
PS:物理スタック交換フォーラムで同じ質問をしたので、コミュニティに質問することを提案しました:-)(/physics/224671/remove-measured-distribution-from-別の配布)