msk_sim.cpp
上传用户：jtjnyq9001
上传日期：2014-11-21
资源大小：3974k
文件大小：8k
- //
- // File = msk_sim.cpp
- //
- #define SIM_NAME "MskSim "
- #define SIM_TITLE "Complex Baseband Simulation for MSK "
- #include "global_stuff.h"
- #include "gausrand.h"
- #include "bitgen.h"
- #include "bitwav.h"
- #include "mskmod.h"
- #include "add_gaus_noise.h"
- #include "spec_analyzer.h"
- #include "phase_genie.h"
- #include "quaddem.h"
- #include "qpskdem.h"
- #include "msk_genie.h"
- #include "integ_dump_slice.h"
- #include "discrete_delay_T.h"
- #include "ber_ctr.h"
- #include "siganchr.h"
- #include "msk_theory.h"
- #include "sigstuff.h"
- #include "parmfile.h"
- #include "model_graph.h"
- #include "sigplot.h"
- #include "reports.h"
- #include "exec.h"
- #include "psstream.h"
-
- //#ifdef _DEBUG
- // ofstream DebugFile("MskModulator.dbg", ios::out);
- //#endif
- //=========================================================
- // Global Stuff
- //ParmFile ParmInput("MskModulator.dat");
- //ofstream LongReport;
- //ofstream ShortReport;
- //PracSimStream DetailedResults;
- //PracSimStream BasicResults;
- //PracSimStream ErrorStream;
- //PracSimModel *CommSystem;
- //ActiveSystemGraph ActSystGraph;
- //Executive Exec;
- //PracSimModel *PrevModelConstr;
- //PracSimModel *ActiveModel;
- //SignalPlotter SigPlot;
- //int PassNumber;
- //int MaxPassNumber;
- //int EnclaveNumber;
- //=========================================================
- main()
- {
- #include "sim_preamble.cpp"
- //=========================================================================
- // Misc special processing
- //-----------------------------------------------------
- // create plot file of theoretical PSD for MSK signal
- MskPsd( 0.0, // carrier_freq
- 1.0, // symb_energy,
- 1.0, // symb_intvl,
- -4.0, // beg_freq,
- 4.0, // end_freq,
- 801, // num_pts,
- true, // plot_in_db,
- "msk_psd_theor.txt "); //out_filename
- //-----------------------------------------------------
- // create plot file of BER for ideal MSK
- MskBer( -10.0, // beg_ebno,
- 18.0, // end_ebno,
- 401, // num_pts,
- "theor_msk_ber.txt "); // out_filename
- //=========================================================================
- // Allocate signals
- BIT_SIGNAL(i_bits);
- BIT_SIGNAL(q_bits);
- BIT_SIGNAL(q_bits_delayed);
- FLOAT_SIGNAL(i_baseband_wave);
- FLOAT_SIGNAL(q_baseband_wave);
- BIT_SIGNAL(i_bit_clock);
- BIT_SIGNAL(q_bit_clock);
- FLOAT_SIGNAL(i_filtered_wave);
- FLOAT_SIGNAL(q_filtered_wave);
- COMPLEX_SIGNAL(modulated_signal);
- COMPLEX_SIGNAL(filtered_sig);
- COMPLEX_SIGNAL(noisy_sig);
- COMPLEX_SIGNAL(phase_ref_sig);
- FLOAT_SIGNAL(power_meas_sig);
- FLOAT_SIGNAL(mag_signal);
- FLOAT_SIGNAL(phase_signal);
- FLOAT_SIGNAL(i_demod_wave);
- FLOAT_SIGNAL(q_demod_wave);
- FLOAT_SIGNAL(i_corr_sig);
- FLOAT_SIGNAL(q_corr_sig);
- FLOAT_SIGNAL(i_filt_demod_wave);
- FLOAT_SIGNAL(q_filt_demod_wave);
- FLOAT_SIGNAL(i_samp_demod_wave);
- FLOAT_SIGNAL(q_samp_demod_wave);
- BIT_SIGNAL(i_bit_decis);
- BIT_SIGNAL(q_bit_decis);
- //============================================================
- // Construct, initialize and connect models
- EnclaveNumber = 0;
- BitGener* i_bit_gen = new BitGener( "i_bit_gen ",
- CommSystem,
- i_bits );
- BitGener* q_bit_gen = new BitGener( "q_bit_gen ",
- CommSystem,
- q_bits );
- DiscreteDelay<bit_t>* q_bit_delay = new DiscreteDelay<bit_t>(
- "q_bit_delay ",
- CommSystem,
- q_bits,
- q_bits_delayed);
- BitsToWave* i_wave_gen = new BitsToWave( "i_wave_gen ",
- CommSystem,
- i_bits,
- i_baseband_wave,
- i_bit_clock );
- BitsToWave* q_wave_gen = new BitsToWave( "q_wave_gen ",
- CommSystem,
- q_bits,
- q_baseband_wave,
- q_bit_clock );
- MskModulator* msk_mod = new MskModulator( "msk_mod ",
- CommSystem,
- i_baseband_wave,
- q_baseband_wave,
- modulated_signal,
- mag_signal,
- phase_signal );
-
- AdditiveGaussianNoise<std::complex<float> >* agn_source =
- new AdditiveGaussianNoise<std::complex<float> >(
- "agn_source ",
- CommSystem,
- modulated_signal,
- noisy_sig,
- power_meas_sig);
- SpectrumAnalyzer<std::complex<float> >* modulator_spec_an =
- new SpectrumAnalyzer<std::complex<float> >(
- "modulator_spec_an ",
- CommSystem,
- modulated_signal );
- PhaseRecoveryGenie* phase_genie = new PhaseRecoveryGenie( "phase_genie ",
- CommSystem,
- phase_ref_sig);
- QuadratureDemod* quad_demod = new QuadratureDemod( "quad_demod ",
- CommSystem,
- noisy_sig,
- phase_ref_sig,
- i_demod_wave,
- q_demod_wave);
- // QpskDemodulator* qpsk_dem = new QpskDemodulator( "qpsk_dem ",
- // CommSystem,
- // noisy_sig,
- // i_demod_wave,
- // q_demod_wave);
- SignalAnchor* temp_anchor = new SignalAnchor( "temp_anchor_1 ",
- CommSystem,
- i_baseband_wave,
- q_baseband_wave,
- 0.0625, //samp_intvl
- 4096 ); //block_size
- MskShapeGenie* msk_genie = new MskShapeGenie( "msk_genie ",
- CommSystem,
- i_corr_sig,
- q_corr_sig);
- IntegrateDumpAndSlice* i_bit_slicer = new IntegrateDumpAndSlice( "i_bit_slicer ",
- CommSystem,
- i_demod_wave,
- //i_corr_sig,
- i_bit_clock,
- i_samp_demod_wave,
- i_bit_decis );
- BerCounter* i_ber_ctr = new BerCounter( "i_ber_ctr ",
- CommSystem,
- i_bits,
- i_bit_decis );
- IntegrateDumpAndSlice* q_bit_slicer = new IntegrateDumpAndSlice( "q_bit_slicer ",
- CommSystem,
- q_demod_wave,
- //q_corr_sig,
- q_bit_clock,
- q_samp_demod_wave,
- q_bit_decis );
- BerCounter* q_ber_ctr = new BerCounter( "q_ber_ctr ",
- CommSystem,
- q_bits_delayed,
- q_bit_decis );
- //=============================================================
- #include "sim_postamble.cpp"
- return 0;
- }
English
