《Matlab中求功率谱密度的代码》由会员分享，可在线阅读，更多相关《Matlab中求功率谱密度的代码（5页珍藏版）》请在金锄头文库上搜索。
1、codeclear;%随机产生长度大于1000的0信号序列，Xt其进行QPS颁制%定义待仿真序列的维数NglobalNN=2000;%定义产生的概率为1pglobalpp=0.5;%产生随机二进制序列s_qpsk=randsrc(1,N,1,0;p,1-p);%画出生成的随机序列图figure(1);stem(s_qpsk);axis(050-0.51.5）;xlabel（维数N）ylabel（信号强度）title（0/1等概分布的信号）%*QPSK信号的数字调制*m_qpsk1,m_qpsk2=qpsk_modulation(s_qpsk);figure(2);plot(m_qpsk1,m_
2、qpsk2,r*);axis(-22-22);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(0100-1.51.5);holdon;plot(insert_qpsk1(1:90);xlabel(实部信号);axis(0100-1.51.5);title(QPSK插值后序列)；subplot(2,1,2);pl
3、ot(insert_qpsk2(1:90),yo);axis(0100-1.51.5);holdon;plot(insert_qpsk2(1:90);xlabel(虚部信号);axis(0100-1.51.5);%*升余弦滤波器滤波*out_qpsk1=rise_cos(insert_qpsk1,N,8*N);out_qpsk2=rise_cos(insert_qpsk2,N,8*N);%画出滤波后的信号figure(5);subplot(2,1,1);n=1:100;plot(n,out_qpsk1(1:100),.-r);holdon;m=25:104;stem(m,insert_qpsk
4、1(1:80),o);legend(滤波输出信号,输入信号)；title(通过平方根升余弦滤波器滤波得到QPS5部输出信号(10个周期);subplot(2,1,2);plot(n,out_qpsk2(1:100),.-r);holdon;stem(m,insert_qpsk2(1:80),y);legend(滤波输出信号,输入信号)；title(通过平方根升余弦滤波器滤波得到QPSKg部输出信号(10个周期);%*输出信号眼图*%滤波后两路信号合并，表达成复数形式eyediagram(out_qpsk1,5*8);title(QPSK实部眼图);eyediagram(out_qpsk2,5*8);title(QPSK虚部眼图)；%*输出信号功率谱密度*out_qpsk=out_qpsk1+i*out_qpsk2;R_I=xcorr(out_qpsk);power_qpsk=fft(R_I);figure(8);plot(10*log10(abs(power_qpsk(1:(length(power_qpsk)+1)/2)-max(10*log10(abs(power_qpsk(1:(length(power_qpsk)+1)/2);gridon;xlabel(频率)；ylabel(dB);title(QPSK功率谱密度);/code此程序来源于网络蜗_牛
