edb_data_calculate_cjjx.go 15 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452
  1. package models
  2. import (
  3. "errors"
  4. "eta_gn/eta_index_lib/global"
  5. "eta_gn/eta_index_lib/utils"
  6. "fmt"
  7. "github.com/nosixtools/solarlunar"
  8. "github.com/shopspring/decimal"
  9. "gorm.io/gorm"
  10. "strconv"
  11. "strings"
  12. "time"
  13. )
  14. // AddCalculateCjjx 超季节性
  15. func AddCalculateCjjx(req *EdbInfoCalculateBatchSaveReq, fromEdbInfo *EdbInfo, edbCode, uniqueCode string, sysUserId int, sysUserRealName string, formulaInt int) (edbInfo *EdbInfo, err error) {
  16. to := global.DEFAULT_DmSQL.Begin()
  17. defer func() {
  18. if err != nil {
  19. to.Rollback()
  20. } else {
  21. to.Commit()
  22. }
  23. }()
  24. fmt.Println("req.EdbInfoId:", req.EdbInfoId)
  25. if req.EdbInfoId <= 0 {
  26. edbInfo = new(EdbInfo)
  27. edbInfo.Source = utils.DATA_SOURCE_CALCULATE_CJJX
  28. edbInfo.SourceName = "超季节性"
  29. edbInfo.EdbCode = edbCode
  30. edbInfo.EdbName = req.EdbName
  31. edbInfo.EdbNameSource = req.EdbName
  32. edbInfo.Frequency = req.Frequency
  33. edbInfo.Unit = req.Unit
  34. edbInfo.ClassifyId = req.ClassifyId
  35. edbInfo.SysUserId = sysUserId
  36. edbInfo.SysUserRealName = sysUserRealName
  37. edbInfo.CreateTime = time.Now()
  38. edbInfo.ModifyTime = time.Now()
  39. edbInfo.UniqueCode = uniqueCode
  40. edbInfo.CalculateFormula = req.Formula
  41. edbInfo.Calendar = req.Calendar
  42. edbInfo.EdbNameEn = req.EdbName
  43. edbInfo.UnitEn = req.Unit
  44. edbInfo.EdbType = 2
  45. edbInfo.Sort = GetAddEdbMaxSortByClassifyId(req.ClassifyId, utils.EDB_INFO_TYPE)
  46. tmpErr := to.Create(edbInfo).Error
  47. if tmpErr != nil {
  48. err = tmpErr
  49. return
  50. }
  51. //关联关系
  52. {
  53. calculateMappingItem := new(EdbInfoCalculateMapping)
  54. calculateMappingItem.CreateTime = time.Now()
  55. calculateMappingItem.ModifyTime = time.Now()
  56. calculateMappingItem.Sort = 1
  57. calculateMappingItem.EdbCode = edbCode
  58. calculateMappingItem.EdbInfoId = edbInfo.EdbInfoId
  59. calculateMappingItem.FromEdbInfoId = fromEdbInfo.EdbInfoId
  60. calculateMappingItem.FromEdbCode = fromEdbInfo.EdbCode
  61. calculateMappingItem.FromEdbName = fromEdbInfo.EdbName
  62. calculateMappingItem.FromSource = fromEdbInfo.Source
  63. calculateMappingItem.FromSourceName = fromEdbInfo.SourceName
  64. calculateMappingItem.FromTag = ""
  65. calculateMappingItem.Source = edbInfo.Source
  66. calculateMappingItem.SourceName = edbInfo.SourceName
  67. calculateMappingItem.FromSubSource = edbInfo.SubSource
  68. err = to.Create(calculateMappingItem).Error
  69. if err != nil {
  70. return
  71. }
  72. }
  73. } else {
  74. edbInfo, err = GetEdbInfoById(req.EdbInfoId)
  75. if err != nil {
  76. return
  77. }
  78. dataTableName := GetEdbDataTableName(utils.DATA_SOURCE_CALCULATE_CJJX, utils.DATA_SUB_SOURCE_EDB)
  79. fmt.Println("dataTableName:", dataTableName)
  80. deleteSql := ` DELETE FROM %s WHERE edb_info_id=? `
  81. deleteSql = fmt.Sprintf(deleteSql, dataTableName)
  82. err = to.Exec(deleteSql, req.EdbInfoId).Error
  83. if err != nil {
  84. return
  85. }
  86. }
  87. //计算数据
  88. err = refreshAllCalculateCjjx(to, edbInfo.EdbInfoId, edbInfo.Source, edbInfo.SubSource, fromEdbInfo, edbInfo.EdbCode, "", "", edbInfo.Calendar, formulaInt)
  89. return
  90. }
  91. // EditCalculateCjjx 超季节性
  92. func EditCalculateCjjx(req *EdbInfoCalculateBatchEditReq, edbInfo, fromEdbInfo *EdbInfo, formulaInt int) (err error) {
  93. to := global.DEFAULT_DmSQL.Begin()
  94. defer func() {
  95. if err != nil {
  96. to.Rollback()
  97. } else {
  98. to.Commit()
  99. }
  100. }()
  101. oldCalculateFormula := edbInfo.CalculateFormula //原先的n值
  102. oldCalendar := edbInfo.Calendar //原先的公历、农历
  103. edbInfo, err = GetEdbInfoById(req.EdbInfoId)
  104. if err != nil {
  105. return
  106. }
  107. //修改指标信息
  108. edbInfo.EdbName = req.EdbName
  109. edbInfo.EdbNameSource = req.EdbName
  110. edbInfo.Frequency = req.Frequency
  111. edbInfo.Unit = req.Unit
  112. edbInfo.ClassifyId = req.ClassifyId
  113. edbInfo.CalculateFormula = req.Formula
  114. edbInfo.Calendar = req.Calendar
  115. edbInfo.EdbNameEn = req.EdbNameEn
  116. edbInfo.UnitEn = req.UnitEn
  117. edbInfo.ModifyTime = time.Now()
  118. err = to.Model(edbInfo).Select([]string{"EdbName", "EdbNameSource", "Frequency", "Unit", "ClassifyId", "CalculateFormula", "Calendar", "ModifyTime", "EdbNameEn", "UnitEn"}).Updates(edbInfo).Error
  119. if err != nil {
  120. return
  121. }
  122. //判断计算指标是否被更换
  123. var existCondition string
  124. var existPars []interface{}
  125. existCondition += " AND edb_info_id=? AND from_edb_info_id=? "
  126. existPars = append(existPars, edbInfo.EdbInfoId, req.FromEdbInfoId)
  127. count, err := GetEdbInfoCalculateCountByCondition(existCondition, existPars)
  128. if err != nil {
  129. err = errors.New("判断指标是否改变失败,Err:" + err.Error())
  130. return
  131. }
  132. if count > 0 && oldCalculateFormula == req.Formula && oldCalendar == req.Calendar { // 指标未被替换,同时N值未修改,同时公历/农历未变更,无需重新计算
  133. return
  134. }
  135. // 指标被替换,或者N值未修改,那么需要重新计算数据
  136. //基础指标被替换了,需要删除原先的 计算指标关联的,基础指标的关联关系
  137. if count <= 0 {
  138. // 需要删除原先的 计算指标关联的,基础指标的关联关系
  139. sql := ` DELETE FROM edb_info_calculate_mapping WHERE edb_info_id = ? `
  140. err = to.Exec(sql, edbInfo.EdbInfoId).Error
  141. if err != nil {
  142. return
  143. }
  144. // 添加新的关联关系
  145. {
  146. calculateMappingItem := &EdbInfoCalculateMapping{
  147. EdbInfoCalculateMappingId: 0,
  148. EdbInfoId: edbInfo.EdbInfoId,
  149. Source: utils.DATA_SOURCE_CALCULATE_CJJX,
  150. SourceName: "超季节性",
  151. EdbCode: edbInfo.EdbCode,
  152. FromEdbInfoId: fromEdbInfo.EdbInfoId,
  153. FromEdbCode: fromEdbInfo.EdbCode,
  154. FromEdbName: fromEdbInfo.EdbName,
  155. FromSource: fromEdbInfo.Source,
  156. FromSourceName: fromEdbInfo.SourceName,
  157. FromTag: "",
  158. Sort: 1,
  159. CreateTime: time.Now(),
  160. ModifyTime: time.Now(),
  161. FromSubSource: fromEdbInfo.SubSource,
  162. }
  163. err = to.Create(calculateMappingItem).Error
  164. if err != nil {
  165. return
  166. }
  167. }
  168. }
  169. //清空原有数据
  170. tableName := GetEdbDataTableName(edbInfo.Source, edbInfo.SubSource)
  171. sql := fmt.Sprintf(` DELETE FROM %s WHERE edb_info_id = ? `, tableName)
  172. err = to.Exec(sql, edbInfo.EdbInfoId).Error
  173. if err != nil {
  174. return
  175. }
  176. //计算数据
  177. err = refreshAllCalculateCjjx(to, edbInfo.EdbInfoId, edbInfo.Source, edbInfo.SubSource, fromEdbInfo, edbInfo.EdbCode, "", "", edbInfo.Calendar, formulaInt)
  178. return
  179. }
  180. // RefreshAllCalculateCjjx 刷新全部超季节性数据
  181. func RefreshAllCalculateCjjx(edbInfoId, source, subSource int, fromEdbInfo *EdbInfo, edbCode, startDate, endDate, calendar string, formulaInt int) (err error) {
  182. to := global.DEFAULT_DmSQL.Begin()
  183. defer func() {
  184. if err != nil {
  185. to.Rollback()
  186. } else {
  187. to.Commit()
  188. }
  189. }()
  190. // 重新计算
  191. err = refreshAllCalculateCjjx(to, edbInfoId, source, subSource, fromEdbInfo, edbCode, startDate, endDate, calendar, formulaInt)
  192. return
  193. }
  194. // refreshAllCalculateCjjx 刷新全部超季节性数据
  195. func refreshAllCalculateCjjx(to *gorm.DB, edbInfoId, source, subSource int, fromEdbInfo *EdbInfo, edbCode, startDate, endDate, calendar string, formulaInt int) (err error) {
  196. if err != nil {
  197. return
  198. }
  199. edbInfoIdStr := strconv.Itoa(edbInfoId)
  200. //计算数据
  201. dataList, err := GetEdbDataListAllByTo(to, fromEdbInfo.Source, fromEdbInfo.SubSource, FindEdbDataListAllCond{
  202. EdbInfoId: fromEdbInfo.EdbInfoId,
  203. }, 1)
  204. if err != nil {
  205. return err
  206. }
  207. var dateArr []string
  208. dataMap := make(map[string]*EdbInfoSearchData)
  209. for _, v := range dataList {
  210. dateArr = append(dateArr, v.DataTime)
  211. dataMap[v.DataTime] = v
  212. }
  213. // 通过插值法补全所有数据(包含周末)
  214. handleDataMap := make(map[string]float64)
  215. _, err = HandleDataByLinearRegression(dataList, handleDataMap)
  216. if err != nil {
  217. return
  218. }
  219. //获取指标所有数据
  220. existDataList := make([]*EdbData, 0)
  221. dataTableName := GetEdbDataTableName(source, subSource)
  222. sql := `SELECT * FROM %s WHERE edb_info_id=? `
  223. sql = fmt.Sprintf(sql, dataTableName)
  224. err = to.Raw(sql, edbInfoId).Scan(&existDataList).Error
  225. if err != nil {
  226. return err
  227. }
  228. existDataMap := make(map[string]string)
  229. removeDataTimeMap := make(map[string]int) //需要移除的日期数据
  230. for _, v := range existDataList {
  231. existDataMap[edbCode+v.DataTime] = v.Value
  232. removeDataTimeMap[v.DataTime] = 1
  233. }
  234. addSql := ` INSERT INTO edb_data_calculate_cjjx(edb_info_id,edb_code,data_time,value,create_time,modify_time,data_timestamp) values `
  235. var isAdd bool
  236. //日度/周度/季度/月度
  237. //isCompatibility := false //是否向上下兼容35天
  238. //if utils.InArrayByStr([]string{"日度", "周度", "季度", "月度"}, fromEdbInfo.Frequency) {
  239. // isCompatibility = true
  240. //}
  241. // 每个年份的日期数据需要平移的天数
  242. moveDayMap := make(map[int]int, 0) // 每个年份的春节公历
  243. var lastDataDay time.Time
  244. if len(dataList) > 0 {
  245. lastDataDay, _ = time.ParseInLocation(utils.FormatDate, dataList[0].DataTime, time.Local)
  246. }
  247. for _, av := range dateArr {
  248. // 如果遇到闰二月,如2.29,去掉该天数据
  249. if strings.Contains(av, "02-29") {
  250. continue
  251. }
  252. currentDate, tmpErr := time.ParseInLocation(utils.FormatDate, av, time.Local)
  253. if tmpErr != nil {
  254. err = tmpErr
  255. return
  256. }
  257. //农历的超季节性运算,只计算11月--次年5月,分段计算,与数据区间和N数值有关
  258. if calendar == "农历" && currentDate.Month() > 5 && currentDate.Month() < 11 {
  259. continue
  260. }
  261. currentItem, ok := dataMap[av]
  262. // 找不到数据就退出当前循环,进入下一循环
  263. if !ok {
  264. continue
  265. }
  266. pastValueList := make([]float64, 0) // 过去几期的数据
  267. //当前日期
  268. pastValueList = append(pastValueList, currentItem.Value)
  269. for i := 1; i < formulaInt; i++ {
  270. //前几年当天公历的日期
  271. historyPreDate := currentDate.AddDate(-i, 0, 0)
  272. moveDay := 0
  273. if calendar == "农历" {
  274. if tmpMoveDay, ok := moveDayMap[historyPreDate.Year()]; !ok {
  275. moveDay, err = getMoveDay(lastDataDay, historyPreDate)
  276. if err != nil {
  277. return
  278. }
  279. } else {
  280. moveDay = tmpMoveDay
  281. }
  282. // 移动天数到对应农历 的 公历 日期
  283. historyPreDate = historyPreDate.AddDate(0, 0, moveDay)
  284. }
  285. historyPreDateStr := historyPreDate.Format(utils.FormatDate)
  286. if tmpValue, ok := handleDataMap[historyPreDateStr]; ok { //上一年同期找到
  287. pastValueList = append(pastValueList, tmpValue)
  288. }
  289. //if findItem, ok := dataMap[hisoryPreDateStr]; ok { //上一年同期找到
  290. // pastValueList = append(pastValueList, findItem.Value)
  291. //} else if isCompatibility { // 如果需要兼容上下35天
  292. // nextDateDay := hisoryPreDate
  293. // preDateDay := hisoryPreDate
  294. // for i := 0; i < 35; i++ {
  295. // nextDateDayStr := nextDateDay.Format(utils.FormatDate)
  296. // if findItem, ok := dataMap[nextDateDayStr]; ok { //上一年同期->下一个月找到
  297. // pastValueList = append(pastValueList, findItem.Value)
  298. // break
  299. // } else {
  300. // preDateDayStr := preDateDay.Format(utils.FormatDate)
  301. // if findItem, ok := dataMap[preDateDayStr]; ok { //上一年同期->上一个月找到
  302. // pastValueList = append(pastValueList, findItem.Value)
  303. // break
  304. // }
  305. // }
  306. // nextDateDay = nextDateDay.AddDate(0, 0, 1)
  307. // preDateDay = preDateDay.AddDate(0, 0, -1)
  308. // }
  309. //}
  310. }
  311. if len(pastValueList) == formulaInt {
  312. delete(removeDataTimeMap, av) //将待删除的日期给移除
  313. val := CjjxSub(currentItem.Value, pastValueList)
  314. if existVal, ok := existDataMap[edbCode+av]; !ok {
  315. timestamp := currentDate.UnixNano() / 1e6
  316. timestampStr := fmt.Sprintf("%d", timestamp)
  317. addSql += GetAddSql(edbInfoIdStr, edbCode, av, timestampStr, val)
  318. isAdd = true
  319. } else {
  320. existValDecimal, err := decimal.NewFromString(existVal)
  321. existStr := existValDecimal.String()
  322. if existStr != val {
  323. sql := ` UPDATE %s SET value=?,modify_time=NOW() WHERE edb_info_id=? AND data_time=? `
  324. sql = fmt.Sprintf(sql, dataTableName)
  325. err = to.Exec(sql, val, edbInfoId, av).Error
  326. if err != nil {
  327. return err
  328. }
  329. }
  330. }
  331. }
  332. existDataMap[edbCode+av] = av
  333. }
  334. //删除已经不存在的超季节性指标数据(由于该指标当日的数据删除了)
  335. {
  336. removeDateList := make([]string, 0)
  337. for dateTime := range removeDataTimeMap {
  338. removeDateList = append(removeDateList, dateTime)
  339. }
  340. removeNum := len(removeDateList)
  341. if removeNum > 0 {
  342. //如果拼接指标变更了,那么需要删除所有的指标数据
  343. tableName := GetEdbDataTableName(source, subSource)
  344. sql := fmt.Sprintf(` DELETE FROM %s WHERE edb_info_id = ? and data_time in (`+utils.GetOrmInReplace(removeNum)+`) `, tableName)
  345. err = to.Exec(sql, edbInfoId, removeDateList).Error
  346. if err != nil {
  347. err = fmt.Errorf("删除不存在的超季节性指标数据失败,Err:" + err.Error())
  348. return
  349. }
  350. }
  351. }
  352. if isAdd {
  353. addSql = strings.TrimRight(addSql, ",")
  354. err = to.Exec(addSql).Error
  355. if err != nil {
  356. return err
  357. }
  358. }
  359. return
  360. }
  361. // getMoveDay 获取两个日期的平移天数
  362. func getMoveDay(lastDataDay, currentDataDay time.Time) (moveDay int, err error) {
  363. if lastDataDay.Month() >= 11 { //最新数据的日期如果大于等于11月份,那么用的是下一年的春节
  364. lastDataDay = lastDataDay.AddDate(1, 0, 0)
  365. }
  366. currentYear := lastDataDay.Year()
  367. currentYearCjnl := fmt.Sprintf("%d-01-01", currentYear) //当年的春节农历
  368. currentYearCjgl := solarlunar.LunarToSolar(currentYearCjnl, false) //当年的春节公历
  369. currentYearCjglTime, tmpErr := time.ParseInLocation(utils.FormatDate, currentYearCjgl, time.Local)
  370. if tmpErr != nil {
  371. err = errors.New("当前春节公历日期转换失败:" + tmpErr.Error())
  372. return
  373. }
  374. //指定年的春节农历
  375. tmpYearCjnl := fmt.Sprintf("%d-01-01", currentDataDay.Year())
  376. //指定年的春节公历
  377. tmpYearCjgl := solarlunar.LunarToSolar(tmpYearCjnl, false)
  378. tmpYearCjglTime, tmpErr := time.ParseInLocation(utils.FormatDate, tmpYearCjgl, time.Local)
  379. if tmpErr != nil {
  380. err = errors.New(fmt.Sprintf("%d公历日期转换失败:%s", currentDataDay.Year(), tmpErr.Error()))
  381. return
  382. }
  383. // 将两个日期同步到同一年,然后计算两个日期相差的天数
  384. tmpCurrentYearCjglTime := currentYearCjglTime.AddDate(currentDataDay.Year()-currentYear, 0, 0)
  385. moveDay = utils.GetTimeSubDay(tmpYearCjglTime, tmpCurrentYearCjglTime)
  386. return
  387. }
  388. // CjjxSub 计算超季节性值
  389. // 计算公式=现值-过去n年(包括今年)均值,n为取数个数,需大于等于1;
  390. // 举例:A指标 2022-10-13值100,2021-10-13值120,2020-10-13值110,设置n=3,则“超季节性”指标计算值为100-(100+120+110)/3=-10。
  391. func CjjxSub(currValue float64, pastValue []float64) (value string) {
  392. num := len(pastValue)
  393. if num == 0 {
  394. return
  395. }
  396. numDecimal := decimal.NewFromInt(int64(num))
  397. af := decimal.NewFromFloat(currValue)
  398. //fmt.Println(af)
  399. bf := decimal.NewFromFloat(pastValue[0])
  400. for k := 1; k < num; k++ {
  401. tmpVal := decimal.NewFromFloat(pastValue[k])
  402. bf = bf.Add(tmpVal)
  403. }
  404. val, _ := af.Sub(bf.Div(numDecimal)).Float64()
  405. //valStr := utils.SubFloatToString(val, 4)
  406. valStr := decimal.NewFromFloat(val).Round(4).String()
  407. return valStr
  408. }