matlab qpsk调制程序,QPSK调制解调完整程序（配有自己的注释）

clc;

clear all;

%假定接收端已经实现载波同步，位同步(盲信号解调重点要解决的问题：载波同步(costas环(未见到相关代码))，位同步(Gardner算法(未见相关代码))，帧同步)

% carrier frequency  for modulation and demodulation

fc=5e6;

%QPSK transmitter

data=5000  ;   %码数率为5MHZ     %原码个数

rand_data=randn(1,5000);

for  i=1:data

if rand_data(i)>=0.5

rand_data(i)=1;

else

rand_data(i)=0;

end

end

%seriel to parallel        %同时单极性码转为双极性码

for  i=1:data

if rem(i,2)==1

if  rand_data(i)==1

I(i)=1;

I(i+1)=1;

else

I(i)=-1;

I(i+1)=-1;

end

else

if rand_data(i)==1

Q(i-1)=1;

Q(i)=1;

else

Q(i-1)=-1;

Q(i)=-1;

end

end

end

% zero insertion   ，此过程称为成形。成形的意思就是实现由消息到波形的转换，以便发射，脉冲成形应该是在基带调制之后。

zero=5;         %sampling  rate  25M HZ  ,明白了，zero为过采样率。它等于 采样率fs/码速率。

for  i=1:zero*data     % 采样点数目=过采样率*原码数目

if rem(i,zero)==1

Izero(i)=I(fix((i-1)/zero)+1);

Qzero(i)=Q(fix((i-1)/zero)+1);

else

Izero(i)=0;

Qzero(i)=0;

end

end

%pulse shape filter， 接着，将进行低通滤波，因为 随着传输速率的增大，基带脉冲的频谱将变宽

%如果不滤波(如升余弦滤波)进行低通滤波，后面加载频的时候可能会出现困难。

%平方根升余弦滤波器

% psf=rcosfir(rf,n_t,rate,