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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 3G开发 > 查看源码
costas_loop_sim.cpp
资源名称：simulation_wireless_c++.rar　[点击查看]
上传用户：jtjnyq9001
上传日期：2014-11-21
资源大小：3974k
文件大小：9k
源码类别：

3G开发

开发平台：

Visual C++

costas_loop_sim.cpp：源码内容
  1. //
  2. //  File = costas_loop_sim.cpp
  3. //
  5. #define SIM_NAME "CostasLoop"
  6. #define SIM_TITLE "Simulation of Costas Loop"
  8. #include "global_stuff.h" 
  10. #include "qpsk_theory.h"
  11. #include "bitgen.h"
  12. #include "bitwav.h"
  13. #include "siganchr.h"
  14. #include "bpsk_mod_bp.h"
  15. #include "phase_rotate.h"
  16. #include "add_gaus_noise.h"
  17. #include "butt_filt_iir.h"
  18. #include "spec_analyzer.h"
  19. #include "quad_carr_genie.h"
  20. #include "costas_loop.h"
  21. #include "contin_delay_T.h"
  22. #include "mixer_bp.h"
  23. #include "quad_mixer_bp.h"
  24. #include "integ_dump_slice.h"
  25. #include "qpskoptbitdem.h"
  26. #include "ber_ctr.h"
  27. #include "qpsk_err_ctr.h"
  29. main()
  30. {
  31. #include "sim_preamble.cpp"
  32.   //=========================================================================
  33.   //  Misc special processing
  35.   //-----------------------------------------------------
  36.   // create plot file of theoretical PSD for QPSK signal
  38.   QpskPsd(  0.0,  // carrier_freq
  39.             1.0,  // symb_energy,
  40.             1.0,  // symb_intvl,
  41.             -4.0, // beg_freq,
  42.             4.0,  // end_freq,
  43.             801,  // num_pts,
  44.             true, // plot_in_db,
  45.             "qpsk_psd_theor.txt"); //out_filename
  47.   //-----------------------------------------------------
  48.   // create plot file of BER for ideal QPSK
  50.   QpskBer(  -10.0,  // beg_ebno,
  51.             18.0, // end_ebno,
  52.             401,  // num_pts,
  53.             "theor_qpsk_ber.txt");  // out_filename
  54.   //=========================================================================
  55.   //  Allocate signals
  57.   BIT_SIGNAL(bit_vals);
  58.   FLOAT_SIGNAL(baseband_wave);
  59.   FLOAT_SIGNAL(modulated_signal);
  60.   FLOAT_SIGNAL(pll_err_sig);
  61.   FLOAT_SIGNAL(pll_filt_err_sig);
  62.   FLOAT_SIGNAL(pll_ref_sig);
  63.   FLOAT_SIGNAL(pll_vco_phase_sig);
  64.   FLOAT_SIGNAL(pll_vco_freq_sig);
  65.   FLOAT_SIGNAL(squared_signal);
  66.   FLOAT_SIGNAL(error_signal);
  67.   FLOAT_SIGNAL(filt_squared_signal);
  68.   FLOAT_SIGNAL(recov_carrier_sig);
  69.   FLOAT_SIGNAL(i_recov_carrier_sig);
  70.   FLOAT_SIGNAL(q_recov_carrier_sig);
  71.   FLOAT_SIGNAL(noisy_sig);
  72.   //FLOAT_SIGNAL(sqrd_noisy_sig);
  73.   BIT_SIGNAL(symb_clock);
  74.   FLOAT_SIGNAL(power_meas_sig);
  76.   //============================================================
  77.   //  Construct, initialize and connect models
  80.   BitGener* bit_gen = new BitGener( "bit_gen",
  81.                                       CommSystem,
  82.                                       bit_vals );
  84.   BitsToWave* wave_gen = new BitsToWave( "wave_gen",
  85.                                            CommSystem,
  86.                                            bit_vals,
  87.                                            baseband_wave,
  88.                                            symb_clock );
  90.   SignalAnchor* temp_anchor_1 = new SignalAnchor( "temp_anchor_1",
  91.                                                 CommSystem,
  92.                                                 baseband_wave,
  93.                                                 //q_baseband_wave,
  94.                                                 0.0078125, //samp_intvl
  95.                                                 4096 ); //block_size
  97.   BpskBandpassModulator* bpsk_mod = new BpskBandpassModulator( "bpsk_mod",
  98.                                                 CommSystem,
  99.                                                 baseband_wave,
  100.                                                 modulated_signal);
  101.   
  102.   AdditiveGaussianNoise<float>* agn_source = 
  103.                   new AdditiveGaussianNoise<float>( 
  104.                                                 "agn_source",
  105.                                                 CommSystem,
  106.                                                 modulated_signal,
  107.                                                 noisy_sig,
  108.                                                 power_meas_sig);
  110.   SpectrumAnalyzer<float>* spec_analyzer = 
  111.                   new SpectrumAnalyzer<float>( 
  112.                                                   "spec_analyzer",
  113.                                                   CommSystem,
  114.                                                   modulated_signal );
  116. //   BandpassSquaringPLL* carr_recov_loop = new BandpassSquaringPLL( "carr_recov_loop",
  117. //                                                                     CommSystem,
  118. //                                                                     noisy_sig,
  119. //                                                                     sqrd_noisy_sig,
  120. //                                                                     pll_err_sig,
  121. //                                                                     pll_filt_err_sig,
  122. //                                                                     pll_ref_sig,
  123. //                                                                     pll_vco_freq_sig,
  124. //                                                                     pll_phase_sig,
  125. //                                                                     recov_carrier_sig);
  127. //   ContinuousDelay< float >* loop_delay_elem = new ContinuousDelay<float>( "loop_delay_elem",
  128. //                                   CommSystem,
  129. //                                   squared_signal,
  130. //                                   delay_squared_signal);
  133.    CostasLoop* carr_recov_loop = new CostasLoop( "carr_recov_loop",
  134.                                                 CommSystem,
  135.                                                 modulated_signal,
  136.                                                 pll_err_sig,
  137.                                                 pll_filt_err_sig,
  138.                                                 pll_ref_sig,
  139.                                                 pll_vco_freq_sig,
  140.                                                 pll_vco_phase_sig,
  141.                                                 recov_carrier_sig);
  143. //  ButterworthFilterByIir<float>* double_freq_filter = 
  144. //            new ButterworthFilterByIir<float>( "double_freq_filter",
  145. //                                        CommSystem,
  146. //                                        squared_signal,
  147. //                                        filt_squared_signal );
  149. //  SpectrumAnalyzer<float>* spec_anal_2 = 
  150. //                  new SpectrumAnalyzer<float>( 
  151. //                                                  "spec_anal_2",
  152. //                                                  CommSystem,
  153. //                                                  sqrd_noisy_sig );
  155.   SpectrumAnalyzer<float>* spec_anal_3 = 
  156.                   new SpectrumAnalyzer<float>( 
  157.                                                   "spec_anal_3",
  158.                                                   CommSystem,
  159.                                                   pll_err_sig );
  161.   SpectrumAnalyzer<float>* spec_anal_4 = 
  162.                   new SpectrumAnalyzer<float>( 
  163.                                                   "spec_anal_4",
  164.                                                   CommSystem,
  165.                                                   pll_filt_err_sig );
  167. //  QuadCarrierRecovGenie* carrier_recovery = new QuadCarrierRecovGenie( "carrier_recovery",
  168. //                                                                   CommSystem,
  169. //                                                                   i_recov_carrier_sig,
  170. //                                                                   q_recov_carrier_sig);
  171.   
  172. //  QuadBandpassMixer* quad_dem = new QuadBandpassMixer( "quad_dem",
  173. //                                                   CommSystem,
  174. //                                                   noisy_sig,
  175. //                                                   i_recov_carrier_sig,
  176. //                                                   q_recov_carrier_sig,
  177. //                                                   i_demod_wave,
  178. //                                                   q_demod_wave );
  180. //  IntegrateDumpAndSlice* i_bit_slicer = new IntegrateDumpAndSlice( "i_bit_slicer",
  181. //                                           CommSystem,
  182. //                                           i_demod_wave,
  183. //                                           symb_clock,
  184. //                                           i_samp_wave,
  185. //                                           i_bit_decis);
  187. //  IntegrateDumpAndSlice* q_bit_slicer = new IntegrateDumpAndSlice( "q_bit_slicer",
  188. //                                           CommSystem,
  189. //                                           q_demod_wave,
  190. //                                           q_symb_clock,
  191. //                                           q_samp_wave,
  192. //                                           q_bit_decis);
  194. //  BerCounter* i_ber_ctr = new BerCounter( "i_ber_ctr",
  195. //                                          CommSystem,
  196. //                                          i_bits,
  197. //                                          i_bit_decis );
  199. //  BerCounter* q_ber_ctr = new BerCounter( "q_ber_ctr",
  200. //                                          CommSystem,
  201. //                                          q_bits,
  202. //                                          q_bit_decis );
  204. //  QpskErrorCounter* ser_ctr = new QpskErrorCounter( "ser_ctr",
  205. //                                                    CommSystem,
  206. //                                                    i_bits,
  207. //                                                    q_bits,
  208. //                                                    i_bit_decis,
  209. //                                                    q_bit_decis );
  211.   //=============================================================
  212.   #include "sim_postamble.cpp"
  213.   return 0;
  214. }