QPSK调制方式频带分析

clc
clear all
fs=1000;
Rs=100;
fc=100;
oversample=fix（fs/Rs）;
numsym=3000;
symbol=randi（413000）-1;
data_mod=pskmod（symbol4pi/4）;
i_data=real（data_mod）;
q_data=imag（data_mod）;
sample_i=[];
sample_q=[];
for i=1:numsym
    sample_i=[sample_ii_data（i）*ones（1oversample）];
    sample_q=[sample_qq_data（i）*ones（1oversample）];
end
  num_sample=length（sample_i）;  
ts=（1:num_sample）/fs; 
i_carrier=cos（2*pi*fc*ts）;
q_carrier=sin（2*pi*fc*ts）;
i_re=i_carrier.*sample_i;
q_re=q_carrier.*sample_q;
data_re=