GMSK MSK 各种数字调制MATLAB源码

% % README Chapter 3 % % by Hiroshi Harada % % If you have any bugs and questions in our simulation programs, please e-mail % to harada@ieee.org. We try to do our best to answer your questions. % In this directory, we can find the twenty-eight files. The relationship between file name and the number of program written in the book is shown in as follows. Program3-1 bpsk.m Program3-2 bpsk_fading.m Program3-3 hrollfcoef.m Program3-4 oversamp.m Program3-5 qpsk.m Program3-6 qpsk_fading.m Program3-7 compconv.m Program3-8 compoversamp.m Program3-9 qpskmod.m Program3-10 qpskdemod.m Program3-11 oqpsk.m Program3-12 oqpsk_fading.m Program3-13 msk.m Program3-14 msk_fading.m Program3-15 msk2.m Program3-16 msk2_fading.m Program3-17 oversamp2.m Program3-18 gmsk.m Program3-19 gmsk_fading.m Program3-20 gaussf.m Program3-21 qam16.m Program3-22 qam16_fading.m Program3-23 qammod.m Program3-24 qamdemod.m Program2-4 comb.m Program2-5 fade.m Program2-6 sefade.m Program2-7 delay.m If you would like to try to use the above programs by using MATLAB.First of all, please copy all of files to your created adequate directory. Then, you start to run MATLAB and you can see the following command prompt in the command window. >> Next, you can go to the directory that have all of programs in this section by using change directory (cd) command. If you copy all of files to /matlabR12/work/chapter3, you only type the following command. >> cd /matlabR12/work/chapter3 In this directory, we can find fourteen main functions, bpsk.m, bpsk_fading.m, qpsk.m, qpsk_fading.m, oqpsk.m, oqpsk_fading.m, msk.m, msk_fading.m, msk2.m, msk2_fading.m, gmsk.m, gmsk_fading.m, qam16.m and qam16_fading.m ######################################################### (1) Simulation of "bpsk.m" ######################################################### This program simulates the transmission performance of BPSK under Additive White Gausian Noise (AWGN) environment. (a) Set paremeters First of all, we set simulation parameters in "bpsk.m". %******************** Preparation part ********************** sr=256000.0; % Symbol rate 256 ksymbol/s ml=1; % Number of modulation levels br=sr.*ml; % Bit rate (=symbol rate in this case) nd = 1000; % Number of symbols that simulates in each loop ebn0=3; % Eb/N0 IPOINT=8; % Number of oversamples %******************** START CALCULATION ********************* nloop=100; % Number of simulation loops (b) Type just the following command >> clear >> bpsk (c) Then, you can see the following simulation result on your command window. (example) 3 2275 100000 2.275000e-002 where first number 3 is Eb/No, second number 2275 is the number of error data, third number 100000 is the number of transmitted data, and fourth number 2.275000e-002 is bit error rate (BER) performance. And, the simulation result is stored in the file (BERbpsk.dat). ######################################################### (2) Simulation of "bpsk_fading.m" ######################################################### This program simulates the transmission performance of BPSK under Rayleigh fading environment. (a) Set paremeters First of all, we set simulation parameters in "bpsk_fading.m". %******************** Preparation part ********************** sr=256000.0; % Symbol rate 256 ksymbol/s ml=1; % Number of modulation levels br=sr.*ml; % Bit rate (=symbol rate in this case) nd = 100; % Number of symbols that simulates in each loop ebn0=10; % Eb/N0 IPOINT=8; % Number of oversamples %******************* Fading initialization ******************** % If you use fading function "sefade", you can initialize all of parameters. % Otherwise you can comment out the following initialization. % The detailed explanation of all of valiables are mentioned in Program 2-8. % Time resolution tstp=1/sr/IPOINT; % Arrival time for each multipath normalized by tstp % If you would like to simulate under one path fading model, you have only to set % direct wave. itau = [0]; % Mean power for each multipath normalized by direct wave. % If you would like to simulate under one path fading model, you have only to set % direct wave. dlvl = [0]; % Number of waves to generate fading for each multipath. % In normal case, more than six waves are needed to generate Rayleigh fading n0=[6]; % Initial Phase of delayed wave % In this simulation four-path Rayleigh fading are considered. th1=[0.0]; % Number of fading counter to skip itnd0=nd*IPOINT*100; % Initial value of fading counter % In this simulation one-path Rayleigh fading are considered. % Therefore one fading counter are needed. itnd1=[1000]; % Number of directwave + Number of delayed wave % In this simulation one-path Rayleigh fading are considered now1=1; % Maximum Doppler frequency [Hz] % You can insert your favorite value fd=160; % You can decide two mode to simulate fading by changing the variable flat % flat : flat fading or not % (1->flat (only amplitude is fluctuated),0->nomal(phase and amplitude are fluctutated) flat =1; %******************** START CALCULATION ********************* nloop=1000; % Number of simulation loops (b) Type just the following command >> clear >> bpsk_fading (c) Then, you can see the following simulation result on your command window. (example) 10 2143 100000 2.143000e-002 The meaning of each value is the same of the result from "bpsk.m". The simulation result is stored in the file (BERbpskfad.dat). ######################################################### (3) Simulation of "qpsk.m" ######################################################### This program simulates the transmission performance of QPSK under Additive White Gausian Noise (AWGN) environment. (a) Set paremeters First of all, we set simulation parameters in "qpsk.m". %******************** Preparation part ********************** sr=256000.0; % Symbol rate 256 ksymbol/s ml=2; % Number of modulation levels br=sr.*ml; % Bit rate (=symbol rate in this case) nd = 1000; % Number of symbols that simulates in each loop ebn0=3; % Eb/N0 IPOINT=8; % Number of oversamples %******************** START CALCULATION ********************* nloop=100; % Number of simulation loops (b) Type just the following command >> clear >> qpsk (c) Then, you can see the following simulation result on your command window. (example) 3 4475 200000 2.237500e-002 The meaning of each value is the same of the result from "bpsk.m". The simulation result is stored in the file (BERqpsk.dat). ######################################################### (4) Simulation of "qpsk_fading.m" ######################################################### This program simulates the transmission performance of QPSK under Rayleigh fading environment. (a) Set paremeters First of all, we set simulation parameters in "qpsk_fading.m". %******************** Preparation part ********************** sr=256000.0; % Symbol rate 256 ksymbol/s ml=2; % Number of modulation levels br=sr.*ml; % Bit rate (=symbol rate in this case) nd = 1000; % Number of symbols that simulates in each loop ebn0=10; % Eb/N0 IPOINT=8; % Number of oversamples %******************* Fading initialization ******************** % If you use fading function "sefade", you can initialize all of parameters. % Otherwise you can comment out the following initialization. % The detailed explanation of all of valiables are mentioned in Program 2-8. % Time resolution tstp=1/sr/IPOINT; % Arrival time for each multipath normalized by tstp % If you would like to simulate under one path fading model, you have only to set % direct wave. itau = [0]; % Mean power for each multipath normalized by direct