2020-11-14

**

Matlab——Hilbert变换之包络解调过程中的一点思考

今天在做包络解调的过程中遇到了一个问题，困扰了一下午，在此记录下来，供大家交流（请大神指导）。

**

两种包络解调计算方法

数字信号的包络解调在matlab中可以用Hilbert变换来实现。常用到的方法如下：

clc;clear all;close all;
fs = 2000;
f0 = 80;
t = 0:1/fs:.1;
sig = sin(2*pi*f0*t);
sig_ref = -cos(2*pi*f0*t);
sig_hilbert = hilbert(sig);
figure;
plot(t,sig,'k','linewidth',2);hold on;
plot(t,imag(sig_hilbert),'b','linewidth',2);hold on;  %matlab中Hilbert变换的结果为复数形式，而教科书中的Hilbert结果为复数的虚部。
legend('原信号','Hilbert变换信号','localization','best');
ylim([-2,2]);
aabs = abs(sig_hilbert);  
figure
plot(t,aabs,'k'); hold on
plot(t,sig,'r');

运行结果如下图

![可以看到原始信号的Hilbert变换结果，如蓝色线条所示](https://img-blog.csdnimg.cn/20201114192629621.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L0ZmMTg4OA==,size_16,color_FFFFFF,t_70#pic_center)

提取到的包络效果还算可以。
然而，此时如果采用如下这种方法求取包络的话：
依据博客中的理论：复数x=a+bi，abs(x)=sqrt(a^2 +b^2)*

https://blog.csdn.net/you23hai45/article/details/13093859

env=sqrt(sig.^2+imag(sig_hilbert).^2); %sig的
figure
plot(t,env,t,sig);

图中的蓝色线条为包络，红色线条为原始信号

至此，这两种计算方法理论上应该是相同的，但是实际情况是什么样的呢？

问题进一步分析

采用程序调试的方法来看的话，在matlab的hilbert变换环节就产生了疑问。

a=real(sig_hilbert);
aa=a-sig;
figure
plot(t,aa);

理论上，计算Hilbert的实部结果应该和原始信号的值保持一致，然而实际情况是二者有一定的误差，虽然误差很小，在e-15级别，但是确实不一样。如下图所示

进一步，采用sig和sig_hilbert的实部分别计算信号的包络

env=sqrt(sig.^2+imag(sig_hilbert).^2); %sig的
env1=sqrt(real(sig_hilbert).^2+imag(sig_hilbert).^2); %sig的
figure
plot(t,env,'--*r',t,sig); hold on
plot(t,env1,'--b','linewidth',2);
legend('方法一','原信号','方法二');

可以看到，基本上两种方法计算结果是一致的。
最后，我感觉应该是由于matlab中内置的实部函数代码存在一定的误差，不过影响很小而已。
如果有大神知道原因，欢迎讨论