开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > matlab例程 > 查看源码
OFDM_LS_directly.m
资源名称：信道估计算法.rar　[点击查看]
上传用户：dgwtdq
上传日期：2013-04-19
资源大小：2k
文件大小：4k
源码类别：

matlab例程

开发平台：

Matlab

OFDM_LS_directly.m：源码内容
  1. %16QAM-OFDM system;
  2. clear all;
  3. clear all;
  5. N = 1024;                                                        %the total number of all subcarriers
  6. Np = 128;                                                         %the number of pilot carriers
  7. B = 2e6;                                                         %bandwidth 2MHz
  8. T = 1/B;
  9. Ng = 256;                                                        %the length of Guard Interval
  10. Xp = 3 + 1.0*j ;                                                      %pilot value
  11. t= 0;  % erro code number
  13. for u = 1:20                                             % run 20 times
  14. %-------------------------------------Data generator--------------------------------------------------
  15. xr1 = rand(1,N-Np);
  16. Xre = (-3) * (xr1 < 0.25) + (-1) * (0.25 <= xr1 & xr1 < 0.5)...
  17.     + (0.5 <= xr1 & xr1 < 0.75) + 3 * (0.75 <= xr1 & xr1 < 1);
  18. xr2 = rand(1,N-Np);
  19. Xim = (-3) * (xr2< 0.25) + (-1) * (0.25 <= xr2 & xr2 < 0.5)...
  20.     + (0.5 <= xr2 & xr2 < 0.75) + 3 * (0.75 <= xr2 & xr2 < 1);
  21. X = complex(Xre,Xim);                                            %16QAM
  22. % Xre = (-1) * (xr1 < 0.5) + 1 * (0.5 <= xr1 & xr1 < 1);
  23. % xr2 = rand(1,N-Np);
  24. % Xim = (-1) * (xr2< 0.5) + 1 * (0.5 <= xr2 & xr2 < 1);
  25. % X = complex(Xre,Xim);                                            %QPSK
  26. %----------------------------------  insert pilot in carriers    -------------------------------------
  27. % pilots locate at 8*i + 3;
  28. Xc = zeros(1,N);                                                 %the carrier value
  29. m = 0;
  30. for i = 1 : N
  31.     if(mod(i,8) == 3)
  32.          Xc(i) = Xp;
  33.          m = m + 1;
  34.     else
  35.         Xc(i) = X(i - m);
  36.     end
  37. end
  38. %----------------------------------OFDM modulation-------------------------------------------------------
  39. Ym = ifft(Xc,N);
  40. %----------------------------------Insert GI  CP---------------------------------------------------------
  41. Y = [Ym(769:N),Ym];
  42. %----------------------------------pass channel-----------------------------------------------------------
  43. Yc = ray_channel(Y,T);
  44. %---------------------------------add noice -------------------------------------------------------------
  45. SNR = 40;
  46. % for S=1:21
  47.   %SNR= 10log10(Es/N0)
  48. % N0 = var(Yc)/(10^(SNR/10));
  49. % sgma = sqrt(N0/2);
  50. % Yawgn=Yc+sgma*randn(1,length(Yc))+sgma*randn(1,length(Yc))*j;
  51. Yawgn = awgn(Yc,SNR,'measured');
  52. %-------------------------------------------------------------------------------------------------------     
  53. % ----------------------------------------------OFDM RECEIVER-------------------------------------------  
  54. %----------------------------------------------OFDM demodulation----------------------------------------
  55. Yr = fft(Yawgn,N);
  56. %-----------------------------------------------Pilot signal Extraction---------------------------------
  57. c = 1;
  58. for i = 3:8:N
  59.     Yp(c) = Yr(i);
  60.     c = c + 1;
  61. end
  62. %----------------------------------- Pilot signal Estimation with LS Algorithm -------------------------   
  63. Hp = Yp ./ Xp;
  64. % figure(1)
  65. % plot(abs(Hp))
  67. %---------------------------------Interpolation Appoach-----------------------------------------------------
  68. pilot_loca = 3:8:N;
  69. Hk = interp1(pilot_loca,Hp,1:N);
  70. H =  fft(Yc,N) ./ fft(Ym,N);
  72. figure(3)
  73. plot(abs(Hk))
  74. hold on
  75. plot(abs(H),'r')
  76. %-------------------------------------data correction ------------------------------------------------------
  77. Xcc = Yr ;
  78. Xcc_line = Yr ./ Hk;
  79. figure(4)
  80. plot(Xcc_line,'+')
  81. hold on
  82. plot(X,'r*')
  83. %----------------------------------------16QAM demodulation--------------
  84. Xdre = zeros(1,N);
  85. Xdim = zeros(1,N);
  86. for i = 1 :N
  87.     if(real(Xcc_line(i)) < -2)
  88.         Xdre(i) = -3;
  89.     elseif((real(Xcc_line(i)) >= -2)&(real(Xcc_line(i)) < 0))
  90.         Xdre(i) = -1;
  91.     elseif((real(Xcc_line(i)) >= 0)&(real(Xcc_line(i)) < 2))
  92.         Xdre(i) = 1;    
  93.     else
  94.         Xdre(i) = 3;
  95.     end
  96. end
  98. for i = 1 :N
  99.     if(imag(Xcc_line(i)) < -2)
  100.         Xdim(i) = -3;
  101.     elseif((imag(Xcc_line(i)) >= -2)&(imag(Xcc_line(i)) < 0))
  102.         Xdim(i) = -1;
  103.     elseif((imag(Xcc_line(i)) >= 0)&(imag(Xcc_line(i)) < 2))
  104.         Xdim(i) = 1;    
  105.     else
  106.         Xdim(i) = 3;
  107.     end
  108. end
  110.  Xd = complex(Xdre,Xdim); 
  111. Xq = zeros(1,N-Np);
  112. m = 0;
  113. for i = 1 : N;
  114.     if(mod(i,8) == 3)
  115.         m = m + 1;
  116.    else
  117.     Xq(i - m) = Xd(i) ;
  118.     end
  119. end
  122. for i= 1:Np
  123.     if(Xq(i) ~= X(i))
  124.         t = t + 1;
  125.     end
  126. end
  128. end
  129. BER = t/(u * (N-Np));
  131. % Xq(1:10)'
  132. % X(1:10)'
  133. %Xcc_lowpass(1:10)
  134. BER
  135. SNR
  136. % t = t/u
  138. % SNR
  139. % figure(10)
  140. % plot(SNR,abs(BER))