untitled2.m

来自「本程序实现产生几种调制方式的眼图并且对其进行比较」· M 代码 · 共 40 行

%——以下求出三种调制方式的功率谱密度
A0=fft(sym0_fc,8000);
Paa0=A0.*conj(A0)/8000;            %－－BPSK调制信号序列的功率谱密度函数

A1=fft(symi_fc+symq_fc,8000);
Paa1=A1.*conj(A1)/8000;            %－－QPSK调制信号序列的功率谱密度函数

A2=fft(symii_fc+symqq_fc,8000);
Paa2=A2.*conj(A2)/8000;            %－－16QAM调制信号序列的功率谱密度函数 

w=zeros(1,16000);
f=8*(-4000:4000)/8000;             %－－频域采样

figure(15);
w(1:8000)=Paa0(8000:-1:1);
w(8001:16000)=Paa0(1:8000);
a0=10*log(w(4000:12000));         
plot(f,a0,'r');                     %——作出－4HZ到4HZ的BPSK调制信号序列的功率谱密度函数图
grid on;
hold on;

w(1:8000)=Paa1(8000:-1:1);
w(8001:16000)=Paa1(1:8000);
a1=10*log(w(4000:12000));
plot(f,a1,'b');                      %——作出－4HZ到4HZ的QPSK调制信号序列的功率谱密度函数图
hold on;

w(1:8000)=Paa2(8000:-1:1);
w(8001:16000)=Paa2(1:8000);
a2=10*log(w(4000:12000));
plot(f,a2,'y');                       %——作出－4HZ到4HZ的16QAM调制信号序列的功率谱密度函数图
hold off;
legend('BPSK','QPSK','16QAM');
title('三种调制方式功率谱密度函数');
xlabel('信号频率f(HZ)');
ylabel('信号功率密度(db)');