zl程序教程

您现在的位置是:首页 >  后端

当前栏目

第十二周-总结符号速率估计的主要方法以及循环谱估计参数的方法

方法循环 总结 以及 参数 主要 符号 估计
2023-09-11 14:20:47 时间

主要工作

1.进行PSK的仿真和QAM的仿真。
2. 了解符号速率的估计方式。基于循环谱估计、基于小波变换、基于瞬时幅度。
3. 着重了解了基于循环谱的符号速率估计。


一、PSK/DPSK

链接

二、QAM

链接

三、符号速率估计

1.瞬时频率法

下面给出了符号速率估计几种算法的优缺点。其中前两种方法上周已有了解,这周用瞬时频率法仿真了一下,但也出现一些问题。
在这里插入图片描述
下面是两种方法的具体过程,留在仿真完PSK\DPSK后,适用这两种方法看能不能实现符号速率的估计。
在这里插入图片描述
我主要是把公式带进去仿真的PSK中,看结果,只不过很奇怪的是峰值一直是多进制M的倍数,而和符号速率无关。由于时间关系,只找到了引用文献,没有来得及公式推导和修正。

N=1024;Nfft=1024;
window=hamming(N); noverlap=N/2;
seta=angle(x1);
ft=diff(seta)/(2*pi)+fc;
[Pxx2,f2]=pwelch(ft,window,noverlap,Nfft);
figure;
subplot(211);
plot(f2*fs/(2*pi),10*log10(Pxx2));title('瞬时频率的功率谱估计');
subplot(212);
f=(0:N-1)*fs/N;
f3=abs(fft(ft,N));
plot(f,10*log10(f3));title('瞬时频率的傅里叶变换');
xlabel('真实频率');

在这里插入图片描述
此时峰值位置是M=8的倍数的位置,当M=16 32也都是M倍数的位置。

2.循环谱法

  1. 循环谱定义
    在这里插入图片描述
  2. 通过循环谱找到载波频偏和符号速率
    先求出根据上述推导过程求出循环谱,如下图所示。在令f=0,此时峰值位置为二倍的频偏;得到频偏fc后,再另f=fc,得到的截面图离y轴最近的峰位置就表示符号速率。
    在这里插入图片描述
    下图是我用PSK仿真的信号做的循环谱图(M=2)
    在这里插入图片描述

四、学习计划

1.验证一下瞬时幅度法的符号速率估计。
2.验证一下循环谱估计的方法。
3.OFDM仿真

参考文献

1.非合作通信信号侦察算法研究与实现_靳鑫
2.一种基于瞬时幅度的PSK信号符号速率估计方法_张晓勇
3.基于循环谱的BPSK信号符号速率盲估计_杨蛟龙