%用m語言編制主程序如下：?
clear?all;?
clc;?
?
%----------------------preparation?part?-----------------------------?
sr=256000.0;????????????????????????????%?symboll?rate?
ml=2;???????????????????????????????????????%?number?of?modulation?levels?
br=sr*ml;??%?bit?rate?
nd=100;????%?number?of?symbol?
ebno=8;????%?Ebno?
?
%-----------------------------------------------Filter?initialization-----------------?
irfn=21;?
IPOINT=8;?
alfs=0.5;?
[xh]=hrollfcoef（irfnIPOINTsralfs1）;???????%?Transmitting?Filter?
[xh2]=hrollfcoef（irfnIPOINTsralfs0）;??????%?Received?Filter?
?
%------------------------------?Spreading?code?initialization--------------------?
user=1;???%?number?of?users?
seq=1;??%?1:m-sequence?????2:Gold????3:orthogonal?Gold?
stage=3;?
ptap1=[1?3];?
ptap2=[2?3];?
regi1=[1?1?1];?
regi2=[1?1?1];?
?
%?-----------------------------?Generation?of?the?spreading?code?--------?
switch?seq?
????case?1?
????????code=mseq（stageptap1regi1user）;?
????case?2?
????????m1=mseq（stageptap1regi1）;?
????????m2=mseq（stageptap2regi2）;?
????????code=goldseq（m1m2user）;?
????case?3?
????????m1=mseq（stageptap1regi1）;?
????????m2=mseq（stageptap2regi2）;?
????????code=[goldseq（m1m2user）zeros（user1）];?
end?
code=code*2-1;?
clen=length（code）;?
?
%?----------------------------------?Fading?initialization?---------------?
rfade=0;???%?Rayleigh?fading??0:nothing??1:consider?
itau=[08];???%?delay?time?
dlvl1=[0.040.0];??%?attenuation?level?
n0=[67];?
th1=[0.00.0];?
itnd1=[30014004];?
now1=2;?
tstp=1/sr/IPOINT/clen;??%?time?resolution?
fd=160;?
flat=1;?
itnde1=nd*IPOINT*clen*30;?
?
%?-----------------------------?Start?calculation?-------------?
nloop=1000;?
noe=0;?
nod=0;?
?
for?ii=1:nloop?
?????
????%?--------------------------?Transmitter?---------------------?
????data=rand（usernd*ml）>0.5;?
????[ichqch]=qpskmod（datauserndml）;?
?????
????[ich1qch1]=spread（ichqchcode）;?
?????
????[ich2qch2]=compoversamp2（ich1qch1IPOINT）;?
????[ich3qch3]=compconv2（ich2qch2xh）;?
?????
????if?user==1?
????????ich4=ich3;?
????????qch4=qch3;?
????else?
????????ich4=sum（ich3）;?
????????qch4=sum（qch3）;?
????end?
?
?%?--------------------------?Fading?channel?-------------------?
?if?rfade==0?
?????ich5=ich4;?
?????qch5=qch4;?
?else?
?????[ich5qch5]=sefade（ich4qch4itaudlvl1th1n0itnd1...?
?????????now1length（ich4）tstpfdflat）;?
?????itnd1=itnd1+itnde1;?
?end?
??
?%?----------------------------------?Receiver??-----------------------------?
?spow=sum（rot90（ich3.^2+qch3.^2））/nd;?
?attn=sqrt（0.5*spow*sr/br*10^（-ebno/10））;?
??
?[ich6qch6]=comb2（ich5qch5attn）;?
??
?[ich7qch7]=compconv2（ich6qch6xh2）;?%?filter?
??
?samp1=irfn*IPOINT+1;?
?ich8=ich7（:samp1:IPOINT:IPOINT*nd*clen+samp1-1）;?
?qch8=qch7（:samp1:IPOINT:IPOINT*nd*clen+samp1-1）;?
??
?[ich9?qch9]=despread（ich8qch8code）;?
??
?demodata=qpskdemod（ich9qch9userndml）;???%?QPSK?demodulation?
??
?%--