Matlab中求功率谱密度的代码.docVIP

[code] clear;%随机产生长度大于1000的‘0’、‘1’信号序列，对其进行QPSK调制%定义待仿真序列的维数 Nglobal NN=2000;%定义产生‘1’的概率为 pglobal pp=0.5;%产生随机二进制序列s_qpsk=randsrc(1,N,[1,0;p,1-p]);%画出生成的随机序列图figure(1);stem(s_qpsk);axis([0 50 -0.5 1.5]);xlabel(维数N)ylabel(信号强度)title(0/1等概分布的信号)%********QPSK信号的数字调制********[m_qpsk1,m_qpsk2]=qpsk_modulation(s_qpsk);figure(2);plot(m_qpsk1,m_qpsk2,r*);axis([-2 2 -2 2]);title(QPSK的信号空间图);%********插值，相邻信号间插入7个零点********insert_qpsk1=upsample(m_qpsk1,8);insert_qpsk2=upsample(m_qpsk2,8);%画出插值后的序列figure(3);subplot(2,1,1);plot(insert_qpsk1(1:90),ro);axis([0 100 -1.5 1.5]);