3

最初の列は 1 から 5000 までのシリアル番号で、2 番目の列は -180 から 180 度の間で変化する 5000 の二面角で構成されます。コマンドを使用して

xmgrace sample1.dat

二面角対シリアル番号のプロットを取得しました。

プロットからヒストグラムを取得するために、Grace のポップアップ ウィンドウで次の操作を行いました。Data -> Trasnformation -> Histograms

と呼ばれる新しくポップアップしたウィンドウでGrace: Histograms、グラフを選択しましたが、グレース チュートリアルのどこかで説明されてsourceいるように、グラフ内で何も選択しませんでした。destination開始を -180、停止を 180、ビン数を 5 に選択しAccept、そのウィンドウでコマンドを実行すると、次のエラー メッセージが表示されました。

[Error] Please select single source and destination graphs

プロットからヒストグラムを取得するのを手伝ってください。xmgraceヒストグラムをプロットするために使用したい。gnuplot などの広く使用されているプロット ツールを使用してプロットする方法を知っています。

4

1 に答える 1

8

あなたはほとんどそれを持っていました!段階的な解決策は次のとおりです。データを 2 列のファイル「sample.dat」に保存します。このようなファイルの内容は次のとおりです。

#n  angle
1   13
2   -162
3   -63
4   53
5   -170
6   80
7   173
8   39
9   113
[ etc ...] 

サンプル ファイルには、ランダムに生成された角度を持つ 40 のエントリがあります。コマンド ラインから xmgrace を呼び出し、データ ファイルをロードします。

xmgrace sample.dat

最初に、角度対 n の折れ線グラフが表示されます。次のように、データをビンに入れ、ヒストグラムをプロットできます。

Data > Transformations > Histograms

これにより、「Grace: Histograms」というタイトルの新しいウィンドウが開きます。データセット「G0.S0」を右クリックして複製します。左側のパネルで、元のデータ セットを 1 回左クリックして選択します。右側のパネルで、新しいデータ セット「G0.S1」を選択します。

ウィンドウの下部で、ビニング パラメータを選択します。

Start at: -180

Stop at: 180 

# of bins: 5. 

をクリックしApplyます。ヒストグラムがプロットされましたが、元のライン プロットは非表示にする必要があります。左側のパネルで、元のデータ セット「G0.S0」を右クリックし、「非表示」を選択します。

これらのすべての手順の後、[ヒストグラム] ウィンドウは次のようになります。

グレース ヒストグラム ウィンドウ

最後に、メインの xmgrace ウィンドウで、プロットを再スケーリングします。角度が 5 つのビンに分割されたヒストグラムが表示されます。角度範囲は x 軸に表示され、各ビン範囲の頻度は y 軸に表示されます。軸にラベルを付けたデフォルトの出力は次のとおりです。

出力画像

すべてのビンからの頻度の合計は、sample.dat のデータ行の数まで加算する必要があります。

于 2015-02-16T15:11:06.007 に答える