Category

Development Platform
Home > Document > Communication Document > View Code
gmsk.m
Package: GMSK-simulink.rar　[view]
Upload User: jy88362166
Upload Date: 2015-04-26
Package Size: 3k
Code Size: 2k
Category:

Communication Document

Development Platform:

Matlab

gmsk.m：Code Content
  1. %++++++++++++++++++++++++++++++++++++++++++++++++++++++
  2. % Filename: gmsk.m
  3. %
  4. % Description: Generate baseband GMSK signal
  5. %
  6. % Author: Jeff Laster, MPRG, Virginia Tech 
  7. %++++++++++++++++++++++++++++++++++++++++++++++++++++++
  9. t0 = clock;
  11. %++++++++++++++++++++++++++++++++
  12. %  Constant Definition
  13. %++++++++++++++++++++++++++++++++
  15. DRate = 1; % data rate
  16. M = 18; % no. of sample per bit
  18. %++++++++++++++++++++++++++++++++
  19. % Variable Definition
  20. %++++++++++++++++++++++++++++++++
  22. N = 12753; % no. of bits for simulation
  23. %EbNodB = 10; % Eb/No (dB scale)
  24. BT = 0.3; % Bandwidth*Period (cannot change )
  26. %++++++++++++++++++++++++++++++++++
  27. %  Program Control Variables
  28. %++++++++++++++++++++++++++++++++++
  30. fig_no = 1; % initialize fig no.
  32. %+++++++++++++++++++++++++++++++++++
  33. % Derived variables
  34. %+++++++++++++++++++++++++++++++++++
  36. T = 1/DRate; % data period
  37. Ts = T/M; % sampling period ( discrete-time )
  38. %EbNo = 10^(EbNodB/10); % Eb/No (linear scale )
  40. %+++++++++++++++++++++++++++++++++++++
  41. % Generates Input Random Data
  42. %+++++++++++++++++++++++++++++++++++++
  43. [data, RandData] = pn_gen(N,M);
  45. %dodata
  47. %+++++++++++++++++++++++++++++++++++++
  48. % Generate GMSK Pulse Response
  49. %+++++++++++++++++++++++++++++++++++++
  50. [g,fig_no] = pulse_re(BT,T,Ts,M,fig_no);
  52. %+++++++++++++++++++++++++++++++++++++++++
  53. % Generate Gaussian Filtered Pulses
  54. %+++++++++++++++++++++++++++++++++++++++++
  55. [b,delay_GF,fig_no] = gauss_pu(g,data,1,fig_no);
  57. %+++++++++++++++++++++++++++++++++++++++++
  58. % Generate GMSK signal
  59. %+++++++++++++++++++++++++++++++++++++++++
  60. [s,fig_no] = gmsk_sig(b,Ts,M,fig_no);
  62. clear T DRate BT
  64. %data=zeros(size(data));
  65. %for i=1:N
  66. % for k=1:M
  67. % data(i*M+M/2-1+k)=Rdiff_mse(i);
  68. % end
  69. %end
  71. Time = etime(clock,t0)