MATLAB仿真资料

MATLAB仿真资料 sr=256000.0; % Symbol rate ml=4; % ml:Number of modulation levels (BPSK:ml=1, QPSK:ml=2, 16QAM:ml=4) br=sr .* ml; % Bit rate nd = 1000; % Number of symbols that simulates in each loop ebn0=6; % Eb/N0 IPOINT=8; % Number of oversamples

% Program 3-1
% bpsk.m
%
% Simulation program to realize BPSK transmission system
%
% Programmed by H.Harada and T.Yamamura
%

%******************** Preparation part **********************

sr=256000.0; % Symbol rate
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

%******************* Filter initialization ********************

irfn=21;     % Number of filter taps          
alfs=0.5;    % Rolloff factor
[xh] = hrollfcoef（irfnIPOINTsralfs1）;   %Transmitter filter coefficients  
[xh2] = hrollfcoef（irfnIPOINTsralfs0）;  %Receiver filter coefficients


%******************** START CALCULATION *********************
nloop=100;  % Number of simulation loops

noe = 0;    % Number of error data
nod = 0;    % Number of transmitted data

for iii=1:nloop
    
%******************** Data generation ********************************  

	data=rand（1nd）>0.5;  % rand: built in function

%******************** BPSK Modulation ***********************  

    data1=data.*2-1;
	[data2] = oversamp（ data1 nd  IPOINT） ;
	data3 = conv（data2xh）;  % conv: built in function

%****************** Attenuation Calculation *****************
	
    spow=sum（data3.*data3）/nd;
	attn=0.5*spow*sr/br*10.^（-ebn0/10）;
	attn=sqrt（attn）;
   
%********************** Fading channel **********************

  % Generated data are fed into a fading simulator
  % In the case of BPSK only Ich data are fed into fading counter
  % [ifadeqfade]=sefade（data3zeros（1length（data3））itaudlvlth1n0itnd1now1length（data3）tstpfdflat）;
  
  % Updata fading counter
  %itnd1 = itnd1+ itnd0;

%************ Add White Gaussian Noise （AWGN） ***************
	
    inoise=randn（1length（data3））.*attn;  % randn: built in function
	data4=data3+inoise;
	data5=conv（data4xh2）;  % conv: built in function

	sampl=irfn*IPOINT+1;
	data6 = data5（sampl:8:8*nd+sampl-1）;
    
%******************** BPSK Demodulation *********************

    demodata=data6 > 0;

%******************** Bit Error Rate （BER） ******************

    noe2=sum（abs（data-demodata））;  % sum: built in function
	nod2=length（data）;  % length: built in function
	noe=noe+noe2;
	nod=nod+nod2;

	fprintf（‘%d\t%e\n‘iiinoe2/nod2）;
end % for iii=1:nloop    

%********************** Output result ***************************

ber = noe/nod;
fprintf（‘%d\t%d\t%d\t%e\n‘ebn0noenodnoe/nod）;
fid = fopen（‘BERbpsk.dat‘‘a‘）;
fprintf（fid‘%d\t%e\t%f\t%f\t\n‘ebn0noe/nodnoenod）;
fclose（fid）;

%******************** end of file ***************************
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     目录          0  2008-10-24 21:26  modulation

     文件        375  2005-04-24 20:39  modulation\BERbpsk.dat

     文件         44  2005-04-24 21:25  modulation\BERbpskfad.dat

     文件         43  2005-04-24 21:27  modulation\BERqpsk.dat

     文件         44  2005-04-24 21:28  modulation\BERqpskfad.dat

     文件       2753  2001-11-30 06:35  modulation\bpsk.m

     文件       4484  2001-11-30 06:30  modulation\bpsk_fading.m

     文件       2838  2005-04-24 21:24  modulation\bpsk_test.m

     文件        658  2001-11-26 18:37  modulation\comb.m

     文件        562  2001-11-30 07:56  modulation\compconv.m

     文件        684  2001-11-30 07:56  modulation\compoversamp.m

     文件        743  2001-11-30 07:55  modulation\delay.m

     文件       2167  2001-11-27 21:49  modulation\fade.m

     文件       1052  2001-11-30 07:55  modulation\gaussf.m

     文件       3387  2001-11-30 07:42  modulation\gmsk.m

     文件       5109  2001-11-30 07:42  modulation\gmsk_fading.m

     文件       1652  2001-11-27 01:37  modulation\hrollfcoef.m

     文件       2933  2001-11-30 06:31  modulation\msk.m

     文件       2967  2001-11-30 06:31  modulation\msk2.m

     文件       4691  2001-11-30 06:31  modulation\msk2_fading.m

     文件       4657  2001-11-30 06:31  modulation\msk_fading.m

     文件       3182  2001-11-30 06:31  modulation\oqpsk.m

     文件       4790  2001-11-30 07:09  modulation\oqpsk_fading.m

     文件        497  2001-11-27 09:39  modulation\oversamp.m

     文件        507  2001-11-30 07:54  modulation\oversamp2.m

     文件       2780  2001-11-30 06:42  modulation\qam16.m

     文件       4592  2001-11-30 06:44  modulation\qam16_fading.m

     文件       1729  2001-11-30 07:54  modulation\qamdemod.m

     文件       1156  2001-11-30 07:53  modulation\qammod.m

     文件       3055  2001-11-30 06:30  modulation\qpsk.m

............此处省略8个文件信息