关于EMD分解的边界效应处理

sssbbbhhh1

1  左端处理
avg=(sum(abs(extrx1))+sum(abs(rightvalue)))/(length(extrx1)*2-2);
       error = rightvalue - extrx1;
中rightvalue应为leftvalue
2.error 不能作为变量名

peterjiasheng

请问那个积分中值求均值的程序谁有，能发我一份吗？我的邮箱是peterjiasheng@163.com 不胜感激

peterjiasheng

[mvalue,indt] = compose(tmin,tmax,mmin,mmax);
      [env,tindex] = mvtd(m,mvalue,indt,0);    % 用中值定理法求均值
请问这个引用compose 函数是不是有问题？我仿真怎么说是:Function 'compose' is not defined for values of class 'double'. 请问这是什么意思？谁能指导一下不？

peterjiasheng

tmax=[0,65,209,358,482,629,775,899,1047,1196,1356,1513,1645,1822,2004,2129,2291,2428,2676,2822,2961,3098,3226,3351,3484,3617,3738,3882,4026,4163,4304,4424,4545,4679,4816,4944,5086,5229,5364,5622,5745,5880,6019,6149,6283,6424,6549,6687,6821,6948,7084,7219,7350,7482,7622,7747,7883,8021,8146,8284,8419,8551,8684,8822,8948,9082,9221,9342,9482,9619,9752,9886,10023,10151,10283,10421,10537,10679,10819,10951,11087,11223,11354,11484,11621,11736,11877,12019,12149,12286,12422,12555,12686,12821,12941,13077,13219,13345,13484,13622,13756,13888,14021,14145,14279,14419,14542,14681,14821,14954,15089,15221,15350,15481,15620,15741,15879,16021,16151,16287,16421,16552,16684,16820,16943,17079,17220,17347,17485,17621,17753,17886,18021,18146,18280,18420,18543,18682,18821,18952,19087,19221,19350,19482,19620,19742,19880,20018];
mmax=[2.782579600351173,2.782579600351173,-0.426844403850422,2.886786415954008,3.484763389713516,0.446486693886786,4.325043100573084,3.766461117786259,0.482757103585305,3.297601637178036,1.005039563227511,1.430418654422075,2.833739172548888,-1.648137294740896,4.179375404232564,3.218405069583238,1.893535689075064,3.22941066615504,-1.70996869698034,2.754155766670606,-1.537104109666618,-0.925471405232276,4.205683624502091,-0.369465811238296,-0.485566052659371,4.170982309646925,-0.433026996341911,-1.332691788030235,2.720009488267894,-2.179629474726049,-1.687719662875929,3.571684736135952,-0.503083733408314,-0.440241963953304,4.554997126622922,0.214049782786196,-0.581384015193431,3.025075976604779,-2.387647262464196,2.935204195789079,-0.773490561037781,-0.553252553837693,4.601528585634375,0.141110037121533,-0.370814280026209,3.642372283174192,-1.776733207932438,-2.217426835510402,2.640132004927345,-1.531594690365861,-1.233021343549119,4.206822651326115,0.121815183670281,-0.069931389302601,4.243617024531574,-1.086756802959381,-1.781331442648525,2.656983999956489,-1.951508340895996,-1.858588530663026,3.603633484225298,-0.222349078421149,-0.075214927213715,4.607002875444909,-0.494171724556054,-1.179816628285236,3.025844147594565,-1.925925103982449,-2.208068547508429,3.018952981930304,-0.819397795485348,-0.454999608548294,4.605183335811754,-0.137792774484172,-0.595926461923857,3.616542546675948,-1.526749047884966,-2.179362334905427,2.664887583085978,-1.414156021684091,-1.056734999702269,4.236841415442747,-0.137206021263745,-0.250062262212344,4.204060992474476,-0.90193260737402,-1.769820884561938,2.680878749754056,-1.801795181575203,-1.669673570387267,3.643069540494533,-0.488720261404116,-0.262896020995269,4.571962372537198,-0.333180610566721,-1.153353008466652,3.050968942753999,-1.824097251149141,-2.067757796257956,3.054060356141056,-1.040316304519916,-0.606410104660073,4.572992989053503,-0.031982864904952,-0.590849399754456,3.639837515875682,-1.441925694384925,-2.058834429532975,2.688947217934875,-1.614409663962974,-1.195010809803914,4.209062803475073,-0.047813649259152,-0.236987915334907,4.226262800894274,-0.868215121052265,-1.693302001774528,2.690652200692535,-1.956325739410066,-1.779377999049231,3.620719555548029,-0.416470988578616,-0.249632319713015,4.58449882038032,-0.306340746929481,-1.104644271329297,3.056026197571724,-1.931199852600907,-2.135261058210451,3.034744040131333,-1.002367250321309,-0.617823643731552,4.57802490229877,0.022038859087893,-0.531468141089857,3.645206022979021,-1.539930067245036,-2.118144296697719,2.675488010040728,-1.583943848946179,-1.20416326728913,4.211880931736641,-0.00580178149512,-0.192814472492541,4.22994263622407,-0.937163449732998,-1.734652982807107,2.460112931886871];

tmin=[1,140,275,417,567,688,838,983,1119,1272,1433,1578,1767,1859,2066,2213,2362,2615,2736,2906,3023,3131,3320,3420,3521,3708,3820,3942,4117,4225,4321,4511,4614,4720,4902,5022,5149,5328,5429,5704,5814,5927,6105,6221,6339,6522,6623,6727,6907,7015,7123,7303,7419,7535,7719,7823,7932,8112,8217,8322,8501,8618,8733,8917,9023,9133,9316,9418,9524,9701,9817,9931,10113,10222,10333,10520,10619,10727,10903,11017,11129,11310,11421,11531,11720,11819,11930,12106,12218,12329,12508,12619,12728,12916,13020,13132,13311,13419,13529,13706,13818,13926,14112,14219,14331,14514,14620,14731,14907,15018,15125,15308,15419,15530,15715,15820,15933,16110,16219,16325,16506,16618,16728,16913,17020,17133,17313,17419,17527,17705,17817,17926,18110,18220,18333,18516,18620,18730,18907,19017,19124,19307,19419,19532,19717,19821,19932,20117];
mmin=[-2.296846525428775,-4.141623383125126,-4.084400242074468,0.165722184124079,-2.581249662433192,-2.116472992139788,1.136570601261942,-3.045081330767765,-2.972607334128952,-1.925543895436031,-4.871520644665601,-0.521316458710067,-2.192876415860649,-1.861963514986633,1.50560196796016,-2.914292656518384,0.519140276285282,-3.883123848963708,-2.806325396027761,-2.237631600793283,-3.221226039291187,-1.071224806274728,-0.463477168999168,-1.660299945403862,-0.664397103160116,-0.526777625109889,-3.360709334153492,-2.280793256440808,-2.599278536768881,-3.583498810153754,-1.704885104834053,-0.629334466285749,-1.968077926144889,-0.755271962500725,-0.039147942604617,-2.810828312810431,-1.515612372912226,-2.587602470818654,-3.686622880335742,-0.992898634030865,-2.650747503176382,-1.059857342489654,-0.17242152090046,-2.133342172288208,-1.122432327791427,-1.849779184199289,-3.567453121163892,-2.467638879805008,-1.773115313930688,-3.191370307708022,-1.532393381456038,-0.239410434908048,-1.93102391459998,-0.725532981877858,-1.169369353871392,-3.110854081074671,-2.23262198890181,-2.097533990735767,-3.644475842453983,-2.120377601127208,-0.66299325932594,-2.143704482476587,-0.614079032953674,-0.614031647854722,-2.497890976565047,-1.701615968272349,-1.986760566989064,-3.817983650228092,-2.553637686498472,-1.28058664316398,-2.64853996479398,-0.873588261974826,-0.343924125988255,-2.000157010512718,-1.0932630027658,-1.544083183954068,-3.616255743283878,-2.655737881200884,-1.808347179470624,-3.26063830051846,-1.416267211320117,-0.448382168487173,-1.812790735921519,-0.656267804581154,-0.915896862890417,-3.123076577064205,-2.38988688574425,-2.078428391319418,-3.735177928240775,-2.02869869489549,-0.881850028828932,-2.016003576784173,-0.562557670202894,-0.388343471438504,-2.516540075926342,-1.835828563506895,-1.970311318325455,-3.890841372179616,-2.493810947919833,-1.474870604670952,-2.540261707017983,-0.838820338896716,-0.166740104751433,-2.012138332084731,-1.204415019661657,-1.516998345942193,-3.683751537147896,-2.613449487021327,-1.981328346299294,-3.164114276511798,-1.399836502872536,-0.303680664333675,-1.837874599006318,-0.746069480804387,-0.926225414265623,-3.168330811298507,-2.353297308007195,-2.206431940274037,-3.661880101337778,-2.020687925490941,-0.771780254833345,-2.044695323077893,-0.629092413112705,-0.397906901877186,-2.545752266250568,-1.813895590765952,-2.069884659645433,-3.844333038767063,-2.470203233444309,-1.387488094867794,-2.555824767738202,-0.897152642123674,-0.153787065247497,-2.047205964498041,-1.196905816530512,-1.611901670682193,-3.646196591169452,-2.586880495506465,-1.914720647936704,-3.175446974512419,-1.442574391930074,-0.295889165773106,-1.862858113645833,-0.738553757370997,-0.993466440201802,-3.143090804173905,-2.330644234655106,-2.700512798743145];
我产生的数据是这样的，用[mvalue,indt] = compose(tmin,tmax,mmin,mmax);好像不可以？通过不了，谁能帮我看看，是不是这个函数用错啦，请问谁知用哪个函数不？:loveliness:

nkdtxf

请问哪位知道神经网络，或者支持向量机的延拓方法程序呢。谢谢了

aishuishou

回复 9 # cleverblue 的帖子

为什么我写的延拓程序放到emd里面，只迭代一次就停止了，另外请问你怎么看GRILLING镜像法的：
function [tmin,tmax,zmin,zmax] = boundary_conditions(indmin,indmax,t,x,z,nbsym)
lx = length(x);
if (length(indmin) + length(indmax) < 3) //这个你怎么看？在emd理论里没有啊？！
  error('not enough extrema')
end

lanceinmf

先学习学习了

younger1985

回复 20 # nkdtxf 的帖子

其实你用神经网络或者支持向量机做的效果不见的好，你试试就知道了 ！哎，感觉好多理论科学都是实际都脱节！

younger1985

回复 21 # aishuishou 的帖子

有吧 你参考一下浙江大学胡劲松博士论文里面就有这个

aishuishou

回复 24 # younger1985 的帖子

哈哈，谢谢提醒啊！！！我当初看这篇论文感觉有点难，没细看，刚才有回去看了一遍确实，他在博士毕业论文里2.3.2节说了，至少两个极值点（和主题里论文的要求极值点大于2没有本质区别，因为因为emd工具箱处理是按至少三个极值点处理的，结果就是：两个极值点时最后的趋势想为y=0，而且程序未予显示，并且三个极值点的话就会产生一个单调趋势项），终于明白了！！！

nkdtxf

谢谢，楼上，提醒我了，明白了，趋势项的理论依据了

hitdely

回复 20 # nkdtxf 的帖子

请问神经网络，或者支持向量机的延拓方法程序  你那有吗 可以分享一下吗 谢谢

yuwenamy

去掉emd.m中的边界延拓部分的程序,怎样去啊 ！

yueemd

恩啊，我也有跟ls一样的问题啊，不知道该怎么去掉emd中的延拓部分，希望大家帮忙解释一下。。。谢谢了

peterzitong

回复 15 # cleverblue 的帖子

不知道楼主还在此论坛里面吗