jens/matlab
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

- auto compensate cyclic prefix

Browse Source 
git-svn-id: http://moon:8086/svn/matlab/trunk@25 801c6759-fa7c-4059-a304-17956f83a07c
This commit is contained in:
jens
2015-04-05 14:07:25 +00:00
parent 81650fd21b
commit bddaaee2d5
2 changed files with 29 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files
ofdm/fdet_wav.m
+25 -5
View File
@@ -28,21 +28,41 @@ else
x = xwav.*exp(-j*(2*pi*fc/fs*(0:len-1)'));
end
fprintf('Perform frequency correction\n');
x = x.*exp(j*(2*pi*foff/fs*(0:len-1)'));
[omega, err, dmod, omega_coarse, omega_est] = fdet(x, N, Ng, omega_min, omega_step, omega_max);
% Remove cyclic prefix
fprintf('Perform CP calculation\n');
x0 = x(1:2*(N+Ng));
x1 = x(N+Ng+1:3*(N+Ng));
c = xcorr(x0, x1);
[v, ii] = max(abs(c));
cp_off = mod(ii, Ng+N)
% Detect CFO
fprintf('Perform frequency detection\n');
[omega, err, dmod, omega_coarse, omega_est] = fdet(x, N, Ng, cp_off, omega_min, omega_step, omega_max);
f_coarse = omega_coarse*fs
f_est = omega_est*fs
close all;
subplot(4,1,1)
plot(abs(c)); grid;
figure;
subplot(5,1,1)
plot(1:lge(omega), real(dmod)/N, 1:lge(omega), imag(dmod)/N, 1:lge(omega), abs(dmod)/N); grid; legend('Re', 'Im', 'abs()');
subplot(4,1,2)
subplot(5,1,2)
plot(1:lge(omega), err*fs, '-'); grid; legend('error');
subplot(4,1,3)
subplot(5,1,3)
plot(1:lge(omega), omega*fs); grid; legend('f');
subplot(4,1,4)
subplot(5,1,4)
plot(1:lge(omega), angle(dmod), '-', 1:lge(omega), diff(angle([0 dmod])), 'r-'); grid; legend('phi', 'dphi');
subplot(5,1,5)
plot(1:lge(omega), abs(dmod)/N - real(dmod)/N, '-'); grid; legend('diff');