Matlab中求功率谱密度的代码

[code]clear;%随机产生长度大于1000的‘0’、‘1’信号序列，对其进行QPSK调制%定义待仿真序列的维数NglobalNN=2000;%定义产生‘1’的概率为pglobalpp=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_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);plot(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_qpsk1(1:80),'o');legend('滤波输出信号','输入信号');title('通过平方根升余弦滤波器滤波得到QPSK实部输出信号（10个周期）');subplot(2,1,2);plot(n,out_qpsk2(1:100),'.-r');holdon;stem(m,insert_qpsk2(1:80),'y');legend('滤波输出信号','输入信号');title('通过平方根升余弦滤波器滤波得到QPSK虚部输出信号（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]此程序来源于网络作者：蜗_牛