%------------------------------S&C算法----------------------------
%兩個(gè)序列，使用已知的差分序列在接收端循環(huán)位移，
clc;
clear?all;
%_-----------------發(fā)送------------
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?;%發(fā)送頻域序列
C=ifft（A）;
len=N;
%?%--------------信道-加頻偏和噪聲-----------------
C=C.*exp（sqrt（-1）*2*pi*norm_fd.*[0:len-1]/N）;%變成256數(shù)據(jù)做fft的點(diǎn)數(shù)，
%???plot（abs（C））;%時(shí)域里前后兩個(gè)部分是相同的

%?????????????%———————————————加噪聲———————
%?????????????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;

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

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







?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????290218??2013-03-18?13:26??頻偏估計(jì)S&C方法及實(shí)現(xiàn)matlab\Robust?Frequency?and?Timing?同步OFDM系統(tǒng).pdf

?????文件???????2112??2013-03-18?15:14??頻偏估計(jì)S&C方法及實(shí)現(xiàn)matlab\SandC.m

?????目錄??????????0??2013-03-18?21:30??頻偏估計(jì)S&C方法及實(shí)現(xiàn)matlab

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

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