作者a34021501 (CARD)
看板CodeJob
标题[情报] 单维度空间频率分析程式码!
时间Sat Mar 11 18:25:09 2017
大家好
我今天又要来分享一下我写过的程式码
这个程式码我只用於分析我自己的心跳
我作完之後觉得没有很想要再继续分析
http://cargon.net/HHSA/
这个程式码除了用Empirical Mode Decomposition(EMD)分析出多组Mode
还继续分析每组Intrinsic Mode Function(IMF)中隐藏的Amplitude资讯
藉此可以看出各种频谱中的组成份的震幅讯息等资讯在时间轴上的现象!
今天是个特别的日子
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.171.26.170
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/CodeJob/M.1489227913.A.A87.html
※ 编辑: a34021501 (1.171.26.170), 03/11/2017 18:26:53
1F:→ a34021501: 我的ECG是量测三个点,是含有空间资讯叠加的时序讯号! 03/11 19:17
2F:→ laject: 你怎麽还在,没被抓走吗 03/11 19:47
除非我用这种技术搭配多维空间频率分析,透过多点能量波相交於某个点的叠加高能震波
这种分析方法可用於当代所有领域的研究,就可以知道多维度科技早就已经应用於生活中
基本上我只有写程式分析我个人的心跳讯号,没有法律问题!那时候我的心跳没有很平顺
※ 编辑: a34021501 (1.171.26.170), 03/11/2017 20:21:24
3F:推 shadowjohn: 脚踏车打到快55ps的感觉,活塞就出来了 03/11 22:35
4F:→ shadowjohn: 抱歉,要回讯息按到推文了 03/11 22:36
5F:推 qxxrbull: 这个跟GMEMD的有什麽不同呢 03/11 22:48
主要用於分析时序资料,HSA(Time,IF,AM) && HHS(IF,HoloIF,HoloAM)
这是从不同面相看资料,时间、瞬时频率、振福 瞬时频率、波峰瞬时频率、波峰瞬时振幅
这样就可以看到被感应器侦测到的讯号中有不同空间频率的能量叠加在整个 RAW 资料之中
Data -> [CEEMD] -> IMF
IMF -> [AM-FM] -> AM,FM
FM -> [NDQ] -> IF
IMF波峰 -> [CEEMD] -> Holo_IMF
Holo_IMF-> [AM-FM] -> Holo_AM,Holo_FM
Holo_FM -> [NDQ] -> Holo_IF
其中[CEEMD]为[EEMD]的双Noise版本,加入正向与反向会抵销的Noise
其中[EEMD] 系以加Noise补足个频带避免mode moxing, 计算Ensemble次的[EMD]後作平均
其中[EMD] 系以反覆扣除平均包络线产生出modg数量组的高频到低频的IMF讯号
其中[AM-FM]是将数据调变为振福AM与频率FM
其中[NDQ] 使用数值分析技巧将频率FM转换成瞬时频率IF
其中[HHSA] 为使用[EMD]演算法分析IMF之振幅波峰的方法
总之[HSA] 是将原始Data分析後投影到Time-IF-AM
而 [HHS] 是将nIMF组分析之後投影到IF-HoloIF-HoloAM
总而言之.. [HHS] 之投影时序我没有写完整!需要後辈努力!
反正就是一种透过Amplitude来载入空间波的方式可以被分析!
我的[GMEMD]可以分析出许多组不同空间频率的空间分布影像!
若有错误敬请指证,谢谢大家!
特别感谢:
中央研究院 黄锷 Norden E. Huang 院士与其团队 所发表的多篇着作
让我能在短时间内理解空间频率与时序分析之间的关联性! Thanks!!
※ 编辑: a34021501 (1.171.26.170), 03/12/2017 00:08:56
6F:→ a34021501: 还有The second layer of HHSA可以分析更精微的资讯! 03/12 16:38
7F:→ a34021501: 在HHSA的论文里是分析length-of-day结果包含天文周期! 03/12 16:40
8F:→ a34021501: 已修正HHSA程式输出模式,增加HSA与HHS的资料与图档! 03/14 00:21
9F:→ a34021501: 程式码重新编排、修正HHS底色异常、新增HHSA时序图表! 03/14 10:46
目前最新Build: 20170314_01
http://cargon.net/HHSA/HHSA_CPU_code_20170314_01.7z
因为我新增了 [HHS2] 以输出 HHS 在各个 HSA 时间区块的细节,即 HHSA_TimeCell####
即 HHS 是张 Freq_Cell * Holo_Freq_cell 在每个 Time_cell 区间的 Holo_AM 累积量
(IF) (Holo_IF)
而 HSA 是张 Time_Cell * Freq_Cell 的 AM 累绩量
(Time) (IF)
但我觉得若要看懂 HHSA_TimeCell#### 的资讯,要拿 HSA 来参考比对个时间点的 AM
先简单解释一下的 HSA 资讯包含的能量强度,是 Data 经过拆频後个频带 IMF 的 AM
先来讲解一下何谓 AM , 即 Signal 经过 [FM-AM] 调变後能代表 Signal 之振幅强度
那麽另一个变量即 FM , 即 Signal 在时间上不考虑其振幅的波形,含有频率的资讯
理论上 Signal(t) = FM(t) * AM(t) 即在特定时间(t)之在某空间参考点测量到的能量
而 IF 是经由 FM 的 [瞬时频率-顺时相位] 分析後转换成 IF 与 IP 即与角频率有关系
理论上 HHS 是某段时间 Duration 测量之 AM 的 Fluctuation 即 IMF 上包络线的资讯
而 HHS 资讯包含的能量强度 Holo_AM 是各频带 Fluctuation 在各个频谱中的分析结果
总而言之 Dectector 所测得 Signal 之 IMF 频谱中的 AM 所包含相乘频谱即 Holo_AM
某个空间参考点测量到的能量在时间上的变化可能与感测器之设计与环境之空间频率有关
因为如果环境电浆有疏密的现象流过感测器则可能会影响 Holo_FM 所分析出之 Holo_IF
当然还有可能包含其他因素例如地球自转时白天与晚上接收太阳能量对脑电波的频谱组成
(这是HHSA论文内容)
不过脑电波的测量可能有其他讯号参杂在回路上如其它能量源的功率因素、因子、因数等
※ 编辑: a34021501 (36.229.39.135), 03/14/2017 18:46:24
或许有人担忧环境之背景电磁波影响情绪的情况只要穿上抗UV的衣物基本上就能降低干扰
------------------------------------------------------------------------------
来提一个我观察到的正面研究提供给大家参考,毕竟太阳能量波似乎每19年一次就是今年
(这是HHSA论文内容,使用second layer HHSA)
在分析地球、月球与太阳的互动下所影响的每天时间差(LOD)会在 AM Frequency表现出来
似乎太阳与更遥远的星体间互动的现象会造成更大的能量波之周期,有听闻光子带的现象
因为尚未确认太阳是否和天狼星一样是个双星系统,但至少会和银河系中央的能量源互动
在天文上着名的周期有
Olympiad cycle (4-year) 2^2
https://en.wikipedia.org/wiki/Olympiad
Octaeteris cycle (8-year cycle of antiquity) 2^3
https://en.wikipedia.org/wiki/Octaeteris
Saros cycle of eclipses (18-year) 2 * 3^2
https://en.wikipedia.org/wiki/Saros_(astronomy)
Metonic cycle (19-year) 19
https://en.wikipedia.org/wiki/Metonic_cycle
Callippic cycle (76-year cycle from 330 bc) 2^2 * 19
https://en.wikipedia.org/wiki/Callippic_cycle
Hipparchic cycle (304-year cycle from 2nd century bc) 2^4 * 19 (2-cycle)
https://en.wikipedia.org/wiki/Hipparchic_cycle
GCD(4,8,18,19,76,152)=1
LCM(4,8,18,19,76,152)=1368
1368 = 2^3 * 3^2 * 19
或许1368年是一个更大的周期,毕竟 Period 或 Duration 的共同倍数会让 AM 能量提升
毕竟1368的因数中只有19不是二次方以上,若将其补足到二次方以上就是 1368*19=25992
很接近许多文章中26000年的循环周期,且还有人画出其椭圆形轨道似乎是类似双星系统
因为一个带有能量的物体在空间中移动,如物质、电浆、粒子、光子等,都会造成空间波
如此一来就可预测更大的潜在周期,我们先分析近似26000的 25992 = 2^3 * 3^2 * 19^2
可以发现这些周期皆是由质数组合而成,看来太阳可能是在19这个质数的轨道因子上运行
感觉看起来就像填周期轨道时,会先填在较小的质数,不知道这样的现象是否发生在宇宙
1368 = 2^3 * 3^2 * 19
12996 = 2^2 * 3^2 * 19^2
25992 = 2^3 * 3^2 * 19^2
77976 = 2^3 * 3^3 * 19^2
以我对天文的看法,这有点像椭圆轨道面积 ab * M_PI,而在空间中则是体积 abc * M_PI
而我们的 Sensor 通常的是一个平面的小面积积分,所以应该只要以 HHSA 就能分析通量
※ 编辑: a34021501 (36.229.39.135), 03/15/2017 06:23:23
※ 编辑: a34021501 (36.229.39.135), 03/15/2017 06:29:13
来提一下我以前做过的相关研究,让大家更容易理解环境现象
这样的星际现象在大气上的表徵能看出较小的 Vortex 在其中
http://cargon.net/GMEMD/20140921_Fung-Wong_Typhoon.gif
这样2D影像频谱经过数值分析後应该能转换成我们熟知的光谱
http://cargon.net/GMEMD/lena_GMEMD_Spectrum.xlsx
http://cargon.net/GMEMD/20141204_GMEMD_Spectrum.pptx
我觉得这些研究都要考虑 Dynamic Light Scattering 的现象
http://cargon.net/GMEMD/20100620_Dynamic_Light_Scattering.pdf
毕竟电浆层是一种粒子在空间之 Spin 与移动与 Vortex 现象
※ 编辑: a34021501 (36.229.39.135), 03/15/2017 16:09:55
今天有人问我是否拥有相关专利,我其实还没申请过任何专利
简单说明,不同光谱在空间中转弯幅度不同,但空间仍有变数
如重力透镜不影响大体积之量重,可是影响空间,涡旋方向未知
如果用简单物理解释,UFO内部就是要量重不变身体才不会崩溃
应该可以先模拟出电磁驱动之重力波是否适合让载具飞上天空
但是能量源是个庞大的议题,毕竟在地球表面磁波充能效率差
毕竟地球的磁场有一些歪斜,可能没有和太阳磁场准确地校准
幸好有月球可以稳定地球周围磁场,让太阳能量较容易被吸收
地磁据说已经越来越弱了,或许可以用HHSA分析地磁等待反转
其实磁场可以有光的现象,例如电磁铁拥有特殊频率与强度等
※ 编辑: a34021501 (36.229.39.135), 03/16/2017 06:05:23
这种在生理上的现象就像是人体的血压会有些微小的波动一样
http://imgur.com/a/HTcTO Oo收缩压-舒张压|心律oO
类似的现象也应用於无线网路的QAM, 我觉得Phase也值得研究
※ 编辑: a34021501 (36.229.39.135), 03/16/2017 14:58:56