test.m

国外大学的同志作出来的好程序

clear all, close all, clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%OFDM Blocks (Turbo+QPSK)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sim=construct_sim;
ofdm_sim=construct_ofdm_sim;
sim=set_comm_ofdm_parameters(sim);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Nfft=512; Ncp=200;
sim.packet_time=7.5;      %if 0, then need calculation; if not 0, frame_num will be calculated.
sim.blank_time=300e-3;       %if 0 then no quiesence period.
sim.lfm_type='lfm_doppler50';      %if empty then, no lfm.
sim.lfm_time=300e-3;
sim.frame_time = (Nfft+Ncp)/sim.symbol_rate;
sim.frame_num =   0;                  %if 0, calculation from packet_time.
sim.coding_type='802_16_ctc_random';     %if empty then, no coding.
sim.coding_rate=3/4;
sim.mapping_type='qam64';
[sim, err]=chk_parameters(sim);
if err==-1,
    disp('Errors in paramter-setting.');
    return;
end

ofdm_sim.Nfft=Nfft;
ofdm_sim.Ncp=Ncp;                      %Number of prefix
ofdm_sim.Np=62;                        %Number of pilot tones
ofdm_sim.Nd=510-62;                    %Number of data carriers.
ofdm_sim.pilot_type='mseq1';           %Pilot type
ofdm_sim.preamble_type='mseq1';        %Preamble type.

[ofdm_sim, sim, err]=chk_ofdm_parameters(ofdm_sim, sim);
if err==-1,
    disp('Errors in paramter-setting of OFDM systems.');
    return;
end

%disp_sim(sim, ofdm_sim);
RNG_seed=10816600;
rand('state', RNG_seed);


sim_vec=[];
waves=[];
[packet_rf_data, info, coded_bits, baseband_symbols, ofdm_baseband_symbols]=single_ofdm_packet(sim, ofdm_sim);

packet_rf_data=single(packet_rf_data);

sim_vec=[sim_vec; sim];
waves=[waves; packet_rf_data];

figure
plot([1: length(packet_rf_data)]/sim.sampling_rate, packet_rf_data);
set(gca, 'Xlim', [0 length(packet_rf_data)/sim.sampling_rate])
xlabel('time (s)')