%------------------------------S&C算法----------------------------
%兩個序列，使用已知的差分序列在接收端循環位移，
clc;
clear?all;
%_-----------------發送------------
N=128;

norm_fds=0:1:16;

snr=10;

f=zeros（1length（norm_fds））;
%?for?s=1:length（snrs）
%?????snr=snrs（s）;
for?s=1:length（norm_fds）
????norm_fd=norm_fds（s）;

a1=zeros（1N）;
a1（1:2:N/2）=randint（1N/4[19]）+j*randint（1N/4[19]）;

A=a1?;%發送頻域序列
C=ifft（A）;
len=N;
%?%--------------信道-加頻偏和噪聲-----------------
C=C.*exp（sqrt（-1）*2*pi*norm_fd.*[0:len-1]/N）;%變成256數據做fft的點數，
%???plot（abs（C））;%時域里前后兩個部分是相同的

%?????????????%———————————————加噪聲———————
%?????????????signal_after_fft1?=?fft（C（1:N/2））;??
%?????????????signal_after_fft2?=?fft（C（N/2+1:N））;?
%?????????????signal_power?=?mean（mean（signal_after_fft1.*?conj（signal_after_fft2）））;
%?
%?????????????%-----------------------------niose?power----------------------------------
%?????????????noise_power?=?signal_power?./?10^（snr/10）;
%?????????????noise_voltage?=?sqrt（noise_power?./?2）;
%?
%?????????????%?-------------------------------------------------------------------------
%?????????????len_noise?=?length（N）;
%?????????????noise0?=randn（1N）;%zeros（1N）;%
%?????????????noise?=?noise_voltage?*?（noise0?+?sqrt（-1）?.*?noise0）;
%?????????????C?=C?+?noise;

%--------------小數倍頻偏估計----時域進行的-----------
p=sum（conj（C（1:N/2））.*C（N/2+1:N））;
f1=angle（p）/pi;
%--------------------補償------------------
C=?C.*exp（sqrt（-1）*（-1）*2*pi*f1.*[0:len-1]/N）;%

%----------------接收-頻域的--------------
D=fft（C）;%接收頻域序列
power=sum（abs（D）.^2）;
%------------------整數倍頻偏----------------------------
for?g=0:N/2-1
B=circshift（A[0?g]）;%上下不移動，向右移動g位
r（g+1）=sum（B.*conj（D））;??
F（g+1）=?r（g+1）.^2/power;
end
[F1g1]=max（F）;%g=N/2時，循環移位后的序列與之前的序列相等，有最大值，即沒有頻偏時
g0=mod（g1N/2）;
%-----------------------------------------------------
f（s）=g0-1+f1
end
figureplot（norm_fdsf‘r-o‘）
legend（‘估計值‘）
hold?on
grid?on







?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????290218??2013-03-18?13:26??頻偏估計S&C方法及實現matlab\Robust?Frequency?and?Timing?同步OFDM系統.pdf

?????文件???????2112??2013-03-18?15:14??頻偏估計S&C方法及實現matlab\SandC.m

?????目錄??????????0??2013-03-18?21:30??頻偏估計S&C方法及實現matlab

-----------?---------??----------?-----??----

???????????????292330????????????????????3