MATLAB绘制包络线可以通过hilbert或者envelope实现。Hilbert变换是一个很有用的变换,用它来做包络分析更是一种有效的数据处理方法。 实例1 程序clc;clearall;closeall;t0:0。01:6pi;ysin(t)。sin(9。t);hhilbert(y);Hilbert变换y13cos(t)。sin(9。t);h1hilbert(y1);Hilbert变换figure;subplot(2,1,1);plot(t,y);holdonplot(t,abs(h),r,linewidth,1。5)holdonplot(t,1abs(h),r,linewidth,1。5)gridonsubplot(2,1,2);plot(t,y1);holdonplot(t,abs(h1),g。,linewidth,1。5)holdonplot(t,1abs(h1),g。,linewidth,1。5)gridon 运行结果 实例2 程序Hilbert变换测试clcclearallcloseallts0。001;fs1ts;N200;f50;k0:N1;tkts;信号变换结论:sin信号Hilbert变换后为cos信号ysin(2pift);yhhilbert(y);matlab函数得到信号是合成的复信号yiimag(yh);虚部为书上定义的Hilbert变换figuresubplot(211)plot(t,y)title(原始sin信号)subplot(212)plot(t,yi)title(Hilbert变换信号)检验两次Hilbert变换的结果(理论上为原信号的负值)结论:两次Hilbert变换的结果为原信号的负值yihhilbert(yi);yiiimag(yih);max(yyii)信号与其Hilbert变换的正交性结论:Hilbert变换后的信号与原信号正交sum(y。yi)谱分析结论:Hilbert变换后合成的复信号的谱没有大于奈氏频率的频谱,即其谱为单边的NFFT2nextpow2(N);ffslinspace(0,1,NFFT);Yfft(y,NFFT)N;YHfft(yh,NFFT)N;figuresubplot(211)plot(f,abs(Y))title(原信号的双边谱)xlabel(频率f(Hz))ylabel(Y(f))subplot(212)plot(f,abs(YH))title(信号Hilbert变换后组成的复信号的双边谱)xlabel(频率f(Hz))ylabel(YH(f)) 运行结果 实例3 程序包络分析(高中心频率的窄带信号分析)基于:两个信号乘积的Hilbert变换取决于高频信号的Hilbert变换clcclearallcloseallts0。001;fs1ts;N200;k0:N1;tkts;原始信号f110;f270;acos(2pif1t);包络1a2cos(2pif1t);包络2a1。(1t。250);包络3msin(2pif2t);调制信号ya。m;信号调制figuresubplot(241)plot(t,a)title(包络)subplot(242)plot(t,m)title(调制信号)subplot(243)plot(t,y)title(调制结果)包络分析结论:Hilbert变换可以有效提取包络、高频调制信号的频率等yhhilbert(y);aabsabs(yh);包络的绝对值aangleunwrap(angle(yh));包络的相位afdiff(aangle)2pi;包络的瞬时频率,差分代替微分计算NFFT2nextpow2(N);NFFT2nextpow2(10244);改善栅栏效应ffslinspace(0,1,NFFT);YHfft(yh,NFFT)N;Hilbert变换复信号的频谱Afft(aabs,NFFT)N;包络的频谱subplot(245)plot(t,aabs,t,a,。)title(包络的绝对值)legend(包络分析结果,真实包络)subplot(246)plot(t,aangle)title(调制信号的相位)subplot(247)plot(t(1:end1),affs)title(调制信号的瞬时频率)subplot(244)plot(f,abs(YH))title(原始信号的Hilbert谱)xlabel(频率f(Hz))ylabel(YH(f))subplot(248)plot(f,abs(A))title(包络的频谱)xlabel(频率f(Hz))ylabel(A(f)) 运行结果 envelope绘制包络线 程序clc;clearall;closeall;t0:0。01:6pi;ysin(t)。sin(9。t);figure;plot(t,y);holdon包络线绘制〔up,lo〕envelope(y);plot(t,up,t,lo,linewidth,1。5)legend(signal,上包络线,下包络线)holdoff 运行结果 本文部分内容来源于长笛人倚楼Gloria的博客,仅供参考学习,如内容、图片有任何版权问题,请联系处理,24小时内删除。 作者郭志龙 编辑郭志龙 校对郭志龙