资源简介
整数和小数倍频偏估计 S&C方法 matlab 实现
代码片段和文件信息
%------------------------------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
- 上一篇:四维超混沌及其图形
- 下一篇:MATLAB中计算psnr值的实现
评论
共有 条评论