前言
只是看到很多地方都有说到单载频脉冲信号的时宽带宽积为1,却没有具体的推导,在这里推导一下试试。
一、单载频脉冲信号
令单载频脉冲信号为f(x)f(x)f(x),那么有f(x)=gτ(x)cos(ω0x)f(x)=g_\tau(x)cos(\omega_0x)f(x)=gτ(x)cos(ω0x),gτ(x)g_\tau(x)gτ(x)为宽度为τ\tauτ的门函数,ω0\omega_0ω0为载频。那么具体波形如下。
代码:
tau = 0.02;
w0 = 2000;
t = -0.02:0.0001:0.02;
f = (heaviside(t+0.01)-heaviside(t-0.01));
plot(t,f);
f1 = (heaviside(t+0.01)-heaviside(t-0.01)).*cos(w0*t);
hold on
plot(t,f1,'--');
title('w0 = 2000, 门宽 = 0.2');
xlabel('t');
ylabel('f');
时宽带宽积推导
对f(x)f(x)f(x)利用频域卷积定理,有
F(jω)=12πGτ(jω)∗π(δ(ω−ω0)+δ(ω+ω0))F(j\omega)=\frac{1}{2\pi} G_\tau(j\omega)*\pi(\delta(\omega-\omega_0)+\delta(\omega+\omega_0))F(jω)=2π1Gτ(jω)∗π(δ(ω−ω0)+δ(ω+ω0))
即
F(jω)=12{Gτ(j(ω−ω0))+Gτ(j(ω+ω0))}F(j\omega)=\frac{1}{2}\{G_\tau(j(\omega-\omega_0))+G_\tau(j(\omega+\omega_0))\}F(jω)=21{Gτ(j(ω−ω0))+Gτ(j(ω+ω0))}
对于门函数gτ(t)g_\tau(t)gτ(t)而言,它的傅里叶变换结果Gτ(jω)G_\tau(j\omega)Gτ(jω)为
Gτ(jω)=τSa(ωτ2)G_\tau(j\omega)=\tau Sa(\frac{\omega\tau}{2})Gτ(jω)=τSa(2ωτ)
从而
F(jω)=τ2{Sa(τ2(ω−ω0))+Sa(τ2(ω+ω0))}F(j\omega)=\frac{\tau}{2}\{Sa(\frac{\tau}{2}(\omega-\omega_0))+Sa(\frac{\tau}{2}(\omega+\omega_0))\}F(jω)=2τ{Sa(2τ(ω−ω0))+Sa(2τ(ω+ω0))}
现在结果很明显了,单载频脉冲函数的频谱就是两个SaSaSa函数偏移了ω0\omega_0ω0个单位而已。
SaSaSa函数结果如下所示
代码:
tau = 0.02;
w0 = 2000;
w = -20000:0.01:20000;
%%F = tau/2*( ((sin(tau/2*(w-w0)))./(tau/2*(w-w0))) + ((sin(tau/2*(w+w0)))./(tau/2*(w+w0))));
F = tau/2*( ((sin(tau/2*(w-w0)))./(tau/2*(w-w0))));
plot(w,F);
title('单边Sa函数');
xlabel('频率w');
ylabel('F');
带宽仅仅看正频率区即可,所以在上图我只画了右半频率区。带宽计算如下所示。
在上课时,我记得老师讲过,SaSaSa函数的−3dB-3dB−3dB带宽和主瓣宽度的一般是差不多相等的,也就是SaSaSa函数取到峰值和取到000之间的距离,虽然优点误差,但是大致误差不超过10%(279与304)10\%(279与304)10%(279与304)。同时载频对带宽并不起作用,于是单载频脉冲信号的带宽就是门函数的带宽。这样就可以利用公式计算,如下,
Sa(ωτ2)=0Sa(\frac{\omega\tau}{2})=0Sa(2ωτ)=0,即ωτ2=π\frac{\omega\tau}{2}=\pi2ωτ=π,则ωτ=2π\omega\tau=2\piωτ=2π换算成频率,左右同时除以2π2\pi2π,得fτ=1f\tau=1fτ=1,此时B=f−0=fB=f-0=fB=f−0=f,所以Bτ=1B\tau=1Bτ=1,大致如此。
总结
仅仅是自行推导,有问题的话欢迎大家提出。