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=