hilpsk.asv

MATLAB用于AWGN信道中产生512点算法!

%compute hilbert transform of xn


%************************ Preparation part *****************************
NUM = 1024;            % The number of FFT point
sr=1000000.0;          % Symbol rate
fc = 1000000;          % The frequency of the carrier
fs = 32000000;         % The sample rate
ml=1;                  % Number of modulation levels
br=sr.*ml;             % Bit rate (=symbol rate in this case)
nd=NUM*fc/fs;         % The number of symbols
n =fc/sr;              % The number of the carrier in a symbol
ebno=10;               % Eb/N0
IPOINT=32;             % Number of oversamples
Ts = 1/sr;             % The period of the sourse signal

t=0:1/fs:1/fc*nd*n-1/fs;    % The sample point in one carrier period
fcos = cos(2*pi*fc*t);      % cosin carrier

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

irfn=21;                                 % Number of filter taps
alfs=0.33;                                % Rolloff factor
[xh]=hrollfcoef(irfn,IPOINT,sr,alfs,1);  % Transmitter filter coefficients
%fprintf('%f\n',xh);

%*********************** Data generation **********************************
data=rand(1,nd)>0.5;                     % rand: built in function
fprintf('%d',data);

%*********************** Bpsk modulation **********************************

data1=data.*2-1;
[data2]=oversample(data1,nd,IPOINT);
data3=conv(data2,xh);                                %(conv: build in function)
%fprintf('%f\n',data3);
data4=data3(irfn*IPOINT/2-2:NUM+irfn*IPOINT/2-3);    % date4 is the baseband singal of BPSK 
modulation = data4.*fcos;                            % The modulation is the BPSK singal
hxn=hilbert(modulation);
yn=imag(hxn);
%y=ceil(1000.*yn);
fprintf('%d,\n',y);
fid=fopen('hilpsk.h','a');
fprintf(fid,'%f\n',yn);
fclose(fid);

a=