Skip to content

Instantly share code, notes, and snippets.

@leliel12
Created September 13, 2017 04:18
Show Gist options
  • Save leliel12/a3ee1c0de5a7f2cf72100c99c220f05c to your computer and use it in GitHub Desktop.
Save leliel12/a3ee1c0de5a7f2cf72100c99c220f05c to your computer and use it in GitHub Desktop.
#!/usr/bin/env python
# -*- coding: utf-8 -*-
import numpy as np
from FATS.FeatureFunctionLib import StructureFunction_index_21
time = np.array([
55309.23823026, 55309.23853661, 55309.23884425, 55451.07538252,
55451.07585872, 55451.07632649, 55780.20269214, 55780.20306809,
55780.20344416, 55806.09274343, 55806.09315782, 55806.09358313,
55824.08408433, 55824.08447836, 55824.08485471, 55839.02564351,
55839.02604943, 55839.02644263, 55847.04765615, 55847.04812098,
55847.04861453, 56004.29941575, 56004.2997771 , 56004.3001366 ,
56020.28833589, 56020.28874771, 56020.28916008, 56095.11881866,
56095.11924039, 56095.11965876, 56102.07178317, 56102.07220034,
56102.07264671, 56114.08291562, 56114.08329834, 56114.08368282,
56115.09852871, 56115.09893635, 56115.09935218, 56126.09510226,
56126.0954939 , 56126.09590994, 56127.03591857, 56127.03629479,
56127.03668146, 56128.17309234, 56128.17350684, 56128.17389667,
56129.16427653, 56129.1646718 , 56129.16505505, 56130.19315306,
56130.19355928, 56130.19397699, 56140.12788675, 56140.12828595,
56140.12869109, 56141.10950439, 56141.10992254, 56141.1103592 ,
56144.01186111, 56144.01223689, 56144.0126127 , 56158.18208834,
56158.18247022, 56158.18284837, 56195.04119073, 56195.04160086,
56195.04199084, 56196.02600544, 56196.02639798, 56196.0268009 ,
56197.03482072, 56197.03523538, 56197.03562637, 56198.06336708,
56198.06376135, 56198.06414412, 56366.32925902, 56366.32964511,
56366.33002107, 56367.31829269, 56367.31867585, 56367.3190516 ,
56373.37865335, 56373.37906072, 56373.37948436, 56422.40322124,
56422.4036387 , 56422.40404936, 56461.16498389, 56461.16537816,
56461.1657366 , 56484.17994532, 56484.18032831, 56484.18072449,
56487.20889724, 56487.20928786, 56487.20968486, 56506.99674115,
56506.99716518, 56506.99761709, 56507.06055399, 56507.06094397,
56507.06135206, 56509.04455691, 56509.04497718, 56509.04536527,
56527.13820874, 56527.13863534, 56527.13906212, 56534.10144673,
56534.10183848, 56534.10222238, 56559.00205499, 56559.00248304,
56559.00287138, 56560.04968177, 56560.05008353, 56560.05046433,
56572.03685178, 56572.03720617, 56572.03756689, 56582.98766069,
56582.98801412, 56582.98839076, 56717.39100094, 56717.39137697,
56717.39175329, 56770.31076709, 56770.31116901, 56770.3115758 ,
56779.42289437, 56779.42328744, 56779.42367959, 56798.31350297,
56798.31390545, 56798.31432163, 56799.12873954, 56799.12913413,
56799.12953574, 56799.15304962, 56799.15342448, 56799.1537996 ,
56810.2125785 , 56810.21472037, 56810.21515343, 56817.07577736,
56817.07616034, 56817.07656084, 57113.40952299, 57113.409948 ,
57113.41037473, 57226.1351118 , 57226.13554475, 57226.1373534 ,
57241.17917508, 57241.17948152, 57241.17978775, 57253.10787655,
57253.10825253, 57253.10861606, 57268.01650433, 57268.01687509,
57268.01725025, 57273.05390254, 57273.05429794, 57273.05469298,
57274.02238872, 57274.0227351 , 57274.02305764])
magnitude = np.array([
12.16 , 12.131, 12.151, 12.151, 12.117, 12.129, 12.158,
12.142, 12.151, 12.187, 12.194, 12.178, 12.185, 12.163,
12.154, 12.223, 12.191, 12.265, 12.16 , 12.102, 12.133,
12.237, 12.246, 12.158, 12.248, 12.205, 12.21 , 12.174,
12.152, 12.16 , 12.239, 12.168, 12.161, 12.163, 12.138,
12.143, 12.176, 12.152, 12.145, 12.159, 12.177, 12.14 ,
12.315, 12.203, 12.209, 12.195, 12.137, 12.15 , 12.15 ,
12.144, 12.14 , 12.183, 12.149, 12.135, 12.159, 12.148,
12.143, 12.178, 12.155, 12.146, 12.168, 12.157, 12.158,
12.167, 12.16 , 12.151, 12.186, 12.168, 12.168, 12.156,
12.17 , 12.159, 12.201, 12.183, 12.156, 12.14 , 12.127,
12.16 , 12.155, 12.136, 12.155, 12.163, 12.16 , 12.143,
12.161, 12.166, 12.19 , 12.169, 12.169, 12.163, 12.161,
12.212, 12.152, 12.256, 12.191, 12.182, 12.217, 12.255,
12.143, 12.297, 12.264, 12.21 , 12.156, 12.141, 12.127,
12.244, 12.212, 12.215, 12.165, 12.159, 12.135, 12.149,
12.143, 12.14 , 12.156, 12.152, 12.153, 12.188, 12.202,
12.195, 12.171, 12.164, 12.208, 12.18 , 12.136, 12.166,
12.19 , 12.212, 12.165, 12.193, 12.224, 12.143, 12.158,
12.153, 12.172, 12.149, 12.159, 12.157, 12.162, 12.151,
12.143, 12.166, 12.152, 12.166, 12.249, 12.182, 12.235,
12.173, 12.168, 12.163, 12.326, 12.172, 12.144, 12.152,
12.148, 12.146, 12.151, 12.152, 12.156, 12.145, 12.14 ,
12.16 , 12.16 , 12.154, 12.155, 12.176, 12.141, 12.159,
12.15 , 12.145, 12.151])
sfunction = StructureFunction_index_21()
sfunction.fit((magnitude, time))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment