Статья Автор: Брагина Екатерина

Самматив Лебедь X-1

import matplotlib.pyplot as plt
import csv

X = []
Y = []

with open('cygX-1.csv', 'r') as datafile:
    plotting = csv.reader(datafile, delimiter=',')

    for ROWS in plotting:
        X.append(float(ROWS[0]))
        Y.append(float(ROWS[1]))

plt.xlabel('Юлианская дата')
plt.ylabel('Видимая звездная величина')
plt.scatter(X, Y)
plt.show()

import matplotlib.pyplot as plt
import csv
import pandas as pd
import numpy as np


X = []
Y = []

with open('cygX-1.csv', 'r') as datafile:
    plotting = csv.reader(datafile, delimiter=',')

    for ROWS in plotting:
        X.append(float(ROWS[0]))
        Y.append(float(ROWS[1]))

for i in range(len(X)):
   X[i-1] = ((X[i-1]-41162.851)/5.5993863)%1

plt.xlabel('Фаза')
plt.ylabel('Видимая звездная величина')
plt.scatter(X, Y)
plt.show()

import matplotlib.pyplot as plt
import csv
import pandas as pd
import numpy as np


X = []
Y = []

with open('cygX-1_2.csv', 'r') as datafile:
    plotting = csv.reader(datafile, delimiter=';')

    for ROWS in plotting:
        X.append(float(ROWS[0]))
        Y.append(float(ROWS[1]))

for i in range(len(X)):
   X[i-1] = ((X[i-1]-41162.851)/5.5993863)%1
for t in range(len(Y)):
    Y[t-1] = ((Y[t-1]-6678.15)/6678.15)*299792.458

plt.xlabel('Фаза')
plt.ylabel('Лучевая скорость')
plt.scatter(X, Y)
plt.show()

maxY = Y[0]
for k in range(len(Y)):
    if Y[k]>maxY:
        maxY=Y[k]
#print(maxY)

minY = Y[0]
for n in range(len(Y)):
    if Y[n]<minY:
        minY=Y[n]
#print(minY)

amplitude = (maxY-minY)/2
print(amplitude)

Прикрепленные файлы
cygX-1.csv
cygX-1_2.csv
Задача по анализу данных, используя библиотеки Pandas и Matplotlib.pdf
самматив лебедь х1.pdf

Загрузка...

Чтобы оставить комментарий, необходимо авторизоваться

💬

Пока нет комментариев. Будьте первым!

Печать