プログラム一覧に戻る

import numpy as np  # データ読み込み機能などを使えるようにする
import matplotlib.pyplot as plt  # グラフプロット機能を使えるようにする

data = np.loadtxt("d2r2.txt")  # データファイルの読み込み

plt.xlabel("time t")  # 横軸の名前
plt.ylabel("position x")  # 縦軸の名前

# 0列目 vs 1列目 を線なし・黒色・丸印でプロット
plt.plot(data[:, 0], data[:, 1], linestyle="", color="black", marker=".")

# 解析解データの作成とプロット
g = 1.5  # 重力加速度
x0 = data[0, 1]  # 初期位置 ＝ ファイルの0行1列目のデータ
t = np.linspace(0, 10, 200)  # 0から10までを200等分した値の入った配列を用意
x_analytic = x0 + 0.5 * g * (t**2)  # 各tに対する解析解の値が入った配列を用意
plt.plot(t, x_analytic, color="red")  # t vs x_analytic を赤色の線でプロット

plt.show()  # グラフを表示