Matlab中求功率谱密度的代码

[code]

clear;

%随机产生长度大于1000的‘0’、‘1’信号序列，对其进行QPSK调制%定义待仿真序列的维数N

global N

N=2000;

%定义产生‘1’的概率为p

global p

p=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]);

hold on;

plot(insert_qpsk1(1:90));

xlabel('实部信号');

axis([0 100 -1.5 1.5]);

title('QPSK插值后序列');

subplot(2,1,2);

plot(insert_qpsk2(1:90),'yo');

axis([0 100 -1.5 1.5]);

hold on;

plot(insert_qpsk2(1:90));

xlabel('虚部信号');

axis([0 100 -1.5 1.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');

hold on;

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');

hold on;

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)))));

grid on;

xlabel('频率');

ylabel('dB');

title('QPSK功率谱密度');

[/code]

此程序来源于网络作者：蜗_牛