ВНИМАНИЕ! На форуме начался конкурс - астрофотография месяца ЯНВАРЬ!
0 Пользователей и 1 Гость просматривают эту тему.
Потому, что я не знаю, что должно произойти, чтобы она не вспыхнула в этом году, а к примеру в 2027 году или ещё позже.
Может уже накопила массу и вспыхнет в 2026 году, но может ещё не накопила, а накопит к 2046 году? Кто ж сейчас это знает? Интрига.
а накопит к 2046 году? Кто ж сейчас это знает? Интрига.
Было бы накоплено периодов пять, тогда можно было бы примерно рассчитать.
import numpy as npt = np.array([ 1217.8, 1787.9, 1866.4, 1946.1 ])n = np.array([ -9., -2., -1., 0.])a = np.array([len(n)*[1.], n]).TQ, R = np.linalg.qr(a)invR = np.linalg.inv(R)sol = invR@(Q.T@t)epoch, period = solprint(f"{epoch, period=}")invAtA = invR@invR.T # np.linalg.inv(a.T@a)cov = (np.linalg.norm(t - a@sol)**2/(a.shape[0] - a.shape[1]))*invAtAprint(f"{cov=}")print(f"{epoch + period*(n[-1] + 1)=}")print(f"{np.sqrt(cov[0][0] + cov[1][1])=}")
epoch, period=(1947.6820000000005, 81.04400000000017)cov=array([[1.744338, 0.243396], [0.243396, 0.081132]])epoch + period*(n[-1] + 1)=2028.7260000000006np.sqrt(cov[0][0] + cov[1][1])=1.3510995522167293
t[-1] + period, np.sqrt(cov[1][1])
(2027.144, 0.2848367953758976)
Но ведь так можно сказать, что она вообще только в следующем столетии вспыхнет.
Тем более, уже был в 2023 году провал её яркости, который сигнализирует о приближающейся вспышки.