%——三種調(diào)制方式的功率譜密度比較
  BPSK=out_sig;
  QPSK=out_symi+j*out_symq;
  QAM16=out_symii+j*out_symqq;
  F_BPSK=fftshift(fft(BPSK));     %——頻譜對稱平移
  F_QPSK=fftshift(fft(QPSK));     %——頻譜對稱平移
  F_QAM16=fftshift(fft(QAM16));   %——頻譜對稱平移
 
  %——三種功率譜密度圖，橫軸采樣點數(shù)，縱軸功率譜密度  
  figure(15);
  f=18*(-5000:4999)/40000;
  semilogy(f,abs(F_BPSK(15001:25000))/max(abs(F_BPSK)),'r');
  hold on;
  f1=18*(-5000:2:4999)/40000;
  semilogy(f1,abs(F_QPSK(7501:12500))/max(abs(F_QPSK)),'b');
  hold on;
  f2=18*(-5000:4:4999)/40000;
  semilogy(f2,abs(F_QAM16(3751:6250))/max(abs(F_QAM16)),'y');
  grid on;
  title('BPSK,QPSK,16QAM調(diào)制功率譜密度比較');
  ylabel('功率譜密度(db)');
  xlabel('頻率(HZ)');
  legend('F_BPSK','F_QPSK','F_QAM16');