qpsk频谱图matlab,QPSK调制得出的调制信号频谱波形多了一个点频信号

clc;

clear all;

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

%carrier frequency  for modulation and demodulation

fc=2.5e6;

f_code=2.4e6;

fs=24e6;

%QPSK transmitter

data=2500  ;                %原码个数

rand_data=round(rand(1,2500));    %产生5000个随机数 均值为0  方差为1

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

for  i=1:data

if rem(i,2)==1    %如果i是奇数

if  rand_data(i)==1   %如果第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   ，此过程称为成形。成形的意思就是实现由消