teager.zip_HRY_Teager 能量解调_energy operator_teager 解调_teager能量算子；

% function testPDESA clc clear close all fs=5120; Ns=2048; t=(0:Ns-1)/fs; x2=exp(-5*t).*cos(2*pi*20*t); % 实际包络 x=x2.*sin(2*pi*500*t); [y,f]=DESA(x); % 能量算子解包络 yhilbert=abs(hilbert(x)); % Hilbert解包络 figure subplot(2,1,1), plot(t,[y; yhilbert; x]) title('解包络') legend('能量算子包络','Hilbert包络','原始数据'); grid on subplot(2,1,2), plot(t,[abs(x2)-y;abs(x2)-yhilbert]) grid on axis([0 0.4 -0.02 0.02]); legend('能量算子包络误差','Hilbert包络误差'); title('解包络误差') figure plot(abs(f)/2/pi*fs) title('能量算子解调得到的频率')