emd分解话出的包络线好奇怪

yanling66520

我是这几天才开始学emd，在网上下载了plot_hht的程序，运行了下自己的信号，得出的包络线很奇怪，在这里贴出来给大家看看，希望知道的人，能指导一下我，我苦于周围没有人交流，能问的人都问了，自己学校这边也没有人搞这个的，希望大家能帮帮我O(∩_∩)O~
如图，其中红色的是上包络线，黑色的为下包络线，绿色的为均值序列，蓝色为原信号

xray

这个问题是典型的包络拟合失败！

正常情况下，上包络线应该在原始信号的上方，而下包络线在原始信号的下方。而你的信号由于极值点个数很少（极大值点3个，极小值点2个），而且信号幅度变化剧烈，在利用三次样条对极值点插值的时候，出现了“过冲“的现象，也就是说下包络跑到了信号的上方，上包络跑到了信号的下方。

这个例子用事实说明了三次样条插值的缺陷！

如果可以的话，建议你把原始数据也贴上来，大家可以用来研究更好的包络插值算法。

yanling66520

本来是一系列信号的，因为全部贴上来，根本看不清楚，识别不了，所以就剪取了一个周期的。全部信号是采样频率为200hz，采样点数2048，具体如下:p=[71 70 69 68 68 68 66 65 64 63 61 62 61 61 61 62 60 59 60 60 59 59 60 58 58 58 57 57 56 56 56 55 56 56 57 59 63 69 77 83 95 103 113 123 132 141 148 156 164 168 172 176 175 176 175 174 172 170 166 163 160 156 151 147 143 139 135 131 128 124 121 117 115 111 109 108 104 103 101 99 99 96 94 93 91 89 88 86 85 83 81 81 81 81 81 81 81 83 83 85 86 87 89 90 92 94 95 96 98 99 101 101 100 101 101 101 100 99 98 97 96 96 93 92 90 89 87 85 84 82 81 80 77 76 75 74 72 71 70 70 68 67 66 65 64 62 62 61 61 60 60 60 59 58 58 59 58 58 58 58 58 57 58 58 57 58 57 57 56 56 57 56 56 57 58 61 63 67 74 80 89 97 107 117 127 137 145 152 158 163 167 170 172 173 174 174 171 169 166 162 159 156 152 148 144 141 136 132 128 125 121 118 115 113 111 108 107 104 103 102 100 99 97 96 95 93 91 90 88 87 86 84 84 84 85 85 86 87 88 89 90 92 93 94 96 98 100 101 102 103 104 104 105 105 105 105 104 104 103 102 100 99 98 96 94 92 91 89 88 86 85 83 81 81 79 79 77 76 75 74 73 72 72 71 70 69 69 69 68 67 66 66 66 65 65 64 63 63 63 64 64 63 63 62 62 61 61 60 61 61 62 63 64 67 72 79 87 96 105 114 123 133 142 150 158 163 169 174 177 179 180 180 178 176 174 171 168 164 160 156 152 147 146 141 136 134 131 128 125 122 120 117 116 114 112 111 109 108 106 105 103 102 100 99 97 96 95 94 92 91 91 91 91 91 92 93 94 95 97 98 100 101 102 103 104 106 105 108 107 107 107 107 107 106 106 104 104 103 101 100 99 98 95 94 92 91 89 88 87 85 84 82 81 80 78 78 77 75 74 74 73 72 71 71 70 70 69 68 67 67 66 67 66 65 65 65 65 64 64 64 64 65 64 64 63 63 63 63 62 62 60 61 61 61 62 64 68 73 80 88 97 107 116 126 135 145 153 161 167 172 176 179 181 182 182 181 179 176 173 171 167 163 159 155 152 149 145 141 138 135 132 130 127 124 122 120 118 117 116 115 113 112 111 110 108 107 106 103 103 101 100 99 97 95 93 94 94 94 95 96 96 97 98 98 99 100 101 103 104 105 105 105 105 105 105 105 105 104 103 103 103 102 100 99 97 96 95 94 92 92 90 88 87 85 84 83 81 80 79 78 77 76 75 73 73 72 71 71 70 69 68 69 67 67 66 66 66 66 65 65 65 64 64 63 64 64 64 64 64 65 64 64 64 64 63 63 62 61 61 61 61 61 61 62 65 68 74 82 88 97 107 117 125 135 144 152 158 164 169 173 176 177 179 179 178 175 173 170 167 163 159 156 152 148 145 142 138 135 132 130 127 125 122 121 119 117 116 115 113 112 111 109 109 107 106 105 104 103 101 100 98 97 96 96 94 95 94 95 95 95 95 96 98 98 99 101 101 102 103 103 102 103 104 103 103 102 102 101 100 98 98 97 95 94 93 91 90 89 88 85 85 84 83 81 80 79 77 76 75 74 73 72 71 69 68 66 67 66 66 65 64 63 61 61 59 58 58 57 56 55 55 55 54 54 55 53 53 53 54 53 54 55 55 54 53 54 52 51 51 51 51 50 51 51 51 53 54 59 65 72 81 91 100 110 119 128 136 144 150 156 161 165 167 168 167 167 165 163 160 158 154 151 147 143 139 134 133 129 126 123 120 117 114 112 110 108 106 105 104 103 101 99 98 97 96 95 94 92 90 89 88 86 83 85 82 83 82 83 83 84 84 86 86 88 89 90 90 91 92 93 93 94 95 95 95 95 94 93 93 92 90 90 88 87 87 85 83 82 81 79 78 76 75 73 72 71 70 69 67 66 65 64 64 64 63 61 61 60 59 58 58 57 57 57 56 55 55 54 53 52 53 52 52 53 52 52 52 51 51 51 52 51 51 51 50 52 51 51 51 51 52 51 51 51 50 51 50 52 49 49 49 49 48 48 50 52 57 63 70 78 86 96 105 114 124 133 141 148 153 158 162 164 165 166 166 165 163 161 158 154 151 146 142 138 135 132 128 124 121 118 116 112 109 108 105 103 102 100 99 98 96 94 94 93 92 91 89 88 85 85 84 83 81 81 80 81 81 82 83 84 85 86 87 89 90 92 93 95 95 96 97 97 98 98 97 96 96 95 94 94 94 91 91 88 87 85 83 82 80 78 77 75 74 72 71 70 69 67 66 67 65 65 64 64 61 62 61 59 59 59 59 58 58 56 56 56 55 54 55 55 55 55 55 54 55 54 53 54 53 54 55 54 54 55 54 55 54 54 54 54 54 53 53 52 51 50 51 51 50 51 52 55 59 65 72 81 89 99 108 117 126 133 141 148 154 157 161 165 166 167 166 165 162 160 157 154 150 147 142 138 135 131 127 124 121 118 115 112 110 107 106 104 103 101 100 98 96 95 94 93 91 90 88 88 87 85 84 82 82 82 82 82 83 83 86 86 88 88 90 91 93 95 97 97 98 99 99 100 101 100 101 101 100 101 99 98 96 94 94 93 91 89 88 86 84 83 81 80 78 77 76 75 73 72 72 70 69 67 67 66 65 65 64 64 63 62 61 61 61 60 60 59 59 58 58 59 57 58 57 58 58 57 57 56 56 56 56 57 57 58 58 58 57 56 57 57 57 57 56 56 56 56 54 55 55 54 55 56 57 60 65 73 81 90 100 109 119 129 139 146 154 160 166 171 174 176 178 179 179 177 175 172 170 166 162 159 156 151 148 144 141 137 134 131 128 126 123 121 119 117 115 114 113 112 111 109 108 106 104 102 102 99 99 98 96 96 95 93 91 91 91 93 94 95 95 96 97 97 98 100 102 103 104 105 106 105 105 105 105 106 107 107 106 105 103 101 100 98 98 97 95 94 92 92 89 88 86 84 83 82 81 79 80 77 75 74 73 72 72 71 69 68 68 68 66 66 64 65 64 64 64 63 63 62 63 62 61 61 61 61 61 60 61 61 61 61 60 60 60 61 61 61 62 62 62 61 61 60 59 59 59 58 58 58 58 57 58 59 59 63 68 74 82 91 101 111 120 131 140 149 157 163 170 175 178 181 183 183 184 183 181 179 176 173 170 165 161 158 154 151 148 144 141 139 136 133 129 128 125 123 122 121 119 118 117 115 114 114 112 111 110 109 107 105 104 103 101 101 100 99 99 99 99 99 99 100 101 102 103 104 105 107 106 108 108 108 109 109 110 110 110 109 109 108 107 106 106 105 103 101 99 97 95 94 93 91 90 88 87 85 84 81 81 79 79 79 77 76 75 74 72 70 70 69 68 68 67 67 66 65 65 63 63 62 63 63 63 63 63 62 62 62 61 61 60 62 60 60 60 61 61 61 61 61 61 60 60 60 60 60 59 59 59 58 57 57 56 57 58 61 64 71 78 86 95 105 115 125 135 144 153 161 166 172 176 180 181 183 183 182 181 179 175 172 169 165 161 157 154 150 146 143 139 136 133 130 128 124 122 120 119 116 115 113 112 110 108 107 107 104 105 102 101 99 96 96 95 93 93 92 91 92 91 92 92 93 94 94 97 98 99 100 101 102 103 104 104 105 106 107 106 106 105 104 102 102 101 99 97 96 95 92 90 88 86 86 85 83 81 79 78 76 74 73 72 71 70 68 67 66 65 64 64 63 63 61 61 60 59 58 58 58 57 58 57 57 57 56 56 56 56 56 56 56 55 55 55 56 57 57 56 56 57 57 57 57 57 56 56 55 55 56 55 55 56 58 60 65 70 77 85 94 104 114 122 132 141 149 156 162 166 170 174 176 176 176 175 175 171 169 165 161 158 154 149 145 142 139 136 132 129 125 122 120 117 115 114 113 111 110 109 106 105 104 103 102 101 100 98 96 95 93 92 91 90 90 90 90 92 92 93 94 96 97 99 101 102 104 105 106 107 108 109 110 111 111 110 110 109 109 108 107 106 105 104 102 101 98 97 94 94 92 91 90 88 87 85 84 83 82 81 80 79 78 78 76 76 74 74 73 74 72 72 71 71 70 70 69 69 68 68 68 68 68 68 68 68 68 69 68 67 68 68 67 66 66 66 67 66 66 65 65 67 66 68 72 78 84 92 99 109 118 127 136 144 153 160 165 172 176 179 181 182 182 181 179 176 174 171 168 165 160 156 153 147 144 142 139 134 132 129 127 124 123 120 118 117 116 114 113 111 110 109 108 107 106 104 102 101 99 100 97 97 97 97 97 97 98 99 100 102 104 106 107 108 109 110 111 111 112 113 113 113 113 113 113 114 112 111 110 109 107 106 104 104 101 100 98 96 93 92 91 90 89 88 87 85 84 82 81 79 78 78 79 77 76 75 75 73 73 72 71 71 70 71 70 68 68 67 67 67 68 68 69 69 68 68 68 68 68 68 69 68 68 67 66 66 66 66 66 67 67 71 74 78 84 93 103 111 121 131 140 149 157 165 171 176 181 184 185 186 186 185 185 182 180 178 173 170 166 162 158 155 151 148 145 142 140 137 134 131 129 127 127 124 123 122 121];

gzslbs2003

是不是你下的emd程序有问题？
如果需要，我这里有emd的原程序

xray

pp = conv(ones(1,20),p);
pp = pp(20:end-20);
稍微处理一下，进行平滑滤波，减少了部分噪声，包络效果好了一些。

yanling66520

我试做了下，效果果真好很多，谢谢你的帮忙，顺便问下，平滑处理之后，原信号变了，那我接下来再对平滑过的信号进行hht，对结果会不会有影响，还是我赋值变量改一下，之后的处理都要用我最原初的信号呢？

yanling66520

谢谢你的热心回复，我想问下你的EMD程序是哪个版本的？如果方便的话，也希望传我一份，我看一下，我这里也有好几个版本，不知道怎么区分，有些还是做不出来

xray

我用的是G. Rilling的那个版本，其实不同版本差别不大，关键是对于EMD算法本身的理解。

至于平滑以后是否会对处理结果产生影响，关键是看你需要获得的信息是什么。平滑肯定是会损失信息，但是只要你需要的关键信息还保留就可以了。

宁静致远

这世界太美好了，谢谢！！

yanling66520

你好，我用plot_hht的程序分析自己的信号，emd部分已经完成了，但是后面的hilbert变换不怎么会，hilbert的原理我知道的，imf的频率图不会画（幅值——频谱），边际谱也不知道怎么得到的。本来想自己编程序把那些谱图画出来，结果昨天研究了一天也没出结果，最后把IMF分量用FFT变换了下，得到一些图，大概也算是做了无用功吧。另外，我在网站上下载的G. Rilling版本的程序，不怎么会用，我昨天也在论坛里看到你对这个版本程序的注解，但是我似乎根本就不知道这个程序下载下来是如何来运用的，所以即使看到了你的注解版本，我依旧是云里雾里。谢谢你的帮助，希望得到您的回复O(∩_∩)O~
下面我附上我用FFT对IMF做出的频谱图

小笨蛋0516

谢谢 正在学习中