残余频率估计

现有这样一个通信系统的载波同步，在载波同步后，还有残余50HZ左右的频率，用M次幂的方法将残余的频率估计出来。具体MATLAB仿真是这样做：
1.如果将基带信号乘以EXP(1i*2*pi*50*t)（假设残余频率为50HZ），通过莱斯信道并且将Fd设为0,那么采用M次幂去除调制信息后，做FFT，求最大峰值对应的频率，那么可以得到其频率为50HZ；但是但将Fd设为100，估计的频率就不对了。
  -------- 问题：当莱斯信道的Fd小于100HZ估计得残余频率还是准的，但实际Fd很可能大于100HZ，也就是说Fd对残余频率估计影响很大，还有什么好的算法能减小Fd带来的影响。
2. matlab中莱斯信道函数ricianchan(Ts,Fd,K)中的Fd说是最大多普勒频率，那为什么我将基带信号通过莱斯信道时（假设频率为50HZ）检测不出来这个频率。（我是这么做的：基带调制方式为QPSK，复基带信号经过莱斯信道，接收端采用M次幂去除调制信息，然后再做FFT，求最大峰值对应的频率）。
  -------问题：按matlab自带的描述，Fd为最大多普勒频率，应该可以检测通过M次幂方法检测出来的,为什么不行尼。

M次幂方法 是什么方法？

比如QPSK调制，4个星座点，4次幂后，调制信息就被去掉了，剩下的就是载波信息。

传说中有个Gardner算法，试试能不能解决你的问题
另外，有多普勒的话什么方法都不准了吧 因为多普勒可以理解成无限多个频偏的叠加

这种FFT后通过峰值求频偏的方法精度很差的。

回复haojia：
传说中有个Gardner算法，试试能不能解决你的问题
  ----------Gardner算法是做符号同步用的，并且是一般是数据连续传输场景，不太适合短突发的符号同步情况。
另外，有多普勒的话什么方法都不准了吧 因为多普勒可以理解成无限多个频偏的叠加
------------感觉是这样的。
回复：acarpente
    这种FFT后通过峰值求频偏的方法精度很差的。
----------------FFT点数比较大，基带速率不大，所以频率分辨率很细，我试过几种算法，这种事最好的。
问题：不知道问什么在MATLAB仿真过程中，经过莱斯信道后，还要再乘以多普勒频率，才可以通过算法给检测出来。  

莱斯信道是有个主径的吧，会不会matlab默认主径是没有多普勒的？你测到的准的时候是主径足够强的时候？