希尔伯特变换Hilbert Transform

希尔伯特变换公式

x^(t)=H[x(t)]=x(t)∗1πt\hat x(t)=\mathcal{H}[x(t)]=x(t)*\frac{1}{\pi t}x^(t)=H[x(t)]=x(t)πt1
也就是说,信号的希尔伯特变换就是把它经过一个冲激响应h(t)=1πth(t)=\frac{1}{\pi t}h(t)=πt1的LTI线性时不变系统得到的输出,此系统的频率响应:
H(jw)=−jsgn(w)={ −j,w&gt;0j,w&lt;0H(jw)=-jsgn(w)=\left\{ \begin{aligned} -j,w&gt;0\\ j,w&lt;0 \end{aligned} \right.H(jw)=jsgn(w)={ j,w>0j,w<0

so,显然,如下图,它的幅频特性全为1,是一个全通滤波器。
而对于相位,对于正频,它相移−π2-\frac \pi 22π,(e−jπ2=−je^{-j\frac \pi 2}=-jej2π=j);
对于负频率部分,相移π2\frac \pi 22π
在这里插入图片描述
画图代码:

clear all,clc;
w=-2*pi:0.01:2*pi;
i=sqrt(-1);
h=-i*sign(w);
figure,
subplot(211)
plot(w,abs(h));
axis([-2*pi 2*pi 0 2])
xlabel('w'),ylabel('|H(jw)|'),title('Hilbert滤波器的幅频响应'),grid on
subplot(212)
plot(w,angle(h));
axis([-2*pi 2*pi -pi pi])
xlabel('w'),ylabel('\psi(w)'),title('Hilbert滤波器的相频响应'),grid on

对比无失真系统的频率响应,可知这个系统是有失真的,会通过相移引入新的频率成分。

无失真系统的输入输出关系必须是r(t)=Ks(t−t0)r(t)=Ks(t-t0)r(t)=Ks(tt0),其中s(t)s(t)s(t是输入,r(t)r(t)r(t)是输出,无失真即只能对输入有时延和幅度缩放,根据傅里叶变换求得系统响应H(jw)=Ke−jwt0H(jw)=Ke^{-jwt_0}H(jw)=Kejwt0

在这里插入图片描述
画图代码:

clear all,clc;
w=-2*pi:0.01:2*pi;
j=sqrt(-1);K=1;t0=1;
H=K*exp(-j*w*t0);
subplot(211)
plot(w,abs(H));
axis([-2*pi 2*pi 0 2])
xlabel('w'),ylabel('|H(jw)|'),title('无失真传输系统的幅频响应 K=1;t0=1;'),grid on
subplot(212)
plot(w,angle(H));
axis([-pi pi -pi pi])
xlabel('w'),ylabel('\psi(w)'),title('无失真传输系统的相频响应 K=1;t0=1;'),grid on

一些有趣的性质

  1. 逆变换
    因为−H(jw)H(jw)=−j2sgn2(w)=1-H(jw)H(jw)=-j^2sgn^2(w)=1H(jw)H(jw)=j2sgn2(w)=1
    所以
    H−1[]=−H[]\mathcal H^{-1}[ \quad]=-\mathcal H[\quad]H1[]

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值