寻找包络线

yangshangbin

close all
fs=30;
t=0:1/fs:200;
x6=sin(2*pi*2*t)+sin(2*pi*4*t);
x6 = x6 + randn(1, length(x6));
% 包络算法，未考虑边界条件
d = diff(x6);
n = length(d);
d1 = d(1:n-1);
d2 = d(2:n);
indmin = find(d1.*d2<0 & d1<0)+1;
indmax = find(d1.*d2<0 & d1>0)+1;
envmin = spline(t(indmin),x6(indmin),t);
envmax = spline(t(indmax),x6(indmax),t);
figure
hold on
plot(t,x6);
plot(t,envmin,'r');
plot(t,envmax,'m');
hold off
xlim([0 9])

matlab的小调

对于有毛刺的情况效果不好啊,怎么去噪声呢???

yangshangbin

matlab的小调 发表于 2015-5-27 14:27
对于有毛刺的情况效果不好啊,怎么去噪声呢???

先对原始数据做滤波，然后再寻找包络。

风子师兄

楼主MATLAB哪个版本的，我按照你的代码运行不行哦，envmin = interp1(t(indmin),x6(indmin),t,'spline'); 用interp1插值才能出包络的效果

admin

风子师兄 发表于 2015-6-11 15:26
楼主MATLAB哪个版本的，我按照你的代码运行不行哦，envmin = interp1(t(indmin),x6(indmin),t,'spline');  ...

不会吧,你的版本spline都没有啊