>> EMD2.java - Hilbert-Huang (HHT) 轉換方法是黃
点击查看更多 ▼
EMD2.java - Hilbert-Huang (HHT) 轉換方法是黃
import java.io.F
import java.io.FileInputS
import java.io.FileNotFoundE
import java.io.InputS
import java.util.ArrayL
import java.util.S
public class EMD2 {
// 利用序列预测法延拓序列，目标是要在断点处获得接近真实的包络并且算法耗时尽量短。
//序列预测法有很多种：神经网络、时间序列模型、多项式拟合等等。根据延拓对象又可以分为
//数据波形延拓、极值延拓等。个人目前认为直接延拓波形可能计算量大导致耗时较长
//类EMD采用镜像延拓，获得了一定的效果，但处理还是过于简单，端点效应仍然明显。
// 采用镜像延拓处理端点，但极值点没有改进，新建文件“jiandan1中将用多项式拟合来处里极值问题”,HOWEVER,jiandan1伤了老夫，多项
// 式拟合貌似有问题啊。。。。。。。。详情间jiandan1
//function imf = emd(x,n);%%
public static int itmax=14;
//public static double data[]=new double[datalength];
public static Minmax extr(double[] a)
//变量定义
int length=a.
int[] indmax0=new int[a.length];//用来存放极大值的indice（坐标）
int[] indmin0=new int[a.length];//用来存放极小值的indice（坐标）
int lenmax=0;//极大值个数
int lenmin=0;//极小值个数
//判断是否是极大或极小值
///第一步，求微分
double[] d=new double[a.length-1];
d=diff(a);
///如果d[n-1]*d[n]&0,并且d[n-1]&0,则a[n]为极大值，反之，若d[n-1]&0,a[n]则为极小
for(int i=1;i&length-1;i++)
if(d[i-1]*d[i]&0)
if(d[i-1]&0)
indmax0[lenmax]=i;
lenmax=lenmax+1;
if(d[i-1]&0)
indmin0[lenmin]=i;
lenmin=lenmin+1;
///如果出现一个或连续几个d[i]=0，首先看这些等值点是否是在序列的首尾出现，如果是在首尾
///则不能判断是否是极值点，否则应算作一个极值点，坐标取这一段的平均值；
double[] d1=new double[a.length-1];
for(int i=0;i&length-1;i++)
if(d[i]==0)d1[i]=1;
else d1[i]=0;
//找出等值序列和非等值序列的转折点,将等值序列的起始点分别存在begins和ends中
int[] begins=new int[a.length-1];
int[] ends=new int[a.length-1];
int lbegin=0;//起始点的个数
int lend=0;//终止点的个数
for(int i=1;i&length-1;i++)
if((d1[i]-d1[i-1])==1)
begins[lbegin]=i;
lbegin=lbegin+1;
if((d1[i]-d1[i-1])==-1)
ends[lend]=i;
lend=lend+1;
//首尾的等值序列不能判断是否是极值点，剔除
if(begins[0]&ends[0])//首端出现等值序列
for(int i=0;i&lend-1;i++)
ends[i]=ends[i+1];
lend=lend-1;//剔除过后，ends个数减一
if(lbegin&lend)lbegin=lbegin-1;//如果末端出现等值序列则会导致begins的个数比ends个数多1，也要剔除；
if((lbegin==lend)&&lbegin!=0)
for(int i=0;i&i++)
if(d[begins[i]-1]*d[ends[i]]&0)
if(d[begins[i]-1]&0)
indmax0[lenmax]=(begins[i]+ends[i])/2;
lenmax=lenmax+1;
indmin0[lenmin]=(begins[i]+ends[i])/2;
lenmin=lenmin+1;
//以下分别对得到的极大值和极小值的坐标进行排序
for(int j=0;j&j++)
if((indmax0[lenmax-lbegin+j]&indmax0[0]))
int temp=indmax0[lenmax-lbegin+j];
for(int k=lenmax-lbegin+j;k&0;k--)indmax0[k]=indmax0[k-1];
indmax0[0]=
else if((indmax0[lenmax-lbegin+j]&indmax0[lenmax-lbegin+j-1]));
for(int i=0;i&lenmax-lbegin+j-1;i++)
if((indmax0[lenmax-lbegin+j]&indmax0[i])&&(indmax0[lenmax-lbegin+j]&indmax0[i+1]))
int temp=indmax0[lenmax-lbegin+j];
for(int k=lenmax-lbegin+j;k&i+1;k--)indmax0[k]=indmax0[k-1];
indmax0[i+1]=
for(int j=0;j&j++)
if((indmin0[lenmin-lbegin+j]&indmin0[0]))
int temp=indmin0[lenmin-lbegin+j];
for(int k=lenmin-lbegin+j;k&0;k--)indmin0[k]=indmin0[k-1];
indmin0[0]=
else if((indmin0[lenmin-lbegin+j]&indmin0[lenmin-lbegin+j-1]));
for(int i=0;i&lenmin-lbegin+j-1;i++)
if((indmin0[lenmin-lbegin+j]&indmin0[i])&&(indmin0[lenmin-lbegin+j]&indmin0[i+1]))
int temp=indmin0[lenmin-lbegin+j];
for(int k=lenmin-lbegin+j;k&i+1;k--)indmin0[k]=indmin0[k-1];
indmin0[i+1]=
//将结果存到minmax中
int[] indmax=new int[lenmax];//用来存放极大值的indice（坐标）
int[] indmin=new int[lenmin];//用来存放极小值的indice（坐标）
for(int i=0;i&i++)indmax[i]=indmax0[i];
for(int i=0;i&i++)indmin[i]=indmin0[i];
Minmax minmax=new Minmax();//用来存放极大值、极小值的个数和坐标
minmax.setIndmax(indmax);
minmax.setIndmin(indmin);
minmax.setLenmin(lenmin);
minmax.setLenmax(lenmax);
return(minmax);
//端点处理，利用镜像法，输入是利用extr抽取出的极值坐标和原始信号，输出时经过端点延伸后
//的极值坐标和值
public static Minmax1 boundary_extend(int[] indmax,int[] indmin,double[] x,int num)
//num为在每端延拓的点数
int length=x.
int lenmax=indmax.
int lenmin=indmin.
//左镜像点坐标
展开＞ ＜收缩
以上是部分代码预览下载源码到电脑，阅读使用更方便
Sponsored links
源码文件列表
温馨提示： 点击源码文件名可预览文件内容哦 ^_^
11.17 kB 11:23
6.03 kB 20:08
814.00 B 14:41
17.06 kB 11:00
25.94 kB 20:36
40.95 kB 11:28
2.41 kB 22:09
10.06 kB 14:44
378.00 B 21:58
615.00 B 17:09
711.00 B 11:28
10.84 kB 15:12
4.06 kB 16:18
4.06 kB 16:18
256.00 B 19:23
Sponsored links
23 篇源代码 21 篇源代码 18 篇源代码 13 篇源代码 9 篇源代码
285 篇源代码 173 篇源代码 48 篇源代码 42 篇源代码 36 篇源代码
评价成功，多谢！
您的CF币不足
支付宝优惠套餐快速获取 10 CF币
支付宝订单完成后，CF币将自动加入到您的账号。以下是优惠期的人民币价格，优惠期过后将恢复美元价格。
您本次下载所消耗的CF币将转交上传作者。
同一源码，30天内重复下载，只扣除一次CF币。
登录 CodeForge
还没有CodeForge账号？
Switch to the English version?
^_^"呃 ...
Sorry!这位大神很神秘，未开通博客呢，请浏览一下其他的吧
请按 Ctrl+D 键添加到收藏夹。关于孩子的文章- 美文故事-散文随笔
关于/有关孩子的文章
点击:14610
　　在电影院排队买票，我前面是一对年轻恋人，刚排到他们，一位妈妈领着孩子急匆匆挤过来，直接冲售票小姐说：我们的已经开场了，先给我们出票吧。 我前面的姑娘不乐意了，说您排一下队好吗？ 那位妈妈完全不理，直接递钱给售票小姐：孩子急着看，麻烦你先给我们出吧。 ()
点击:14400
　　无论光阴之外，或光阴之里，我们都是时间的孩子。它牵着我们，深一脚，浅一脚地往前走，带给我们欢欣喜悦，也带给我们逼仄疼痛。 1 打开衣柜，着实吃了一惊。不知什么时候，衣服基本都变成深色系了，以一片肃穆寡欲为主，即便带点颜色，也是偏向暗沉的浊色。 以前的我()
点击:19118
　　作为一名教育和保健护理工作者,我曾经和数不清的感染上艾滋病病毒的孩子打过交道.我和这些特殊的孩子之间的关系是生活赋予的恩赐.他们教会我许多东西,我尤其懂得了即使在最弱小的人物身上也能发现其所蕴含的巨大勇气.让我告诉你泰勒的故事()
点击:16321
　　请放下手机，就一会儿……。你年迈的父亲盼星星盼月亮，终于盼到你放假回家，但好几天了，你总埋头于手机，对于父亲的问话，你只是简单地“嗯嗯”，再没有多出的话语。你很忙吗？你在忙什么？以前那个亲热的孩子哪去了？年迈的父亲可怜巴巴地望着你，苦苦地等待着……()
点击:20155
　　一场考试，触及灵魂――题记 期中考试后第三天， 儿9岁，三年级，他告之我说，“明天，要开家长会，老师说，必须要带一支笔。” 其实不是家长会，没有往常的通告、总结、家长代表讲话。确切的说，这是一场考试，面对家长的考试，或者是一次调查。 题目简单，直接，都是()
　　网上看到一篇文章，确切地说是文章题目，说好多父母让自己的孩子参观别墅。我想当然的认为，这样做的目的显而易见，就是希望孩子现在要好好努力，刻苦学习，将来过上富足的生活。这些家长们真可谓‘用心良苦’。 想让孩子过上丰衣足食，无忧无虑的好生活，这是所有家长()
点击:18372
　　每天学生上学的时间，学校门口总会聚集着许多家长。这当中，想把孩子亲自送到教室的低年级家长络绎不绝。尽管值班教师告诉他们：学生自己会走进教室的。但还是有许多家长总是找各种各样的借口，不是怕孩子找不到“作业本”，就是说怕孩子把该拿的东西弄掉了，或者说帮()
点击:24141
　　文/秋色 在央企干了一辈子的我，退休后去了一家审计公司，还是干老本行，私企的工作节奏与央企相比有着天壤之别，紧张有序工作经常忙得团团转， 但感觉还是很充实，心情特别舒畅，与眼前这些九零后的孩子们同台济事，心中有说不出的滋味，感觉时光如梭转眼间就成了这些()
点击:19457
　　再见，迷茫 你走在路上，迷茫的你时常会问自己这样一个问题：为什么一棵树过得比一个人好？一棵树从一支树苗长成一棵参天大树，恰如一个人从一个幼稚的孩子长成一个成熟的大人。是树经历的磨难比人少吗？不是。一棵树在它整个人生阶段要每时每刻都面临风吹雨打、日晒雨()
点击:36528
　　如果你知道 ――致孩子 如果你知道 在米色的阳光后面也有一片看不见的黑暗 你是否还会一如既往地憧憬那种温暖 如果你知道 在温顺的外表下猫咪也会有锋利的爪子 你是否还会毫不犹豫地抚摸它 如果你知道 在美丽的森林里也会有一些不为人知的野兽 你是否还会坚持走近它 如()
共53页/530篇短文
本站会员关于孩子的日志随笔：
Copyright &
版权所有.,,故事,优秀作文在线阅读欣赏当前位置： > 领导留言
日晚，我弟弟被人枪杀在海南省五指山市国税局附近（位于五指山市河南派出所不足五十米远处）。案件定性为聚众斗殴，咨询警方斗殴原因避而不答说是案件正在侦查中不可透露，我们家属表示有些不解但也只能任由，后尸检报告得出结论“气枪”所致。我不知所谓的气枪是否有那么大的威力也不知哪位神枪手手法如此精准能够正中太阳穴穿脑而过。案发过程（目击者口诉）我弟弟同朋友坐着三轮车前面跑后面的人坐着皮卡车开着枪扫射追沿路追。警方说双方都持有危险性武，可为什么没有受伤的人员只有我弟弟一人死亡，不是双方斗殴吗？这条路也是五指山市的主线路都有监控，可是公安部门却告知说监控已受损。监控受损、治安不力、查枪不严难道说公安局如同摆设吗？后家属不信多次上公安局咨询就给看一部分视频可是不足以证明什么，那完整的视频监控又是怎么样的为什么不敢公开。事发至今家属对案情一无所知，去公安局咨询案情，一会说凶手已经归案一会说还在逃，人民警察都是如此纠结的吗？看着老爸老妈一天天的憔悴下去，心理满是难过却无能为力，我不知道弟弟到底做错了什么事还是得罪了什么人只是希望公安能尽快查清案情也能够公平的处理，可是有关部门除了来催促尽快火化并无其它了并且要求自己处理还下文件要求家属签字。家属疑惑能够目无法纪公然持枪当街杀人幕后是否存在保护伞不然为何至今案件毫无进展一拖再拖的忽悠受害者家属。
日晚，我弟弟被人枪杀在海南省五指山市国税局附近（位于五指山市河南派出所不足五十米远处）。案件定性为聚众斗殴，咨询警方斗殴原因避而不答说是案件正在侦查中不可透露，我们家属表示有些不解但也只能任由，后尸检报告得出结论“气枪”所致。我不知所谓的气枪是否有那么大的威力也不知哪位神枪手手法如此精准能够正中太阳穴穿脑而过。案发过程（目击者口诉）我弟弟同朋友坐着三轮车前面跑后面的人坐着皮卡车开着枪扫射追沿路追。警方说双方都持有危险性武，可为什么没有受伤的人员只有我弟弟一人死亡，不是双方斗殴吗？这条路也是五指山市的主线路都有监控，可是公安部门却告知说监控已受损。监控受损、治安不力、查枪不严难道说公安局如同摆设吗？后家属不信多次上公安局咨询就给看一部分视频可是不足以证明什么，那完整的视频监控又是怎么样的为什么不敢公开。事发至今家属对案情一无所知，去公安局咨询案情，一会说凶手已经归案一会说还在逃，人民警察都是如此纠结的吗？看着老爸老妈一天天的憔悴下去，心理满是难过却无能为力，我不知道弟弟到底做错了什么事还是得罪了什么人只是希望公安能尽快查清案情也能够公平的处理，可是有关部门除了来催促尽快火化并无其它了并且要求自己处理还下文件要求家属签字。家属疑惑能够目无法纪公然持枪当街杀人幕后是否存在保护伞不然为何至今案件毫无进展一拖再拖的忽悠受害者家属。
用户名：(空为匿名)
?在发布信息时，请遵守中华人民共和国有关法律、法规，并尊重网上道德；
?管理人员有权对留言内容进行删改；
?因您的言论直接或间接导致的民事或刑事法律责任由您个人承担。
1、报料以提供线索为主。
2、报料内容为突发事件、不法现象、不文
明行为、奇闻趣事等。
3、报料内容要求真实、客观。
4、请不要转载已刊发新闻。
报料热线：966123
电子邮箱：
传真电话：5
官方QQ群：
本网为拓展本土新闻资源，广开新闻渠道，对提供新闻线索并经本网采写成新闻的报料人进行奖励。最高奖励现金5000元。[]
消费者投诉举报咨询　　　　12315
省旅游投诉中心　　　　　　12301
商品价格投诉　　　　　　　12358
文化市场投诉电话　　　　　12318
质量投诉举报热线　　　　　12365
法律维权咨询　　　　　　　12348
海南日报报业集团南海网 版权所有　
电话：(86)5　 传真：5
地址：海南省海口市金盘路30号新闻大厦9楼
本网法律顾问：海南方圆律师事务所　李君律师}